Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,908,514 A→T 9.8% intergenic (‑148/+35) pstB ← / ← pstA phosphate ABC transporter ATPase/phosphate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,908,5140AT 9.8% 1311.4 / 22.7 454intergenic (‑148/+35)pstB/pstAphosphate ABC transporter ATPase/phosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (269/118);  new base T (2/57);  total (271/183)
Fisher's exact test for biased strand distribution p-value = 8.28e-24
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 470 total aligned reads displayed.

ATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAACGGG  >  NC_000913/3908365‑3908663
                                                                                                                                                     |                                                                                                                                                     
 ttGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                                      <  1:3101495/150‑1 (MQ=255)
   gCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAggg                                                                                                                                                    <  2:3898096/150‑1 (MQ=255)
    cAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGt                                                                                                                                                   >  2:2265371/1‑150 (MQ=255)
    cAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGt                                                                                                                                                   >  2:614783/1‑150 (MQ=255)
     aaTCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGtt                                                                                                                                                  <  1:5713789/150‑1 (MQ=255)
     aaTCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGtt                                                                                                                                                  >  2:6204882/1‑150 (MQ=255)
      aTCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTa                                                                                                                                                 >  2:1547601/1‑150 (MQ=255)
       tCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAg                                                                                                                                                <  1:5660433/150‑1 (MQ=255)
       tCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTACCCTCTCCCCAGAGGGGCGAggggacc                                                                                                                                                <  1:1620790/150‑5 (MQ=39)
        cTTGTCCTCCTCGGGCCTGATGCACTTCGCTAATCCGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAgg                                                                                                                                               <  2:3377934/150‑1 (MQ=255)
         ttttCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAggg                                                                                                                                              >  2:1637279/1‑150 (MQ=255)
             ccTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAg                                                                                                                                          >  1:2079605/1‑150 (MQ=255)
             ccTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGATAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAg                                                                                                                                          <  2:787178/150‑1 (MQ=255)
              cTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgg                                                                                                                                         <  1:5885854/150‑1 (MQ=255)
              cTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgg                                                                                                                                         <  2:1366257/150‑1 (MQ=255)
                       ccTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCg                                                                                                                                >  2:4781337/1‑150 (MQ=255)
                       ccTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCg                                                                                                                                >  1:2527680/1‑150 (MQ=255)
                          gATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCgtgt                                                                                                                             >  2:2802925/1‑150 (MQ=255)
                            tGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCCCGCCGCATCCGGCAAAGGGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCGCGCCCCTCTGGGGAGCGGGTTCGGGTGCGGGGAAAACCGTGTCa                                                                                                                           <  2:2938659/150‑1 (MQ=255)
                              cACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGc                                                                                                                         >  1:4960828/1‑150 (MQ=255)
                               aCTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCATAATGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGct                                                                                                                        >  2:2699250/1‑149 (MQ=255)
                               aCTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCa                                                                                                                        >  1:2297999/1‑150 (MQ=255)
                                 ttCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAat                                                                                                                      >  2:2656625/1‑150 (MQ=255)
                                 ttCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAat                                                                                                                      >  2:1161931/1‑150 (MQ=255)
                                 ttCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAat                                                                                                                      >  2:3602915/1‑150 (MQ=255)
                                 ttCGATAATCTGGATTACCGTTGGCCGGATCTGGCGTTCACGCCGCATCCGGCAATGTGTTAACCGTGTACAGAGTCAAAAGGGGGCTGGGTGGGTCCCCTCGCCCCTGTGGGGTGTGGGTTTGGTTGAGGGGAAAACCGTGTCAGCAat                                                                                                                      <  1:2221699/150‑1 (MQ=255)
                                   cGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAatat                                                                                                                    >  1:4182776/1‑150 (MQ=255)
                                      tAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCaa                                                                                                                 >  1:1814113/1‑150 (MQ=255)
                                      tAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCaa                                                                                                                 >  1:2954220/1‑150 (MQ=255)
                                      tAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCaa                                                                                                                 >  1:6510109/1‑150 (MQ=255)
                                        aTCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAAcc                                                                                                               >  1:5141861/1‑150 (MQ=255)
                                         tCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCg                                                                                                              >  1:6192691/1‑150 (MQ=255)
                                         tCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCg                                                                                                              >  1:3887676/1‑150 (MQ=255)
                                             gCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGtt                                                                                                          >  2:4186135/1‑150 (MQ=255)
                                             gCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGtt                                                                                                          >  2:5649666/1‑150 (MQ=255)
                                             gCTGACCGTCGGCCGGATCTGGCGGTCACGCCGCATCCGGCAACGTGTTAACCGTGTACAGAGTCACAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGTGTGGGTTTGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGtt                                                                                                          <  2:2007786/150‑1 (MQ=38)
                                                  ccGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGTGTGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttc                                                                                                     <  1:1449467/150‑1 (MQ=255)
                                                    gTAGGCCGTACCAGGCGTTCACGCCGCCTCCGGCAATGTGTTAACCGTGTACCGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttctt                                                                                                   <  2:2430859/150‑1 (MQ=38)
                                                    gTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttctt                                                                                                   <  2:4960828/150‑1 (MQ=255)
                                                    gTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttctt                                                                                                   >  1:3179950/1‑150 (MQ=255)
                                                     tAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttcttc                                                                                                  >  2:1471139/1‑150 (MQ=255)
                                                     tAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttcttc                                                                                                  >  2:379619/1‑150 (MQ=255)
                                                      aGGCCGGATCTGGCGTTCACGCCGCATCGGGCAAAGTTTTAACCGTTTACAGTGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCg                                                                                                 <  2:3172389/150‑1 (MQ=38)
                                                       ggCAGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAAACGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCACGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGc                                                                                                >  2:5884294/1‑150 (MQ=255)
                                                        gCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCa                                                                                               >  1:1551329/1‑150 (MQ=255)
                                                         ccGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCaa                                                                                              >  1:1927742/1‑150 (MQ=255)
                                                           ggATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGCAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCaaaa                                                                                            >  2:1834614/1‑150 (MQ=255)
                                                           ggATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCaaaa                                                                                            >  1:3592192/1‑150 (MQ=255)
                                                            gATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCaaaaa                                                                                           >  2:4356443/1‑150 (MQ=255)
                                                                  gCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgc                                                                                     >  2:1454781/1‑150 (MQ=255)
                                                                  gCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGACAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGATAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgc                                                                                     >  2:3730463/1‑150 (MQ=255)
                                                                     ttCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgcgcg                                                                                  >  1:1193197/1‑150 (MQ=255)
                                                                     ttCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgcgcg                                                                                  >  1:6589984/1‑150 (MQ=255)
                                                                        aCGCCGCATCCGGCACAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCACAAAAACAACGCGCGCCa                                                                               >  1:5263768/1‑150 (MQ=255)
                                                                          gccgcATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGa                                                                             >  1:3278370/1‑150 (MQ=255)
                                                                          gccgcATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGa                                                                             >  2:2577687/1‑150 (MQ=255)
                                                                            cgcATCCGGCAATTTGTTAACCGTGAACAGAGACAAAATGGCGCTCGGTCGGCCCCCTCGCCCCTCTGGGGAGAGGGTTTGGGTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAAt                                                                           <  1:3705355/150‑1 (MQ=38)
                                                                            cgcATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCATCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAAt                                                                           <  2:3971613/150‑1 (MQ=255)
                                                                             gcATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATg                                                                          >  2:5822076/1‑150 (MQ=255)
                                                                                tCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAATGGTGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTc                                                                       <  2:5349417/150‑1 (MQ=255)
                                                                                  cGGAAATGTGTTAACCGTGTACAGAGTCAAAATGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGTTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcag                                                                     <  1:4877174/150‑1 (MQ=255)
                                                                                   ggCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagc                                                                    >  1:5556293/1‑150 (MQ=255)
                                                                                    gCAAAGTGATAACCGTGTACAGAGTAAAAAAGGCGCTCGGTCGGTCCCCTCGCACCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACAGTCTCAGCAATATCAACCTTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagaa                                                                   >  2:3128044/1‑148 (MQ=38)
                                                                                     cAATGTGTTACCCGTGTACAGAGTCAAAACGGCGCTGGGTGGGGCCCCTCGCCCCGCTGGGGTGTGGTTTAGGGTGTGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagcag                                                                  <  1:4396751/150‑1 (MQ=255)
                                                                                      aaaGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagcagc                                                                 >  2:162344/1‑150 (MQ=255)
                                                                                      aaaGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagcagc                                                                 >  2:4443503/1‑150 (MQ=255)
                                                                                      aaaGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTAAGAAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCTCGCAAGAATGTTcagaaga                                                                 >  2:914120/1‑149 (MQ=38)
                                                                                            ttAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACg                                                           <  2:4082231/150‑1 (MQ=255)
                                                                                              aaCCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGca                                                         >  2:3053125/1‑150 (MQ=255)
                                                                                              aaCCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGca                                                         >  2:757854/1‑150 (MQ=255)
                                                                                                  gtgtACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAgg                                                     >  2:5166986/1‑150 (MQ=255)
                                                                                                  gtgtACAGAGAAAAAACGGGGCTCGGTCGGCCCCCTCGCCCCTGTGGGTTGTGGGTTTGGGTGTGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAgg                                                     <  2:1120957/150‑1 (MQ=255)
                                                                                                       cAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAAt                                                <  2:5508862/150‑1 (MQ=255)
                                                                                                        agagTCAAACAGGCGCTCGGTCGGTCCCCTCGCCCCGCGGGGGAGAGGTTTTGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATg                                               <  1:3431262/150‑1 (MQ=255)
                                                                                                        agagTCAAAAAGGCGCTCGGTCGGTCCCCGCGCCCCTCTGGGGAGAGGGTTTGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATg                                               <  1:4099074/150‑1 (MQ=255)
                                                                                                          agTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGAt                                             >  1:4777778/1‑150 (MQ=255)
                                                                                                            taaaaaaGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCa                                           <  1:2616143/148‑1 (MQ=255)
                                                                                                            tCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTCATGATCa                                           >  1:4519172/1‑150 (MQ=255)
                                                                                                              aaaaaGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAAt                                         >  1:6029996/1‑150 (MQ=255)
                                                                                                              aaaaaGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAAt                                         >  2:5952134/1‑150 (MQ=255)
                                                                                                                aatggCTCTGGGTCGGACCCCTCGCCCGTCTGGGTTGTGGTTTAGTGTGAGGGCAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATAc                                       <  1:4691271/147‑1 (MQ=255)
                                                                                                                aaaGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATAc                                       >  2:5263418/1‑150 (MQ=255)
                                                                                                                    gcgcTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCg                                   >  2:6463518/1‑150 (MQ=255)
                                                                                                                       cTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTCGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGcc                                <  2:6392605/150‑1 (MQ=255)
                                                                                                                       cTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGcc                                <  1:3023470/150‑1 (MQ=255)
                                                                                                                        tcggtcggTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGccc                               >  2:5124494/1‑150 (MQ=255)
                                                                                                                          ggtcggtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAg                             <  2:2079605/150‑1 (MQ=255)
                                                                                                                               gtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTTAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCaa                        <  2:4542055/150‑1 (MQ=255)
                                                                                                                                    cTCGCCCCTCTGGGGAGAGGGTTCGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgct                   <  1:744749/150‑1 (MQ=255)
                                                                                                                                    cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgct                   <  1:5421233/150‑1 (MQ=255)
                                                                                                                                     tCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgctg                  >  2:3884611/1‑150 (MQ=255)
                                                                                                                                     tCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgctg                  >  2:1224620/1‑150 (MQ=255)
                                                                                                                                     ctaacccTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgctg                  >  2:23319/5‑150 (MQ=255)
                                                                                                                                      cGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgctgc                 <  1:5937829/150‑1 (MQ=255)
                                                                                                                                         cccTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCAt              >  1:1865097/1‑150 (MQ=255)
                                                                                                                                         cccTCTGGGGAGAGGGATAGGGTGAGGGGAAAACCGTTTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTCCCAt              >  1:3973787/1‑150 (MQ=255)
                                                                                                                                          ccTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCAtt             <  2:1004554/150‑1 (MQ=255)
                                                                                                                                                  gagaGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAAc     >  2:1084843/1‑150 (MQ=255)
                                                                                                                                                  gagaGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAAc     <  2:6031544/150‑1 (MQ=255)
                                                                                                                                                  gagaGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAAc     >  1:1013355/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAACGGG  >  NC_000913/3908365‑3908663

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: