Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,708,429 T→C 5.5% intergenic (‑724/+187) 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,4290TC 5.5% 719.0 / 22.1 220intergenic (‑724/+187)dppA/proKdipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor/tRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base T (94/112);  new base C (7/5);  total (102/118)
Fisher's exact test for biased strand distribution p-value = 5.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

Only 100 of 224 total aligned reads displayed.

AATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAT  >  NC_000913/3708280‑3708578
                                                                                                                                                     |                                                                                                                                                     
aaTAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGtt                                                                                                                                                       <  1:1537174/150‑1 (MQ=255)
aaTAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGtt                                                                                                                                                       <  1:1718975/150‑1 (MQ=255)
aaTAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGtt                                                                                                                                                       <  2:1445064/150‑1 (MQ=255)
  tAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCACAAAAGCCTGCTCGTTGa                                                                                                                                                     >  2:2882093/1‑150 (MQ=255)
    aGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGc                                                                                                                                                   <  1:2375671/150‑1 (MQ=255)
       cGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAgg                                                                                                                                                <  1:2817935/150‑1 (MQ=255)
          tCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCtt                                                                                                                                             >  2:2426706/1‑150 (MQ=255)
          tCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCtt                                                                                                                                             >  2:1359223/1‑150 (MQ=255)
            gcgcCGCATCCGACACTGGCCGGCTCCGTGGGCGGTCGCGACGCTGACGCGACTTATCAGGCCTAGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCtttt                                                                                                                                           <  2:1049953/150‑1 (MQ=17)
             cgccgcATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTc                                                                                                                                          <  1:3089448/150‑1 (MQ=255)
             cgccgcATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTc                                                                                                                                          <  2:3354813/150‑1 (MQ=255)
             cgccgcATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCAACGAGCAGGCttttt                                                                                                                                          <  2:3870183/150‑2 (MQ=255)
             cgccgcATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCAACGAGCAGGCttttt                                                                                                                                          <  1:4226112/150‑2 (MQ=255)
               tcgcATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGAGGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGa                                                                                                                                        >  2:3937290/2‑150 (MQ=255)
               ccgcATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGa                                                                                                                                        >  2:3023590/1‑150 (MQ=255)
                 gcATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAAt                                                                                                                                      <  1:3645952/150‑1 (MQ=255)
                 gcATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAAt                                                                                                                                      <  1:3475803/150‑1 (MQ=255)
                  cATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAAtt                                                                                                                                     >  2:3503827/1‑150 (MQ=255)
                    tCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTg                                                                                                                                   >  1:1783218/1‑150 (MQ=255)
                    tCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTg                                                                                                                                   <  1:2098899/150‑1 (MQ=255)
                       gACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTc                                                                                                                                >  2:2376932/1‑150 (MQ=255)
                         cactGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCgg                                                                                                                              <  2:850914/147‑1 (MQ=255)
                         cAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCgg                                                                                                                              <  2:4080243/150‑1 (MQ=255)
                         cAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCgg                                                                                                                              <  2:115990/150‑1 (MQ=255)
                         cAATGGCCGGCACATTTGCCTGATGCGACGCTGACGAGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCgg                                                                                                                              <  2:4037286/150‑1 (MQ=255)
                            tGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCAACCTCCTGAGCTGCAGATGCAAAAAAGCCTGCGCGTTGAGCAGGCTTTTCGAAGTTGGTCGGTGa                                                                                                                           >  1:653772/1‑150 (MQ=255)
                              gCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATa                                                                                                                         <  1:4067104/150‑1 (MQ=255)
                               ccGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATag                                                                                                                        >  1:129558/1‑150 (MQ=255)
                                cGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATaga                                                                                                                       <  2:1632692/150‑1 (MQ=255)
                                 ggCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGACTTTGGTCGGTGATagag                                                                                                                      >  1:3944031/1‑150 (MQ=255)
                                  gCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAgg                                                                                                                     >  2:4216939/1‑150 (MQ=255)
                                  gCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCAGCTCCTTGGTCAGGCATTTCGAATTTGGTCGGGGGTTGAgg                                                                                                                     >  2:1352792/1‑150 (MQ=18)
                                   cacaTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGa                                                                                                                    <  2:3414559/150‑1 (MQ=255)
                                           ccTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGc                                                                                                            <  1:2748138/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)
                                                     cGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTgg                                                                                                  <  1:2781236/150‑1 (MQ=255)
                                                      gCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAATATTGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGTTAGTGGATTAGTAGCTACTTCGTggg                                                                                                 >  2:2233338/1‑149 (MQ=17)
                                                      gCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGc                                                                                                 <  1:3971324/150‑1 (MQ=255)
                                                          aCGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGc                                                                                             <  2:2046076/150‑1 (MQ=255)
                                                            gcgTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTAGTTGAGCAGGCTGTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGccc                                                                                           >  1:452784/1‑150 (MQ=255)
                                                            gcgTCTTATCAGGCCAACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGccc                                                                                           >  1:1474534/1‑150 (MQ=255)
                                                                   aTCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGc                                                                                    <  1:2944207/150‑1 (MQ=255)
                                                                   aTCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGc                                                                                    <  1:1564107/150‑1 (MQ=255)
                                                                   aTCAGGACTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCATTCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGc                                                                                    <  2:608449/150‑1 (MQ=255)
                                                                     cAGGCCTACGTCGAGGCGAACCTCCCTGCCCCCGGAGGTTCCCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCg                                                                                  <  1:823881/150‑1 (MQ=255)
                                                                        gCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTg                                                                               <  1:294156/150‑1 (MQ=255)
                                                                        gCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTg                                                                               >  2:2250832/1‑150 (MQ=255)
                                                                         ccTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGc                                                                              <  1:3420438/150‑1 (MQ=255)
                                                                          cTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATTCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCt                                                                             <  2:3205081/150‑1 (MQ=255)
                                                                          cTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCt                                                                             <  2:4215287/150‑1 (MQ=255)
                                                                          cTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCt                                                                             <  1:2499722/150‑1 (MQ=255)
                                                                          cTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCt                                                                             <  1:3496103/150‑1 (MQ=255)
                                                                          cTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGACGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCt                                                                             <  1:3434525/150‑1 (MQ=255)
                                                                            aCGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCa                                                                                                                                                  <  1:1181416/79‑1 (MQ=34)
                                                                             cGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGgc                                                                          <  1:2858563/150‑1 (MQ=255)
                                                                             cGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGgc                                                                          <  1:2908151/150‑1 (MQ=255)
                                                                                  agtcgaACCTCCCTGCCCCCGGAGGGTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAAc                                                                     <  1:1936373/150‑1 (MQ=255)
                                                                                    tcgaACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGGTGCTGGCGCAACGt                                                                   >  1:2069719/1‑150 (MQ=255)
                                                                                        aCCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTtctc                                                               <  1:742887/150‑1 (MQ=255)
                                                                                            cccTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGTTGCAAAATAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCCCTTCGTGGCTCGCCCGGCGGGCCGGTGCTGGCGCAACGTGCTCtcgc                                                           >  2:2662839/1‑150 (MQ=18)
                                                                                            cccTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTAGCTCGCCCACCGTACTCtcgc                                                           >  2:2622714/1‑150 (MQ=18)
                                                                                               tGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAATGCCTGCTCGGTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACGTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCTACGGTCTCtcgcttg                                                        >  2:2481740/1‑149 (MQ=21)
                                                                                                 cccccGGAGGTTCTCATCCTCATGATCTGCCGGTGCAAAAAAGACTTCTAGTTGAGCAGGCATTTCGAATTTGGTCGGTGCTAGAGGATTACTCGCCACATCGTGGCTCGCCCTGCGGGCCGGTGCTGGCGCGAAGCTCTCtcgcttcgc                                                      >  1:1605750/1‑150 (MQ=18)
                                                                                                     cGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGa                                                  <  1:971072/150‑1 (MQ=255)
                                                                                                           ttCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGaa                                            >  1:986519/1‑150 (MQ=255)
                                                                                                           ttCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGaa                                            >  1:4116407/1‑150 (MQ=255)
                                                                                                           ttCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGaa                                            <  1:2426706/150‑1 (MQ=255)
                                                                                                              tcATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCt                                         >  1:3922669/1‑150 (MQ=255)
                                                                                                               cATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTc                                        >  1:1746387/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)
                                                                                                                    tCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGc                                   >  1:2686677/1‑150 (MQ=255)
                                                                                                                      aTGAGCTGCAGATGCTCACAAGCCGGCTCGTTGAGCTGGCTTTTCGACTTTGGGCGGTGAGGGAGGAGTACGGGGCTCTGGGGGGCTCGCGGGGCGGGCGGGTGGGGGCGGATCGTTCTCTGGCGTCGGTCGAGTCGAACCTCCCTGccc                                 <  2:653772/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:3125090/1‑150 (MQ=255)
                                                                                                                                     aaaaaaGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATc                  >  1:184409/1‑150 (MQ=255)
                                                                                                                                      aaaaaGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATcc                 >  2:213648/1‑150 (MQ=255)
                                                                                                                                       aaaaGCCTGCTCAACGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGTGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCt                >  1:4146641/1‑150 (MQ=255)
                                                                                                                                        aaaGCCTGCTCGTTGATCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTc               <  1:4051381/150‑1 (MQ=255)
                                                                                                                                        aaaGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTc               <  1:3136935/150‑1 (MQ=255)
                                                                                                                                        aaaGCCTGCTCAACGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATccac               >  1:3303649/1‑148 (MQ=255)
                                                                                                                                         aaGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCa              <  1:3000191/150‑1 (MQ=255)
                                                                                                                                         aaGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTGCTCATCCTCa              >  1:3515345/1‑150 (MQ=255)
                                                                                                                                         aaGCCTGCTCAACGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCa              >  1:2585855/1‑150 (MQ=255)
                                                                                                                                           gCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATg            <  2:1133210/150‑1 (MQ=255)
                                                                                                                                             cTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCGTGAg          >  2:3700051/1‑150 (MQ=255)
                                                                                                                                               gCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCGGCCCCCGGAGGTTCTCATCCTCATGAGCt        >  2:1842968/1‑150 (MQ=255)
                                                                                                                                               gCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCGCCCAGCCCCCCGAGGTTCTCATCCTCAAGAGCt        >  1:4013813/1‑150 (MQ=255)
                                                                                                                                               gCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTAGCTGGCGCAACGGTCTCTCGCTTCGCTCGAGTCGAACCTCCCAGCCCCCGGAGGTTCTCATCCGCATGAGCt        >  2:148301/1‑150 (MQ=255)
                                                                                                                                                cTCGTAGAGCAGGCTTTTCGAAATTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTg       <  1:2447443/150‑1 (MQ=255)
                                                                                                                                                 tCGTTGAGCAGTCTTTTCGCAGGTGGTCGGCGCTAGAGGATTACTCGCCACTGCGTGGCTCTCCCTGCGGGCCGTTGCTGGCGCAACGTTATCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGc      <  2:1398839/150‑1 (MQ=37)
                                                                                                                                                     tGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAt  <  2:2175499/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAT  >  NC_000913/3708280‑3708578

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: