Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,708,427 G→A 5.8% intergenic (‑722/+189) dppA ← / ← proK dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor/tRNA‑Pro

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,708,4270GA 5.8% 705.6 / 20.8 208intergenic (‑722/+189)dppA/proKdipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor/tRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base G (80/115);  new base A (8/4);  total (88/120)
Fisher's exact test for biased strand distribution p-value = 1.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 209 total aligned reads displayed.

AATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAG  >  NC_000913/3708280‑3708576
                                                                                                                                                    |                                                                                                                                                     
   aaGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAg                                                                                                                                                  >  1:25207/1‑150 (MQ=255)
     ggCGGTCGCGCCGCCTCCGACAATGGCCGGCACCTTTGCCTGATGCGACGCGGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCa                                                                                                                                                <  2:468103/150‑1 (MQ=255)
          tCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCtt                                                                                                                                           <  2:1798704/150‑1 (MQ=255)
          tCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCAACGAGCAGGcat                                                                                                                                           <  2:1672872/150‑3 (MQ=255)
          tCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCAACGAGCAGGcat                                                                                                                                           <  1:1776420/150‑3 (MQ=255)
          tCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCAACGAGCAGGCtt                                                                                                                                           <  2:1126440/150‑1 (MQ=255)
            gcgcCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTCCGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGcgggg                                                                                                                                         >  2:1523413/1‑146 (MQ=255)
            gcgcCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAATAGCCTGCTCGTTGAGCAGGCtttt                                                                                                                                         >  2:2713636/1‑150 (MQ=255)
               ccgcATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCg                                                                                                                                       >  1:1379354/1‑150 (MQ=255)
               ccgcATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGa                                                                                                                                      >  2:1499616/1‑150 (MQ=255)
               ccgcATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGa                                                                                                                                      >  1:1814596/1‑150 (MQ=255)
                  cATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAAtt                                                                                                                                   <  1:2733359/150‑1 (MQ=255)
                  cATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAAtt                                                                                                                                   <  1:2541434/150‑1 (MQ=255)
                  cATCCGACAAGGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAGGC‑AAACAAGCCTGCTCGTTGAGCAGGCTTTTCGAAtt                                                                                                                                   <  2:568783/150‑1 (MQ=255)
                    tCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTg                                                                                                                                 <  2:1295384/150‑1 (MQ=255)
                    gccGACAATGGCCGGCACTTTTGGCTGTTGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTg                                                                                                                                 <  2:810279/149‑1 (MQ=255)
                          aaTGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGt                                                                                                                           <  1:443520/150‑1 (MQ=255)
                          aaTGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGt                                                                                                                           <  1:2105988/150‑1 (MQ=255)
                              gCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATa                                                                                                                       >  1:2614455/1‑150 (MQ=255)
                                   cacaTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGa                                                                                                                  <  2:2441976/150‑1 (MQ=255)
                                   cacaTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGa                                                                                                                  <  2:177473/150‑1 (MQ=255)
                                   cacaTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGa                                                                                                                  <  2:2262857/150‑1 (MQ=255)
                                      aTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTa                                                                                                               <  1:152309/150‑1 (MQ=255)
                                       tttGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTAc                                                                                                              >  1:1839097/1‑150 (MQ=255)
                                          gCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCg                                                                                                           <  1:1687624/150‑1 (MQ=255)
                                          gCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCg                                                                                                           <  1:2705345/150‑1 (MQ=255)
                                           ccTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTAGc                                                                                                          >  2:709249/1‑150 (MQ=255)
                                           ccTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGAGATAGAGGATAACTCGc                                                                                                          >  1:548918/1‑150 (MQ=255)
                                            cTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTagaacca                                                                                                         >  1:69291/1‑144 (MQ=255)
                                              gATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCAc                                                                                                       >  1:1366843/1‑150 (MQ=255)
                                                  cGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTGTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCg                                                                                                   >  1:2121231/1‑150 (MQ=255)
                                                    aCGCTGACGCGTCTTATCCGGCCTACGTCGAGTCGAACCTCCCGGCCCCCGGAGGTTCTCATCCTCCTGAGCGGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTg                                                                                                 <  2:1834945/150‑1 (MQ=255)
                                                    aCGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTg                                                                                                 <  2:632818/150‑1 (MQ=255)
                                                        tGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTc                                                                                             <  2:2337443/150‑1 (MQ=255)
                                                        tGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTc                                                                                             <  1:1193338/150‑1 (MQ=255)
                                                        tGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTc                                                                                             <  1:1302565/150‑1 (MQ=255)
                                                        tGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTc                                                                                             <  1:1324/150‑1 (MQ=255)
                                                          aCGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCAACTTCGTGGCTCGc                                                                                           >  2:1268895/1‑150 (MQ=255)
                                                                 ttATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCgg                                                                                    <  1:1935929/150‑1 (MQ=255)
                                                                 ttATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCgg                                                                                    <  1:512216/150‑1 (MQ=255)
                                                                 ttATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCgg                                                                                    <  1:630165/150‑1 (MQ=255)
                                                                 ttATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCATAATTAGCTGAAGATGA‑AAAAAAGCCTGCTAGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCgg                                                                                    <  1:1391515/150‑1 (MQ=18)
                                                                        gCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGTCCTGCGGGCCGTTg                                                                             <  1:1450766/150‑1 (MQ=255)
                                                                        gCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTg                                                                             <  2:2119706/150‑1 (MQ=255)
                                                                        gCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTg                                                                             <  2:2119721/150‑1 (MQ=255)
                                                                        gCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCATTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTg                                                                             <  2:2214464/150‑1 (MQ=255)
                                                                         ccTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCAACTTCCTGGCTAGCCCTGCGGGCCTTTGc                                                                            >  2:1159939/1‑150 (MQ=255)
                                                                          cTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCACATGC‑‑AAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTg                                                                          <  1:817992/150‑1 (MQ=255)
                                                                           tACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCACATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTg                                                                          <  1:418199/150‑1 (MQ=255)
                                                                             cGTCGAGTCGAACGTCCCTGCCCCCGGAGGTTCTCATCCTCTTGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGttgcttgc                                                                        <  2:1767114/150‑1 (MQ=255)
                                                                             cGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGgc                                                                        <  1:366770/150‑1 (MQ=255)
                                                                             cGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGgc                                                                        <  1:1643400/150‑1 (MQ=255)
                                                                                 gagtcgaACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCaa                                                                    <  1:891961/150‑1 (MQ=255)
                                                                                  agtcgaACCTCCCTGCCCCCGGAGGTTCTCCTCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGCTTACTCGCCACGTCGTGGCTCGCCCTGCGGGCCGGAGCTGGCGCacc                                                                   >  2:1982189/1‑148 (MQ=17)
                                                                                   gtcgaACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACg                                                                  <  1:1099625/150‑1 (MQ=255)
                                                                                   gtcgaACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACg                                                                  <  1:2466716/150‑1 (MQ=255)
                                                                                   gtcgaACCTCCCTGCCCCCGGAGGGTCTCATCCTCATGAGCGGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGACTTTGGTCGGTGATAGAGGATTACTCGCCACGTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACg                                                                  <  2:2167171/150‑1 (MQ=255)
                                                                                     cgaACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGtt                                                                <  1:2082694/150‑1 (MQ=255)
                                                                                       aaCCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTtct                                                              >  1:1541008/1‑150 (MQ=255)
                                                                                         ccTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTtctct                                                            <  1:1145357/150‑1 (MQ=255)
                                                                                          cTCCCTGCCCCCGGCGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTtctctc                                                           <  2:2599/150‑1 (MQ=255)
                                                                                           tCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGGTGCTGGCGCAACGGTCTCtcg                                                          >  1:286925/1‑150 (MQ=255)
                                                                                             ccTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTGTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCACCGGTcgctcgct                                                        >  1:2404990/1‑150 (MQ=255)
                                                                                                gCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCtcgcttcg                                                     <  1:2356573/150‑1 (MQ=255)
                                                                                                  ccccGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCtcgcttcgct                                                   <  1:1036076/150‑1 (MQ=255)
                                                                                                    ccGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAATCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCg                                                 >  1:913114/1‑150 (MQ=255)
                                                                                                        aGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGTCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTc                                             <  1:1637709/150‑1 (MQ=255)
                                                                                                        aGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTc                                             <  1:824926/150‑1 (MQ=255)
                                                                                                          gTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGa                                           >  2:2318010/1‑150 (MQ=255)
                                                                                                           ttCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGaa                                          >  1:998818/1‑150 (MQ=255)
                                                                                                           ttCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGaa                                          <  2:25207/150‑1 (MQ=255)
                                                                                                            tctcATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAAc                                         >  1:1883128/1‑150 (MQ=255)
                                                                                                            tctcATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCCTGGCACGCCCTGCGGGCCGTAGCTGGCGCACAGTTCACACGCATCGATCGAGTCGAAc                                         >  1:1168947/1‑150 (MQ=18)
                                                                                                                aTCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTcc                                     <  2:1821088/150‑1 (MQ=255)
                                                                                                                aTCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTcc                                     >  2:392402/1‑150 (MQ=255)
                                                                                                                aTCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCGCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTcc                                     >  2:1579829/1‑150 (MQ=255)
                                                                                                                        gAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCTCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGccccc                             >  2:1313370/1‑150 (MQ=255)
                                                                                                                        gAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGccccc                             <  2:1609109/150‑1 (MQ=255)
                                                                                                                        gAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGccccc                             >  2:1921092/1‑150 (MQ=255)
                                                                                                                         aGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCg                            >  2:1480144/1‑150 (MQ=255)
                                                                                                                           cTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGACCCTCCCTGCCCCCGGa                          >  1:2167248/1‑150 (MQ=255)
                                                                                                                           cTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGa                          >  1:1754943/1‑150 (MQ=255)
                                                                                                                           cTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGGTCTCTCGCTTCGCTCGAGTCGAACCGCCCTGCCCCCGGa                          >  1:474821/1‑150 (MQ=255)
                                                                                                                            tGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAg                         >  2:434764/1‑150 (MQ=255)
                                                                                                                               aGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGtt                      <  2:2370978/150‑1 (MQ=255)
                                                                                                                               aGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGtt                      <  2:307841/150‑1 (MQ=255)
                                                                                                                                gcggC‑AACACAGCCTGCGCGGGGCGCAGGCTTTTCGAATTTGGTCGGTGCTAGAGGATGACTCGCCACGTCGGGGCGCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCGGCCCCCGGAGGTtc                     <  2:1321801/147‑1 (MQ=255)
                                                                                                                                      aaaaaaGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATc                >  2:1954103/1‑150 (MQ=255)
                                                                                                                                       aaaaaGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACTTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATcc               >  2:1358694/1‑150 (MQ=255)
                                                                                                                                       aaaaaGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGAGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATcc               >  1:2436687/1‑150 (MQ=255)
                                                                                                                                         aaaGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTc             >  1:1406082/1‑150 (MQ=255)
                                                                                                                                         aaaGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTc             >  1:2212922/1‑150 (MQ=255)
                                                                                                                                          aaGCCTGCTCAACGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCAGTTCGTGGCTCGCCCTGCGGGCCGTTGGTGGCGCAACGTTCTCTCGCTTCGCTCGTGTCGAACCTCCCTG‑CCCCGGAGGTTCTCATGCTCag           >  2:708150/1‑149 (MQ=38)
                                                                                                                                          aaGCCTGCTCAACGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCCTCCTCa            >  1:2152942/1‑150 (MQ=255)
                                                                                                                                          aaGCCTGCTCAACGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCa            >  2:1426359/1‑150 (MQ=255)
                                                                                                                                          aaGCCTGCTCAACGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCGCATCCTCa            >  2:620346/1‑150 (MQ=255)
                                                                                                                                             gcagCTCTTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGACACTTCGTGGCTCGCCCTGCGGGCCGTTGCATGCGCATCGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGa         <  2:204102/147‑1 (MQ=38)
                                                                                                                                              cTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAg        <  2:1544446/150‑1 (MQ=255)
                                                                                                                                                gCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCt      >  1:446579/1‑150 (MQ=255)
                                                                                                                                                gCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCAGCCCCCGGAGGTTCTCATCCTCATGAGCt      >  1:381344/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGC‑AAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAG  >  NC_000913/3708280‑3708576

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: ‑