Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,908,507 T→G 6.8% 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 6.8% 629.8 / 6.4 236intergenic (‑141/+42)pstB/pstAphosphate ABC transporter ATPase/phosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (118/95);  new base G (4/17);  total (124/112)
Fisher's exact test for biased strand distribution p-value = 2.19e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 249 total aligned reads displayed.

CATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGC  >  NC_000913/3908358‑3908656
                                                                                                                                                     |                                                                                                                                                     
cATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCctct                                                                                                                                                       <  1:1254291/150‑1 (MQ=255)
     tCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGga                                                                                                                                                  <  1:275209/150‑1 (MQ=255)
        ttGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                               <  1:811530/150‑1 (MQ=255)
        ttGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                               >  1:51765/1‑150 (MQ=255)
        ttGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                               >  2:477750/1‑150 (MQ=255)
         tGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                              <  1:684029/150‑1 (MQ=255)
         tGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                              >  1:1647211/1‑150 (MQ=255)
         tGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                              >  1:1779635/1‑150 (MQ=255)
         tGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                              <  1:3205351/150‑1 (MQ=255)
             aTCTTTTCCTCCTCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTa                                                                                                                                          <  1:3043701/150‑1 (MQ=255)
                ttttCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGGTAggg                                                                                                                                       >  1:2222443/1‑150 (MQ=255)
                 tttCCTCCTCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGt                                                                                                                                      <  2:2950214/150‑1 (MQ=255)
                 tttCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGt                                                                                                                                      >  2:168174/1‑150 (MQ=255)
                 tttCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGt                                                                                                                                      >  2:1818277/1‑150 (MQ=255)
                    ccTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAg                                                                                                                                   >  2:2727492/1‑150 (MQ=255)
                     ctcctcGTGCCTGCTGCACGTCGCGAAGCAGGCGTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTACCCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgg                                                                                                                                  <  2:2510864/150‑1 (MQ=38)
                     cTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgg                                                                                                                                  <  2:1212212/150‑1 (MQ=255)
                              ccTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGACAACCg                                                                                                                         >  2:1982633/1‑150 (MQ=255)
                                 gttGCACTTCGATAATCTGGCTTACGGTTGGCCGGATCTGGCGTTCACGCCGCATCGGGAAAAGTGTTAACCGTGTACTGAGTAAAAATGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGTTTTGGGTTTGGGTGTGGGGAAAACCgtgt                                                                                                                      <  1:2250569/148‑1 (MQ=255)
                                 gATGCACTGCGCTAATCCGGCTTACCGTCGGCCGGATCTGGCGTTCCCGCCGCAGCCGGCAATGTGTGAACCGGGTACCGAGTCAAAAAGGCGCTCGGTCGGTCCCCGCGCCCCTCGGGGGTGCGGGTTTGGGTGCGGGGAAAACCgtgt                                                                                                                      <  2:257233/150‑1 (MQ=255)
                                   tGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCa                                                                                                                    >  2:707447/1‑150 (MQ=255)
                                        ttCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAat                                                                                                               >  1:1842432/1‑150 (MQ=255)
                                        ttCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAat                                                                                                               >  1:682654/1‑150 (MQ=255)
                                        ttCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAat                                                                                                               >  1:1376677/1‑150 (MQ=255)
                                        ttCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAat                                                                                                               >  1:2587212/1‑150 (MQ=255)
                                        ttCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAat                                                                                                               <  2:1064338/150‑1 (MQ=255)
                                         tCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAata                                                                                                              >  2:743842/1‑150 (MQ=255)
                                         tCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAata                                                                                                              >  1:2458501/1‑150 (MQ=255)
                                           gCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATc                                                                                                            >  1:2019473/1‑150 (MQ=255)
                                           gCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATc                                                                                                            >  1:3077167/1‑150 (MQ=255)
                                           gCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATc                                                                                                            >  1:2159987/1‑150 (MQ=255)
                                           gCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCGGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATc                                                                                                            >  1:1467404/1‑150 (MQ=255)
                                            cTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCa                                                                                                           <  2:2084721/150‑1 (MQ=255)
                                             tAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTCCCGGGGCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCGCGGGGGCGTGGGTGAGGGGGAGGGGAAGACCGTGTCAGCaagctcac                                                                                                          >  2:2103027/1‑144 (MQ=37)
                                              aaTCTGGCTTACCGTCGGCCGGATCTGGCGTTCCCGCCGCATCCGGCAAAGTGTTAACCGTTTACTGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTTGGGTGTGGGGAAAACCGTGTCAGCAATATCAAc                                                                                                         <  1:750641/150‑1 (MQ=38)
                                              aaTCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAAc                                                                                                         <  2:1281892/150‑1 (MQ=255)
                                               aTCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAAcc                                                                                                        >  2:1577766/1‑150 (MQ=255)
                                                 ctggCTTACCGTCGGCGGGATCAGGCGTTCCCGCCGCATCGGGCAAAGTGTTAACCGTGTACCGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCGGGGGTGTGGGTTGGGGTGTGGGGAAAACCGTGTCAGCAATATCAACCgt                                                                                                      <  2:2224919/148‑1 (MQ=255)
                                                    gCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGtt                                                                                                   >  1:1437537/1‑150 (MQ=255)
                                                     cTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGttt                                                                                                  >  1:1923295/1‑150 (MQ=255)
                                                     cTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGttt                                                                                                  >  1:1327200/1‑150 (MQ=255)
                                                     cTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGttt                                                                                                  >  1:2832345/1‑150 (MQ=255)
                                                     cTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGttt                                                                                                  >  1:238523/1‑150 (MQ=255)
                                                      ttACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGGTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTa                                                                                                 >  2:388892/1‑150 (MQ=255)
                                                        aCCTTAGGCCGGATCTGGCGTACACGCCGCATCCGGCAAAGTTTTAACCGTGTAAAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGTTGAGGTTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAtt                                                                                               <  2:133453/150‑1 (MQ=38)
                                                          cGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttct                                                                                             >  1:468655/1‑150 (MQ=255)
                                                            tAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAATGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttcttc                                                                                           <  2:1083440/150‑1 (MQ=255)
                                                            tAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttcttc                                                                                           <  2:2396121/150‑1 (MQ=255)
                                                            gcggCCGGAGCCGGCGTTCACGCCGCAGCGGGCAATTGGTTAACCGTGTACAGATTCAAAATGGCGGGCGGTGGGTCCCCTCGCCCCTCTGGGGTTTGGGTTTGGGTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttcttc                                                                                           <  2:1145316/148‑1 (MQ=255)
                                                             cggCCGGATCAGGCGTTCACGCCGCATCCGGAAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGTGGTTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCg                                                                                          <  2:334861/149‑1 (MQ=255)
                                                               gCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCa                                                                                        >  1:392645/1‑150 (MQ=255)
                                                                   gATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCaaaaa                                                                                    >  1:1143755/1‑150 (MQ=255)
                                                                       aGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAaacaac                                                                                >  2:3087750/1‑150 (MQ=255)
                                                                         gCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgc                                                                              >  2:2264748/1‑150 (MQ=255)
                                                                         gCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATAGCAACCGTGTTTATTCTTCGCAAAAACAAcgc                                                                              >  2:2374045/1‑150 (MQ=255)
                                                                          cGTTCACCCCGCAACCGCAAAAGTGTTAACCGTCTACAGACACAAAAAGGCGCTGGGTCGGCCCCCTCGCCCCTCGGGGTAGAGGGTCCGGGTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgcg                                                                             <  1:614866/150‑1 (MQ=255)
                                                                            ttCCCGCCGCAGCCGGCCCGGGGTTACCCGGGCACCGTGTCAAACCGGGGCGGGGGGGGGCCCCTCGCCCCGGGGGGGCGCGGGTTGGGGTGGGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgcgcg                                                                           <  2:1177890/150‑1 (MQ=255)
                                                                            ttCACGCCGCAACCGGCAAAGTGTTAACCGTGAACAGAGTCAACAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgcgcg                                                                           <  2:2134120/150‑1 (MQ=255)
                                                                            tccaggccgcATCCGGCATTTTTTTAACGGTGTACTGAGTAAAATTGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGTTGCGGTTTTGGGTGTGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgcgcg                                                                           <  2:1467402/145‑1 (MQ=255)
                                                                             tCACGCCGCATCGGGCAAAGTGTTAACCGTGTACAGAGTCAAATAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTTGGTTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACgcgcgc                                                                          <  2:149284/150‑1 (MQ=255)
                                                                             tCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTTGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACgcgcgc                                                                          <  2:2961100/150‑1 (MQ=255)
                                                                                cgccgcATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAg                                                                       >  2:2620362/1‑150 (MQ=255)
                                                                                cgccgcATCCGGAAAAGTGTTAACCGTGTACAGAGTCAAAACGGCGCTCGGTCGGACCCCTCGCCCCTCTGGGGAGTGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAg                                                                       <  2:3056589/150‑1 (MQ=255)
                                                                                 gccgcATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGCCCCCTCGCCCCTCTGGGGTGTGGGTTTGGGTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGa                                                                      <  1:431168/150‑1 (MQ=255)
                                                                                       cccGGCAAAGTGTTAACCGTGTACAGAGACAAAATGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGTTGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTc                                                                <  2:2764334/149‑1 (MQ=255)
                                                                                         cGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGACCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAAGATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGACTGTTcag                                                              >  1:2131980/1‑150 (MQ=255)
                                                                                           gCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagca                                                            >  1:262536/1‑150 (MQ=255)
                                                                                              atgtgtAACCCGCGTACAGACTCACACCGGCGGTCGGGCGGCCCCCTCGCCCCGGGGGGGATTGGGTTCGGGTGTGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAATAACTACGCGCGCCTGAATGTTcagcagca                                                         <  1:2548844/148‑1 (MQ=255)
                                                                                              aattttaaaCCGTGTACAGACAAAAAAAGCCGCTCGGTCGGTCCCCTCCCCCCTCTGGTGAGTGGTTTAGGTTGAGGGCAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagcagca                                                         <  1:926053/143‑1 (MQ=255)
                                                                                              aaGTGTTACCCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCCCGCCCCTCTGGGGAGTGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagcagca                                                         <  1:1077089/150‑1 (MQ=255)
                                                                                                  gtTAACGGTGTACAGAGTCAAAATGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGTGTGGGTTAGGGTGTGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCTCCAGAATGTTCAGCAGCAGTAc                                                     <  1:1992467/150‑1 (MQ=39)
                                                                                                     aaCCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGca                                                  >  2:1873856/1‑150 (MQ=255)
                                                                                                        cGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAg                                               <  2:2025509/150‑1 (MQ=255)
                                                                                                         gtgtACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAgg                                              >  2:609064/1‑150 (MQ=255)
                                                                                                              cAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAAt                                         <  2:1596634/150‑1 (MQ=255)
                                                                                                                 agTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGTGGTTTTGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGAt                                      <  2:1243425/150‑1 (MQ=255)
                                                                                                                 agTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGAt                                      >  1:1609071/1‑150 (MQ=255)
                                                                                                                 agTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGAAGTACGCACAGGGTAATGAt                                      >  1:2916497/1‑150 (MQ=255)
                                                                                                                   tCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCa                                    >  1:2967109/1‑150 (MQ=255)
                                                                                                                        aaGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGTAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATAcc                               <  2:2874726/150‑1 (MQ=255)
                                                                                                                        aaGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATAcc                               >  2:3237873/1‑150 (MQ=255)
                                                                                                                         tggCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGTGTGGTTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATAccc                              <  2:20542/149‑1 (MQ=255)
                                                                                                                         cggCGCTCGGTCGGGCCCCTCGCCCCTGTGGGGTGTGGGTTCGGGTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATAccc                              <  2:1405576/149‑1 (MQ=255)
                                                                                                                         aGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGTTTAGGGTGTGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATAccc                              <  2:250727/150‑1 (MQ=255)
                                                                                                                            cgcTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAATCAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTCATGATCAATACCCCgg                           >  2:2785563/1‑150 (MQ=255)
                                                                                                                              ctgggtcggccCCCGCGCCCGTGTGGGGAGTGGTTTTGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGcc                         <  2:904787/147‑1 (MQ=255)
                                                                                                                                cggtcggtCCCCGCGCCCCTCTGGGGAGAGGGTTAGGGTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCa                       <  1:288365/150‑1 (MQ=255)
                                                                                                                                 ggtcggtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAg                      >  2:971394/1‑150 (MQ=255)
                                                                                                                                 ggtcggtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAg                      >  2:970604/1‑150 (MQ=255)
                                                                                                                                   tcggtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGc                    <  2:2161196/150‑1 (MQ=255)
                                                                                                                                    cctaacccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGcc                   >  1:2698094/6‑150 (MQ=255)
                                                                                                                                       tcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAAt                >  2:2426137/1‑150 (MQ=255)
                                                                                                                                       tcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAAt                <  1:2374071/150‑1 (MQ=255)
                                                                                                                                         cccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtg              >  1:2050102/1‑150 (MQ=255)
                                                                                                                                         cccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtg              >  1:361033/1‑150 (MQ=255)
                                                                                                                                           cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgct            <  2:1058497/150‑1 (MQ=255)
                                                                                                                                                 ccTGTGGGGTGAGGGTTCGGGTGCGGGCAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCAtt      <  1:2727468/150‑1 (MQ=255)
                                                                                                                                                 ccTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCAtt      <  2:2443316/150‑1 (MQ=255)
                                                                                                                                                  ctctGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTc     >  1:2434116/1‑150 (MQ=255)
                                                                                                                                                   tctGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTcc    >  1:3241684/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGC  >  NC_000913/3908358‑3908656

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: