Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,908,505 T→A 7.9% intergenic (‑139/+44) 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,5050TA 7.9% 236.1 / 7.9 95intergenic (‑139/+44)pstB/pstAphosphate ABC transporter ATPase/phosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (37/49);  new base A (4/4);  total (41/54)
Fisher's exact test for biased strand distribution p-value = 7.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

TACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCC  >  NC_000913/3908360‑3908654
                                                                                                                                                 |                                                                                                                                                     
tACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTgg                                                                                                                                                   <  2:1478747/150‑1 (MQ=255)
tACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTgg                                                                                                                                                   <  1:1328513/150‑1 (MQ=255)
 aCTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTggg                                                                                                                                                  <  2:127869/150‑1 (MQ=255)
 aCTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTggg                                                                                                                                                  <  2:1410223/150‑1 (MQ=255)
  cTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAAACAGGCTTACCGTAGGCCGGATCAGGCGATCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCACGGTCGGTCCCCTCGCCCCACTgggg                                                                                                                                                 >  1:388772/1‑150 (MQ=255)
  cTCATTGCAATCTTTTCCGCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGAGCAGGCGTGCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTgggg                                                                                                                                                 <  1:154471/150‑1 (MQ=255)
    cATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTAGGTCCACTCGCCCCTCTGGGgag                                                                                                                                               >  2:412568/1‑150 (MQ=255)
       tGCATTCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCTCCCCAGAGGGGCGAgg                                                                                                                                            <  2:1184828/150‑1 (MQ=39)
       tGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTACTCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAACGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                            <  2:111049/150‑1 (MQ=255)
       tGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                            <  2:1349317/150‑1 (MQ=255)
       tGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCACTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                            <  2:385218/150‑1 (MQ=255)
        gCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAggg                                                                                                                                           >  2:650074/1‑150 (MQ=255)
         cAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCTCCCCAGAGGGGCGAgggg                                                                                                                                          <  1:571514/150‑2 (MQ=39)
         cAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGGCGGTCCCCTCTCCCCAGAGGGGCGAgggg                                                                                                                                          <  1:1400068/150‑2 (MQ=39)
          aaTCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGTATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGtt                                                                                                                                         <  2:871123/150‑1 (MQ=255)
          aaTCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGtt                                                                                                                                         <  2:404832/150‑1 (MQ=255)
            tCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACCGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAg                                                                                                                                       <  2:775096/150‑1 (MQ=255)
             cTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAgg                                                                                                                                      <  1:1275190/150‑1 (MQ=255)
                  cctcctCGTGCCTGATGCACTTCGCTAATCAGGCTTACGGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTCACCGTGTTCAGAGTCTAATAGGGGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAg                                                                                                                                 <  1:14994/150‑1 (MQ=38)
                  ccTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAg                                                                                                                                 >  1:1021506/1‑150 (MQ=255)
                     cATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgggg                                                                                                                              >  1:884872/1‑150 (MQ=255)
                     cATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgggg                                                                                                                              >  1:34237/1‑150 (MQ=255)
                      aTCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTTGGGTGAGGGGa                                                                                                                             >  2:792360/1‑150 (MQ=255)
                       tCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGaa                                                                                                                            >  1:405238/1‑150 (MQ=255)
                                aTGCACTTCGCTAATCAGGCTTACCGTAGGCCGGAGCCGGCGTTCACGCCGCATCGGGCAATGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGTTCCCCTCGCCCCTCTGGGGTGAGGGTTAGGGTGAGGGGAAAACCGTGTc                                                                                                                   <  1:1692992/150‑1 (MQ=39)
                                 tGCACTTCGCTAATCAGGCTTACCGTCGGCCGGATCCGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCGCGGTCGGGCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGTGGGGAAAACCGTGTCa                                                                                                                  <  1:437793/150‑1 (MQ=39)
                                  gCACTTCGCTAATCAGGCGTACCTTAGGCGGGATCTGGCGTTCAGGCCGCATCCGGCAAAGTGTTAACCGTGTTAAGTGTCAAAAAGGCGCGCGGGCGGCCCCCTCGCCCCTCTGGGGTTTGGGTTAGGGTGAGGGTAAAACCGTGTCAg                                                                                                                 <  1:1068570/150‑1 (MQ=255)
                                      ctcccTAATCCGGCTAACCGTAGGCCGGCTCCGGCGTTCACGCCGCACCCGGCAATGTGTTAACCGTGTACAGATTCAAAAAGGCGCTCGGTCGGGCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAat                                                                                                             <  1:174862/146‑1 (MQ=37)
                                       tCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGGTAGGGGGAGGGAAAACCGTGTCAGCAatat                                                                                                           >  2:399112/1‑150 (MQ=255)
                                       tCGCAAATCAGGCTTCCGGTAGGCCGGCTCTGGCGTTCACGCCGCATCCGGAAAAGTGTTAACCGTGTCCAGCGGCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGTGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAata                                                                                                            <  2:1253523/150‑1 (MQ=37)
                                        cGCAAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTACCCGTGTACAGAGTCAAAAAGCCGCTCGGGCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAACAatat                                                                                                           <  1:916789/150‑1 (MQ=38)
                                             aTCTGGCTTAGCGTAGGAGGGTTCAGGCGTTCACGCGGCTTCCGGCAAAGTGTTAACCGTGTACAGATTAAAAAAGGCGGTCGGTGGGTCCCCTCCCCCATCTGGGGAGTGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAAcc                                                                                                      <  2:1423436/150‑1 (MQ=255)
                                             aTCAGGCTTACCTTAGGCGGGATCAGGCGTTCACGCCGCATCCGGCAATGTGTTAACCGTGTACAGACTCAAAAAGGCGCTGGGTCGTTCCCCTCGCCCCTCTGGGGAGTGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAAcc                                                                                                      <  1:1614529/150‑1 (MQ=38)
                                              tccggCTTACCGTTGGCCGGATCTGCCGTTCCCGCCGCATCCGGCAAAGTGTTAACCGTGAACAGAGTCAAAATGGCGCTCGGGCGGTCCCCTCGCCCCTCGGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCg                                                                                                     <  2:83845/147‑1 (MQ=38)
                                                 ggCTTACCGTAGGCAGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCCAGTAAAGAGTCAAAAAGGCGCTCGGTCAGTCCCCTCGCCCCTCTCAAGAGAGGTTTAGTGTGAGGGGAAAACAGTGTCATCAATATCAACCttgt                                                                                                  >  2:931532/1‑150 (MQ=255)
                                                  gCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGtt                                                                                                 >  1:889258/1‑150 (MQ=255)
                                                    ttACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGTTATGTCGGTCCCCTCGCCACTCTGTGGAGATGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTa                                                                                               <  1:276093/150‑1 (MQ=39)
                                                      aCCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTAtt                                                                                             >  1:83766/1‑150 (MQ=255)
                                                           cggCGGGACCAGGCGTGCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGTGTCAAAACGCCGCGCGGTCGGCCCCCTCGCCCCGCGGGGGTGAGGGTTAGGGGGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCg                                                                                        <  1:553146/149‑1 (MQ=255)
                                                             gCCGGATGCGGCGTGAACGCCTGATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                            <  2:151654/96‑1 (MQ=37)
                                                             gCCGGATGCGGCGTGAACGCCTGATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAgg                                                                                                                                            >  1:151654/1‑96 (MQ=37)
                                                             gCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCa                                                                                      >  1:1603291/1‑150 (MQ=255)
                                                             gCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCa                                                                                      >  1:1222189/1‑150 (MQ=255)
                                                              ccGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCaa                                                                                     <  2:256818/150‑1 (MQ=255)
                                                               gggATCAGGCGTTAACGCCGCATCCGGAAAAGTGTTAACCGTGTAAAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTGGCCCCTCTGGGGTGAGGGTTTGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCaaa                                                                                    <  2:1404345/149‑1 (MQ=255)
                                                               cGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAATGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCaaa                                                                                    <  2:167212/150‑1 (MQ=255)
                                                               cGGATCAGGCGGTCTCGCCGAATCCGTCAAAGTGTTATCCGTGTACAGAGTCAAAAACGCGCGCGGTCGGTCCCCTCGCCCCTCTGGTGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCaaa                                                                                    <  2:274334/150‑1 (MQ=38)
                                                                ggATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCaaaa                                                                                   >  1:852703/1‑150 (MQ=255)
                                                                ggATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCaaaa                                                                                   >  1:475495/1‑150 (MQ=255)
                                                                  ttaaGGCGTTCACGCGGCATGCGGCAAAGTGTTAACGGTGTACAGATTCAAAAAGGCGCTCGGGCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGTGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAaac                                                                                 <  2:388772/147‑1 (MQ=38)
                                                                  aTCAGGGGTTCACGCCTCATCCGGCAAAGTGTTAACCGTGTAATGCGTCAAATAGCCGCGAGGGCGGCCCGCGCCCCCCTGTTGGTAGAGGGTTAGGGGGCGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAaac                                                                                 <  2:556376/150‑1 (MQ=255)
                                                                      ggCGTTCAGGCCGCATCCGGCAAAGTATTAACCGTGTACAGAGGCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcg                                                                             <  2:251689/150‑1 (MQ=255)
                                                                        cGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTGCGCAACAACAAcgcg                                                                           >  1:84253/1‑150 (MQ=255)
                                                                        cGTGCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGGCAAAAAGGCTCTCGGTCGGTCCCCTCGCCCCTCTGGGGTGAGGGTTTGGGTGAGGGGAAAACCGTGGCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgcg                                                                           <  2:319566/150‑1 (MQ=39)
                                                                        agTTCACGCCGCATCCGGCAAAGTGTTAACCGTTTACAGAGTCAAACAGGCGCGCGGGCGGTCCCCGCGCCCCTCGGGGGAGAGGGTTTGGGGGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAAcgcg                                                                           <  2:314888/149‑1 (MQ=255)
                                                                          gtCACGCCGCATCCGGCAAAGTGGTAACCGGTTACAGAGTCAAAAAGGCGCGCGGGCGGTACCCTCGCCCCTCTGGGGAGAGGGTGAGGGGGAGGGGAAACCCGTGTCCGCACTATCAACCGTGGTTATTCTTCGCAAAAACAAcgcgcg                                                                         <  2:262636/149‑1 (MQ=37)
                                                                            cACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGcc                                                                       >  2:357581/1‑150 (MQ=255)
                                                                            cACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGcc                                                                       >  1:656737/1‑150 (MQ=255)
                                                                              cgccgcATCCGGCAAAGTGTTAACCGTGTACAGATTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGAGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAg                                                                     >  1:311218/1‑150 (MQ=255)
                                                                              cgccgcATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGGGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAg                                                                     <  1:650074/150‑1 (MQ=255)
                                                                                    aTCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGtt                                                               <  2:389333/150‑1 (MQ=255)
                                                                                     tCCTGCTAAGTGTTAACCGTGTACTGAGTCAAAAAGGCGCTGGGGCGGGCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGATAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTc                                                              <  2:1064135/150‑1 (MQ=38)
                                                                                     tCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTc                                                              >  2:1718228/1‑150 (MQ=255)
                                                                                       cGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcag                                                            >  2:379638/1‑150 (MQ=255)
                                                                                         gCAAAGTGATAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagca                                                          <  1:53380/150‑1 (MQ=255)
                                                                                          cAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagcag                                                         >  2:711/1‑150 (MQ=255)
                                                                                           aaaGTTTTAACCGTGAACAGCGTCACAACGGCGCTCGGGGGGTCCCCGCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTATATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagcagc                                                        <  2:1067292/150‑1 (MQ=38)
                                                                                             aGTGTTAACCGTGTACAGAGTCAAAACGGCGCTCGGGCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTTAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTcagcagcag                                                      <  2:1375506/150‑1 (MQ=255)
                                                                                                    aCCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTAGGcac                                               >  1:821534/1‑150 (MQ=255)
                                                                                                      cGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCGTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGGGTTGATTCTTCGCAAACCAACGCGCGCCAGAATGTTCAGCCGCAGTACGCACAgg                                            >  1:1150199/1‑150 (MQ=38)
                                                                                                         gtACAGAGTCAAACAGGCGCTCGGGCGGTCCCGTCGCCCCTCTGGGGCGAGGGTTAGGTTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGt                                          <  2:1713777/150‑1 (MQ=255)
                                                                                                           aCAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCACAATGTTCAGCAGCAGTACGCACAGGGaaa                                        >  2:1536262/1‑150 (MQ=255)
                                                                                                               agTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGAt                                    >  1:171572/1‑150 (MQ=255)
                                                                                                                gTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACACCGCGCGCCAGAATGTTCAGCAGCAGTCCGCACAGGGTAATGATc                                   >  2:1454675/1‑150 (MQ=255)
                                                                                                                gTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATc                                   >  1:1313132/1‑150 (MQ=255)
                                                                                                                  cAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCaa                                 >  1:376704/1‑150 (MQ=255)
                                                                                                                   aaaaaGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGCTccag                                >  1:304422/1‑147 (MQ=255)
                                                                                                                       cggCGCTCGGTCGGGCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATAccc                            <  1:1238284/149‑1 (MQ=255)
                                                                                                                       aGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATAccc                            >  2:226382/1‑150 (MQ=255)
                                                                                                                            cTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGcc                       <  2:405238/150‑1 (MQ=255)
                                                                                                                             tcggtcggTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGccc                      <  2:1148405/150‑1 (MQ=255)
                                                                                                                               ggtcggtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAACAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTTATGTGCAATACCCCGGCCCAg                    >  2:263737/1‑150 (MQ=255)
                                                                                                                               ggtcggtcCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAAGGATCAATACCCCGGCCCAg                    >  2:1440278/1‑150 (MQ=255)
                                                                                                                                ggcggccCCCTCGCCCCTCTGGGGTGAGGGTTAGGGGGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAgg                   <  2:753678/148‑1 (MQ=255)
                                                                                                                                 ccctaacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCa                                  >  2:1283406/7‑134 (MQ=255)
                                                                                                                                 ccctaacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCa                                  <  1:1283406/128‑1 (MQ=255)
                                                                                                                                 ccctaacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGc                  >  2:1291207/7‑150 (MQ=255)
                                                                                                                                    taacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCaa               >  1:879468/4‑150 (MQ=255)
                                                                                                                                      ccccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTCCGCAAAAACAACGCGCGCCAGAATGTTCCGCAGCAGTACGCACAGGGTAAAGAGCACTACCCCGGCCCACGCCACtt             >  2:1342511/1‑150 (MQ=38)
                                                                                                                                      acccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAACCCGAGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACCGGGTAATGATCACTACCCCGGCCCAGGCCCAtt             >  1:237555/2‑150 (MQ=39)
                                                                                                                                         cTCGCCCCTCTGGGGTGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgct          >  1:361479/1‑150 (MQ=255)
                                                                                                                                          tCGCCCCTCTGGGGAGAGGGTTAGGGTTAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgctg         <  1:1337682/150‑1 (MQ=255)
                                                                                                                                          tCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgctg         <  2:34237/150‑1 (MQ=255)
                                                                                                                                          gcGCCCCTCGGGGGTGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgctg         <  2:884872/149‑1 (MQ=255)
                                                                                                                                           cGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATtgctgc        <  2:83766/150‑1 (MQ=255)
                                                                                                                                                ctctGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTc   >  2:1596594/1‑150 (MQ=255)
                                                                                                                                                 tctGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTcc  >  1:694449/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
TACTCATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCC  >  NC_000913/3908360‑3908654

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: