Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,622,8160GA41.0% 31.8 / 35.3 40intergenic (+75/+144)yneM/dgcZinner membrane‑associated protein/diguanylate cyclase, zinc‑sensing
Reads supporting (aligned to +/- strand):  ref base G (22/1);  new base A (0/16);  total (22/18)
Fisher's exact test for biased strand distribution p-value = 4.51e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAA  >  NC_000913/1622668‑1622964
                                                                                                                                                      |                                                                                                                                                    
ttATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGCTAATCCCTCACCTAACCGGCCCCTTGTTAAAGTTGGGTACAAGGGGC‑C‑TGATTTTTATAACGGCGa                                                                                                                                                     >  1:915119/1‑150 (MQ=255)
ttATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGa                                                                                                                                                     >  2:573020/1‑150 (MQ=255)
  aTGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGaaa                                                                                                                                                   >  2:158264/1‑150 (MQ=255)
    ggCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGaaaaa                                                                                                                                                 >  2:426570/1‑150 (MQ=255)
     gCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGaaaaaa                                                                                                                                                >  1:459843/1‑150 (MQ=255)
          aCTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGAAAAAAAACCg                                                                                                                                           <  1:136018/150‑1 (MQ=255)
                        ttttttttgtttttctggCCGCGTTTTCCACCCAAAAAGGGAATGACAAAGAAGCGGATAAAACCCCCCACCAACCCGGCCTTTTTTCATTGTAGTGTA‑AAGGGGCTG‑TATTTTTTGTGCGGGCAAAAAAAAACCGCCAGTAAACcggcg                                                                                                                             <  2:607327/148‑1 (MQ=255)
                        tattttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACcggcgg                                                                                                                            >  1:774454/1‑151 (MQ=255)
                         attttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACcggcgg                                                                                                                            >  1:202241/1‑150 (MQ=255)
                             tctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAAt                                                                                                                       >  1:132208/1‑151 (MQ=255)
                               tggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATg                                                                                                                      >  1:607291/1‑150 (MQ=255)
                                 ttttcctggCCGCGTATTTCACCCAAAAGTGGAATGATTAATGAACGGATATACTCCCTCGCCTACCCGCCACCTTGAAACATTTATTCAAAAGCGCC‑C‑TGTTTTTTATGCCGGAAAAAAAAACCCGCCAGTAAACCGGCGGTGAAt                                                                                                                       <  1:403872/146‑1 (MQ=255)
                                    ttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGc                                                                                                                 >  1:879419/1‑150 (MQ=255)
                                     tctggCCGCTTTTTTCACCCAAAAAGGGGATGACTAAGGAAGGGAAAAACCCCCTCCCCTACCCGCCCCTTTTTTACATTTGT‑AAAAAGGGGC‑CTTTTTTTTTAGGCCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCa                                                                                                                <  2:283967/150‑1 (MQ=255)
                                         gCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGAt                                                                                                           >  2:964597/1‑151 (MQ=255)
                                              tttttacggCCAAAAAGGGGATAAATAATGAACGAAAAAAACCCCTCCCCTACCCGCCCCTTTTTTACTGTTGTGAAAAGGGGGC‑T‑TGTTTTTTACGGGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATa                                                                                                          <  1:12653/139‑1 (MQ=255)
                                                 ttccggCCAAAAAGGGGAGAACAAAGGACGGGAAATACCCCCCCACCAAACCGCCCCCTTGTAAGAGTGGTGAAAGAGGGTC‑T‑TGTTTTTGGTGCAGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAttt                                                                                                     <  1:935595/144‑1 (MQ=255)
                                                  ttCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTtgtg                                                                                                  >  2:590588/1‑151 (MQ=255)
                                                          aaaaggggtaaaaaaaTGAAGGGAGAAACCCCCTCACCTACCCGCCCCTTTTTTAAGGTTGTGTAGAGGGGGT‑T‑TTTTTTTTGGGGGGGAAAAAAAAATCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCtt                                                                                          <  1:762800/137‑1 (MQ=255)
                                                           aaaGGGGATGAAAAGTAAACGGATATAACCCCTCCCAAACCCGCCCCCTTTTCACAGTTGTG‑AAAAGGGGC‑CTTGTTTTTTAAAAGGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCt                                                                                           <  1:203612/149‑1 (MQ=255)
                                                              tGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTAc                                                                                      >  1:116466/1‑151 (MQ=255)
                                                                   tGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCG‑AAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAAc                                                                                >  2:277484/1‑151 (MQ=255)
                                                                   tGACTAATGAACGGAGATAATCCCTCACCTAACCGGCACCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACg                                                                                     >  1:463996/1‑147 (MQ=255)
                                                                                    tcccccctcccctcCCCGGCCTTTTGTTACGTTTGTGAAAAAGGGCC‑T‑TTTTTTTTATGAGGAAAAAAAAAACCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGc                                                                  <  1:639408/145‑1 (MQ=255)
                                                                                        cccTCCCCTACCCGCCCCTTTTTTACATTGGTGAACAGGGGGC‑T‑TTTTTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACt                                                               <  1:136961/148‑1 (MQ=255)
                                                                                              ccTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGa                                                       >  1:698773/1‑150 (MQ=255)
                                                                                              ccTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTAAGCTAACAGGCATTTTCCTGCACTGATAACGa                                                       >  1:456335/1‑150 (MQ=255)
                                                                                                    cGGCCCCTTTTTTCGGTTGTGAACAAGGGCT‑C‑TTTTTTTTAGGAGGGC‑AAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGtt                                                 <  1:783831/149‑1 (MQ=255)
                                                                                                    cGGCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGGTAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTg                                                >  2:284288/1‑151 (MQ=255)
                                                                                                     ggCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGGGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGATAACAGGCATTTTCCTGCACTGATAACGAATCg                                                   >  1:134191/1‑147 (MQ=255)
                                                                                                     gcccctTTTTTAAATTTGTGAACAGGGGGC‑T‑TGATTTTTATGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGtt                                                 <  1:716480/147‑1 (MQ=255)
                                                                                                           tttttACATTGGTTTACAAGGCGG‑C‑TTTTTTTTATGACGACAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGAcaca                                           <  1:421100/146‑1 (MQ=255)
                                                                                                                   gTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATc                                   >  1:769005/1‑149 (MQ=255)
                                                                                                                    ttGTGTAAAGGGGGC‑T‑TGTTTTTAAGGACGAAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCAt                                    <  1:944114/147‑1 (MQ=255)
                                                                                                                       tgtACAGGGGCC‑C‑TGTTTTTTATGAGGAAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGttt                              <  2:710878/150‑1 (MQ=255)
                                                                                                                       tgtACAAGGGGC‑C‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGtttt                             >  1:484911/1‑151 (MQ=255)
                                                                                                                                ggggctttttttttAGGCCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaat                       <  2:132478/141‑1 (MQ=255)
                                                                                                                                      tgttttttATGCGGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCaatgaa                    <  1:966390/144‑1 (MQ=255)
                                                                                                                                         tttttAGGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGt                 <  1:665221/147‑1 (MQ=255)
                                                                                                                                                     cGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTaaa  >  1:204985/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
TTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGC‑C‑TGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAA  >  NC_000913/1622668‑1622964

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