Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 285,887 C→A 7.5% L165M (CTG→ATG)  yagG → CP4‑6 prophage; putative sugar transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913285,8870CA 7.5% 607.9 / 31.3 173L165M (CTG→ATG) yagGCP4‑6 prophage; putative sugar transporter
Reads supporting (aligned to +/- strand):  ref base C (82/78);  new base A (8/5);  total (90/83)
Fisher's exact test for biased strand distribution p-value = 5.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.49e-01

Only 100 of 176 total aligned reads displayed.

TCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGG  >  NC_000913/285739‑286032
                                                                                                                                                    |                                                                                                                                                 
tctcCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCgcgc                                                                                                                                                   <  1:655583/149‑1 (MQ=255)
 ctcCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCt                                                                                                                                                  <  3:1450531/149‑1 (MQ=255)
  tcCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTg                                                                                                                                                 >  3:990216/1‑149 (MQ=255)
  tcCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTg                                                                                                                                                 >  1:1701885/1‑149 (MQ=255)
      gACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCAGCGCGATGCCGc                                                                                                                                             <  1:1183102/149‑1 (MQ=255)
         ccTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGATAGAGTGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTgg                                                                                                                                          >  4:1097860/1‑149 (MQ=255)
         ccTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTTGCTCGCTATCAGCGGCATCGCGCTGCCGCTgg                                                                                                                                          <  3:1218719/149‑1 (MQ=255)
          cTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGt                                                                                                                                         <  1:979504/149‑1 (MQ=255)
           tGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTg                                                                                                                                        >  1:518549/1‑149 (MQ=255)
            ggTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGa                                                                                                                                       <  2:1442207/149‑1 (MQ=255)
            ggTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGa                                                                                                                                       <  4:1408098/149‑1 (MQ=255)
            ggTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGa                                                                                                                                       <  4:1497351/149‑1 (MQ=255)
            ggTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGAGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGa                                                                                                                                       <  4:1167033/149‑1 (MQ=255)
                   acCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGcatcat                                                                                                                                >  3:1515716/1‑149 (MQ=255)
                    ccTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCatcatc                                                                                                                               <  3:1270054/149‑1 (MQ=255)
                    ccTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCatcatc                                                                                                                               <  3:922754/149‑1 (MQ=255)
                     cTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCg                                                                                                                              >  1:1269124/1‑149 (MQ=255)
                      ttCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCgg                                                                                                                             <  3:1772333/149‑1 (MQ=255)
                      ttCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCgg                                                                                                                             >  4:1621641/1‑149 (MQ=255)
                              cGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAgggg                                                                                                                     <  1:1616586/149‑1 (MQ=255)
                               gTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGa                                                                                                                    <  1:1445060/149‑1 (MQ=255)
                                tGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGAc                                                                                                                   >  4:1680122/1‑149 (MQ=255)
                                        tGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGAc                                                                                                                   <  4:1835817/141‑1 (MQ=255)
                                         gcgcCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTg                                                                                                          <  4:1865356/149‑1 (MQ=255)
                                          cgcCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTgg                                                                                                         <  2:1593154/149‑1 (MQ=255)
                                            ccATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGc                                                                                                       <  3:463161/149‑1 (MQ=255)
                                               tGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTAc                                                                                                    >  3:658997/1‑149 (MQ=255)
                                                        tcatcaCCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGcc                                                                                           >  4:492027/1‑149 (MQ=255)
                                                              ccgccgACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTgg                                                                                                                                          >  1:582705/1‑96 (MQ=255)
                                                                   gACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCt                                                                                >  4:232176/1‑149 (MQ=255)
                                                                     cccGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTgg                                                                              <  4:629909/149‑1 (MQ=255)
                                                                       cGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTgggg                                                                            <  3:1680122/149‑1 (MQ=255)
                                                                        gAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGTATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTggg                                                                             <  3:1493473/147‑1 (MQ=255)
                                                                        gAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGc                                                                           >  4:180048/1‑149 (MQ=255)
                                                                          aaGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTg                                                                         <  2:1502389/149‑1 (MQ=255)
                                                                           agagCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATTTGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGa                                                                        >  3:539028/1‑149 (MQ=255)
                                                                           agagCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGa                                                                        <  4:658997/149‑1 (MQ=255)
                                                                           agagCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGa                                                                        <  4:539028/149‑1 (MQ=255)
                                                                           agagCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGa                                                                        <  3:1854797/149‑1 (MQ=255)
                                                                            gagCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAg                                                                       >  1:1470449/1‑149 (MQ=255)
                                                                             agCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGc                                                                                                                                             >  1:1340015/1‑78 (MQ=255)
                                                                                    cGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGc                                                                                                                                     <  2:1695413/79‑1 (MQ=255)
                                                                                     gCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGa                                                                        <  2:964063/139‑1 (MQ=255)
                                                                                     gCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCgtggt                                                              >  2:206728/1‑149 (MQ=255)
                                                                                       ccTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGtgc                                                            <  3:232176/149‑1 (MQ=255)
                                                                                       ccTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGtgc                                                            <  2:1154824/149‑1 (MQ=255)
                                                                                          gCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAg                                                                                                             >  1:183671/1‑97 (MQ=255)
                                                                                          gCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAg                                                                                                             <  2:183671/97‑1 (MQ=255)
                                                                                              tCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTg                                                                               <  1:1063667/123‑1 (MQ=255)
                                                                                              tCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTg                                                                               >  2:1063667/1‑123 (MQ=255)
                                                                                                 tGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGt                                                  >  2:164960/1‑149 (MQ=255)
                                                                                                  ggCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTc                                                 <  2:1269124/149‑1 (MQ=255)
                                                                                                   gcgcTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCt                                                >  1:233611/1‑149 (MQ=255)
                                                                                                   gcgcTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCt                                                >  1:1568814/1‑149 (MQ=255)
                                                                                                           tcCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCttcttc                                        <  3:1888738/149‑1 (MQ=255)
                                                                                                            ccTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACCAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCa                                       <  1:1003966/149‑1 (MQ=255)
                                                                                                             cTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCAc                                      <  3:1418458/149‑1 (MQ=255)
                                                                                                                  ggcggcggGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGtgctg                                                          >  3:1241515/1‑124 (MQ=255)
                                                                                                                  ggcggcggGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGtgctg                                                          <  4:1241515/124‑1 (MQ=255)
                                                                                                                    cggcggGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCaaa                               >  3:352967/1‑149 (MQ=255)
                                                                                                                     ggcggGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGtgctgc                                                         <  4:1479504/122‑1 (MQ=255)
                                                                                                                        gggCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTgg                                                                                                         <  1:1688572/71‑1 (MQ=255)
                                                                                                                        gggCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAgc                           >  1:607852/1‑149 (MQ=255)
                                                                                                                         ggCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTTTACGTCTGCTTCTTCACGACCAAAGAgcg                          <  3:1207553/149‑1 (MQ=255)
                                                                                                                         ggCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAgcg                          <  4:471191/149‑1 (MQ=255)
                                                                                                                           ctcgctcgCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGc                                              >  4:1380662/1‑127 (MQ=255)
                                                                                                                           ctcgctcgCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGc                                              <  3:1380662/127‑1 (MQ=255)
                                                                                                                            tcgctcgcTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCa                                                                                          >  1:1612289/1‑82 (MQ=255)
                                                                                                                            tcgctcgcTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCa                                                                                          <  2:1612289/82‑1 (MQ=255)
                                                                                                                            tcgctcgcTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAaga                             <  1:772849/143‑1 (MQ=255)
                                                                                                                            tcgctcgcTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAaga                             >  2:772849/1‑143 (MQ=255)
                                                                                                                            tcgctcgcTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTa                       >  3:490849/1‑149 (MQ=255)
                                                                                                                             cgctcgctATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAgcgc                         >  3:620706/1‑146 (MQ=255)
                                                                                                                              gctcgctATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAgcgc                         <  2:365966/145‑1 (MQ=255)
                                                                                                                               ctcgctATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATgtgcgt                                                                                   >  2:432766/1‑86 (MQ=255)
                                                                                                                               ctcgctATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATgtgcgt                                                                                   <  1:432766/86‑1 (MQ=255)
                                                                                                                                  gctATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCaaa                               >  3:1081841/1‑135 (MQ=255)
                                                                                                                                   ctATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGtgctg                                                          <  2:1679594/107‑1 (MQ=255)
                                                                                                                                    tATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTg               >  4:1314854/1‑149 (MQ=255)
                                                                                                                                    tATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTg               >  2:200893/1‑149 (MQ=255)
                                                                                                                                    tATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGAGGGCTACTTCGGCGCCAt                                                                                         >  2:554127/1‑75 (MQ=255)
                                                                                                                                     accAGCGGCAGCGCGATGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAgcgg                                                                    <  4:1055906/93‑1 (MQ=255)
                                                                                                                                     accAGCGGCAGCGCGATGCCGCTGATGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGa              >  1:357524/3‑149 (MQ=255)
                                                                                                                                     aTCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGctct                                                      <  3:746730/109‑1 (MQ=255)
                                                                                                                                      tCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAg             <  2:665113/149‑1 (MQ=255)
                                                                                                                                      tCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAg             >  2:499714/1‑149 (MQ=255)
                                                                                                                                       cAGCGGCAGCGCGATGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGcttc                                           <  3:1824834/118‑1 (MQ=255)
                                                                                                                                       cAGCGGCAGCGCGATGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCt                  >  1:1857264/1‑143 (MQ=255)
                                                                                                                                       cAGCGGCAGCGCGATGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGAGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAgg            >  2:234186/1‑149 (MQ=255)
                                                                                                                                        aGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGt           <  4:990216/149‑1 (MQ=255)
                                                                                                                                         gCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAg             >  4:957015/1‑146 (MQ=255)
                                                                                                                                         gCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTg          >  2:47686/1‑149 (MQ=255)
                                                                                                                                             cATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGc      >  4:563492/1‑149 (MQ=255)
                                                                                                                                              aTCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGcttc                                           <  4:1884714/111‑1 (MQ=255)
                                                                                                                                              aTCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAg             <  3:829287/141‑1 (MQ=255)
                                                                                                                                              aTCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAg             >  4:829287/1‑141 (MQ=255)
                                                                                                                                              aTCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGcc     >  1:939407/1‑149 (MQ=255)
                                                                                                                                               tCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCg    <  4:352967/149‑1 (MQ=255)
                                                                                                                                               tCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCg    >  2:1547713/1‑149 (MQ=255)
                                                                                                                                                cgcgCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCgg   <  1:405749/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TCTCCTGACCCTGGTCTACACCTTCGTTAACGTGCCGTACTGCGCCATGCCGGGCGTCATCACCGCCGACCCGAAAGAGCGTCACGCCCTGCAGTCCTGGCGCTTCTTCCTGGCGGCGGCGGGCTCGCTCGCTATCAGCGGCATCGCGCTGCCGCTGGTGAGCATCATCGGCAAAGGGGACGAGCAGGTGGGCTACTTCGGCGCCATGTGCGTGCTGGGGCTGAGCGGCGTGGTGCTGCTCTACGTCTGCTTCTTCACGACCAAAGAGCGCTACACCTTTGAGGTGCAGCCGGG  >  NC_000913/285739‑286032

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