Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,908,512 A→C 6.8% intergenic (‑146/+37) 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,5120AC 6.8% 289.7 / 10.7 114intergenic (‑146/+37)pstB/pstAphosphate ABC transporter ATPase/phosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (55/46);  new base C (4/4);  total (61/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-01

Only 100 of 115 total aligned reads displayed.

CATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAA  >  NC_000913/3908364‑3908659
                                                                                                                                                    |                                                                                                                                                   
cATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGgag                                                                                                                                                    >  2:1174027/1‑150 (MQ=255)
   tGCACTCTTTTCCTCATCGTGCCTGATGCACTTCGCTCATCACGCTAACCGTAGGCCGGATCAGGCGATCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCTCCCCAGAGGGGCGAgg                                                                                                                                                 <  2:698723/150‑1 (MQ=255)
    ggattCTTTTCCTCATCGTGCCCGTTGCACCTCGCTTATCACGCCTACCGTAGGCCGGAGCAGGCGTGCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGGCAAAAAGGCGCTCGGTCGGTCCCCTCTCCCCAGAGGGGCGAggg                                                                                                                                                <  1:658442/146‑1 (MQ=255)
    gCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAggg                                                                                                                                                <  1:365056/150‑1 (MQ=255)
    gCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAggg                                                                                                                                                <  1:1098380/150‑1 (MQ=255)
       aTCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCTCCCCAGAGGGGCGAggggac                                                                                                                                             <  1:1219623/150‑4 (MQ=39)
       aTCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTa                                                                                                                                             <  2:1335307/150‑1 (MQ=255)
        tCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCTCCCCTCTGGGGAGAGGGTTAg                                                                                                                                            >  2:1395423/1‑150 (MQ=255)
          ttttCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAggg                                                                                                                                          <  2:712462/150‑1 (MQ=255)
          ttttCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAggg                                                                                                                                          >  1:1445780/1‑150 (MQ=255)
            ttCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTg                                                                                                                                        <  2:1683471/150‑1 (MQ=255)
             tCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGa                                                                                                                                       >  1:673759/1‑150 (MQ=255)
              ccTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGTGTTAGGGTGAg                                                                                                                                      <  1:1335290/150‑1 (MQ=255)
               cTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgg                                                                                                                                     <  1:625067/150‑1 (MQ=255)
                tCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAggg                                                                                                                                    <  2:1263621/150‑1 (MQ=255)
                  aTCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTACCCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGCGGGTTAGGGTGAGGGGa                                                                                                                                  >  2:1054542/1‑150 (MQ=255)
                   tCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaa                                                                                                                                 >  2:745105/1‑150 (MQ=255)
                    cGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaa                                                                                                                                >  1:1382884/1‑150 (MQ=255)
                      tGCCTGATGCACTTCGCTATTGAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGAAAAGTGTTAACCGTGTACAGAGTCAAAAAGGGGCGCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAc                                                                                                                              <  1:1314315/150‑1 (MQ=38)
                       gCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTTTACAGAGTCAAAAAGGCGCTCGGTCGGGCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAcc                                                                                                                             <  1:826687/150‑1 (MQ=255)
                        ccTGATTCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGGGTTAACCGGGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCGCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCg                                                                                                                            <  1:1606675/150‑1 (MQ=255)
                        ccTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGAGATCGGaagagta                                                                                                                            >  1:1389641/1‑145 (MQ=255)
                           gATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAACGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCgtgt                                                                                                                         <  1:722977/150‑1 (MQ=255)
                                               cTTACCGTAGGCCGGATCAGGCGTTCACCCCGCATCCGGAAATGTGTTAACCGTGTACAGAGTCAAAACGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGTGTGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGttt                                                                                                     <  2:1445780/150‑1 (MQ=39)
                                                 tACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAACCCGTGTCAGCAAGATCAACCGTGTTTAt                                                                                                   >  2:897344/1‑150 (MQ=255)
                                                      tAGGCCGGAGCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttcttc                                                                                              <  1:25798/150‑1 (MQ=255)
                                                       ttgCCGGATGCGGCGTGAACGCCTGATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCg                                                                                             >  2:893265/3‑150 (MQ=38)
                                                       tggCCGGAGCAGGCGTTCACGCCGCATCCGGAAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGTGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCg                                                                                             <  2:1356801/149‑1 (MQ=255)
                                                       aGGCGGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTAAAGATTCAAAATGGGGCGCGGGCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTTGGGTGTGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCg                                                                                             <  2:1341880/150‑1 (MQ=38)
                                                       aGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCg                                                                                             <  1:632373/150‑1 (MQ=255)
                                                       aGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCg                                                                                             >  2:922094/1‑150 (MQ=255)
                                                        ggCCGGAGCAGGGGTGCACGCCGCATCGGGAAAAGTGTTAACCGTGTACAGAGTCAAAACGGCGCTGGGGCGGGCCCCTCGCCCCTCTGGGGTTAGGGTTAGGGGGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGc                                                                                            <  2:1253564/150‑1 (MQ=255)
                                                          ccGGATGCGGCGTGAACGCCTGATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTTGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCaa                                                                                          >  1:747257/1‑150 (MQ=38)
                                                          ccGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGCGTGTCAGCAATATCACCCGTGTTTATTCTTCGCag                                                                                          >  1:117320/1‑149 (MQ=255)
                                                           cGGATCAGGCGTTCACGCCGCATCCGGAAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCaaa                                                                                         <  2:388198/150‑1 (MQ=255)
                                                             gATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGATGTGAAGACCGTGTCAGCAATATCAACAGGGTATATGCTTCGcctaac                                                                                       >  1:94778/1‑145 (MQ=38)
                                                               tCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAaaca                                                                                     <  2:956856/150‑1 (MQ=255)
                                                                cAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAaacaa                                                                                    >  1:1294678/1‑150 (MQ=255)
                                                                cAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAaacaa                                                                                    >  2:740336/1‑150 (MQ=255)
                                                                cAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGATTATTCTTCGCAAAaacac                                                                                    >  2:934093/1‑149 (MQ=255)
                                                                  ggCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcg                                                                                  >  2:1466517/1‑150 (MQ=255)
                                                                  ggCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCATCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcg                                                                                  >  2:251087/1‑150 (MQ=255)
                                                                  ggCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCATCGCCCCTCTGGGGAGAGGGTGAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAacaacacag                                                                                  >  2:282038/1‑147 (MQ=255)
                                                                    cGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGATGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgcg                                                                                <  1:113794/150‑1 (MQ=255)
                                                                        cACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGcc                                                                            >  2:191483/1‑150 (MQ=255)
                                                                          cgccgcATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAACCAACGCGCGCCAg                                                                          >  1:842014/1‑150 (MQ=255)
                                                                          cgccgcATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAg                                                                          >  2:882153/1‑150 (MQ=255)
                                                                               cATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGa                                                                         >  2:268823/1‑146 (MQ=255)
                                                                               cATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGa                                                                         <  1:268823/146‑1 (MQ=255)
                                                                                 tCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCGGGGGAGAGGGTTAGGGTGAGGGGAAATCCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTc                                                                   >  1:320662/1‑150 (MQ=255)
                                                                                  ccGGCAAAGTGTTAACCGTGTACAGCGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCACTGGGGATAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTca                                                                  <  1:84077/150‑1 (MQ=255)
                                                                                  ccGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTca                                                                  <  1:1031219/150‑1 (MQ=255)
                                                                                   cGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcag                                                                 >  2:259544/1‑150 (MQ=255)
                                                                                   cGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcag                                                                 >  1:1041895/1‑150 (MQ=255)
                                                                                    ggCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagc                                                                >  1:458978/1‑150 (MQ=255)
                                                                                     gCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATAAACAGTGTTTATTCTTAGCAAAAACAACGCTCGCCAGAATGTTcacca                                                               >  2:1460921/1‑150 (MQ=38)
                                                                                        aaGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagcagca                                                            <  2:228680/150‑1 (MQ=255)
                                                                                              tAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGc                                                      <  2:404352/150‑1 (MQ=255)
                                                                                                  cGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAg                                                  >  1:1151560/1‑150 (MQ=255)
                                                                                                      tACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTa                                              >  1:419125/1‑150 (MQ=255)
                                                                                                       aCAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTaa                                             >  1:1288208/1‑150 (MQ=255)
                                                                                                           agTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGAt                                         >  2:393387/1‑150 (MQ=255)
                                                                                                            gTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATc                                        >  2:1496969/1‑150 (MQ=255)
                                                                                                            gTCAAAAAGGCGCTCGGTCGGGCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATc                                        <  2:427403/150‑1 (MQ=255)
                                                                                                                      ggcTCGGTCGGTCCCCTCGCCCCTGTGGGGAGAGGGTTAGGGTGAGGGGAAAACAGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCgg                              <  1:934093/149‑1 (MQ=255)
                                                                                                                      cgcTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAACACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACACCGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCgg                              >  2:340108/1‑150 (MQ=255)
                                                                                                                       gcTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGc                             >  1:1638702/1‑150 (MQ=255)
                                                                                                                        cgcggtcggcCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGcc                            <  2:320662/148‑1 (MQ=255)
                                                                                                                        cTCGGTCGGTCCCCTCGCCCCTCTGGGGAGTGGGTTATGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTATTGATCATTAccctgtgc                            <  2:457414/150‑6 (MQ=39)
                                                                                                                        cTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGcc                            >  2:1413185/1‑150 (MQ=255)
                                                                                                                            gtcggtaCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCTGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAgg                        >  2:623176/1‑150 (MQ=255)
                                                                                                                            accctaacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAAGATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAgg                        >  1:491117/8‑150 (MQ=255)
                                                                                                                             ccctaacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGACAACCGGGTCAGCAAGATCAACCGTGTTTATTCTTCGCAAACGCAGCGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTACTGATCTTTCCCCCGGCCCAGGc                       >  1:857631/7‑150 (MQ=255)
                                                                                                                              cggtcCCCTCGCCCCTCTGGGGTGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCCACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGAACGCACAGGGTAATGATCAATACCCCGGCCCAGGcc                      >  1:1335820/1‑150 (MQ=255)
                                                                                                                              cggtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATAcccc                                <  2:166129/140‑1 (MQ=255)
                                                                                                                              cggtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATAcccc                                <  2:171720/140‑1 (MQ=255)
                                                                                                                              cggtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATAcccc                                >  1:166129/1‑140 (MQ=255)
                                                                                                                              cggtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATAcccc                                >  1:171720/1‑140 (MQ=255)
                                                                                                                              cggtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGcc                      <  2:1638702/150‑1 (MQ=255)
                                                                                                                              cctaacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTCCGCACAGGGTAATGATCACTACCCCGGCCCAGGcc                      >  1:537807/6‑150 (MQ=255)
                                                                                                                               ggtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAAt                                     >  1:30849/1‑134 (MQ=255)
                                                                                                                               ggtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAAt                                     <  2:30849/134‑1 (MQ=255)
                                                                                                                                gtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCaa                    >  1:654438/1‑150 (MQ=255)
                                                                                                                                gtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCaa                    >  2:580214/1‑150 (MQ=255)
                                                                                                                                  ccccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAAtt                  >  2:1090917/1‑150 (MQ=255)
                                                                                                                                   cccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAATCTACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtg                 >  1:266578/1‑150 (MQ=255)
                                                                                                                                   cccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGTGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtg                 >  1:1574229/1‑150 (MQ=255)
                                                                                                                                     cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgct               >  2:1484454/1‑150 (MQ=255)
                                                                                                                                     cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgct               <  1:1670485/150‑1 (MQ=255)
                                                                                                                                      tCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgctg              >  2:1571898/1‑150 (MQ=255)
                                                                                                                                        gCCCCTGTGGGGAGAGGGTTAGGGTTAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGcc            <  1:1054542/150‑1 (MQ=255)
                                                                                                                                        gCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGcc            <  2:1294678/150‑1 (MQ=255)
                                                                                                                                        gCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGGCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGcc            <  1:893265/150‑1 (MQ=255)
                                                                                                                                         ccccTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCa           >  1:907671/1‑150 (MQ=255)
                                                                                                                                         ccccTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCa           >  2:24930/1‑150 (MQ=255)
                                                                                                                                          cccTCTGGGGAGAGGTTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACATGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCAt          <  1:970235/150‑1 (MQ=255)
                                                                                                                                            ctctGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTc        >  1:1262192/1‑150 (MQ=255)
                                                                                                                                              ctGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCg      <  2:57949/150‑1 (MQ=255)
                                                                                                                                                ggggAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCa    >  2:1217538/1‑150 (MQ=255)
                                                                                                                                                  ggAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCaaa  >  1:1409688/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAA  >  NC_000913/3908364‑3908659

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: