Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,908,506 C→G 5.8% intergenic (‑140/+43) 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,5060CG 5.8% 1685.0 / 29.8 503intergenic (‑140/+43)pstB/pstAphosphate ABC transporter ATPase/phosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (276/191);  new base G (5/30);  total (281/222)
Fisher's exact test for biased strand distribution p-value = 1.87e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 531 total aligned reads displayed.

CATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCG  >  NC_000913/3908358‑3908655
                                                                                                                                                    |                                                                                                                                                     
cATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCctct                                                                                                                                                      <  2:3985304/150‑1 (MQ=255)
cATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCctct                                                                                                                                                      <  1:2722272/150‑1 (MQ=255)
 aTACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTg                                                                                                                                                     >  2:5317221/1‑150 (MQ=255)
  tACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTgg                                                                                                                                                    <  1:1139659/150‑1 (MQ=255)
  tACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTgg                                                                                                                                                    <  2:2495236/150‑1 (MQ=255)
  tACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTgg                                                                                                                                                    <  2:5030049/150‑1 (MQ=255)
     tCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGga                                                                                                                                                 <  1:5379560/150‑1 (MQ=255)
       aTTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGgaga                                                                                                                                               <  2:115931/150‑1 (MQ=255)
            aaTCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCTCCTCTCCCCAGAGGGGCGAgggga                                                                                                                                          <  1:5663374/150‑3 (MQ=39)
            aaTCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGtt                                                                                                                                          >  1:2290484/1‑150 (MQ=255)
             aTCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTa                                                                                                                                         <  2:4545882/150‑1 (MQ=255)
              tcatttCCTCCTCGTGCCTGATGCTCTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCACAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCGGGGGAGAGGGTTAg                                                                                                                                        <  2:6235729/147‑1 (MQ=39)
              tCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAg                                                                                                                                        >  1:4281248/1‑150 (MQ=255)
                ttttCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCTGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAggg                                                                                                                                      >  1:1836375/1‑150 (MQ=255)
                   tCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGa                                                                                                                                   <  2:4881311/150‑1 (MQ=255)
                    ccTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgg                                                                                                                                 >  2:2862486/1‑150 (MQ=255)
                     cTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgg                                                                                                                                 <  1:2939287/150‑1 (MQ=255)
                     cTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgg                                                                                                                                 <  2:1243708/150‑1 (MQ=255)
                     cTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgg                                                                                                                                 <  2:4053092/150‑1 (MQ=255)
                      tCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAggg                                                                                                                                >  2:2021341/1‑150 (MQ=255)
                             gCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGCGTCAAAAAGGCGCTCGGTCGGGCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAcc                                                                                                                         <  2:3662109/150‑1 (MQ=255)
                               cTGATGCACTTCGCTAATCCGGCTTACGTTTGGCCGGATCTGGCGTTCACGCCGCATCCGGAATTGTTTTAACCGTGTAAAGATTAAAAATGGCGCTGGGTCGGTCCCCTCGCCCCTCTGGGTTGTGGGTTTGGGTGAGGGGAAAACCgt                                                                                                                       <  2:2421010/150‑1 (MQ=255)
                               cTGATGCACTTCGATAATCAGGCTTACCTTAGGCCGGTTAAGGCGTTCACGCCGCATCCGGCAAAGTTTTAACCGTGTACAGAGTAAAAATGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGTGAGGTTTAGGGTGAGGGGAAAACCgt                                                                                                                       <  2:5480177/150‑1 (MQ=38)
                                 gATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCgtgt                                                                                                                     >  1:905214/1‑150 (MQ=255)
                                    gCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGTGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAg                                                                                                                  >  1:3429221/1‑150 (MQ=255)
                                      aCTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCa                                                                                                                >  1:2610279/1‑150 (MQ=255)
                                      aCTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCa                                                                                                                >  1:4751464/1‑150 (MQ=255)
                                      aCTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCa                                                                                                                >  1:4546115/1‑150 (MQ=255)
                                      aCTTCGATAATCTGGCTTACCGTAGGCCGTATCTGGCGTTCACGCCGCATCCGGAAATGTTTTAACCGTTTACAGTGTAAAAATGGCGCTCGGTGGGTCCCCTCGCCCCTCTGGGTTGTGGTTTTGGTTGTGGGGAAAACCGTGTCAGCa                                                                                                                <  1:6195308/150‑1 (MQ=255)
                                        ttCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAat                                                                                                              >  1:1087771/1‑150 (MQ=255)
                                        ttCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAat                                                                                                              >  2:6622562/1‑150 (MQ=255)
                                           gCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                              >  1:2108568/1‑115 (MQ=255)
                                           gCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                              <  2:2108568/115‑1 (MQ=255)
                                           gCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCCGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATc                                                                                                           <  2:1503313/150‑1 (MQ=255)
                                            cTAATCCGGCTTACCGTCGGCCGGATCAGGCGTTCCCGCCGCATCCGGCAAAGTTTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCGGGGGAGTGGGTTCGGGTGTGGGGAAAACCGTGTCAGCAATATCa                                                                                                          <  2:2610464/150‑1 (MQ=255)
                                            cTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCa                                                                                                          >  2:5961994/1‑150 (MQ=255)
                                            cTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCa                                                                                                          >  1:2362319/1‑150 (MQ=255)
                                             tAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCaa                                                                                                         >  1:5764565/1‑150 (MQ=255)
                                             tAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGACAACCGTGTCAGCAATATCaa                                                                                                         >  1:5206646/1‑150 (MQ=255)
                                              aaTCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAAc                                                                                                        >  1:3527218/1‑150 (MQ=255)
                                               aTCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAAcc                                                                                                       >  1:4638837/1‑150 (MQ=255)
                                                 cAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCgt                                                                                                     >  2:737286/1‑150 (MQ=255)
                                                  aGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCgtg                                                                                                    >  1:3185012/1‑150 (MQ=255)
                                                    gCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGTGGGTTAGGGTGTGGGGAAAACCGTGTCAGCAATATCAACCGTGtt                                                                                                  <  1:6105069/150‑1 (MQ=255)
                                                    gCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGtt                                                                                                  >  1:3824729/1‑150 (MQ=255)
                                                     cTTCCCGTAGGCCGGATCAGGCTTTCACGCCGCATCCGGCAAAGTGTTACCCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTTGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGttt                                                                                                 <  2:539889/150‑1 (MQ=255)
                                                         ccGTAGGCCGTATCTGGCGTTCACGCCGCATCCGGCAAAGTTTTAACCGTGTACAGATTCAAAAAGGCGCTCGGTGGGTCCCCTCGCCCCTGTGGGTAGTGGTTTAGGGTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttc                                                                                             <  2:13705/150‑1 (MQ=255)
                                                           gTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttctt                                                                                           <  1:6643647/150‑1 (MQ=255)
                                                            tttgCCGGATGCGGCGTGAACGCCTGATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTa                                                                                                                                         >  2:2229842/4‑103 (MQ=37)
                                                            tAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttcttc                                                                                          >  1:1344328/1‑150 (MQ=255)
                                                            tAGGCCGGAGCTGGCGGTCAGGCCGCTCCCGGCAATGTGTTAAGGGTGTACAGTGAAAAAACGGCGCTCGGGCGGTCCCCTGGCCCCTGTGGGGTGTGGGTTAGGGTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttcttc                                                                                          <  2:5113990/150‑1 (MQ=255)
                                                            gcggCCGGATCTGGCGTTCACGCCGCAGCCGGCATTGTGTTAACCGTGTACAGTGGCAAAAAGGCGCGGGGTCGGTCCCCTCGCCCCTGTGGGGAGGGGGTTCGGGTGGGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttcttc                                                                                          <  2:1304265/148‑1 (MQ=255)
                                                             aGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCg                                                                                         <  1:6530412/150‑1 (MQ=255)
                                                              ggCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGc                                                                                        <  2:1087771/150‑1 (MQ=255)
                                                               gCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCa                                                                                       >  1:2043313/1‑150 (MQ=255)
                                                               gCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCa                                                                                       >  1:2150252/1‑150 (MQ=255)
                                                               gCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCa                                                                                       >  1:694871/1‑150 (MQ=255)
                                                                 cGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCaaa                                                                                     >  2:1687998/1‑150 (MQ=255)
                                                                     tCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGGTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAaaca                                                                                 >  2:4224184/1‑150 (MQ=255)
                                                                       aGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAaacaac                                                                               >  2:184938/1‑150 (MQ=255)
                                                                       aGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGGCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAaacaac                                                                               <  2:6137362/150‑1 (MQ=255)
                                                                        ggCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAATCAAcg                                                                              >  1:4880085/1‑150 (MQ=255)
                                                                        ggCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcg                                                                              >  2:522424/1‑150 (MQ=255)
                                                                        ggCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcg                                                                              >  1:5873462/1‑150 (MQ=255)
                                                                        ggCGTTCACGCCGCATCCGGCAAAGTGTTAACAGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcg                                                                              >  2:5557693/1‑150 (MQ=255)
                                                                         gCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgc                                                                             <  2:5648402/150‑1 (MQ=255)
                                                                         gCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGCCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAAAATCAACCGTGTTTATTCTTCGCAAAAACAAcgc                                                                             >  2:4768739/1‑150 (MQ=255)
                                                                            ttCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgcgcg                                                                          >  2:4455763/1‑150 (MQ=255)
                                                                               aCGCCGCATCCGTATAAGTTTTAACCGTGTACATAGAAAAAACGGCGCTGGGTCGGTCCCCTCGCCCCTCTGGGGTGAGGTTTAGGTTGTGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCa                                                                       <  2:6348507/150‑1 (MQ=255)
                                                                                 gccgcATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGa                                                                     >  2:4623603/1‑150 (MQ=255)
                                                                                  ccgcATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGaa                                                                    >  2:1883377/1‑150 (MQ=255)
                                                                                   cgcATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAAt                                                                   >  2:2990107/1‑150 (MQ=255)
                                                                                      aTCCGGCAAAGTTTTAACCGTGTACAGAGTCAAAACGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTCGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGtt                                                                <  2:6399213/150‑1 (MQ=255)
                                                                                      aTCCGGCAAAGTGTTAACCGTGTACTGAGTCAAAATGGCGCTCGGTCGGCCCCCTCGCCCCTCGGGGTAGAGGTTTCGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGtt                                                                <  2:268004/150‑1 (MQ=255)
                                                                                         cGGCAACGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTTGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcag                                                             <  1:5964157/150‑1 (MQ=255)
                                                                                          ggCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagc                                                            >  1:1008657/1‑150 (MQ=255)
                                                                                            cAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagcag                                                          >  2:4132790/1‑150 (MQ=255)
                                                                                              attttttAACGGTGTACTTATTAAAAACGGCGCTGGGTGGGTCCCCTCGCCCCTCTGGGGTGAGGGTTTGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagcagca                                                        <  1:169129/145‑1 (MQ=37)
                                                                                              aaGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagcagca                                                        >  2:3799209/1‑150 (MQ=255)
                                                                                                     aaCCGTGTACAGAGTCAAAACGGCGCTCGGTCGGCCCCCTCGCCCCTCTGGGGAGTGGGTTTGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGca                                                 <  1:4381165/150‑1 (MQ=255)
                                                                                                     aaCCGTGTAAAGAGTCAAAATGGCGCTCGGTCGGTCCCCTCGCCCGTCTGGGGTGTGGTTTCGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGca                                                 <  1:5480141/150‑1 (MQ=255)
                                                                                                       ccGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTa                                                                                                <  2:6703825/101‑1 (MQ=255)
                                                                                                         gtttACTGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGTATTGGGTTTGGGTTAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAgg                                             <  1:3691179/147‑1 (MQ=255)
                                                                                                         gtgtACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGTAGTACGCACAgg                                             >  2:3313159/1‑150 (MQ=255)
                                                                                                         gtgtACAGAGTCAAAAAGGCGCTCGGTCGGTCACCTCGCCCCTCTGGGGAGATGGTTAGGGTGAGTGGAAAACCGTGTCAGCAATATCAACCGTTTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGcacatg                                             >  1:5113007/1‑148 (MQ=255)
                                                                                                         gtgtACAGAGACAAATCGGCGCGCGGTCGGCCCCCTCGCCCCTCTGGGGAGTGGGTTCGGGTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAgg                                             <  1:3077583/150‑1 (MQ=38)
                                                                                                                 agTCAAAATGGCGCTCGGTGGGCCCCCGCGCCCCTCTGGGGAGTGGTTTCGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGAt                                     <  1:5915758/150‑1 (MQ=255)
                                                                                                                     aaaaaGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAAt                                 >  2:857210/1‑150 (MQ=255)
                                                                                                                     aaaaaGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAAt                                 >  2:350292/1‑150 (MQ=255)
                                                                                                                           gcgcTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGTGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCg                           <  2:3891208/150‑1 (MQ=255)
                                                                                                                            cgcTCGGTCGGTCCCCTCGCCCCTCTGGGGAGTGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCgg                          <  2:196941/150‑1 (MQ=255)
                                                                                                                                         cccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtg             >  1:6406775/1‑150 (MQ=255)
                                                                                                                                          ccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTa                                          >  1:3361042/1‑120 (MQ=255)
                                                                                                                                           cTCGCCCCTCTGGGGTGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgct           <  2:1817144/150‑1 (MQ=255)
                                                                                                                                            tCGCCCCTCTGGGGAGAGGGTTCGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgctg          <  1:3896804/150‑1 (MQ=255)
                                                                                                                                                 ccTCGGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCAtt     <  2:5741145/150‑1 (MQ=255)
                                                                                                                                                  ctctGGGGAGTGGGTTCGGGTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTc    <  2:5873400/150‑1 (MQ=255)
                                                                                                                                                  ctctGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTc    >  2:4290095/1‑150 (MQ=255)
                                                                                                                                                   tctGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTcc   <  1:1533726/150‑1 (MQ=255)
                                                                                                                                                   tctGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAAATGCTGCCATTcc   >  1:4097809/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCG  >  NC_000913/3908358‑3908655

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: