Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,781,022 T→G 8.9% intergenic (+5/‑193) lldD → / → trmL L‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,781,0220TG 8.9% 1691.6 / 43.1 619intergenic (+5/‑193)lldD/trmLL‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base T (148/399);  new base G (67/3);  total (216/403)
Fisher's exact test for biased strand distribution p-value = 1.82e-30
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 638 total aligned reads displayed.

TCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTG  >  NC_000913/3780873‑3781170
                                                                                                                                                     |                                                                                                                                                      
tCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGtt                                                                                                                                                        >  2:4692921/1‑150 (MQ=255)
 cTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGttt                                                                                                                                                       <  2:177680/150‑1 (MQ=255)
  tGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGtgtg                                                                                                                                                      >  1:6369306/1‑150 (MQ=255)
  tGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTg                                                                                                                                                      >  1:1103318/1‑150 (MQ=255)
   gATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTgc                                                                                                                                                     <  1:5281165/150‑1 (MQ=255)
   gATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTgc                                                                                                                                                     <  1:709567/150‑1 (MQ=255)
     tCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTgcgc                                                                                                                                                   <  2:4568435/150‑1 (MQ=255)
      cGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGcc                                                                                                                                                  >  1:3284442/1‑150 (MQ=255)
       gAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGccc                                                                                                                                                 >  2:3504880/1‑150 (MQ=255)
       gAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGccc                                                                                                                                                 <  2:5760525/150‑1 (MQ=255)
        aaaaaGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGcccc                                                                                                                                                <  2:6483227/150‑1 (MQ=255)
             gagaTGAAAGTGTCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATATTCGTTTGCGCCCCCTCAc                                                                                                                                           <  2:5279975/150‑1 (MQ=255)
             gagaTGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCAc                                                                                                                                           <  2:3610737/150‑1 (MQ=255)
               gaTGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCAccc                                                                                                                                         <  2:5045643/150‑1 (MQ=255)
                aTGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCccccccccccc                                                                                                                                        >  1:2812115/1‑149 (MQ=255)
                aTGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCGCAccca                                                                                                                                        >  2:301769/1‑149 (MQ=255)
                 tGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTa                                                                                                                                       <  1:4338777/150‑1 (MQ=255)
                  gAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTaa                                                                                                                                      <  1:1768445/150‑1 (MQ=255)
                    aaGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAAcc                                                                                                                                    <  1:1610638/150‑1 (MQ=255)
                     aGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAAccc                                                                                                                                   <  2:3857388/150‑1 (MQ=255)
                      ctGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCct                                                                                                                                  <  2:2271838/149‑1 (MQ=255)
                       tGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCACCCAAAaccac                                                                                                                                 >  2:1561434/1‑148 (MQ=255)
                        ggCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCCCACCCCAAcccaca                                                                                                                                >  1:3580516/1‑147 (MQ=255)
                         gCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCtctc                                                                                                                               <  1:2655535/150‑1 (MQ=255)
                         gCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCtctc                                                                                                                               <  2:367315/150‑1 (MQ=255)
                           gATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCACCCCAACCCGCTccc                                                                                                                             >  2:898981/1‑150 (MQ=255)
                           gATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCCCCCCCGAAACCACTccc                                                                                                                             >  2:631706/1‑150 (MQ=255)
                           gATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCACACCCCAAaccacaacc                                                                                                                             >  2:6134269/1‑144 (MQ=255)
                               aCGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑g                                                                                                                        <  2:1505893/150‑1 (MQ=255)
                               aCGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGACGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑g                                                                                                                        <  1:4542239/150‑1 (MQ=255)
                                cGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑gg                                                                                                                       <  2:6226521/150‑1 (MQ=255)
                                cGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCAAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑gg                                                                                                                       <  2:4105013/150‑1 (MQ=255)
                                  cTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGA‑GGga                                                                                                                     <  2:3267102/150‑1 (MQ=255)
                                    gACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGgaga                                                                                                                   <  1:1013525/150‑1 (MQ=255)
                                     aCTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGA‑GGGagag                                                                                                                  <  2:1655041/150‑1 (MQ=255)
                                      cgggCGCGAAAGCGATCAGCGAAAGTCCGCAAGCTTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCGCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTACCCCTCTCCCTGA‑GGGAGAgg                                                                                                                 <  2:1147831/148‑1 (MQ=255)
                                      cTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTGA‑GGGAGAgg                                                                                                                 <  1:4805967/150‑1 (MQ=255)
                                      cTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAgg                                                                                                                 <  1:2592571/150‑1 (MQ=255)
                                      cTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAgg                                                                                                                 <  2:3448623/150‑1 (MQ=255)
                                         gcgcGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGa                                                                                                              <  1:4959458/150‑1 (MQ=255)
                                         gcgcGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGa                                                                                                              <  1:2318985/150‑1 (MQ=255)
                                         gcgcGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGa                                                                                                              <  1:6116300/150‑1 (MQ=255)
                                          cgagAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGACGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAgggtta                                                                                                             <  2:5938702/147‑4 (MQ=255)
                                           gcgAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCAAAACCACTCCCTCA‑GGGAGAGGGGAcc                                                                                                            >  1:2189748/1‑150 (MQ=255)
                                           gcgAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTGGCGCCCCCACCCCCCAAACCCCACCCACC‑GGGAGGGGGGccc                                                                                                            >  1:4437590/1‑150 (MQ=255)
                                               aaTCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTc                                                                                                        <  2:1368617/150‑1 (MQ=255)
                                                 tCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCgg                                                                                                      <  1:6158330/150‑1 (MQ=255)
                                                    aTCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCAAACCCACACCCACA‑GGGAGAGGGGACCCTTCGgcgc                                                                                                   >  1:6423115/1‑150 (MQ=255)
                                                    aTCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCACACCCAAAACCACTACCACA‑GGGAGAGGGGACCCGTCGgcgc                                                                                                   >  1:3071170/1‑150 (MQ=255)
                                                       aGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCtgt                                                                                                <  1:390118/150‑1 (MQ=255)
                                                       aGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCACCCAAAACCACACCCCCA‑GGGTGAGGGGGCCCGTCGGCGCtgg                                                                                                >  2:5854658/1‑149 (MQ=255)
                                                         cGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCCCCCCCCAAACCACAACCACA‑GGGGGCGGGGACCCGTCGGCGCtggat                                                                                              >  1:487729/1‑150 (MQ=255)
                                                         cGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCAAAACCAAACCCGAAACCACC‑GGGAGAGTGGACCGGTCGGCGCggtct                                                                                              >  1:6520958/1‑148 (MQ=255)
                                                             aTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTAc                                                                                          <  2:4402358/150‑1 (MQ=255)
                                                               tACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTc                                                                                        <  1:3388002/150‑1 (MQ=255)
                                                                  gCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCCCACCCCAAACCCCACCCCCA‑CGGCGGGGGGGCCCGTGGGCGCGGGATTTACTcccc                                                                                     >  2:5650022/1‑149 (MQ=255)
                                                                   cAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCCCACCCGAAACCTCACCCACA‑GGGAGAGGGGACCCTTCGGCGCTGTATGTAcaaccgc                                                                                    >  2:1044710/1‑144 (MQ=255)
                                                                    aaGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCa                                                                                   <  2:4018167/150‑1 (MQ=255)
                                                                      gATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCAct                                                                                 <  1:472119/150‑1 (MQ=255)
                                                                         tCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑Tg                                                                             <  2:6004551/150‑1 (MQ=255)
                                                                          cGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGa                                                                            <  1:1948087/150‑1 (MQ=255)
                                                                                             tAAAGAGTTGCCTGCGGCACTGGCTCCCTTGGCGAAAGGGATTGCGGCATTGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGTGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCtttt                                                         <  1:5554412/150‑1 (MQ=255)
                                                                                              aaaGAGTTGCCTGCGGCACTGGCTCCCTTGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCttttt                                                        <  1:893785/150‑1 (MQ=255)
                                                                                              aaaGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCttttt                                                        <  1:810182/150‑1 (MQ=255)
                                                                                                agagTTGCCTGCGGCACTGACTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCttttttt                                                      <  1:5401229/150‑1 (MQ=255)
                                                                                                       ccTGCGGCACTGGCTCCCATGGCGACAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTc                                               <  2:6560173/150‑1 (MQ=255)
                                                                                                       ccTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTc                                               <  1:1286573/150‑1 (MQ=255)
                                                                                                       ccTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTc                                               <  1:5864900/150‑1 (MQ=255)
                                                                                                       ccTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTc                                               <  1:1582170/150‑1 (MQ=255)
                                                                                                       ccTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCACACCCGAACCCTCACCCTCA‑GGGCGAGGGGCCCGTTCGGCGCGGTTTGTACTCCCTCACAC‑AGAAACCACACCCCACTCTTTTTTTCTCCCTc                                               >  2:3808895/1‑150 (MQ=255)
                                                                                                          gCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCGcc                                            <  2:5503935/150‑1 (MQ=255)
                                                                                                           cGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCCAACCCAAAACCAAAACCAAA‑GGGGGAGGGGCCCCGTAGGCGCGGTATGTAATACCACACAC‑GGAAAAGACAACCCCCTCATTTTTTATACCTCGccc                                           >  1:5952637/1‑150 (MQ=255)
                                                                                                             gCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTTCGCCCCCGCACCCAAAACCACAACCACAGGGGGGGGGGGACCG‑TCGGCGGGGTTTGTAATCCCGCCCTC‑TGAAACGGCACCGCCCTGTTTTTTTATGCCAcccccca                                         >  2:3564775/1‑149 (MQ=255)
                                                                                                             gCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCCCACCCCAAACCCCAACCCCC‑GGGGGTGGGGGCCCGGCGGGGCGGTTTGTACACCCACACACAGGAAACCA‑ACCGCACTCTTTTTTTTTCCCGCGcccca                                         >  1:1486064/1‑149 (MQ=255)
                                                                                                                 tGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCgg                                     <  2:1863442/150‑1 (MQ=255)
                                                                                                                    cTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCGTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGga                                  <  2:2795028/150‑1 (MQ=255)
                                                                                                                    cTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGga                                  <  1:2006242/150‑1 (MQ=255)
                                                                                                                     tCCCATGGCGACAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGgag                                 <  2:2917136/150‑1 (MQ=255)
                                                                                                                      cccATGGCGACAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCC‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGgaga                                <  1:3797230/150‑1 (MQ=255)
                                                                                                                      cccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGgaga                                <  1:3527062/150‑1 (MQ=255)
                                                                                                                       ccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGagag                               <  2:2191922/150‑1 (MQ=255)
                                                                                                                       ccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCGCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCCCTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGagag                               <  2:6098236/150‑1 (MQ=255)
                                                                                                                        cATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCGCACCCCAACCCACTCCCCCC‑CGGGGAGGGGGCCCGTAGGCGCGGTTTTTAATACCACACAA‑AGAAACCACACCCCACACTTTTTTTATACCACCCCCCACCCGGGGGAgg                              >  2:2060425/1‑150 (MQ=255)
                                                                                                                          tGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCTCCCCGGAGGGGCGAggga                            <  1:1648025/150‑2 (MQ=255)
                                                                                                                          tGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTATTTTCTCCCTCTCCCCGGAGGGGCGAggga                            <  1:3188822/150‑2 (MQ=255)
                                                                                                                           ggCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGcc                           <  1:5152583/150‑1 (MQ=255)
                                                                                                                            gCGAAAGGGAATGCGGAATATTCGTTTGCGCCCCCACACACTAACCCTCTCCTTCC‑GGGAGTGGGGACGTTTCGGGGTTGTATTTACACCCTCACTA‑TGAAACGACACCGAACTCTTTTTTTCTCCCTCGCCCCTCCGGGGCGTGGGcgg                          <  2:2658217/150‑3 (MQ=255)
                                                                                                                             cGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCgg                         <  1:4023903/150‑1 (MQ=255)
                                                                                                                              gAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCTCACCCGAACCCGCTCCCGCA‑GGGTGCGGGGGCCGTTCGGCGCGGTTTGTACGCCCTCACTC‑GGGAACGACACCGCACACTTTTTTTCTCCCACGCCCCACCGGGGGGAGGGGCggg                        >  2:3168743/1‑150 (MQ=255)
                                                                                                                              gAAAGGGAATGCGGCATAGTCGTGTGCGCCCCCCCCCCCCAACCCACGACCACC‑GGGTGCGGGGCCCCGTAGGCGCGGGCTGTAAGCCCGCCCGC‑AGGAACGCCCCCGCCCGCGTTTGTTTTCCCGCGCCCCGCCGGGGGGAGGGCCggg                        >  2:1981272/1‑150 (MQ=255)
                                                                                                                                aaGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGt                      <  1:3687481/150‑1 (MQ=255)
                                                                                                                                   ggAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCC‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAg                   <  1:2073124/150‑1 (MQ=255)
                                                                                                                                   ggAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAg                   <  1:5705122/150‑1 (MQ=255)
                                                                                                                                    gAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCCCTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                  <  1:4117798/150‑1 (MQ=255)
                                                                                                                                         cGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaa             <  2:6045607/150‑1 (MQ=255)
                                                                                                                                          ggCATAGTCGTTTGCGCCCCCTCACCCTTACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TTAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGTGGAGAGGGCCGGGGTGAGGGGataa            <  2:487685/150‑1 (MQ=255)
                                                                                                                                          ggCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACGGTTCGGCGGTGTTTGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCGGGGGAGCGGGCCGGGGTGAGGGGaaaa            <  1:1188086/150‑1 (MQ=255)
                                                                                                                                          ggCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa            <  2:4294050/150‑1 (MQ=255)
                                                                                                                                          ggCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa            <  1:4669718/150‑1 (MQ=255)
                                                                                                                                                 tCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCa     >  2:857153/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
TCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCA‑GGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTC‑TGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTG  >  NC_000913/3780873‑3781170

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑