Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,908,507 T→G 6.6% intergenic (‑141/+42) pstB ← / ← pstA phosphate ABC transporter ATPase/phosphate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,908,5070TG 6.6% 473.2 / 5.4 173intergenic (‑141/+42)pstB/pstAphosphate ABC transporter ATPase/phosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (80/77);  new base G (2/12);  total (82/91)
Fisher's exact test for biased strand distribution p-value = 1.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 178 total aligned reads displayed.

ATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCG  >  NC_000913/3908359‑3908655
                                                                                                                                                    |                                                                                                                                                    
aTACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTg                                                                                                                                                     <  2:757955/150‑1 (MQ=255)
aTACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTg                                                                                                                                                     >  2:2294386/1‑150 (MQ=255)
aTACTCATTGCAATCGTGTCCTCATCGTGCCTGATGCACTTCGCTAATCTGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAATGGCGCTCGGTCGGTCCCCTCGCCCCTCTg                                                                                                                                                     <  2:1673933/150‑1 (MQ=255)
 tACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTgg                                                                                                                                                    <  1:531897/150‑1 (MQ=255)
  aCTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTggg                                                                                                                                                   <  2:1524058/150‑1 (MQ=255)
  aCTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTggg                                                                                                                                                   <  1:368237/150‑1 (MQ=255)
    tCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGga                                                                                                                                                 >  1:1269352/1‑150 (MQ=255)
    tCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGga                                                                                                                                                 >  2:1699763/1‑150 (MQ=255)
    tCATTGCAATCTTTTCCTCATCGTGCCTGAGGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGga                                                                                                                                                 >  2:1216549/1‑150 (MQ=255)
      aTTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTCGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGgaga                                                                                                                                               <  2:877229/150‑1 (MQ=255)
      aTTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGgaga                                                                                                                                               <  1:654454/150‑1 (MQ=255)
       ttGCAATCTTTTCCTCCTCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                              <  2:1642859/150‑1 (MQ=255)
       ttGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                              <  2:415643/150‑1 (MQ=255)
         gCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAggg                                                                                                                                            <  2:604929/150‑1 (MQ=255)
             tcgtttCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAg                                                                                                                                        <  1:82049/147‑1 (MQ=255)
             tCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAg                                                                                                                                        <  2:83259/150‑1 (MQ=255)
               ttttCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTCCAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAggg                                                                                                                                      >  2:1414105/1‑150 (MQ=255)
               ttttCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAggg                                                                                                                                      >  1:259034/1‑150 (MQ=255)
               ttttCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAggg                                                                                                                                      >  2:232136/1‑150 (MQ=255)
                  tCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAGAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGa                                                                                                                                   >  1:1669876/1‑150 (MQ=255)
                     tCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAACAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGCggg                                                                                                                                <  2:1186970/150‑1 (MQ=255)
                       aTCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGa                                                                                                                              >  2:2145753/1‑150 (MQ=255)
                           ggCCTGATGCACTTCGCTAATCAGGCTTACGGTAGGCCGGATCAGGCGTTCACGCCGCATCGGGCAATGTGTTAACCGTGTAAAGTGTAAAAATGGCGCTCGGTCGTTCCCCTCGCCCCTCTGGGGTGTGGGTTAGGTTGAGGGGAAAAc                                                                                                                          <  2:460597/149‑1 (MQ=255)
                            gCCTTATGCACTTCGCTAATCTGGCTTACGGTCGGCGGTATCTGGCTTTCACGCCGCATCCGGCAATGTTTTAACCGTGTACTGATTCAAAATGGCGCTGGGTGGGTCCCCTCGCCCCTCTGGGGTTTGGGTTTGGGTGTGGGGAAAAcc                                                                                                                         <  1:1758098/150‑1 (MQ=255)
                               tGATGCCCTACGCTATGCCGGCTTACGGTCGCCCGGATCTGGCTTACCCGCCGCACCGGGCAATGTGTTAACGGTTTACTGAGTCAAAAAGGGGCCGGGTGGGTCCCCTCGCCCCTGGGGGGTGTGGGTTAGGGTGAGGGGAAAACCgtg                                                                                                                      <  2:2042622/150‑1 (MQ=255)
                               tGATGCACTTCGCGAATCAGGCTTCCCGTAGGCCGGAGCAGGCGTTCACGCCGCAGCCGGCAATGTGTTAACGGTTTACAGAGTAAAAAAGGCGCTCGGTGGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCgtg                                                                                                                      <  2:1955836/150‑1 (MQ=37)
                               tGAGGCAGTTCGCTAATCCGGCTTACCGTTGGCCGGATCCGGGTTGCCCGCCGCATCGGGAAACTTGTTAACCGTGGACAGTGTAAAAATGGCGCTGGGGCGGTCCCCTCGCCCCTCTGGGGTGTGGTTTTGGGGGTGGGGAAAACCgtg                                                                                                                      <  2:875913/150‑1 (MQ=255)
                               tGAGGCACTTCGCTACTCCGGCTTCCCGTCGGCCGGATCCGGCGTTCACGCCGCAACCGGCAAAGTGTTAACCGTGTACAGCGGCAAACAGGCGCTCGGGCGGCCCCCTCGCCCCTCTGGGGAGTGGGTTAGGGTGTGGGGAAAACCgtg                                                                                                                      <  2:1729336/150‑1 (MQ=255)
                                 aTGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTCGGGTGAGGGGCAAACCGTGTc                                                                                                                    >  2:1848758/1‑150 (MQ=255)
                                 aTGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTc                                                                                                                    >  2:1446708/1‑150 (MQ=255)
                                 aTGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTc                                                                                                                    >  2:1446611/1‑150 (MQ=255)
                                 aTGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTc                                                                                                                    >  2:207751/1‑150 (MQ=255)
                                 aTGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTc                                                                                                                    >  2:52554/1‑150 (MQ=255)
                                     aCTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCa                                                                                                                >  2:95388/1‑150 (MQ=255)
                                      ctccGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGGCCCCTCGCCCCTCTGGGGAGTGGGTTAGGGTGAGGGGAAAACCGTGTCAGCaa                                                                                                               <  1:1679065/147‑1 (MQ=255)
                                          gCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATc                                                                                                           >  1:198971/1‑150 (MQ=255)
                                             aaTCTGGCTTACCGTAGGCCGTATCAGGCGTTCACGCCGCATCCGTCAAAGTGTTAACCGTGTAAAGAGTCAAAAAGGCGCTCGTTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGTGGGGAAAACCGTGTCAGCAATATCAAc                                                                                                        <  2:171172/150‑1 (MQ=38)
                                             aaTCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAAc                                                                                                        >  2:1831983/1‑150 (MQ=255)
                                             aaTCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCCCCCCTCTGGGGATAGGTTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAAc                                                                                                        <  2:171196/150‑1 (MQ=255)
                                              aTCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAAcc                                                                                                       >  2:545783/1‑150 (MQ=255)
                                                cAGGCTTACCGTTGGCCGGATCTGGCGTTCACGCCGCATCCGGAAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCgt                                                                                                     <  2:198971/150‑1 (MQ=255)
                                                cAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCgt                                                                                                     >  2:2362132/1‑150 (MQ=255)
                                                    cTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGttt                                                                                                 >  2:257189/1‑150 (MQ=255)
                                                    cTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGttt                                                                                                 >  2:642231/1‑150 (MQ=255)
                                                    cTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGttt                                                                                                 >  2:642282/1‑150 (MQ=255)
                                                     ttACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTa                                                                                                >  1:477375/1‑150 (MQ=255)
                                                       aCCGTAGGCCGGATCAGGCTTTCACGCCGCATCCGCCAATGTGTTAACCGTGTACAGAGACAAAACGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGTGGTTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAtt                                                                                              <  2:1091244/150‑1 (MQ=38)
                                                        ccGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttc                                                                                             >  2:210508/1‑150 (MQ=255)
                                                           tAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAttcttc                                                                                          <  2:1965455/150‑1 (MQ=255)
                                                             ggCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAACAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGc                                                                                        <  2:1293789/150‑1 (MQ=255)
                                                             ggCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGTTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGc                                                                                        <  2:800050/150‑1 (MQ=255)
                                                              gCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCa                                                                                       >  2:397151/1‑150 (MQ=255)
                                                                      aGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAaacaac                                                                               >  1:1089476/1‑150 (MQ=255)
                                                                          gtgcACGCCGCCTCCGGCAAAGTGTTAACCGGGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGTGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgcgc                                                                           <  2:1960562/147‑1 (MQ=255)
                                                                          gTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgcgc                                                                           <  2:972074/150‑1 (MQ=255)
                                                                          gTTCACGCCGCATCCGGAAAAGTGTTAACCGGGTACAGCGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTGGGGGGTTTGGTTTAGGGTGAGGGGATAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgcgc                                                                           <  2:1161192/150‑1 (MQ=255)
                                                                                gccgcATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGa                                                                     >  1:235160/1‑150 (MQ=255)
                                                                                 ccgcATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGaa                                                                    >  1:1878680/1‑150 (MQ=255)
                                                                                     aTCCGGCAAATTTTTAACCGTTTACAGAGTCAAACAGGCGCTCGGTCGGTCCCCTCGCCCCTGGGGGGAGAGGGTTAGGGTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGtt                                                                <  2:1627680/150‑1 (MQ=255)
                                                                                     aTCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAATGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGTGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGtt                                                                <  1:777443/150‑1 (MQ=255)
                                                                                      tCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTc                                                               >  1:1054606/1‑150 (MQ=255)
                                                                                       ccGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTca                                                              <  2:593005/150‑1 (MQ=255)
                                                                                         ggCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagc                                                            >  2:1654413/1‑150 (MQ=255)
                                                                                             atgtgtTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTGTGGGGTGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagcagca                                                        <  1:1492627/148‑1 (MQ=255)
                                                                                             aattttaaaCCGTGTACAGACACAAAAAGGCGCTCGGTCGCACCCCTCCCCCCTCTGGTCATAGGTTTAGGGTGAGGGCAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagcagca                                                        <  1:1878308/143‑1 (MQ=255)
                                                                                                 gtTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTAc                                                    >  1:1919849/1‑150 (MQ=255)
                                                                                                       cGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAg                                              >  1:2387610/1‑150 (MQ=255)
                                                                                                        ttgtACAGAGTAAAATTGGCGCTGGGTGGGTCCCCTCGCCCTTTTGGTTTGTGGTTTAGGTTTTGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAgg                                             <  1:1699763/149‑1 (MQ=255)
                                                                                                            aCAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGTGGGTTTGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTaa                                         <  2:2279111/150‑1 (MQ=255)
                                                                                                             cAGAGTCAAACCGGCGCTCGGTCGGTCCCCTCGCCCCTCGGGGGAGTGGGTTTGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAAt                                        <  1:1874604/150‑1 (MQ=255)
                                                                                                             cAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAAt                                        <  1:222278/150‑1 (MQ=255)
                                                                                                             cAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAAt                                        >  1:833266/1‑150 (MQ=255)
                                                                                                                  tCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGTGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCa                                   <  2:1922038/150‑1 (MQ=255)
                                                                                                                    aaaaaGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAAt                                 >  1:1233445/1‑150 (MQ=255)
                                                                                                                    aaaaaGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAAt                                 >  2:1754710/1‑150 (MQ=255)
                                                                                                                         ggCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGTTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATAcccc                            <  2:2092911/150‑1 (MQ=255)
                                                                                                                         ggCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATAcccc                            >  1:1033751/1‑150 (MQ=255)
                                                                                                                         ggCGCTCGGGCGGACCCCTCGCCCCGCGGGGGAGCGGGTGAGGGGGAGGGGAAAACCGTGTCAGCAATATCCCCCGTGTTTATTCTTCGCAACAACAACGCGCGCCTGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATCcccc                            <  2:1478274/150‑1 (MQ=255)
                                                                                                                               cggtcggtCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCa                      >  1:1312759/1‑150 (MQ=255)
                                                                                                                                 gtcggtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAgg                    <  1:1755903/150‑1 (MQ=255)
                                                                                                                                     gtcCCCTCGCCCCTCTGGGGTGTGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCaa                <  1:1962750/150‑1 (MQ=255)
                                                                                                                                     gtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCCACAACACCGCGCGCCAGACTGGTCAGCAGCAGTACGCACCGGGTACTGATCAATACCCCGGCCCAGGCCaa                >  2:1353761/1‑150 (MQ=38)
                                                                                                                                     gtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCaa                >  2:1829306/1‑150 (MQ=255)
                                                                                                                                     gtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCaa                >  1:428290/1‑150 (MQ=255)
                                                                                                                                     gtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCaa                <  1:1188088/150‑1 (MQ=255)
                                                                                                                                     gtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCaa                >  1:1349436/1‑150 (MQ=255)
                                                                                                                                     ggccccTCGCCCCGGTGGGGAGTGGGTTCGGGTGTGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCaa                <  2:705338/148‑1 (MQ=255)
                                                                                                                                         ccTCGCCCCTCTGGGTAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgc            <  1:953804/150‑1 (MQ=255)
                                                                                                                                          cTCGCCCCTCTGGGGAGTGGTTTCGGGTGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgct           <  1:1037881/150‑1 (MQ=255)
                                                                                                                                          cTCGCCCCTCTGGGGAGCGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgct           <  1:1964809/150‑1 (MQ=255)
                                                                                                                                          cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTACAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgct           <  2:1268547/150‑1 (MQ=255)
                                                                                                                                            cGCCCCTCTGGGGTGTGGGTTTGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgctgc         <  2:259034/150‑1 (MQ=255)
                                                                                                                                             gCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGcc        >  2:2185932/1‑150 (MQ=255)
                                                                                                                                             gCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGcc        >  2:1309195/1‑150 (MQ=255)
                                                                                                                                              ccccTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAAAACCCCGGCCCAGGCCAATTGCTGCCa       >  2:1148765/1‑150 (MQ=255)
                                                                                                                                                 ctctGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTc    >  1:1444531/1‑150 (MQ=255)
                                                                                                                                                 ctctGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTc    >  1:1444616/1‑150 (MQ=255)
                                                                                                                                                   ctGGGGTTAGGTTTAGGGTGTGGGGAAATCCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCg  <  1:1935799/150‑1 (MQ=255)
                                                                                                                                                   ctGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCg  >  1:1490691/1‑150 (MQ=255)
                                                                                                                                                   ctGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCg  >  2:1005249/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCG  >  NC_000913/3908359‑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: