Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,622,8020AT42.9% 61.5 / 34.6 49intergenic (+61/+158)yneM/dgcZinner membrane‑associated protein/diguanylate cyclase, zinc‑sensing
Reads supporting (aligned to +/- strand):  ref base A (21/7);  new base T (0/21);  total (21/28)
Fisher's exact test for biased strand distribution p-value = 3.24e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTAATATGAATGTTTTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATG  >  NC_000913/1622652‑1622948
                                                                                                                                                         |                                                                                                                                                     
ggTAATATGAATGTTTTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGa                                                                                                                                                       <  2:214745/151‑1 (MQ=255)
   aaTATGAATGTTTTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGAttt                                                                                                                                                   <  2:935988/151‑1 (MQ=255)
   aaTATGAATGTTTTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGAtt                                                                                                                                                    <  2:639882/150‑1 (MQ=255)
    atatGAATGTTTTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGAtttt                                                                                                                                                  >  2:550454/1‑151 (MQ=255)
    atatGAATGTTTTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGAtttt                                                                                                                                                  <  2:560536/151‑1 (MQ=255)
     tatGAATGTTTTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGAttttt                                                                                                                                                 >  1:667346/1‑151 (MQ=255)
      atGAATGTTTTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGAttttt                                                                                                                                                 >  2:268439/1‑150 (MQ=255)
            gTTTTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACg                                                                                                                                           <  1:104302/150‑1 (MQ=255)
                ttATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGa                                                                                                                                      >  1:594321/1‑150 (MQ=255)
                 tATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGaa                                                                                                                                     >  1:242965/1‑150 (MQ=255)
                    ggCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGaaaaa                                                                                                                                  >  1:679536/1‑150 (MQ=255)
                        gTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAAcc                                                                                                                            >  2:357729/1‑151 (MQ=255)
                        gTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAAcc                                                                                                                            <  2:871019/151‑1 (MQ=255)
                            tGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGcc                                                                                                                         >  2:365767/1‑150 (MQ=255)
                              ggAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGt                                                                                                                      >  1:140035/1‑151 (MQ=255)
                                     ttttattttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACcgg                                                                                                               >  2:914901/1‑151 (MQ=255)
                                      tttattttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACcgg                                                                                                               >  1:339189/1‑150 (MQ=255)
                                       tttttttctgtttttctggCCGCTTTTTCCACCAAAAAAGGGGATGACTAATAAAGGAAAATACCCCCCCCCCTACCCGCCCCTTTTTTACAGTTGTGTACAAGGGGCTTGTTTTTTATGACGGCAAAAAAAAACCGCCAGTAAACcggc                                                                                                              <  1:168591/147‑1 (MQ=255)
                                          ttttctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGt                                                                                                           >  2:305902/1‑150 (MQ=255)
                                          tgttgtgtttttcgggCCGATTATTACACACAAAAGGGGGATGAAAAATGAGGGGAAAAAACCCCTCACAAACCGGCCCCTTGGTTAGATTTGGGAAAGGGGGCCTTGTTTTTTAGGACGAAAAAAAAAAACCGCCAGTAAACCGGCGGt                                                                                                           <  2:113674/142‑1 (MQ=255)
                                           tttttgtttttcgggCCGCTTATTCCACCCAAAAAGGGGATGAATAATGAAGGGAAAAAACCCCCCCCCAACCCGCCCCTTTTTAAATGTTGTGAACAGGGGCCCTGTTTTTTATGCCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTg                                                                                                          <  1:909710/150‑1 (MQ=255)
                                             tctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGaa                                                                                                        >  1:546514/1‑150 (MQ=255)
                                                gtttttcgggCCGGTTTTTTCACCAAAAAAGGGGATAAATAATGAACGGAAATAACCCCTCCCCTACCCGCCCCCTTGTACCGGTGGTGCACAGGGGGCTTGTTTTTTGGGACGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGc                                                                                                     <  1:554708/148‑1 (MQ=255)
                                                 gttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCtt                                                                                                   >  1:732063/1‑151 (MQ=255)
                                                  ttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCtt                                                                                                   >  2:389026/1‑150 (MQ=255)
                                                       tggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATgg                                                                                             >  1:643774/1‑151 (MQ=255)
                                                        gcccGCTTTTTTCACCAAAAAAGGGGATGACAAATGAAGGGAAATACCCCCTCCCCTACCCGCCCCTTTTTTACGTTTGTGTAAAGGGGGCCTGTTTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCAt                                                                                               <  1:99563/146‑1 (MQ=255)
                                                        gcccGCTTATTCCACCCAAAAAGGGGATGTCAAATGAAGGGAGATACCCCCCCCCCAACCCGCCCCTTTTTTAATGTTGTGAAAAAGGGCCTTGTTTTTTATGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGa                                                                                            <  2:89593/149‑1 (MQ=255)
                                                              ttATTCCAGCAAAAAAGGGGAAGACAAAGGAGAGGAGAAATTCCCCCACCTACCCGCCCTTTTTTTAGGGTTGTAAAAAAGGGGCTTGTTTTTTAGGGGGGAAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAg                                                                                         <  1:22816/148‑1 (MQ=255)
                                                                     acccAAAAAGGGGATAACTAAGGAACGGAAATACCCCCCCCCCAACCCGCCCCTTTTTTAAAGTTGTGTACAGGGGGCTTTTTTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGtt                                                                                <  1:119749/148‑1 (MQ=255)
                                                                          aaaaGGGGATGACTAATAAACGAAAATACCCCCTCCCCAACCCGCCCCTTTTTTACGTTTGTGTAAAGGGGGCTTTTTTTTTAGGACGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTg                                                                             <  1:747251/147‑1 (MQ=255)
                                                                           aagtGGGATAACAAATGAACGAGAAAAACCCCTCCCCTACCCCGCCCCTTTTTAAATTTGTTTACAAGGGCCTTTTTTTTTAGCACGGAAAAAAAAATCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCttt                                                                         <  2:758904/148‑1 (MQ=255)
                                                                           aaaTGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCttt                                                                         >  1:157548/1‑151 (MQ=255)
                                                                            aaggggATGAAAAATGACGGGAAAAAACCCCCCCCCTACCCGCCCCTTTTTTACGGTTGTGTACAGGGGGCTTTTTTTTTATGACGGGAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCtttt                                                                        <  1:304861/148‑1 (MQ=255)
                                                                               gggAAGAATAATAAACGGATAAACCCCCCCACCTACCCCCCCCTTTGTTATGGTGGTTAACAAGGGGCTTGTTTTTGATGAGGACAAAAAAAATCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCtttt                                                                        <  2:506328/148‑1 (MQ=255)
                                                                                  aTGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCt                                                                   >  1:756864/1‑150 (MQ=255)
                                                                                        aaTAAAGGGAAATAACCCCTCCCCAACCCGCCCCTTTTTTACATTTTTGAAAAGGGGCCTTTTTTTTTAGGCGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAg                                                              <  1:692412/149‑1 (MQ=255)
                                                                                            aaCGGAAATACCCCCCCCCCTACCCGCCCCTTTGTTAGGTTTGTGTACAGGGGGGGTGTTTTTTAGGAGGGAAAAAAAAATCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACa                                                               <  2:37782/144‑1 (MQ=255)
                                                                                                 aaaaaccccccCACCTACCGGCCCCCTTTTTACAGTGTTTAACAGGGGCCTTGTTTTTTATGGCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAtttt                                                       <  1:937981/140‑1 (MQ=255)
                                                                                                            caccgccccGTCCCTTTGTTGCAGTTGTGTACAAGGGGCTTGTTTTTTTTGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACg                                        <  1:598467/144‑1 (MQ=255)
                                                                                                              ccaacccGGCCCTTTTTTACTGTTGTGTAAAAGGGTCTTTTTTTTGATGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACg                                        <  2:772904/144‑1 (MQ=255)
                                                                                                                   ccGCCCCTTTTTTACATTTGTGTAAAGGGGCCCTTTTTTTTAGGGGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAAt                                     <  1:680077/147‑1 (MQ=255)
                                                                                                                                      ttGTGTAGAAGGGGCTTGTTTTTTAGGACGACAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATc                   <  2:620522/148‑1 (MQ=255)
                                                                                                                                        gtgtACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGttt              >  1:256410/1‑151 (MQ=255)
                                                                                                                                            aaaaGGGGCCTTTTTTTTTAGGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtct           <  1:684662/149‑1 (MQ=255)
                                                                                                                                            aCAGGGGGCTTGTTTTTTAGGGCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtc            <  2:643766/149‑1 (MQ=255)
                                                                                                                                             cAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtctc          >  2:587303/1‑150 (MQ=255)
                                                                                                                                             acaGGGGGCCGATTTTTTAGACGGAGAAAAAAAATCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTtc            <  2:944703/147‑1 (MQ=255)
                                                                                                                                               aGGGGCCTTTATTTTTAGGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGtttt             <  1:608465/146‑1 (MQ=255)
                                                                                                                                                    ccTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatgaatg  >  1:825493/1‑151 (MQ=255)
                                                                                                                                                         |                                                                                                                                                     
GGTAATATGAATGTTTTTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATG  >  NC_000913/1622652‑1622948

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