Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,908,512 A→C 8.3% intergenic (‑146/+37) pstB ← / ← pstA phosphate ABC transporter ATPase/phosphate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,908,5120AC 8.3% 220.4 / 10.2 95intergenic (‑146/+37)pstB/pstAphosphate ABC transporter ATPase/phosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (41/38);  new base C (3/5);  total (45/50)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

CATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAACG  >  NC_000913/3908364‑3908661
                                                                                                                                                    |                                                                                                                                                     
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:475495/1‑150 (MQ=255)
                                                            ggATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCaaaa                                                                                          >  1:852703/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                                                                              >  1:656737/1‑150 (MQ=255)
                                                                        cACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGcc                                                                              >  2:357581/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)
                                                                                                                                              ctGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCg        <  2:1313132/150‑1 (MQ=255)
                                                                                                                                               tGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAAAACCCCGGCCCAGGCCAATTGATGCCATTCCGc       >  1:55852/1‑150 (MQ=255)
                                                                                                                                                ggggAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACACCGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCa      >  1:1017241/1‑150 (MQ=255)
                                                                                                                                                  ggAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCaaa    >  2:543168/1‑150 (MQ=255)
                                                                                                                                                  ggAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCaaa    >  1:465703/1‑150 (MQ=255)
                                                                                                                                                    agagGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAACg  >  1:1575721/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CATTGCAATCTTTTCCTCATCGTGCCTGATGCACTTCGCTAATCAGGCTTACCGTAGGCCGGATCAGGCGTTCACGCCGCATCCGGCAAAGTGTTAACCGTGTACAGAGTCAAAAAGGCGCTCGGTCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAACCGTGTCAGCAATATCAACCGTGTTTATTCTTCGCAAAAACAACGCGCGCCAGAATGTTCAGCAGCAGTACGCACAGGGTAATGATCAATACCCCGGCCCAGGCCAATTGCTGCCATTCCGCAAACG  >  NC_000913/3908364‑3908661

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: