Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 13,356 T→G 16.1% S398S (TCT→TCG dnaK → chaperone Hsp70, co‑chaperone with DnaJ

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_00777913,3560T→G16.1% 158.3 ‑0.1 ‑0.5285S398S (TCT→TCGdnaKchaperone Hsp70, co‑chaperone with DnaJ
Reads supporting (aligned to +/- strand):  new base (22/24):  ref base (119/120):  total (141/144)
Fisher's exact test strand distribution p-value = 8.73e-01
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 3.03e-01

Only 100 of 304 total aligned reads displayed.

TTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACA  >  NC_007779/13256‑13456
                                                                                                    |                                                                                                    
 tAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGtctc                                                                                                     >  1:7621269/1‑101 (MQ=255)
  aaCCCGGACGAGGTTGTCGCAGTTGGTGCTGCTCTGCAGGCCGGCGTTCTGCGCGGTGAGGTCAAGGACGTGCTGCTGCTCGACGTTACCCCGCTGTcgct                                                                                                    >  1:6558789/1‑101 (MQ=255)
  aaCCCGGACGAGGTTGTCGCAGTTGGTGCTGCTCTGCAGGCCGGCGTTCTGCGCGGTGAGGTCAAGGACGTGCTGCTGCTCGACGTTACCCCGCTGTcgct                                                                                                    <  1:3350634/101‑1 (MQ=255)
  aaCCCGGACGAGGTTGTCGCAGTTGGTGCTGCTCTGCAGGCCGGCGTTCTGCGCGGTGAGGTCAAGGACGTGCTGCTGCTCGACGTTACCCCGCTGTcgct                                                                                                    >  1:10783170/1‑101 (MQ=255)
  aaCCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTctct                                                                                                    <  1:2249551/101‑1 (MQ=255)
  aaCCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTctct                                                                                                    <  1:571154/101‑1 (MQ=255)
   aCCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTg                                                                                                   >  1:9544105/1‑101 (MQ=255)
   aCCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTg                                                                                                   <  1:1504022/101‑1 (MQ=255)
    cccGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTgg                                                                                                  <  1:5075433/101‑1 (MQ=255)
    cccGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTgg                                                                                                  >  1:13961272/1‑101 (MQ=255)
     ccGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTggg                                                                                                 >  1:11031843/1‑101 (MQ=255)
       ggACGAGGTTGTCGCAGTTGGTGCTGCTCTGCAGGCCGGCGTTCTGCGCGGTGAGGTCAAGGACGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTa                                                                                               >  1:11689404/1‑101 (MQ=255)
       ggACGAAGCTGTAGCCATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTa                                                                                               <  1:9780674/101‑1 (MQ=255)
       ggACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTa                                                                                               >  1:6516092/1‑101 (MQ=255)
        gACGAGGTTGTCGCAGTTGGTGCTGCTCTGCAGGCCGGCGTTCTGCGCGGTGAGGTCAAGGACGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTAt                                                                                              >  1:10641928/1‑101 (MQ=255)
        gACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTAt                                                                                              >  1:12863118/1‑101 (MQ=255)
        gACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTAt                                                                                              <  1:2370883/101‑1 (MQ=255)
         aCGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATc                                                                                             <  1:3287845/101‑1 (MQ=255)
          cGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCg                                                                                            >  1:5980314/1‑101 (MQ=255)
           gAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGa                                                                                           <  1:2980539/101‑1 (MQ=255)
             aGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGaaa                                                                                         <  1:3575592/101‑1 (MQ=255)
              gCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAAc                                                                                        >  1:4121394/1‑101 (MQ=255)
              gCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAAc                                                                                        <  1:9869435/101‑1 (MQ=255)
               cTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAAcc                                                                                       >  1:14270624/1‑101 (MQ=255)
                tGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCa                                                                                      <  1:8422680/101‑1 (MQ=255)
                 gtcgcagttggTGCTGCTCTGCAGGCCGGCGTTCTGCGCGGTGAGGTCAAGGACGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTATTGAGACCaa                                                                                     <  1:7242989/92‑2 (MQ=255)
                     gagttggTGCTGCTCTGCAGGCCGGCGTTCTGCGCGGTGAGGTCAAGGACGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGc                                                                                 <  1:11740870/96‑1 (MQ=255)
                     cagttggTGCTGCTCTGCAGGCCGGCGTTCTGCGCGGTGAGGTCAAGGACGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGc                                                                                 <  1:709630/96‑1 (MQ=255)
                      aaTCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCg                                                                                <  1:1904336/101‑1 (MQ=255)
                      aaTCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCg                                                                                <  1:9023774/101‑1 (MQ=255)
                       aTCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTGCCCCGCTGTCTCTGGGTATCGAAACCATGGGCgg                                                                               <  1:9420490/101‑1 (MQ=255)
                       aTCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCgg                                                                               >  1:13795525/1‑101 (MQ=255)
                                gctgTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATgac                                                                      <  1:12206666/101‑1 (MQ=255)
                                 ctgTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATgacg                                                                     <  1:13296784/101‑1 (MQ=255)
                                   gTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATgacgac                                                                   <  1:5500787/101‑1 (MQ=255)
                                   gTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATgacgac                                                                   <  1:6871881/101‑1 (MQ=255)
                                   gTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACATACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATgacgac                                                                   <  1:10193052/101‑1 (MQ=255)
                                    ttCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGacgacg                                                                  <  1:1986991/101‑1 (MQ=255)
                                      cAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCt                                                                <  1:14179262/101‑1 (MQ=255)
                                      cAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCt                                                                <  1:3179800/101‑1 (MQ=255)
                                       aGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTg                                                               <  1:9276060/101‑1 (MQ=255)
                                        gggTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGa                                                              <  1:5316935/101‑1 (MQ=255)
                                          gtggtgTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATc                                                            <  1:1611858/101‑1 (MQ=255)
                                           tggtgTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATcg                                                           <  1:12937975/101‑1 (MQ=255)
                                              tgtTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGa                                                        <  1:12019380/101‑1 (MQ=255)
                                                 tCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGaaaa                                                     <  1:1244024/101‑1 (MQ=255)
                                                     aCTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAacac                                                 >  1:7991856/1‑101 (MQ=255)
                                                        ggTGAGGTCAAGGACGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAAcaccac                                              <  1:283575/101‑1 (MQ=255)
                                                         gTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACt                                             <  1:10048871/101‑1 (MQ=255)
                                                         gTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACt                                             <  1:9194095/101‑1 (MQ=255)
                                                          tGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTa                                            <  1:11924110/101‑1 (MQ=255)
                                                          tGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTa                                            >  1:8183539/1‑101 (MQ=255)
                                                           gACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTAt                                           >  1:6390587/1‑101 (MQ=255)
                                                           gACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTAt                                           <  1:11408590/101‑1 (MQ=255)
                                                            aggtcaaggACGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATc                                          >  1:7386556/9‑101 (MQ=255)
                                                            aggtcaaggACGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATc                                          <  1:2649208/93‑1 (MQ=255)
                                                            aCGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACGATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATc                                          >  1:8949511/1‑101 (MQ=255)
                                                            aCGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGGAAAACACCACTATc                                          <  1:11884765/101‑1 (MQ=255)
                                                              gTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATccc                                        <  1:4865504/101‑1 (MQ=255)
                                                              gTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATccc                                        >  1:12619079/1‑101 (MQ=255)
                                                               tcaaggACGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCg                                       >  1:12262521/6‑101 (MQ=255)
                                                               tAATAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACCATCCCg                                       <  1:10118893/101‑1 (MQ=255)
                                                                 aaggACGTGCTGCTGCTCGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCGAc                                     >  1:609567/4‑101 (MQ=255)
                                                                 aaaGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGAc                                     >  1:5462089/1‑101 (MQ=255)
                                                                  aaGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGAcc                                    >  1:9117839/1‑101 (MQ=255)
                                                                  aaGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGAcc                                    <  1:5856260/101‑1 (MQ=255)
                                                                   aGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCa                                   >  1:3764561/1‑101 (MQ=255)
                                                                      cGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGc                                <  1:6207860/101‑1 (MQ=255)
                                                                      cGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGc                                >  1:7512254/1‑101 (MQ=255)
                                                                       gtgctgctgctcGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCGACCAAGcg                               >  1:4004896/4‑100 (MQ=255)
                                                                       gTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGca                               <  1:4826934/101‑1 (MQ=255)
                                                                       gTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGca                               >  1:13084181/1‑101 (MQ=255)
                                                                       gTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCCGATCGCGAAAAACACCACTATCCCGACCAAGca                               <  1:3414473/101‑1 (MQ=255)
                                                                         aCTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCTCTATCCCGACCAAGcaca                             <  1:8829845/101‑1 (MQ=255)
                                                                          ctgctgctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAg                            >  1:13634999/1‑101 (MQ=255)
                                                                          ctgctgctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAg                            >  1:11266680/1‑101 (MQ=255)
                                                                          ctgctgctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAg                            <  1:213903/101‑1 (MQ=255)
                                                                          ctgctgctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAg                            <  1:8334975/101‑1 (MQ=255)
                                                                          ctgctgctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAg                            >  1:456736/1‑101 (MQ=255)
                                                                            gctgctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGcc                          <  1:10867515/101‑1 (MQ=255)
                                                                             ctgctcGACGTTACCCCGCTGTCGCTGGGTATTGAGACCAAGGGCGGTGTCATGACCAAGCTGATCGAGCGCAACACCACTATCCCGACCAAGcgttcgga                         <  1:5584141/101‑8 (MQ=255)
                                                                                ctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGgt                      <  1:5933419/101‑1 (MQ=255)
                                                                                ctgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGgt                      >  1:3523114/1‑101 (MQ=255)
                                                                                 tgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGgtg                     <  1:12552224/101‑1 (MQ=255)
                                                                                 tgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGgtg                     <  1:9499879/101‑1 (MQ=255)
                                                                                 tgGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGgtg                     >  1:5028352/1‑101 (MQ=255)
                                                                                   gACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAACACACCACTATCCCGACCAAGCACAGCCAGGTGtt                   >  1:1745378/1‑101 (MQ=255)
                                                                                     cGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTtct                 >  1:12785404/1‑101 (MQ=255)
                                                                                      gTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTtctc                >  1:1464583/1‑101 (MQ=255)
                                                                                       ttACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTctct               <  1:4448612/101‑1 (MQ=255)
                                                                                        tACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTa              >  1:11138196/1‑101 (MQ=255)
                                                                                        tACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTa              >  1:602994/1‑101 (MQ=255)
                                                                                          ccccGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTAcc            <  1:6786855/101‑1 (MQ=255)
                                                                                           cccGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCACGACCAAGCACAGCCAGGTGTTCTCTACCg           <  1:5559886/101‑1 (MQ=255)
                                                                                              gCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTg        <  1:13550363/101‑1 (MQ=255)
                                                                                              gCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTg        >  1:13258919/1‑101 (MQ=255)
                                                                                              gCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTg        <  1:7822571/101‑1 (MQ=255)
                                                                                                 gTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAg     >  1:13443975/1‑101 (MQ=255)
                                                                                                   ctctGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGAc   >  1:10008204/1‑101 (MQ=255)
                                                                                                    tctGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACa  <  1:4415832/101‑1 (MQ=255)
                                                                                                    |                                                                                                    
TTAACCCGGACGAAGCTGTAGCAATCGGTGCTGCTGTTCAGGGTGGTGTTCTGACTGGTGACGTAAAAGACGTACTGCTGCTGGACGTTACCCCGCTGTCTCTGGGTATCGAAACCATGGGCGGTGTGATGACGACGCTGATCGCGAAAAACACCACTATCCCGACCAAGCACAGCCAGGTGTTCTCTACCGCTGAAGACA  >  NC_007779/13256‑13456

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: 
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.