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

TTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAA  >  NC_000913/1622668‑1622964
                                                                                                                                                      |                                                                                                                                                    
ttATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGCTAATCCCTCACCTAACCGGCCCCTTGTTAAAGTTGGGTACAAGGGGCCTGATTTTTATAACGGCGa                                                                                                                                                     >  1:915119/1‑150 (MQ=255)
ttATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGa                                                                                                                                                     >  2:573020/1‑150 (MQ=255)
  aTGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGaaa                                                                                                                                                   >  2:158264/1‑150 (MQ=255)
    ggCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGaaaaa                                                                                                                                                 >  2:426570/1‑150 (MQ=255)
     gCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGaaaaaa                                                                                                                                                >  1:459843/1‑150 (MQ=255)
          aCTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCg                                                                                                                                           <  1:136018/150‑1 (MQ=255)
                        ttttttttgtttttctggCCGCGTTTTCCACCCAAAAAGGGAATGACAAAGAAGCGGATAAAACCCCCCACCAACCCGGCCTTTTTTCATTGTAGTGTAAAGGGGCTGTATTTTTTGTGCGGGCAAAAAAAAACCGCCAGTAAACcggcg                                                                                                                             <  2:607327/148‑1 (MQ=255)
                        tattttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACcggcgg                                                                                                                            >  1:774454/1‑151 (MQ=255)
                         attttCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACcggcgg                                                                                                                            >  1:202241/1‑150 (MQ=255)
                             tctggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAAt                                                                                                                       >  1:132208/1‑151 (MQ=255)
                               tggttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATg                                                                                                                      >  1:607291/1‑150 (MQ=255)
                                 ttttcctggCCGCGTATTTCACCCAAAAGTGGAATGATTAATGAACGGATATACTCCCTCGCCTACCCGCCACCTTGAAACATTTATTCAAAAGCGCCCTGTTTTTTATGCCGGAAAAAAAAACCCGCCAGTAAACCGGCGGTGAAt                                                                                                                       <  1:403872/146‑1 (MQ=255)
                                    ttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGc                                                                                                                 >  1:879419/1‑150 (MQ=255)
                                     tctggCCGCTTTTTTCACCCAAAAAGGGGATGACTAAGGAAGGGAAAAACCCCCTCCCCTACCCGCCCCTTTTTTACATTTGTAAAAAGGGGCCTTTTTTTTTAGGCCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCa                                                                                                                <  2:283967/150‑1 (MQ=255)
                                         gCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGAt                                                                                                           >  2:964597/1‑151 (MQ=255)
                                              tttttacggCCAAAAAGGGGATAAATAATGAACGAAAAAAACCCCTCCCCTACCCGCCCCTTTTTTACTGTTGTGAAAAGGGGGCTTGTTTTTTACGGGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATa                                                                                                          <  1:12653/139‑1 (MQ=255)
                                                 ttccggCCAAAAAGGGGAGAACAAAGGACGGGAAATACCCCCCCACCAAACCGCCCCCTTGTAAGAGTGGTGAAAGAGGGTCTTGTTTTTGGTGCAGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAttt                                                                                                     <  1:935595/144‑1 (MQ=255)
                                                  ttCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTtgtg                                                                                                  >  2:590588/1‑151 (MQ=255)
                                                          aaaaggggtaaaaaaaTGAAGGGAGAAACCCCCTCACCTACCCGCCCCTTTTTTAAGGTTGTGTAGAGGGGGTTTTTTTTTTGGGGGGGAAAAAAAAATCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCtt                                                                                          <  1:762800/137‑1 (MQ=255)
                                                           aaaGGGGATGAAAAGTAAACGGATATAACCCCTCCCAAACCCGCCCCCTTTTCACAGTTGTGAAAAGGGGCCTTGTTTTTTAAAAGGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCt                                                                                           <  1:203612/149‑1 (MQ=255)
                                                              tGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTAc                                                                                      >  1:116466/1‑151 (MQ=255)
                                                                   tGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAAc                                                                                >  2:277484/1‑151 (MQ=255)
                                                                   tGACTAATGAACGGAGATAATCCCTCACCTAACCGGCACCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACg                                                                                     >  1:463996/1‑147 (MQ=255)
                                                                                    tcccccctcccctcCCCGGCCTTTTGTTACGTTTGTGAAAAAGGGCCTTTTTTTTTATGAGGAAAAAAAAAACCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGc                                                                  <  1:639408/145‑1 (MQ=255)
                                                                                        cccTCCCCTACCCGCCCCTTTTTTACATTGGTGAACAGGGGGCTTTTTTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACt                                                               <  1:136961/148‑1 (MQ=255)
                                                                                              ccTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGa                                                       >  1:698773/1‑150 (MQ=255)
                                                                                              ccTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTAAGCTAACAGGCATTTTCCTGCACTGATAACGa                                                       >  1:456335/1‑150 (MQ=255)
                                                                                                    cGGCCCCTTTTTTCGGTTGTGAACAAGGGCTCTTTTTTTTAGGAGGGCAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGtt                                                 <  1:783831/149‑1 (MQ=255)
                                                                                                    cGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGGTAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTg                                                >  2:284288/1‑151 (MQ=255)
                                                                                                     ggCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGGGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGATAACAGGCATTTTCCTGCACTGATAACGAATCg                                                   >  1:134191/1‑147 (MQ=255)
                                                                                                     gcccctTTTTTAAATTTGTGAACAGGGGGCTTGATTTTTATGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGtt                                                 <  1:716480/147‑1 (MQ=255)
                                                                                                           tttttACATTGGTTTACAAGGCGGCTTTTTTTTATGACGACAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGAcaca                                           <  1:421100/146‑1 (MQ=255)
                                                                                                                   gTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATc                                   >  1:769005/1‑149 (MQ=255)
                                                                                                                    ttGTGTAAAGGGGGCTTGTTTTTAAGGACGAAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCAt                                    <  1:944114/147‑1 (MQ=255)
                                                                                                                       tgtACAGGGGCCCTGTTTTTTATGAGGAAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGttt                              <  2:710878/150‑1 (MQ=255)
                                                                                                                       tgtACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGtttt                             >  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)
                                                                                                                                                      |                                                                                                                                                    
TTATGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGGAGCTTAAA  >  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: