Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,908,507 T→G 5.1% intergenic (‑141/+42) 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,5070TG 5.1% 503.4 / 3.4 173intergenic (‑141/+42)pstB/pstAphosphate ABC transporter ATPase/phosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (97/63);  new base G (3/9);  total (100/73)
Fisher's exact test for biased strand distribution p-value = 2.96e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 178 total aligned reads displayed.

ATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCC  >  NC_000913/3908359‑3908654
                                                                                                                                                    |                                                                                                                                                   
agaCTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTg                                                                                                                                                    <  2:978991/148‑1 (MQ=255)
aTACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTg                                                                                                                                                    <  2:366787/150‑1 (MQ=255)
  aCTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTggg                                                                                                                                                  <  1:310185/150‑1 (MQ=255)
  aCTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTggg                                                                                                                                                  <  2:151893/150‑1 (MQ=255)
  aCTCATTGCAATCTTGTCCTCATCGGGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTggg                                                                                                                                                  <  1:426678/150‑1 (MQ=255)
    tCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGga                                                                                                                                                >  1:691526/1‑150 (MQ=255)
    tCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGga                                                                                                                                                >  1:1105124/1‑150 (MQ=255)
    tCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGga                                                                                                                                                >  1:1107773/1‑150 (MQ=255)
    tCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGga                                                                                                                                                >  1:163839/1‑150 (MQ=255)
    tCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGga                                                                                                                                                <  2:1628123/150‑1 (MQ=255)
        tGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCACAGTGGTAACCGTGTACCGAGTCAAACAGGCGCTCGGTCGGTCCCCTCGCCCCTCGggggcgcgg                                                                                                                                            >  2:1403674/1‑145 (MQ=255)
        tGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                            >  2:656996/1‑150 (MQ=255)
        tGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                            <  2:560617/150‑1 (MQ=255)
            aTCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTa                                                                                                                                        <  2:1700304/150‑1 (MQ=255)
            aTCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCTCTCCCCAGAGGGGCGAggggacc                                                                                                                                       <  2:1039/150‑5 (MQ=38)
             tCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAg                                                                                                                                       >  1:1503348/1‑150 (MQ=255)
              cTTTTCCTCATCGGGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAgg                                                                                                                                      <  1:1496272/150‑1 (MQ=255)
                    cTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgg                                                                                                                                >  1:958321/1‑150 (MQ=255)
                    cTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgg                                                                                                                                >  1:84028/1‑150 (MQ=255)
                     tCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAggg                                                                                                                               >  1:1385912/1‑150 (MQ=255)
                      cATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgggg                                                                                                                              >  1:67430/1‑150 (MQ=255)
                       aTCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGa                                                                                                                             <  2:1098321/150‑1 (MQ=255)
                        tCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaa                                                                                                                            >  1:1059059/1‑150 (MQ=255)
                        tCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaa                                                                                                                            >  2:1146021/1‑150 (MQ=255)
                              cTGATGCAGTTCGATAATCTGGCTTACGGTTGGCCGTATCTGGGGTTCACGCCGCATCGGGAAATGTTTTAACGGTTTACAGAGTAAAATTGGCGCTGGGTGGGTCCCCTCGCCCCTCTGGGTTGTGGTTTTGGGTGTGGGGAAAACCgt                                                                                                                      <  1:1268149/150‑1 (MQ=255)
                                 aTGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTc                                                                                                                   >  1:887668/1‑150 (MQ=255)
                                   gCAGTTCGCTAATCAGGCTTACCGTAGGCCGGATCTGGCGTTCACGCCGCATCCGGATAAGTGTTAACGGTGAACCGAGTCAAATTGGCGCTGGTTCGGTCCCCTCGCCCCTCTGGGGATAGGGTTTGGGTGAGGGGAAAACCGTGTCAg                                                                                                                 <  2:1389567/150‑1 (MQ=37)
                                     aCTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCa                                                                                                               >  2:308601/1‑150 (MQ=255)
                                     aCTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCa                                                                                                               >  2:1284239/1‑150 (MQ=255)
                                       ttCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAat                                                                                                             >  1:537338/1‑150 (MQ=255)
                                        tCGATAATCAGGATTACCTTAGGCCGTATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTATGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAata                                                                                                            <  2:139876/150‑1 (MQ=39)
                                          gCTAATCAGGCTTACCGTAGGCCGGATCCGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCGCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGGGAGGGGAACACCGTGTCCGCCATATc                                                                                                          >  2:1372113/1‑150 (MQ=38)
                                          gCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGTGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATc                                                                                                          >  2:425448/1‑150 (MQ=255)
                                           cTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCa                                                                                                         >  1:1022949/1‑150 (MQ=255)
                                                cAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTCGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCgt                                                                                                    >  2:351375/1‑150 (MQ=255)
                                                cAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCgt                                                                                                    >  1:935247/1‑150 (MQ=255)
                                                cAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCGCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCgt                                                                                                    >  1:429551/1‑150 (MQ=255)
                                                 aGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCgtg                                                                                                   >  1:1295340/1‑150 (MQ=255)
                                                 aGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCgtg                                                                                                   >  2:880767/1‑150 (MQ=255)
                                                 aGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCgtg                                                                                                   >  1:776678/1‑150 (MQ=255)
                                                   gCTTACCGGAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGtt                                                                                                 <  2:562415/150‑1 (MQ=255)
                                                    cTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAATCCGTGTCAGCAATATCAACCGTGttt                                                                                                >  2:867642/1‑150 (MQ=255)
                                                    cTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGttt                                                                                                >  2:867660/1‑150 (MQ=255)
                                                          gTAGGCCGGATCTGGCGTTCACGCCGCATCCGGAAAAGTGTTAACCGTGTACAGATTCAAAATGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTTGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttctt                                                                                          <  2:387832/150‑1 (MQ=39)
                                                          gTAGGCCGGATCTGGCGTTCACGCCGAATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTTGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttctt                                                                                          <  2:1205895/150‑1 (MQ=255)
                                                          gTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttctt                                                                                          <  2:1288469/150‑1 (MQ=255)
                                                                      aGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAaacaac                                                                              >  1:673423/1‑150 (MQ=255)
                                                                      aGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAaacaac                                                                              >  2:140310/1‑150 (MQ=255)
                                                                      aGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAaacaac                                                                              >  1:16674/1‑150 (MQ=255)
                                                                        gCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCATCCGTGTTTATTCTTCGCAAAATCAAcgc                                                                            >  2:143464/1‑150 (MQ=255)
                                                                        gCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgc                                                                            >  1:282003/1‑150 (MQ=255)
                                                                        gCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgc                                                                            >  1:446869/1‑150 (MQ=255)
                                                                           ttCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgcgcg                                                                         >  1:1225812/1‑150 (MQ=255)
                                                                           ttCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgcgcg                                                                         >  1:359926/1‑150 (MQ=255)
                                                                              aCGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCa                                                                      >  2:743572/1‑150 (MQ=255)
                                                                                       ccGGAAAAGTGTTAACCGTGTACAGAGTCAAAATGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTca                                                             <  1:1688707/150‑1 (MQ=255)
                                                                                          gcaacttgtTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGTGTGGTTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagca                                                          <  2:1185248/144‑1 (MQ=39)
                                                                                           cAAAGTGTTACCCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGCCCCCTCGCCCCTCTGGGGAGTGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagcag                                                         <  2:1235304/150‑1 (MQ=255)
                                                                                             atgttttAACCGTGTACTGTGTCAAAATGGCGCTGGTTCGGTCCCCTCGCCCCTGTGGGGAGAGGGTTTGGGTTCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagcagca                                                       <  2:650909/145‑1 (MQ=255)
                                                                                             atgtgtTAACCGTGTACAGAGTCAAACAGGCGCTCGGTCGGTCCCCTCGCCCCTGTGGGGAGTGGGTTAGGGTGTGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagcagca                                                       <  2:1614425/148‑1 (MQ=255)
                                                                                                   tAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGc                                                 >  2:687779/1‑150 (MQ=255)
                                                                                                      ccGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGcaca                                              >  1:1653573/1‑150 (MQ=255)
                                                                                                        gtgtACAGAGTCAAAACGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGTGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAgg                                            <  1:1543892/150‑1 (MQ=255)
                                                                                                        gtgtACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAgg                                            <  2:1059059/150‑1 (MQ=255)
                                                                                                        gtgtACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAgg                                            >  2:305146/1‑150 (MQ=255)
                                                                                                             ccgagTCAAACAGGCGCTGGGTGGGTCCCCCCGCCCCTCTGGGGTTTGGTTTTGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAAt                                       <  1:528928/148‑1 (MQ=255)
                                                                                                             catagTCAAAACGGCGCGGGGTGGGCCCCCGCGCCCCTGTGGGGTGTGGGTTAGGGTGTGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAAt                                       <  1:296414/147‑1 (MQ=255)
                                                                                                              agagTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATg                                      >  1:1099499/1‑150 (MQ=255)
                                                                                                               nagTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTCGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGCATGTTCAGCAGCAGTACGCACAggggcctga                                     >  1:202/2‑144 (MQ=255)
                                                                                                               gagTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCGGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGa                                     >  1:1586999/1‑150 (MQ=255)
                                                                                                                agTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGAt                                    >  1:643049/1‑150 (MQ=255)
                                                                                                                agTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGAt                                    >  1:1323201/1‑150 (MQ=255)
                                                                                                                agTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTCCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGAt                                    >  1:378383/1‑150 (MQ=255)
                                                                                                                agTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTCCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGAt                                    >  1:434244/1‑150 (MQ=255)
                                                                                                                agTAAAAATGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGTGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGAt                                    <  1:917590/150‑1 (MQ=255)
                                                                                                                  tCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCa                                  >  1:1164019/1‑150 (MQ=255)
                                                                                                                  tCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCa                                  <  1:1409216/150‑1 (MQ=255)
                                                                                                                      aaaGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATAc                              >  2:557965/1‑150 (MQ=255)
                                                                                                                        aGGCGCTGGGTCGGTCCCCTCGCCCCTCTGTGGTGTGGTTTTGGGTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATAccc                            <  1:506953/150‑1 (MQ=255)
                                                                                                                           ggcgcgggcgggCCCCGCGCCCCGGGGGGGTGGGGGTTTGGGTGCGGGGAAACCCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCgg                         <  1:1191008/146‑1 (MQ=255)
                                                                                                                            gcTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGc                        >  1:605136/1‑150 (MQ=255)
                                                                                                                              tcggtcggTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGccc                      >  2:1504515/1‑150 (MQ=255)
                                                                                                                               cggtcggtCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCa                     <  2:840656/150‑1 (MQ=255)
                                                                                                                                   cggtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGcc                 <  1:921360/150‑1 (MQ=255)
                                                                                                                                   cggtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGcc                 <  1:869255/150‑1 (MQ=255)
                                                                                                                                    ggtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCa                >  2:394318/1‑150 (MQ=255)
                                                                                                                                        cgggcGCCCCTGTGGGGTGTGGTTTTGGGTGTGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCTGCTGTACGCACAGGGTATTGATCAATACCCCGGCCCAGGCCAATtg            <  1:1594617/146‑1 (MQ=255)
                                                                                                                                         ccTCGCCCATCGTGGGTGTGGGTTAGTGTTCGGTGAAAACCTTGTCAGCAATTTCAACAGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTTCTCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgc           <  1:694920/150‑1 (MQ=255)
                                                                                                                                          cccGCCCCTCTGGGGAGTGGGTTTGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgct          <  2:601792/148‑1 (MQ=255)
                                                                                                                                          cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgct          >  2:1724714/1‑150 (MQ=255)
                                                                                                                                          cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgct          >  2:452587/1‑150 (MQ=255)
                                                                                                                                            cGCCCCTCGGGGGAGAGGTTTCGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgctgc        <  1:779198/150‑1 (MQ=255)
                                                                                                                                            cGCCCCCCTGGGGAGAGGGTTAGGGTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgctgc        <  1:779195/150‑1 (MQ=255)
                                                                                                                                               cccTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtg            >  1:1008036/1‑143 (MQ=255)
                                                                                                                                               cccTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtg            <  2:1008036/143‑1 (MQ=255)
                                                                                                                                               cccTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCAt     >  1:1623091/1‑150 (MQ=255)
                                                                                                                                                 ctctGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTc   >  2:724848/1‑150 (MQ=255)
                                                                                                                                                 ctctGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTc   >  2:898499/1‑150 (MQ=255)
                                                                                                                                                  tctGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTcc  <  1:1284239/150‑1 (MQ=255)
                                                                                                                                                  tctGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAAATGCTGCCATTcc  >  2:1151042/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCC  >  NC_000913/3908359‑3908654

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: