Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,908,512 A→T 6.9% 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,5120AT 6.9% 1861.8 / 20.4 649intergenic (‑146/+37)pstB/pstAphosphate ABC transporter ATPase/phosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (359/216);  new base T (2/58);  total (370/279)
Fisher's exact test for biased strand distribution p-value = 1.04e-20
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 685 total aligned reads displayed.

TCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAACG  >  NC_000913/3908363‑3908661
                                                                                                                                                     |                                                                                                                                                     
tCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGga                                                                                                                                                       >  2:6706014/1‑150 (MQ=255)
tCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTCACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGga                                                                                                                                                       >  2:3162933/1‑150 (MQ=255)
tCATTGCAATCTTTTCATCATCGTGCCTGATGCACTTCGCTAATCAGGCTTAAAGTAGGCCGTATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAAACGTGTACAGAGTCAAAAAGGAGATCGGTCTGTCCCCTCGCCACTCTTGGga                                                                                                                                                       >  2:3402395/1‑150 (MQ=255)
  aTTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGgaga                                                                                                                                                     <  2:3283779/150‑1 (MQ=255)
  aTTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGgaga                                                                                                                                                     <  2:475304/150‑1 (MQ=255)
    tGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                                   <  1:5369727/150‑1 (MQ=255)
    tGCAATCTTGTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                                   <  2:9141154/150‑1 (MQ=255)
       aaTCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGtt                                                                                                                                                >  1:5867111/1‑150 (MQ=255)
       aaTCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTACCCTCTCCCCAGAGGGGCGAgggga                                                                                                                                                <  2:7020008/150‑3 (MQ=39)
         tCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAg                                                                                                                                              >  1:1206247/1‑150 (MQ=255)
         tCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAg                                                                                                                                              >  1:187043/1‑150 (MQ=255)
          cTTTTCCTCCTCGTGCCTGATGCACTTCGCTAATCAGGCTGACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAgg                                                                                                                                             <  1:6424863/150‑1 (MQ=255)
          cTTTGCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAgg                                                                                                                                             <  2:8737042/150‑1 (MQ=255)
            tttCCTCATCGTGCCTGATGCACTTCGCTATTCAGGCTTACCGTTGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTAAAATTGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGTGGGTTAGGGt                                                                                                                                           >  2:2088054/1‑150 (MQ=255)
            tttCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGt                                                                                                                                           >  1:3108272/1‑150 (MQ=255)
                 tCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAggg                                                                                                                                      <  2:7053754/150‑1 (MQ=255)
                 tCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAggg                                                                                                                                      >  2:4621640/1‑150 (MQ=255)
                 tCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAggg                                                                                                                                      >  2:4421704/1‑150 (MQ=255)
                     ggTGCCTGATGCATTTCTATAATCTGGCTTACCGTAGGCCGTATCTGGGGTTCACGCCGCATCCGTAAATGTTTTAACCGTTTACTGAGTAAAAATGGCGGTCGGTGGGTCCCCTCGCCCCTCTGGGTTGTGGGTTTGGGTGTGGGGaaa                                                                                                                                  <  1:4748970/149‑1 (MQ=255)
                      gTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaaaa                                                                                                                                 >  2:1475820/1‑150 (MQ=255)
                          cTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCgt                                                                                                                             <  2:2318666/150‑1 (MQ=255)
                           tGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCgtg                                                                                                                            >  1:1725116/1‑150 (MQ=255)
                            gATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCgtgt                                                                                                                           >  1:6775214/1‑150 (MQ=255)
                               gCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAg                                                                                                                        <  2:5971014/150‑1 (MQ=255)
                                  cTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCaa                                                                                                                     >  2:7385617/1‑150 (MQ=255)
                                   ttCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCGGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGTGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAat                                                                                                                    <  2:1349841/150‑1 (MQ=255)
                                   ttCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAat                                                                                                                    >  2:6594596/1‑150 (MQ=255)
                                   ttCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAat                                                                                                                    >  1:987185/1‑150 (MQ=255)
                                         aaTCTGGCTTACGGTTGGCCGGATCTGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCGGGGGTGTGGGTTTGGGTGAGGGGAAAACCGTGTCAGCAATATCAAc                                                                                                              <  2:279164/150‑1 (MQ=255)
                                            cAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCgt                                                                                                           >  1:2173033/1‑150 (MQ=255)
                                                 ttACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGTGTCAAAATGGCGCTGGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTa                                                                                                      <  2:8853709/150‑1 (MQ=255)
                                                 ttACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTa                                                                                                      >  2:9285968/1‑150 (MQ=255)
                                                    ccGTAGGCCGGAACAGGCGTTCACGCCGCATCCGGCAATGTTTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGTAGAGGGTTAGGTTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttc                                                                                                   <  2:2477431/150‑1 (MQ=39)
                                                     cGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttct                                                                                                  >  1:6186422/1‑150 (MQ=255)
                                                      gttggCCGGCTCAGGCGTTCACGCCGCATCCGTCAAAGTGTTAACCGTGTACAGAGTCAAAACGGCGCTCGGTGGGTCCCCTCGCCCCTCTGGGTAGTGGGTTAGGGTGTGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttctt                                                                                                 <  1:5575088/147‑1 (MQ=255)
                                                      gTAGGGCGGATCAGGCGTTCACGCCGCATCCGGAAATGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGGCCCGTCGCCCCTGTGGGTAGTGGTTTTGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttctt                                                                                                 <  2:3542162/150‑1 (MQ=255)
                                                      gTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttctt                                                                                                 >  1:4016230/1‑150 (MQ=255)
                                                           ccGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCaa                                                                                            >  1:3872132/1‑150 (MQ=255)
                                                           ccGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCaa                                                                                            >  1:6642106/1‑150 (MQ=255)
                                                           ccGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGGTAGGGGGAGGGGAAAACCGTGTCAGCACTATCAACCGTGTTTATTCTTCGCaa                                                                                            >  1:3563554/1‑150 (MQ=255)
                                                              gATCTGGCGTTCCCGCCGCATCCGGCAATGTGTTACCCGTGTACAGAGACAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGTTGTGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCaaaaa                                                                                         <  2:8015657/150‑1 (MQ=38)
                                                              gATCAGGCGTTCACGCCGCATCCGGAAATGTGTTAACCGTGTACAGAGTCAAAACGGCGCTCGGTCGGCCCCCTCGCCCCTCTGGGGTGTGGGTTTGGGTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCaaaaa                                                                                         <  1:9332648/150‑1 (MQ=255)
                                                                 cAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAACTGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGTGGGTTCGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAaacaa                                                                                      <  2:2637865/150‑1 (MQ=255)
                                                                 cAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAaacaa                                                                                      <  1:7809691/150‑1 (MQ=255)
                                                                   ggCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcg                                                                                    >  2:7893271/1‑150 (MQ=255)
                                                                   ggCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcg                                                                                    >  2:6139113/1‑150 (MQ=255)
                                                                    gCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgc                                                                                   >  1:4290663/1‑150 (MQ=255)
                                                                      gcaccccccgcAGCCGGAAATTTTTTAACGGTGTACTGAGAAAAAACGGCGGGGGGTGGGCCCCCTCGCCCGTGTGGGTTGTGGGTTTGGGTGGGGGCAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgcgc                                                                                 <  1:2732603/143‑1 (MQ=255)
                                                                             ccgcATCCGGCAAAGTGTTAACCGTGTACAGATTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGTAGAGGGTTTGGGTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGaa                                                                          <  2:841380/150‑1 (MQ=255)
                                                                             ccgcATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGaa                                                                          <  2:5561880/150‑1 (MQ=255)
                                                                             ccgcATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGaa                                                                          >  1:805379/1‑150 (MQ=255)
                                                                             ccgcATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGaa                                                                          >  1:4875202/1‑150 (MQ=255)
                                                                             ccgcAGCCGCCATTGTGTTAACCGTGTACAGAGACAAACCGGCGCTCGGTCGGCCCCCTCCCCCCTGGGGGGTGGGGTTTTGGGTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGaa                                                                          <  1:2937778/150‑1 (MQ=255)
                                                                              cgcATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAAt                                                                         >  2:2913517/1‑150 (MQ=255)
                                                                                cagcgggCAAAGTTTTACCCGTGTACAGCGACAACACGGCGCTGGGTCGGCCCCCGCGCCCCTCGGGGGTGTGGGTTTGGGTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGt                                                                       <  1:8115961/145‑1 (MQ=255)
                                                                                  tCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTc                                                                     >  2:1754768/1‑150 (MQ=255)
                                                                                   ccGGAAAAGTTTTACCCGTGTACAGCGTCAAAACGGCGCTCGGTCGGCCCCCTCGCCCCTCTGGGGTGAGGGTTAGGGTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTca                                                                    <  1:7491394/150‑1 (MQ=255)
                                                                                    cgtaaaaGTTTTAACCTTGTACAGAGTCAAACAGGCGCTCGGTCGGCCCCCGCGCCCCTCTGGGGTGTGGGTTGGGGTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcag                                                                   <  1:1562385/146‑1 (MQ=255)
                                                                                    cGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcag                                                                   >  2:8535894/1‑150 (MQ=255)
                                                                                    cGGCAAAGTGTTAACCGGTTCCCGAGCCAAAAAGGCGCGGGGGCGGGCCCCTCGCCCCTCGGGGGTGTGGGTTTGGGGGGGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcag                                                                   <  1:8637365/150‑1 (MQ=255)
                                                                                     ggCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagc                                                                  >  2:7552641/1‑150 (MQ=255)
                                                                                     ggCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagc                                                                  >  1:7295691/1‑150 (MQ=255)
                                                                                      gCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagca                                                                 <  2:5089224/150‑1 (MQ=255)
                                                                                       cAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagcag                                                                >  2:6900235/1‑150 (MQ=255)
                                                                                        aaaGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagcagc                                                               >  1:3401358/1‑150 (MQ=255)
                                                                                        aaaGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagcagc                                                               >  1:1337703/1‑150 (MQ=255)
                                                                                            ttttAACGGTGTACAGAGTCAAAATGGCGCGGGGTGGGTCCCCCCGCCCCTGGGGGGTGTGGTTTTGGGTGTGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTa                                                           <  1:8274780/148‑1 (MQ=255)
                                                                                              ttAACCGTTAACAGTGACAAACAGGCGCTGGGTCGGACCCCGCGCCCCTGTGGGTAGTGGTTTTGGGTGTGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACg                                                         <  1:770540/150‑1 (MQ=255)
                                                                                              ttAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACg                                                         >  1:8944089/1‑150 (MQ=255)
                                                                                                        aCAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTaa                                               >  1:6903414/1‑150 (MQ=255)
                                                                                                        aCAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTaa                                               <  1:2313493/150‑1 (MQ=255)
                                                                                                         cAGAGTAAAAACGGCGCTCGGTCGGGCCCCGCCCCCCTGTGGGGAGAGGGTTCGGGTGCGGGCAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAAt                                              <  1:7937362/150‑1 (MQ=255)
                                                                                                          attgTCAAACTGGCGCTCGGTGGGTCCCCTCGCCCCTGTGGGGTGAGGGTTCGGGTGTGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATg                                             <  1:3722776/147‑1 (MQ=255)
                                                                                                           gagTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGa                                            <  2:5335749/150‑1 (MQ=255)
                                                                                                             gTCAAACTGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATc                                          <  2:9014473/150‑1 (MQ=255)
                                                                                                                aaaaaGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAAt                                       >  1:9254999/1‑150 (MQ=255)
                                                                                                                     ggCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATAcccc                                  >  2:8203099/1‑150 (MQ=255)
                                                                                                                       cgcTGGGTCGGTCCCCTCGCCCCTGTGGGGAGTGGGTTTGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCgg                                <  1:4995600/150‑1 (MQ=255)
                                                                                                                        gcTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGc                               >  2:544097/1‑150 (MQ=255)
                                                                                                                        gcTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGc                               <  1:8573163/150‑1 (MQ=255)
                                                                                                                         cTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGTGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGcc                              >  2:5693897/1‑150 (MQ=255)
                                                                                                                          tcggtcggTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGATCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGccc                             >  1:8801494/1‑150 (MQ=255)
                                                                                                                          tcggtcggTCCCCGCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCCGAAGGGTCAGCAGCAGTACGCACAGGGGCATGATCCATACCCCGGccc                             >  1:6000297/1‑150 (MQ=38)
                                                                                                                            ggtcggtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAg                           <  2:121216/150‑1 (MQ=255)
                                                                                                                             gtcggtcCCCGCGCCCCTTTGGGTAGAGGTTTTGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAgg                          <  2:8627751/150‑1 (MQ=255)
                                                                                                                                   ccccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAAtt                    >  2:7614510/1‑150 (MQ=255)
                                                                                                                                    cccTCGCCCCTCTGGGGAGAGGGGTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtg                   >  2:2151427/1‑150 (MQ=255)
                                                                                                                                      cgcGCCCCTGTGGGGAGTGGTTTTGGGGGTGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgct                 <  1:4421704/148‑1 (MQ=255)
                                                                                                                                      cTCGCCCCTGTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgct                 <  2:4993897/150‑1 (MQ=255)
                                                                                                                                      cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgct                 <  2:3108272/150‑1 (MQ=255)
                                                                                                                                       tCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgctg                <  2:5867111/150‑1 (MQ=255)
                                                                                                                                         gCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAACAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGcc              >  2:4060247/1‑150 (MQ=255)
                                                                                                                                           cccagaggggCGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCAt            >  1:1019165/7‑150 (MQ=255)
                                                                                                                                           cccTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCAt            <  1:8986186/150‑1 (MQ=255)
                                                                                                                                             ctctGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTc          >  2:4198583/1‑150 (MQ=255)
                                                                                                                                              tcgggggAGAGGGTTCGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTcc         <  2:6513193/147‑1 (MQ=255)
                                                                                                                                               ggggggTGTGGTTTTGGGGGTGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCg        <  1:7061625/148‑1 (MQ=255)
                                                                                                                                               ctGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCg        >  1:3293893/1‑150 (MQ=255)
                                                                                                                                               ctGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCg        <  1:5097155/150‑1 (MQ=255)
                                                                                                                                                    gagtGGGTTCGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAAc   <  2:3781403/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAACG  >  NC_000913/3908363‑3908661

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: