Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091317,0000CT40.0% 60.1 / 17.0 45intergenic (‑40/‑489)mokC/nhaAregulatory protein for HokC, overlaps CDS of hokC/sodium‑proton antiporter
Reads supporting (aligned to +/- strand):  ref base C (19/8);  new base T (0/18);  total (19/26)
Fisher's exact test for biased strand distribution p-value = 1.10e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGAT  >  NC_000913/16851‑17146
                                                                                                                                                     |                                                                                                                                                   
gcggTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCt                                                                                                                                                    <  2:971359/151‑1 (MQ=255)
     gATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTtct                                                                                                                                               >  1:136157/1‑151 (MQ=255)
        aCAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTa                                                                                                                                            <  2:861288/151‑1 (MQ=255)
                    cAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTc                                                                                                                                  >  2:281560/1‑149 (MQ=255)
                     aGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGc                                                                                                                               >  1:865861/1‑151 (MQ=255)
                       ggCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTa                                                                                                                             <  1:369123/151‑1 (MQ=255)
                           aCAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGc                                                                                                                        >  2:41911/1‑151 (MQ=255)
                             aaTCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCt                                                                                                                      <  2:640746/151‑1 (MQ=255)
                             aaTCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCt                                                                                                                      <  1:374868/151‑1 (MQ=255)
                                        tATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGagacagac                                                                                                             <  2:867642/149‑1 (MQ=255)
                                        tATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGagacagac                                                                                                             >  1:281174/1‑149 (MQ=255)
                                                 tCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGc                                                                                                   >  2:251982/1‑150 (MQ=255)
                                                  cATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCAc                                                                                                 >  1:634179/1‑151 (MQ=255)
                                                       ctctctTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACccgcc                                                                                            >  2:100991/1‑151 (MQ=255)
                                                                   tGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCg                                                                                >  2:896706/1‑151 (MQ=255)
                                                                   tGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCg                                                                                >  2:699281/1‑151 (MQ=255)
                                                                    gACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGATGGTAAAATATCACTCGGGTCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCCCCCGCCGCTATTATATc                                                                                 <  2:633517/149‑1 (MQ=255)
                                                                          gcGGTGAGTAAGAGCCATTCAACTTTTTTTCACAATATAGGACGCGTGGGGTTGATAAAAAAATCCCTCTGGGGTTTTTTCTCTACCTGCCGTTCTGCTAATCCCGAAAACAGACAGCCTAAACCCCCCGCCGCTATTATATCGctct                                                                            <  1:325523/147‑1 (MQ=255)
                                                                                   aaGAGGCACTCTACATGTGTTCAGAATAGAGGAGGCCTGGGGTTGATGGAAAAATACCCCTGGGGTTTTTTCTCTATCTCCCTTTCAGTTAATGCCTGAGACAGACACCCCCAACCCCCCGCCGCTATTATATCGCTCTCTTTAAcc                                                                    <  2:684308/147‑1 (MQ=255)
                                                                                   aaGAGCCCCTCTACATTTTTTCAGCAAATAGGAGCCCTGGGTTTGATGTAAAAATACCCCTGGGGGTTTTTCTCATTCTGCCGTTAAGTTAATGCTGGAGACAGACACCCTCAACCACCCGCCGCTATTATATCGCTCTCTTTaa                                                                      <  2:34363/145‑1 (MQ=255)
                                                                                   aaGAGCCACTCTACATGTGTTCAGAAAATAGGAGGCCTGGGGTTGATGGTAAAATACCACCCGGGTTTTTTCTCTATCGCCCGTTCAGCTAATGCCTGAGACAGACACCCTAAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCa                                                                  <  2:398154/149‑1 (MQ=255)
                                                                                     gagGCACTCTACATTTGTTCAGCATAGAGGGGCCCTCGGGTTGATGATAAAATACCCCGCGGGGTTTTTCTCTCTCCGCCGTTAATCTAATGCCTGAAACAAACACCTACAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCAttt                                                               <  1:941537/150‑1 (MQ=255)
                                                                                     gagGCACTCTACATGTTTTCAGAATAAGGGGGGCCTGGGTTTGATGGAAAAATATCCCTCGGGTTTTTTCTCTATCTGCCTTTCAGCTAATCCCTGAGACAGCCACCCCAAACCCCCCGCCGCTATTATATCGCTCTCTTTAACCCAttt                                                               <  2:325944/150‑1 (MQ=255)
                                                                                         ccctctACATGTGTTCAAAATATAGGAGGCCTGGGGTTGATGGAAAAATTTCCCTCGGGGTTTTTCTTTATCTGCCGTTCAGCTAATGCCTGAGAAAGACAGCCACAACCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGt                                                            <  2:859511/147‑1 (MQ=255)
                                                                                                cATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGAt                                                   >  2:591107/1‑151 (MQ=255)
                                                                                                  tgtgTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTc                                                 >  2:427271/1‑151 (MQ=255)
                                                                                                   gtgtTCATCATATAGGGGGCCTGGGTTTGGTGAAAAAATCTCCCGCGGGTTTTTTCTCTCCCTCCCTTTCAGCTAATCCCTGAGACAGACACCTTAAACCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCg                                                     <  1:528962/146‑1 (MQ=255)
                                                                                                                aggaggCTTGGGGTTAAGGAAAAAATACCACTGGGGGCTTTTCTCTATCCCCCGTTCAGTTAATGCCTGAAACAGCCCCCCTCAACCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCt                                                <  1:641941/138‑1 (MQ=255)
                                                                                                                    ggCCTCGGTTGGATGATAAAATATCCCTGGGGGTTTTTCTCTATCCGCCGTTCATCTAATCCCTGAGACAGACACCCTCAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTc                               <  1:279673/151‑1 (MQ=255)
                                                                                                                    ggCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTc                               >  1:593301/1‑151 (MQ=255)
                                                                                                                      cctgggggtgaagaaaaaaTATCACGCGGGGTTTTTCTCTATCTGCCGTTCACTAATTCCCAGAGACAGACCCCCCAAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGc                             <  2:892216/136‑1 (MQ=255)
                                                                                                                      ccTGGGTTTTATGATAAAATATCCCGGGGCGTTTTTCTCTATCTGCCTTTCATCTATTGCCTGAGACCAACTGCCTCAACCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGAc                                    <  2:237609/144‑1 (MQ=255)
                                                                                                                          gggTTAAGGTTAAAATATCGTTGGGGGTTTTTCTTTTCCTCCCTTTCATCTATTGCCGCAGAAAGACACCCTAAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCAttt                         <  2:739484/151‑1 (MQ=255)
                                                                                                                             tggatgaaaaaaTATCCCTGGGGCTTTTTCTCTACCGGCCTTTCACCTAATCCCTGAGACAGACACCCTCAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCAttttt                       <  1:717027/142‑1 (MQ=255)
                                                                                                                                 tGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGc                  >  2:160169/1‑151 (MQ=255)
                                                                                                                                   gTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGt                >  2:146980/1‑151 (MQ=255)
                                                                                                                                      aaaTACCCCCCGGGGTTTTTCTTTACCTCCCGTTCAGCAAATGCCGGAGACAGACACCCTCAAGCCACCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAAt             <  2:71506/151‑1 (MQ=255)
                                                                                                                                        atatCCCTCGGGCTTTTTCTCTATCTGCCTTTAAGTTAATGCCTGAGACAGACCCCCTAAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAAttt           <  1:236735/151‑1 (MQ=255)
                                                                                                                                             aCTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTa       >  2:349397/1‑150 (MQ=255)
                                                                                                                                              ccgggggTTTTTCTCTTTCTGCCTTTCATCTAATCCTTGAGACAACCACCCTCAACCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAAt     <  1:740415/148‑1 (MQ=255)
                                                                                                                                              cTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAAt     >  1:312677/1‑151 (MQ=255)
                                                                                                                                               tCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATg    >  1:543206/1‑151 (MQ=255)
                                                                                                                                                 ggggTTTTTCTTTATCTGCCGTTCAGCTAATGCCTGAGACAGACACCCTCAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGAt  <  2:281225/151‑1 (MQ=255)
                                                                                                                                                 ggggTTTTTCTCTATCTGCCGTTAAGCAAATGCCGGAGACAGACACCCTCAACCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGAt  <  2:439363/151‑1 (MQ=255)
                                                                                                                                                 ggggCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGAt  >  2:208330/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAATGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGAT  >  NC_000913/16851‑17146

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