Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,708,428 T→A 5.3% intergenic (‑723/+188) 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,4280TA 5.3% 726.3 / 23.1 220intergenic (‑723/+188)dppA/proKdipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor/tRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base T (91/113);  new base A (7/5);  total (102/118)
Fisher's exact test for biased strand distribution p-value = 3.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.92e-01

Only 100 of 225 total aligned reads displayed.

GAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCA  >  NC_000913/3708279‑3708575
                                                                                                                                                     |                                                                                                                                                   
gAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGt                                                                                                                                                     >  2:2535182/1‑150 (MQ=255)
 aaTAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGtt                                                                                                                                                    <  2:1445064/150‑1 (MQ=255)
   tAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCACAAAAGCCTGCTCGTTGa                                                                                                                                                  >  2:2882093/1‑150 (MQ=255)
       gCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAg                                                                                                                                              >  1:4027250/1‑150 (MQ=255)
           tCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCtt                                                                                                                                          >  2:1359223/1‑150 (MQ=255)
           tCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCtt                                                                                                                                          >  2:2427334/1‑150 (MQ=255)
              cgccgcATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTc                                                                                                                                       <  2:3354813/150‑1 (MQ=255)
              cgccgcATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCAACGAGCAGGCttttt                                                                                                                                       <  2:3870183/150‑2 (MQ=255)
               gccgcATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGTGCTGCAGATGCAAAAAAGCCTGCTCAACGAGCAGGCtttttt                                                                                                                                      <  1:2829391/150‑3 (MQ=255)
                tcgcATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGAGGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGa                                                                                                                                     >  2:3937290/2‑150 (MQ=255)
                tcgcATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGa                                                                                                                                     >  2:1331834/2‑150 (MQ=255)
                ccgcATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGa                                                                                                                                     >  2:3023590/1‑150 (MQ=255)
                 cgcATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGaa                                                                                                                                    <  1:3018309/150‑1 (MQ=255)
                 cgcACCCGACACTGGCCGGCACCTGTGCCGGACGCGCCGCTGACGCGTCTTATCAGGCCTACGTCGAGGCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGaa                                                                                                                                    <  1:2391357/150‑1 (MQ=18)
                  gcATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAAt                                                                                                                                   <  1:3475803/150‑1 (MQ=255)
                  gcATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAAt                                                                                                                                   <  1:3645952/150‑1 (MQ=255)
                   cATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAAtt                                                                                                                                  >  2:3503827/1‑150 (MQ=255)
                        gACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTc                                                                                                                             >  2:2376932/1‑150 (MQ=255)
                          cactGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCgg                                                                                                                           <  2:850914/147‑1 (MQ=255)
                          cAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCgg                                                                                                                           <  2:115990/150‑1 (MQ=255)
                          cAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCgg                                                                                                                           <  2:4080243/150‑1 (MQ=255)
                          cAATGGCCGGCACATTTGCCTGATGCGACGCTGACGAGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCgg                                                                                                                           <  2:4037286/150‑1 (MQ=255)
                             tGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGa                                                                                                                        <  2:1472245/150‑1 (MQ=255)
                             tGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCAACCTCCTGAGCTGCAGATGCAAAAAAGCCTGCGCGTTGAGCAGGCTTTTCGAAGTTGGTCGGTGa                                                                                                                        >  1:653772/1‑150 (MQ=255)
                                ccGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATag                                                                                                                     >  1:129558/1‑150 (MQ=255)
                                 cGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATaga                                                                                                                    <  2:1632692/150‑1 (MQ=255)
                                   gCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAgg                                                                                                                  >  2:4216939/1‑150 (MQ=255)
                                   gCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCAGGAGCTGCAGATGCACAAAAGCCTGCTCGGTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAgg                                                                                                                  >  1:1554301/1‑150 (MQ=255)
                                    cacaTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGa                                                                                                                 <  2:3414559/150‑1 (MQ=255)
                                            ccTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGc                                                                                                         <  1:1219297/150‑1 (MQ=255)
                                            ccTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGc                                                                                                         <  1:3746223/150‑1 (MQ=255)
                                            ccTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCAGCTCGTTGATCAGGCTTTTCGAAATTGTTAGGGGATAGAGGATTACTCGc                                                                                                         >  1:545372/1‑150 (MQ=18)
                                                 tGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACtt                                                                                                    >  1:230114/1‑150 (MQ=255)
                                                    gACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGt                                                                                                 <  1:519148/150‑1 (MQ=255)
                                                     aCGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTg                                                                                                <  1:3274632/150‑1 (MQ=255)
                                                     aCGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTg                                                                                                <  1:2081640/150‑1 (MQ=255)
                                                      cGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTgg                                                                                               <  1:2781236/150‑1 (MQ=255)
                                                      cGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTgg                                                                                               <  1:4242787/150‑1 (MQ=255)
                                                      cGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTgg                                                                                               <  1:3906562/150‑1 (MQ=255)
                                                       gCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGc                                                                                              <  1:3971324/150‑1 (MQ=255)
                                                       gCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGc                                                                                              >  2:2233223/1‑150 (MQ=255)
                                                           aCGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGc                                                                                          <  2:2046076/150‑1 (MQ=255)
                                                            cgcgTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGcc                                                                                         <  1:2808232/150‑1 (MQ=255)
                                                             gcgTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGccc                                                                                        >  1:1474451/1‑150 (MQ=255)
                                                                   tATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCggg                                                                                  <  1:1751140/150‑1 (MQ=255)
                                                                   tATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACGTCGAGGCTCGCCCTGCggg                                                                                  >  1:1734235/1‑150 (MQ=255)
                                                                    aTCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGc                                                                                 <  1:2944207/150‑1 (MQ=255)
                                                                    aTCAGGACTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCATTCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGc                                                                                 <  2:608449/150‑1 (MQ=255)
                                                                      cAGGCCTACGTCGAGGCGAACCTCCCTGCCCCCGGAGGTTCCCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCg                                                                               <  1:823881/150‑1 (MQ=255)
                                                                       aGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGt                                                                              <  1:236932/150‑1 (MQ=255)
                                                                       aGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGt                                                                              >  1:3498678/1‑150 (MQ=255)
                                                                         gCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTg                                                                            <  1:1661180/150‑1 (MQ=255)
                                                                         gCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTg                                                                            <  1:294156/150‑1 (MQ=255)
                                                                           cTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATTCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCt                                                                          <  2:3205081/150‑1 (MQ=255)
                                                                           cTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCt                                                                          <  1:3496028/150‑1 (MQ=255)
                                                                           cTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCt                                                                          <  2:4215287/150‑1 (MQ=255)
                                                                           cTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCt                                                                          <  1:2499722/150‑1 (MQ=255)
                                                                           cTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGACGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCt                                                                          <  1:3434525/150‑1 (MQ=255)
                                                                             aCGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCa                                                                                                                                               <  1:1181416/79‑1 (MQ=34)
                                                                              cGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGgc                                                                       <  1:1481839/150‑1 (MQ=255)
                                                                                          ccTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTtctct                                                           <  1:2904017/150‑1 (MQ=255)
                                                                                             cccTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGTTGCAAAATAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCCCTTCGTGGCTCGCCCGGCGGGCCGGTGCTGGCGCAACGTGCTCtcgc                                                        >  2:2662839/1‑150 (MQ=18)
                                                                                             cccTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCtcgc                                                        >  1:3399277/1‑150 (MQ=255)
                                                                                             cccTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTAGCTCGCCCACCGTACTCtcgc                                                        >  2:2622714/1‑150 (MQ=18)
                                                                                                tGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCtcgcttc                                                     <  2:2233674/150‑1 (MQ=255)
                                                                                                  cccccGGAGGTTCTCATCCTCATGATCTGCCGGTGCAAAAAAGACTTCTAGTTGAGCAGGCATTTCGAATTTGGTCGGTGCTAGAGGATTACTCGCCACATCGTGGCTCGCCCTGCGGGCCGGTGCTGGCGCGAAGCTCTCtcgcttcgc                                                   >  1:1605750/1‑150 (MQ=18)
                                                                                                    cccGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTc                                                 <  1:2017097/150‑1 (MQ=255)
                                                                                                      cGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGa                                               <  1:971072/150‑1 (MQ=255)
                                                                                                            ttCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGac                                         >  1:1917389/1‑149 (MQ=255)
                                                                                                            ttCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGaa                                         <  1:2427334/150‑1 (MQ=255)
                                                                                                            ttCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGaa                                         >  1:4116407/1‑150 (MQ=255)
                                                                                                            ttCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGaa                                         >  1:986519/1‑150 (MQ=255)
                                                                                                               tcATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCt                                      >  2:208949/1‑150 (MQ=255)
                                                                                                                 aTCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTcc                                    <  1:1419372/150‑1 (MQ=255)
                                                                                                                  tCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTGTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTccc                                   >  2:1210548/1‑150 (MQ=255)
                                                                                                                     tCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCTCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGc                                >  2:2566910/1‑150 (MQ=255)
                                                                                                                       aTGAGCTGCAGATGCTCACAAGCCGGCTCGTTGAGCTGGCTTTTCGACTTTGGGCGGTGAGGGAGGAGTACGGGGCTCTGGGGGGCTCGCGGGGCGGGCGGGTGGGGGCGGATCGTTCTCTGGCGTCGGTCGAGTCGAACCTCCCTGccc                              <  2:653772/150‑1 (MQ=255)
                                                                                                                         gAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGccccc                            >  1:2667194/1‑150 (MQ=255)
                                                                                                                            cTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGTTTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGa                         <  1:3929591/150‑1 (MQ=255)
                                                                                                                            cTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGa                         >  1:3124768/1‑150 (MQ=255)
                                                                                                                                 gATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTtc                    <  1:1033180/150‑1 (MQ=255)
                                                                                                                                   tGCAAAAAAGCCTGCTCTTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTtctc                  <  2:2975494/150‑1 (MQ=255)
                                                                                                                                   tGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTtctc                  <  2:773804/150‑1 (MQ=255)
                                                                                                                                    gCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCa                 >  2:1036380/1‑150 (MQ=255)
                                                                                                                                    gCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCa                 >  2:3125090/1‑150 (MQ=255)
                                                                                                                                      aaaaaaGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATc               >  1:184409/1‑150 (MQ=255)
                                                                                                                                       aaaaaGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGAATACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATcc              >  1:274008/1‑150 (MQ=255)
                                                                                                                                        aaaaGCCTGCTCAACGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGTGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCt             >  1:4146641/1‑150 (MQ=255)
                                                                                                                                        aaaaGCCTGCTCAACGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACATCGTGGCTCGCCCTGCGGGCCGGTGCTGGCGCCCCGTTCTCTCGCTTCGCGCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATccg             >  1:4146672/1‑149 (MQ=38)
                                                                                                                                         aaaGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTc            <  1:3136935/150‑1 (MQ=255)
                                                                                                                                          aaGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCa           <  1:3000188/150‑1 (MQ=255)
                                                                                                                                            gCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATg         <  2:1133210/150‑1 (MQ=255)
                                                                                                                                              cTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCGTGAg       >  2:3700051/1‑150 (MQ=255)
                                                                                                                                              cTGCTCGTTGAGCAGGCTTTTAGTATTTGGTCGGTGATAGAGGATTACTCGCCACGTCGTGGCTCGCGCTGCGGGCCGTTGCTGGCGCATCGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAg       <  1:4066136/150‑1 (MQ=38)
                                                                                                                                                gCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCGGCCCCCGGAGGTTCTCATCCTCATGAGCt     >  2:1842968/1‑150 (MQ=255)
                                                                                                                                                gCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCGCCCAGCCCCCCGAGGTTCTCATCCTCAAGAGCt     >  1:4013813/1‑150 (MQ=255)
                                                                                                                                                 cTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTg    >  1:2563133/1‑150 (MQ=255)
                                                                                                                                                 cTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCACCCTGCCCCCCGAGGTTCTCATCCTCAAGAGCTg    >  1:721192/1‑150 (MQ=255)
                                                                                                                                                  tCGTTGAGCAGTCTTTTCGCAGGTGGTCGGCGCTAGAGGATTACTCGCCACTGCGTGGCTCTCCCTGCGGGCCGTTGCTGGCGCAACGTTATCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGc   <  2:1398839/150‑1 (MQ=37)
                                                                                                                                                   cGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCa  >  2:2277243/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCA  >  NC_000913/3708279‑3708575

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: