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

GCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAA‑TGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGAT  >  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)
                           aCAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAA‑TGc                                                                                                                        >  2:41911/1‑151 (MQ=255)
                             aaTCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAA‑TGCCt                                                                                                                      <  2:640746/151‑1 (MQ=255)
                             aaTCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAA‑TGCCt                                                                                                                      <  1:374868/151‑1 (MQ=255)
                                        tATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAA‑TGCCTGagacagac                                                                                                             <  2:867642/149‑1 (MQ=255)
                                        tATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAA‑TGCCTGagacagac                                                                                                             >  1:281174/1‑149 (MQ=255)
                                                 tCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAA‑TGCCTGAGACAGACAGCCTCAAGc                                                                                                   >  2:251982/1‑150 (MQ=255)
                                                  cATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAA‑TGCCTGAGACAGACAGCCTCAAGCAc                                                                                                 >  1:634179/1‑151 (MQ=255)
                                                       ctctctTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAA‑TGCCTGAGACAGACAGCCTCAAGCACccgcc                                                                                            >  2:100991/1‑151 (MQ=255)
                                                                   tGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAA‑TGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCg                                                                                >  2:896706/1‑151 (MQ=255)
                                                                   tGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAA‑TGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCg                                                                                >  2:699281/1‑151 (MQ=255)
                                                                    gACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGGGGCCTCGGGTTGATGGTAAAATATCACTCGGGTCTTTTCTCTATCTGCCGTTCAGCTAA‑TGCCTGAGACAGACAGCCTCAAGCCCCCGCCGCTATTATATc                                                                                 <  2:633517/149‑1 (MQ=255)
                                                                          gcGGTGAGTAAGAGCCATTCAACTTTTTTTCACAATATAGGACGCGTGGGGTTGATAAAAAAATCCCTCTGGGGTTTTTTCTCTACCTGCCGTTCTGCTAA‑TCCCGAAAACAGACAGCCTAAACCCCCCGCCGCTATTATATCGctct                                                                            <  1:325523/147‑1 (MQ=255)
                                                                                   aaGAGGCACTCTACATGTGTTCAGAATAGAGGAGGCCTGGGGTTGATGGAAAAATACCCCTGGGGTTTTTTCTCTATCTCCCTTTCAGTTAA‑TGCCTGAGACAGACACCCCCAACCCCCCGCCGCTATTATATCGCTCTCTTTAAcc                                                                    <  2:684308/147‑1 (MQ=255)
                                                                                   aaGAGCCCCTCTACATTTTTTCAGCAAATAGGAGCCCTGGGTTTGATGTAAAAATACCCCTGGGGGTTTTTCTCATTCTGCCGTTAAGTTAA‑TGCTGGAGACAGACACCCTCAACCACCCGCCGCTATTATATCGCTCTCTTTaa                                                                      <  2:34363/145‑1 (MQ=255)
                                                                                   aaGAGCCACTCTACATGTGTTCAGAAAATAGGAGGCCTGGGGTTGATGGTAAAATACCACCCGGGTTTTTTCTCTATCGCCCGTTCAGCTAA‑TGCCTGAGACAGACACCCTAAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCa                                                                  <  2:398154/149‑1 (MQ=255)
                                                                                     gagGCACTCTACATTTGTTCAGCATAGAGGGGCCCTCGGGTTGATGATAAAATACCCCGCGGGGTTTTTCTCTCTCCGCCGTTAATCTAA‑TGCCTGAAACAAACACCTACAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCAttt                                                               <  1:941537/150‑1 (MQ=255)
                                                                                     gagGCACTCTACATGTTTTCAGAATAAGGGGGGCCTGGGTTTGATGGAAAAATATCCCTCGGGTTTTTTCTCTATCTGCCTTTCAGCTAA‑TCCCTGAGACAGCCACCCCAAACCCCCCGCCGCTATTATATCGCTCTCTTTAACCCAttt                                                               <  2:325944/150‑1 (MQ=255)
                                                                                         ccctctACATGTGTTCAAAATATAGGAGGCCTGGGGTTGATGGAAAAATTTCCCTCGGGGTTTTTCTTTATCTGCCGTTCAGCTAA‑TGCCTGAGAAAGACAGCCACAACCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGt                                                            <  2:859511/147‑1 (MQ=255)
                                                                                                cATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAA‑TGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGAt                                                   >  2:591107/1‑151 (MQ=255)
                                                                                                  tgtgTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAA‑TGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTc                                                 >  2:427271/1‑151 (MQ=255)
                                                                                                   gtgtTCATCATATAGGGGGCCTGGGTTTGGTGAAAAAATCTCCCGCGGGTTTTTTCTCTCCCTCCCTTTCAGCTAA‑TCCCTGAGACAGACACCTTAAACCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCg                                                     <  1:528962/146‑1 (MQ=255)
                                                                                                                aggaggCTTGGGGTTAAGGAAAAAATACCACTGGGGGCTTTTCTCTATCCCCCGTTCAGTTAA‑TGCCTGAAACAGCCCCCCTCAACCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCt                                                <  1:641941/138‑1 (MQ=255)
                                                                                                                    ggCCTCGGTTGGATGATAAAATATCCCTGGGGGTTTTTCTCTATCCGCCGTTCATCTAA‑TCCCTGAGACAGACACCCTCAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTc                               <  1:279673/151‑1 (MQ=255)
                                                                                                                    ggCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAA‑TGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTc                               >  1:593301/1‑151 (MQ=255)
                                                                                                                      cctgggggtgaagaaaaaaTATCACGCGGGGTTTTTCTCTATCTGCCGTTCA‑CTAATTCCCAGAGACAGACCCCCCAAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGc                             <  2:892216/136‑1 (MQ=255)
                                                                                                                      ccTGGGTTTTATGATAAAATATCCCGGGGCGTTTTTCTCTATCTGCCTTTCATCTAT‑TGCCTGAGACCAACTGCCTCAACCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGAc                                    <  2:237609/144‑1 (MQ=255)
                                                                                                                          gggTTAAGGTTAAAATATCGTTGGGGGTTTTTCTTTTCCTCCCTTTCATCTAT‑TGCCGCAGAAAGACACCCTAAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCAttt                         <  2:739484/151‑1 (MQ=255)
                                                                                                                             tggatgaaaaaaTATCCCTGGGGCTTTTTCTCTACCGGCCTTTCACCTAA‑TCCCTGAGACAGACACCCTCAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCAttttt                       <  1:717027/142‑1 (MQ=255)
                                                                                                                                 tGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAA‑TGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGc                  >  2:160169/1‑151 (MQ=255)
                                                                                                                                   gTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAA‑TGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGt                >  2:146980/1‑151 (MQ=255)
                                                                                                                                      aaaTACCCCCCGGGGTTTTTCTTTACCTCCCGTTCAGCAAA‑TGCCGGAGACAGACACCCTCAAGCCACCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAAt             <  2:71506/151‑1 (MQ=255)
                                                                                                                                        atatCCCTCGGGCTTTTTCTCTATCTGCCTTTAAGTTAA‑TGCCTGAGACAGACCCCCTAAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAAttt           <  1:236735/151‑1 (MQ=255)
                                                                                                                                             aCTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAA‑TGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTa       >  2:349397/1‑150 (MQ=255)
                                                                                                                                              ccgggggTTTTTCTCTTTCTGCCTTTCATCTAA‑TCCTTGAGACAACCACCCTCAACCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAAt     <  1:740415/148‑1 (MQ=255)
                                                                                                                                              cTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAA‑TGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAAt     >  1:312677/1‑151 (MQ=255)
                                                                                                                                               tCGGGGCTTTTCTCTATCTGCCGTTCAGCTAA‑TGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATg    >  1:543206/1‑151 (MQ=255)
                                                                                                                                                 ggggTTTTTCTTTATCTGCCGTTCAGCTAA‑TGCCTGAGACAGACACCCTCAAGCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGAt  <  2:281225/151‑1 (MQ=255)
                                                                                                                                                 ggggTTTTTCTCTATCTGCCGTTAAGCAAA‑TGCCGGAGACAGACACCCTCAACCCCCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGAt  <  2:439363/151‑1 (MQ=255)
                                                                                                                                                 ggggCTTTTCTCTATCTGCCGTTCAGCTAA‑TGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGAt  >  2:208330/1‑151 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GCGGTGATACAGATGACGATCAGGGCGACAATCATCGCCTTATGCTGCTTCATTGCTCTCTTCTCCTTGACCTTACGGTCAGTAAGAGGCACTCTACATGTGTTCAGCATATAGGAGGCCTCGGGTTGATGGTAAAATATCACTCGGGGCTTTTCTCTATCTGCCGTTCAGCTAA‑TGCCTGAGACAGACAGCCTCAAGCACCCGCCGCTATTATATCGCTCTCTTTAACCCATTTTGTTTTATCGATTCTAATCCTGAAGACGCCTCGCATTTTTGTGGCGTAATTTTTTAATGAT  >  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: ‑