Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 963,761 T→C 7.6% intergenic (+93/‑67) rpsA → / → ihfB 30S ribosomal subunit protein S1/integration host factor (IHF), DNA‑binding protein, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913963,7610TC 7.6% 336.5 / 14.6 118intergenic (+93/‑67)rpsA/ihfB30S ribosomal subunit protein S1/integration host factor (IHF), DNA‑binding protein, beta subunit
Reads supporting (aligned to +/- strand):  ref base T (49/60);  new base C (2/7);  total (51/67)
Fisher's exact test for biased strand distribution p-value = 2.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.71e-01

Only 100 of 118 total aligned reads displayed.

AACTTCTCCAACAACGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTAAA  >  NC_000913/963615‑963907
                                                                                                                                                  |                                                                                                                                                  
aaCTTCTCCAACAACGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACg                                                                                                                                                 <  2:667262/150‑1 (MQ=255)
  cTTCTCCAACAACGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGATGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGgc                                                                                                                                               >  1:534272/1‑150 (MQ=255)
     ctcCAACAACGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTGCGGTATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGgccgc                                                                                                                                            >  2:278048/1‑150 (MQ=255)
      tcCAACAACGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGgccgcc                                                                                                                                           <  1:990289/150‑1 (MQ=255)
       ccAACAACGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGccc                                                                                                                                          >  1:1607120/1‑150 (MQ=255)
       ccAACAACGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGccc                                                                                                                                          >  2:1418759/1‑150 (MQ=255)
       ccAACAACGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGccc                                                                                                                                          >  2:1548074/1‑150 (MQ=255)
       ccAACAACGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGccc                                                                                                                                          >  1:1350291/1‑150 (MQ=255)
       ccAACAACGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTATACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGccc                                                                                                                                          >  1:1627805/1‑150 (MQ=255)
        caacaaCGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCt                                                                                                                                         >  2:1696460/1‑150 (MQ=255)
        caacaaCGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCt                                                                                                                                         >  2:175229/1‑150 (MQ=255)
        caacaaCGCAATGGCTGAAGCATTCAAAGCAGCTAACGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGTAGTTTGTTGAGTTTACTTGACAGATTTCAGGTTTCTTCCTGTAATCAAGCACTAAGGGCGGCCGTAGCCGCCCt                                                                                                                                         <  2:249421/150‑1 (MQ=255)
         aacaacGCAATGGCTGAAGCTTTCAAAGCAGCTAAATGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCCGTAGCCGCCCtt                                                                                                                                        <  1:371193/150‑1 (MQ=255)
          acaacGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTa                                                                                                                                       >  2:1376099/1‑150 (MQ=255)
            aacGCAATGTCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTCATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGGAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAAt                                                                                                                                     <  2:121968/150‑1 (MQ=255)
            aacGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGAGGCTACGGCCGCCCTTAAt                                                                                                                                     >  2:425993/1‑150 (MQ=255)
             acGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATc                                                                                                                                    <  2:1105872/150‑1 (MQ=255)
               gCAATGGCTGAAGCTTTCAAAGCAGGTAAAGTCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCCGTAGCCGCCCTTagtgct                                                                                                                                  <  2:466346/150‑6 (MQ=255)
               gCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCTAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCCGTAGCCGCCCTTagtgct                                                                                                                                  <  2:1197308/150‑6 (MQ=255)
               gCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTTCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCaa                                                                                                                                  <  1:601117/150‑1 (MQ=255)
                 aaTGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATg                                                                                                                                >  1:1199154/1‑150 (MQ=255)
                 aaTGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCCGTAGCCGCCCTTagtgcttg                                                                                                                                <  2:1267403/150‑8 (MQ=255)
                 aaTGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCCGTAGCCGCCCTTagtgcttg                                                                                                                                <  1:459649/150‑8 (MQ=255)
                  aTGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATgc                                                                                                                               <  1:757939/150‑1 (MQ=255)
                    ggCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATgcag                                                                                                                             >  1:1004385/1‑150 (MQ=255)
                    ggCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATgcag                                                                                                                             <  1:411085/150‑1 (MQ=255)
                     gCTGAAGCTTTCAAATCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTTCTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATgcagc                                                                                                                            <  1:960749/150‑1 (MQ=255)
                       tGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATATTTATTCAGAAGTTTGTAGAGTTTACTTGACAGATTGCAGCTATCTTCATGTAATCAATCACTAAGGGAGGCTACGGACGCCCTTAATCAATGCAGCaa                                                                                                                          >  1:131095/1‑150 (MQ=255)
                           gCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCATCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAAcagc                                                                                                                      >  1:278300/1‑150 (MQ=255)
                            cTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAAcagca                                                                                                                     >  1:177562/1‑150 (MQ=255)
                            cTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAAcagca                                                                                                                     <  2:1010189/150‑1 (MQ=255)
                            cTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAAcagca                                                                                                                     <  2:672339/150‑1 (MQ=255)
                             tttCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATTCAGCAAcagcag                                                                                                                    >  2:1624472/1‑150 (MQ=255)
                                 aaaGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGc                                                                                                                <  2:388336/150‑1 (MQ=255)
                                   agcagcTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCtt                                                                                                              <  2:374259/150‑1 (MQ=255)
                                        cTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAAttt                                                                                                         >  2:717932/1‑150 (MQ=255)
                                         tAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTg                                                                                                        <  1:1418759/150‑1 (MQ=255)
                                           aaGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGcc                                                                                                      <  1:845802/150‑1 (MQ=255)
                                               cGAGTAATTCTCTGACTCTTCGTGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCTTCCTGTAATCAAGCACTCAGGGCGTCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTa                                                                                                  <  1:280236/150‑1 (MQ=255)
                                                        ctctGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCgg                                                                                         <  2:1149240/150‑1 (MQ=255)
                                                        ctctGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCgg                                                                                         >  1:636117/1‑150 (MQ=255)
                                                        ctctGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTAGTCCTGTAATCAAGAACTAAGGGCGGCTACGGCCGACCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAAACgg                                                                                         >  1:296184/1‑150 (MQ=255)
                                                            taCTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCggagga                                                                                     <  1:1075638/149‑1 (MQ=255)
                                                            gACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCggagga                                                                                     >  1:1302257/1‑150 (MQ=255)
                                                            gACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCggagga                                                                                     <  1:1376099/150‑1 (MQ=255)
                                                             aCTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGaa                                                                                    <  1:175229/150‑1 (MQ=255)
                                                                ctTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCa                                                                                 <  1:1378713/150‑1 (MQ=255)
                                                                  tCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATg                                                                               >  1:430101/1‑150 (MQ=255)
                                                                     ggATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGAcc                                                                            <  1:987099/150‑1 (MQ=255)
                                                                      gATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCaa                                                                          >  1:1594036/1‑150 (MQ=255)
                                                                        tttttATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAg                                                                         >  1:485286/1‑150 (MQ=255)
                                                                           ttATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCa                                                                      <  1:1106109/150‑1 (MQ=255)
                                                                           ttATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCa                                                                      <  1:58434/150‑1 (MQ=255)
                                                                           ttATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCa                                                                      <  2:1600424/150‑1 (MQ=255)
                                                                            tATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAg                                                                     >  2:1148743/1‑150 (MQ=255)
                                                                            tATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGCAACCGGAGGAATCATGACAAAGTCAg                                                                     >  1:1611137/1‑150 (MQ=255)
                                                                               tCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGaaa                                                                  >  2:1329417/1‑149 (MQ=255)
                                                                               tCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAAt                                                                  <  1:287905/150‑1 (MQ=255)
                                                                                ccGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAAtt                                                                 <  1:632905/150‑1 (MQ=255)
                                                                                ccGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAAtt                                                                 <  1:977487/150‑1 (MQ=255)
                                                                                   aaGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGAt                                                              <  1:1548074/150‑1 (MQ=255)
                                                                                    aGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATa                                                             <  2:1275598/150‑1 (MQ=255)
                                                                                     gTTTGTTGAGTCTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAg                                                            <  2:1627805/150‑1 (MQ=255)
                                                                                      tttGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGa                                                           <  1:1696460/150‑1 (MQ=255)
                                                                                      tttGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGa                                                           <  1:1624472/150‑1 (MQ=255)
                                                                                        tgttgAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCGTGACCAAGTCAGAATTGATAGaaa                                                         <  1:278048/150‑1 (MQ=255)
                                                                                        tgttgAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGaaa                                                         >  2:1108503/1‑150 (MQ=255)
                                                                                        tgttgAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTAGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGaaa                                                         <  2:1609839/150‑1 (MQ=255)
                                                                                           tgAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGAc                                                      <  1:324322/150‑1 (MQ=255)
                                                                                            gAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGTCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACt                                                     <  1:346697/150‑1 (MQ=255)
                                                                                               tttACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGc                                                  <  2:1685145/150‑1 (MQ=255)
                                                                                                ttACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGcc                                                 <  2:518395/150‑1 (MQ=255)
                                                                                                   cTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCAcc                                              >  1:751184/1‑150 (MQ=255)
                                                                                                    ttGACAGATAGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGTCCGCCATTAATCAATGAATCAAAAGCAGCCGCTTAATTTGCCTTTAAGGAACCGCAGGAATCATGTCCAAATCAGAATTGATAGAAAGACTTGCCAccc                                             >  2:1087104/1‑150 (MQ=255)
                                                                                                     tGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCa                                            >  1:461491/1‑150 (MQ=255)
                                                                                                     tGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCa                                            >  2:219686/1‑150 (MQ=255)
                                                                                                      gACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAg                                           <  1:717932/150‑1 (MQ=255)
                                                                                                           aTTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATc                                      <  2:1700091/150‑1 (MQ=255)
                                                                                                             tGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGc                                    <  1:368707/150‑1 (MQ=255)
                                                                                                                  gTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACAtt                               >  2:1470671/1‑150 (MQ=255)
                                                                                                                      cGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCg                           <  2:1607120/150‑1 (MQ=255)
                                                                                                                       gTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCg                                     <  2:247645/139‑1 (MQ=255)
                                                                                                                         ccTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCa                        <  2:278300/150‑1 (MQ=255)
                                                                                                                           tGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGAGTTGCCACCCAGCAATCGCACATTCCCGCCAAg                      <  2:745265/150‑1 (MQ=255)
                                                                                                                           tGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAg                      <  2:537876/150‑1 (MQ=255)
                                                                                                                                tCAAGCACTAATGGCGGCTACTGCCGCCCTTAATCAATGCATCAACAGCAGCCTCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGt                 <  1:934650/150‑1 (MQ=255)
                                                                                                                                tCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGt                 >  2:918707/1‑150 (MQ=255)
                                                                                                                                 cAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGtt                >  1:220795/1‑150 (MQ=255)
                                                                                                                                  aaGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTg               <  1:1148743/150‑1 (MQ=255)
                                                                                                                                   aGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGa              >  2:1313808/1‑150 (MQ=255)
                                                                                                                                      aCTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTTAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGa           <  1:929930/150‑1 (MQ=255)
                                                                                                                                       ttAAGGGCGGCCGTAGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGAt          >  1:56677/2‑150 (MQ=255)
                                                                                                                                         aaGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGc        >  2:1399741/1‑150 (MQ=255)
                                                                                                                                         aaGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGc        >  1:485969/1‑150 (MQ=255)
                                                                                                                                          aGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCa       <  1:886866/150‑1 (MQ=255)
                                                                                                                                          aGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCa       >  1:149520/1‑150 (MQ=255)
                                                                                                                                           gggCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATTACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAg      >  2:1183322/1‑150 (MQ=255)
                                                                                                                                           gggCGGCTACGGCCGACCTTAATCAATGCAGCAACAGTAGCTGCTTAATTTGCCTTTATGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAg      <  1:114818/150‑1 (MQ=255)
                                                                                                                                            ggcggcTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAGt     >  1:604177/1‑150 (MQ=255)
                                                                                                                                               ggcTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTaaa  <  2:149520/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AACTTCTCCAACAACGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTAAA  >  NC_000913/963615‑963907

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: