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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,908,5140AT 8.9% 1897.8 / 29.5 649intergenic (‑148/+35)pstB/pstAphosphate ABC transporter ATPase/phosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (368/195);  new base T (2/75);  total (370/279)
Fisher's exact test for biased strand distribution p-value = 1.03e-28
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 679 total aligned reads displayed.

ATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAACGGG  >  NC_000913/3908365‑3908663
                                                                                                                                                     |                                                                                                                                                     
   gCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCTCCCCAGAGGGGCGAggg                                                                                                                                                    <  2:8800948/150‑1 (MQ=39)
    cAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGt                                                                                                                                                   <  1:981032/150‑1 (MQ=255)
      aTCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTa                                                                                                                                                 <  1:4224049/150‑1 (MQ=255)
       tCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAg                                                                                                                                                >  1:1206247/1‑150 (MQ=255)
       tCTTGTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAg                                                                                                                                                <  2:4648718/150‑1 (MQ=255)
        cTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAgg                                                                                                                                               <  2:3778295/150‑1 (MQ=255)
         ttttCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAggg                                                                                                                                              >  2:8219496/1‑150 (MQ=255)
          tttCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTCGGGt                                                                                                                                             >  1:3229426/1‑150 (MQ=255)
              cTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgg                                                                                                                                         <  1:4815397/150‑1 (MQ=255)
              cTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgg                                                                                                                                         >  1:810750/1‑150 (MQ=255)
               tCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAggg                                                                                                                                        >  2:1874664/1‑150 (MQ=255)
               tCATCGTGCATTATGCACTTCTCTAATCAGGCTTACAGTAGGCCTGATCAGTCGTTCACGCCGCATCCGTCAAAGTGTTAAACGTGTACAGAGTCAAAAAGGCTCTCTGTCTGTCCCCTCTCCCATCTGGTGAGATGGTTAgtgttattt                                                                                                                                        <  2:7054745/150‑6 (MQ=255)
                     tGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGCCAAAGTGTTAACCGAGTACAGAGTCAAAAAGGCGCTCGGTCGCTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAc                                                                                                                                  >  2:8871364/1‑150 (MQ=255)
                          gATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCgtgt                                                                                                                             >  1:304394/1‑150 (MQ=255)
                          gATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCgtgt                                                                                                                             >  1:6775214/1‑150 (MQ=255)
                             gCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAg                                                                                                                          <  2:5971014/150‑1 (MQ=255)
                             gCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAg                                                                                                                          >  1:9245668/1‑150 (MQ=255)
                             gCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAg                                                                                                                          >  1:8268571/1‑150 (MQ=255)
                               aCTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCa                                                                                                                        >  1:5668116/1‑150 (MQ=255)
                                cTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCaa                                                                                                                       >  2:7385617/1‑150 (MQ=255)
                                 ttCGCTAATCAGGCTTCCCGTCGGCCGGATCAGGCGTTCCCGCCGCTCCCGGCAAAGTGTTAACCGGGTACAGCGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAat                                                                                                                      <  2:1349736/150‑1 (MQ=38)
                                 ttCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAat                                                                                                                      >  2:6594596/1‑150 (MQ=255)
                                      tAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCaa                                                                                                                 >  1:3351433/1‑150 (MQ=255)
                                       aaTCAGGCTTACCGTCGGCCGTATCTGGCGTTCCGGCCGCATCCGGAAATTTTTTAACCGTTTACAGAGAAAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGTGGGTTCGGGTGCGGGGAAAACCGTGTCAGCAATATCAAc                                                                                                                <  2:2435963/150‑1 (MQ=25)
                                       aaTCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAAc                                                                                                                >  2:3146238/1‑150 (MQ=255)
                                             gCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGtt                                                                                                          >  2:4442388/1‑150 (MQ=255)
                                             gCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGtt                                                                                                          >  2:4018650/1‑150 (MQ=255)
                                             gCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGtt                                                                                                          >  2:1726306/1‑150 (MQ=255)
                                               ttACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTa                                                                                                        >  2:9286533/1‑150 (MQ=255)
                                                tACCGTCGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTACCCTCGCCCCTCTGGGGAGTGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAt                                                                                                       <  1:4599514/150‑1 (MQ=255)
                                                  ccGTAGGCCGGATCAGGCGTTCACGCCGCTTCCGGCAAATTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTTGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttc                                                                                                     <  2:557468/150‑1 (MQ=255)
                                                    gTAGGCCGGATCTGGCTTTCACGCCGCATCCGGCAAATTGTTAACCGTTTACAGAGTCAAACAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGTGCGGGTTCGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttctt                                                                                                   <  1:8777606/150‑1 (MQ=38)
                                                     tcggCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGGGTACAGAGTTAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCGCTGGGGAGAGGGTCAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttcttc                                                                                                  <  2:4768829/148‑1 (MQ=255)
                                                     tAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGTGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttcttc                                                                                                  <  1:2089484/150‑1 (MQ=255)
                                                     tAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttcttc                                                                                                  >  1:8907100/1‑150 (MQ=255)
                                                      cggCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGTGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCg                                                                                                 <  2:1206054/149‑1 (MQ=255)
                                                         ccGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCaa                                                                                              >  1:6642106/1‑150 (MQ=255)
                                                          cGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAACATCAACCGTGTTTATTCTTCGCaaa                                                                                             >  2:865087/1‑150 (MQ=255)
                                                            gATCTGGGGTCCCCGCCGCATCGGTCAATGTTTTAACGGTGTCCAGAGTCAAACTGGCGCTCGGTGGGTCCCCACGCCCCTCGGGGTATTGGGTTTGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCaaaaa                                                                                           <  2:3518740/150‑1 (MQ=255)
                                                            gATCAGGCGTTCACGCCGCATCCGGAAATGTGTTAACCGTGTACAGAGTCAAAACGGCGCTCGGTCGGCCCCCTCGCCCCTCTGGGGTGTGGGTTTGGGTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCaaaaa                                                                                           <  1:9332648/150‑1 (MQ=255)
                                                               cAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAaacaa                                                                                        <  1:7809691/150‑1 (MQ=255)
                                                                 ggCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcg                                                                                      >  2:6139113/1‑150 (MQ=255)
                                                                  gCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgc                                                                                     >  2:674138/1‑150 (MQ=255)
                                                                  gCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgc                                                                                     >  1:5455165/1‑150 (MQ=255)
                                                                  gCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgc                                                                                     >  2:7935528/1‑150 (MQ=255)
                                                                   cgtacccgccgcATCCGGCATTGTGTTAACCGTGTACAGCGTAAAACCGGCGCTGGGTGGGTCCCCTCGCCCCTCTGGGGTGCGGTTTAGGGTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgcg                                                                                    <  1:5305013/144‑1 (MQ=255)
                                                                    gtgcccgccgcCTCCGGCAACGTGTTAACGGTGTCCAGAGTCAAAACGGCGCTGGGTGGGTCCCCGCGCCCCTCTGGGGTGTGGGTTGGGGGGTGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgcgc                                                                                   <  2:3234814/145‑1 (MQ=255)
                                                                    gcaccccccgcAGCCGGAAATTTTTTAACGGTGTACTGAGAAAAAACGGCGGGGGGTGGGCCCCCTCGCCCGTGTGGGTTGTGGGTTTGGGTGGGGGCAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgcgc                                                                                   <  1:2732603/143‑1 (MQ=255)
                                                                           ccgcATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGTTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGaa                                                                            <  2:841613/150‑1 (MQ=255)
                                                                           ccgcATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGaa                                                                            >  2:9215492/1‑150 (MQ=255)
                                                                           ccgcATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGaa                                                                            >  2:7030936/1‑150 (MQ=255)
                                                                           ccgcATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGaa                                                                            >  2:2278257/1‑150 (MQ=255)
                                                                           ccgcATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGaa                                                                            >  1:805379/1‑150 (MQ=255)
                                                                           ccgcAGCCGCCATTGTGTTAACCGTGTACAGAGACAAACCGGCGCTCGGTCGGCCCCCTCCCCCCTGGGGGGTGGGGTTTTGGGTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGaa                                                                            <  1:2937778/150‑1 (MQ=255)
                                                                              cacccGTCAATGTTTTACCCGTGTACTGAGACAAAACGGCGCGCGGTCGGTCCCCCCGCCCCTGTGGGGTGTGGGTTTGGGTTCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGt                                                                         <  1:396941/147‑1 (MQ=255)
                                                                              cATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGt                                                                         >  2:2857031/1‑150 (MQ=255)
                                                                               aTCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGtt                                                                        >  1:6290444/1‑150 (MQ=255)
                                                                                tCCGGAAAAGTGTTAACCGTTTACAGAGTCAAAAAGGCGCTCGGTCGGCCCCCTCGCCCCTCTGGGGTGAGGGTTAGGGTTAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTc                                                                       <  2:211677/150‑1 (MQ=255)
                                                                                 ccGGAAAAGTTTTACCCGTGTACAGCGTCAAAACGGCGCTCGGTCGGCCCCCTCGCCCCTCTGGGGTGAGGGTTAGGGTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTca                                                                      <  1:7491394/150‑1 (MQ=255)
                                                                                  cGGCAAAGTGTTAACCGGTTCCCGAGCCAAAAAGGCGCGGGGGCGGGCCCCTCGCCCCTCGGGGGTGTGGGTTTGGGGGGGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcag                                                                     <  1:8637365/150‑1 (MQ=255)
                                                                                   ggCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagc                                                                    >  1:7295691/1‑150 (MQ=255)
                                                                                   ggCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagc                                                                    >  2:8040543/1‑150 (MQ=255)
                                                                                    gCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagca                                                                   >  2:570345/1‑150 (MQ=255)
                                                                                    gCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagca                                                                   <  2:5089224/150‑1 (MQ=255)
                                                                                     cAAAGTGATAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCACCTCGCCCCTCTCGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAATCTTCGCAAAAACAACGCGCGCCAGAATGATcagcag                                                                  >  2:6900695/1‑150 (MQ=39)
                                                                                      aaaGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagcagc                                                                 >  2:7654124/1‑150 (MQ=255)
                                                                                      aaaGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagcagc                                                                 >  1:3720028/1‑150 (MQ=255)
                                                                                          tgtTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTa                                                             >  1:1618686/1‑150 (MQ=255)
                                                                                            ttAACCGTTAACAGTGACAAACAGGCGCTGGGTCGGACCCCGCGCCCCTGTGGGTAGTGGTTTTGGGTGTGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACg                                                           <  1:770540/150‑1 (MQ=255)
                                                                                            ttAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACg                                                           >  1:8976959/1‑150 (MQ=255)
                                                                                                  gtgtACAGAGTCAAAAAGGCACTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAgg                                                     >  2:7567747/1‑150 (MQ=255)
                                                                                                    gtACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGt                                                   >  2:9083122/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)
                                                                                                         gagTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGa                                              <  2:5335749/150‑1 (MQ=255)
                                                                                                          agTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTACAGCAGCAGTACGCACAGGGTAATGAt                                             >  1:5062590/1‑150 (MQ=255)
                                                                                                           gTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATc                                            >  2:3028401/1‑150 (MQ=255)
                                                                                                             cAAAACGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGCGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCaa                                          <  2:4206329/150‑1 (MQ=255)
                                                                                                              aaaaaGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAAt                                         >  2:4994562/1‑150 (MQ=255)
                                                                                                              aaaaaGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAAt                                         >  2:1470053/1‑150 (MQ=255)
                                                                                                                 aaGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATAcc                                      >  1:4407016/1‑150 (MQ=255)
                                                                                                                      gcTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGc                                 <  1:8573163/150‑1 (MQ=255)
                                                                                                                        tcggtcggTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGccc                               >  2:2655312/1‑150 (MQ=255)
                                                                                                                        tcggtcggTCCCCGCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCCGAAGGGTCAGCAGCAGTACGCACAGGGGCATGATCCATACCCCGGccc                               >  1:6000297/1‑150 (MQ=38)
                                                                                                                          ggtcggtcCCCTCGCCCCTCTGGGGTGTGGGTTTGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAg                             <  2:8549908/150‑1 (MQ=255)
                                                                                                                            gcggtcCCCTCGCCCCTCTGGGGTGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGc                           <  1:4925946/149‑1 (MQ=255)
                                                                                                                                 ccccTCGCCCCTCTGGGGAGAGGGTTCGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAAtt                      <  1:6440125/150‑1 (MQ=255)
                                                                                                                                    cTCGCCCCTCTGGGGTGTGGGTTCGGGTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgct                   <  1:173218/150‑1 (MQ=255)
                                                                                                                                    cTCGCCCCTCTGGGGAGAGGTTTTGGGTTTGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgct                   <  1:3533554/150‑1 (MQ=255)
                                                                                                                                     tCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgctg                  <  2:5867111/150‑1 (MQ=255)
                                                                                                                                         cccagaggggCGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCAt              >  1:1019165/7‑150 (MQ=255)
                                                                                                                                         cccTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCAt              >  1:3813879/1‑150 (MQ=255)
                                                                                                                                            tctGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATg                                                                          <  1:1355849/87‑1 (MQ=255)
                                                                                                                                            tcgggggAGAGGGTTCGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTcc           <  2:6513193/147‑1 (MQ=255)
                                                                                                                                             ggggggTGTGGTTTTGGGGGTGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCg          <  1:7061625/148‑1 (MQ=255)
                                                                                                                                                 ggAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCaaa      <  2:5455165/150‑1 (MQ=255)
                                                                                                                                                  gagaGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAAc     >  2:3302204/1‑150 (MQ=255)
                                                                                                                                                   agagGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAACg    >  2:2580288/1‑150 (MQ=255)
                                                                                                                                                    gagGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGc                           >  1:1864348/1‑126 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAACGGG  >  NC_000913/3908365‑3908663

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: