Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 963,762 A→G 7.6% intergenic (+94/‑66) 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,7620AG 7.6% 331.5 / 15.0 118intergenic (+94/‑66)rpsA/ihfB30S ribosomal subunit protein S1/integration host factor (IHF), DNA‑binding protein, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (49/60);  new base G (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 = 5.71e-01

Only 100 of 119 total aligned reads displayed.

AACTTCTCCAACAACGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAG  >  NC_000913/963615‑963911
                                                                                                                                                   |                                                                                                                                                     
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:1350291/1‑150 (MQ=255)
       ccAACAACGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGccc                                                                                                                                              >  2:1418759/1‑150 (MQ=255)
       ccAACAACGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGccc                                                                                                                                              >  1:1607120/1‑150 (MQ=255)
       ccAACAACGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTATACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGccc                                                                                                                                              >  1:1627805/1‑150 (MQ=255)
        caacaaCGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCt                                                                                                                                             >  2:175229/1‑150 (MQ=255)
        caacaaCGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCt                                                                                                                                             >  2:1696460/1‑150 (MQ=255)
        caacaaCGCAATGGCTGAAGCATTCAAAGCAGCTAACGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGTAGTTTGTTGAGTTTACTTGACAGATTTCAGGTTTCTTCCTGTAATCAAGCACTAAGGGCGGCCGTAGCCGCCCt                                                                                                                                             <  2:249421/150‑1 (MQ=255)
         aacaacGCAATGGCTGAAGCTTTCAAAGCAGCTAAATGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCCGTAGCCGCCCtt                                                                                                                                            <  1:371193/150‑1 (MQ=255)
            aacGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGAGGCTACGGCCGCCCTTAAt                                                                                                                                         >  2:425993/1‑150 (MQ=255)
               gCAATGGCTGAAGCTTTCAAAGCAGGTAAAGTCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCCGTAGCCGCCCTTagtgct                                                                                                                                      <  2:466346/150‑6 (MQ=255)
               gCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCTAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCCGTAGCCGCCCTTagtgct                                                                                                                                      <  2:1197308/150‑6 (MQ=255)
               gCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTTCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCaa                                                                                                                                      <  1:601117/150‑1 (MQ=255)
               gCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCCGTAGCCGCCCTTagtgct                                                                                                                                      <  2:1155942/150‑6 (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                                                                                                                                 >  2:368707/1‑150 (MQ=255)
                     gCTGAAGCTTTCAAATCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTTCTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATgcagc                                                                                                                                <  1:960749/150‑1 (MQ=255)
                       tGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATATTTATTCAGAAGTTTGTAGAGTTTACTTGACAGATTGCAGCTATCTTCATGTAATCAATCACTAAGGGAGGCTACGGACGCCCTTAATCAATGCAGCaa                                                                                                                              >  1:131095/1‑150 (MQ=255)
                           gCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCATCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAAcagc                                                                                                                          >  1:278300/1‑150 (MQ=255)
                            cTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAAcagca                                                                                                                         <  2:1010189/150‑1 (MQ=255)
                            cTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAAcagca                                                                                                                         >  1:177562/1‑150 (MQ=255)
                            cTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAAcagca                                                                                                                         <  2:672339/150‑1 (MQ=255)
                             tttCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATTCAGCAAcagcag                                                                                                                        >  2:1624472/1‑150 (MQ=255)
                             tttCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAAcagcag                                                                                                                        >  2:1617879/1‑150 (MQ=255)
                                 aaaGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGc                                                                                                                    <  2:388336/150‑1 (MQ=255)
                                      agcTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAAt                                                                                                               >  2:1295325/1‑150 (MQ=255)
                                        cTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAAttt                                                                                                             >  2:717932/1‑150 (MQ=255)
                                          aaaGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGc                                                                                                           >  1:537876/1‑150 (MQ=255)
                                               cGAGTAATTCTCTGACTCTTCGTGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCTTCCTGTAATCAAGCACTCAGGGCGTCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTa                                                                                                      <  1:280236/150‑1 (MQ=255)
                                                 aGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCATCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAg                                                                                                    <  1:438678/150‑1 (MQ=255)
                                                        ctctGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCgg                                                                                             >  1:636117/1‑150 (MQ=255)
                                                        ctctGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCgg                                                                                             <  1:711459/150‑1 (MQ=255)
                                                        ctctGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCgg                                                                                             <  2:1149240/150‑1 (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:58434/150‑1 (MQ=255)
                                                                           ttATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCa                                                                          <  1:1106109/150‑1 (MQ=255)
                                                                           ttATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCa                                                                          <  2:1600424/150‑1 (MQ=255)
                                                                            tATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAg                                                                         >  2:1148743/1‑150 (MQ=255)
                                                                               tCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGaaa                                                                      >  2:1329417/1‑149 (MQ=255)
                                                                               tCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAAt                                                                      <  1:287905/150‑1 (MQ=255)
                                                                                ccGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAAc                                                                                                <  2:1476394/123‑1 (MQ=255)
                                                                                ccGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAAc                                                                                                >  1:1476394/1‑123 (MQ=255)
                                                                                   aaGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGAt                                                                  <  1:1548074/150‑1 (MQ=255)
                                                                                    aGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATa                                                                 <  2:1275598/150‑1 (MQ=255)
                                                                                      tttGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGTCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGa                                                               <  1:1617879/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:177562/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)
                                                                                                      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)
                                                                                                                       gTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCg                                         >  1:247645/1‑139 (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)
                                                                                                                                  aaGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTg                   <  1:1148743/150‑1 (MQ=255)
                                                                                                                                   aGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGa                  >  2:1313808/1‑150 (MQ=255)
                                                                                                                                    gCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGaa                 <  2:1004385/150‑1 (MQ=255)
                                                                                                                                      aCTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTTAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGa               <  1:929930/150‑1 (MQ=255)
                                                                                                                                       ttAAGGGCGGCCGTAGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGAt              >  1:56677/2‑150 (MQ=255)
                                                                                                                                         aaGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGc            >  2:1399741/1‑150 (MQ=255)
                                                                                                                                         aaGGGCGGCCGTAGCCGCCCTTAATCAATGCAGCAACCGCAGCCGCTTAATTTGCCTTTAAGGATCCGGAGGAATCATGCCCAAGTCAGGATTGATAGAAAGACTTGCCACCCAGCCATCGCCCATTCCCGCCAAGACGGTTGAAGATGc            >  1:306041/1‑150 (MQ=255)
                                                                                                                                          aGGGCGGCTACGGCCGCCCTTAATCAATTCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCa           <  1:892875/150‑1 (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)
                                                                                                                                                   aCGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTAAAagag  >  2:1660842/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AACTTCTCCAACAACGCAATGGCTGAAGCTTTCAAAGCAGCTAAAGGCGAGTAATTCTCTGACTCTTCGGGATTTTTATTCCGAAGTTTGTTGAGTTTACTTGACAGATTGCAGGTTTCGTCCTGTAATCAAGCACTAAGGGCGGCTACGGCCGCCCTTAATCAATGCAGCAACAGCAGCCGCTTAATTTGCCTTTAAGGAACCGGAGGAATCATGACCAAGTCAGAATTGATAGAAAGACTTGCCACCCAGCAATCGCACATTCCCGCCAAGACGGTTGAAGATGCAGTAAAAGAG  >  NC_000913/963615‑963911

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: