Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,622,8160GA57.1% ‑6.3 / 12.6 14intergenic (+75/+144)yneM/dgcZinner membrane‑associated protein/diguanylate cyclase, zinc‑sensing
Reads supporting (aligned to +/- strand):  ref base G (6/0);  new base A (0/8);  total (6/8)
Fisher's exact test for biased strand distribution p-value = 3.33e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.05e-03

GTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGG  >  NC_000913/1622676‑1622956
                                                                                                                                              |                                                                                                                                            
gTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAAc                                                                                                                                     >  2:34246/1‑150 (MQ=255)
                                 cccGCTTATTTCACCCAAAAGGGGGTTGACAAATGAAGGGAAATAACCCCTCCCCTAACCGCCCCTTTTTTACATTTGTGTAAAGGGGGCTTGTTTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGAt                                                                                                   <  2:98429/150‑1 (MQ=255)
                                      tttttccAGCCAAAAAGGGGATGACTAATGAACGGAAATAATCCCCCCCCTACCCGCCCCCTTGTTACAGTTGTGTACAGGGGGCCTATTTTTTAGGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAtt                                                                                              <  1:109813/145‑1 (MQ=255)
                                       tttttCAGCCAAAAGGGGAAATCTAATTGAAGGGATATAATCCCTCCCCAACCCGCCCCCTTTTGGAGGTGTGTTACAGGGGGCTTTATTTTTGTGAGGGAAAAAAAAAGCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGAtt                                                                                              <  2:164853/148‑1 (MQ=255)
                                                        ggATGATTAATGAAGGGAGATAACCCCCCCCCAACCCGGCCCCTTTTTACATTTGTGAACAGGGGGCCTGTTTTTTATGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACg                                                                             <  1:239872/150‑1 (MQ=255)
                                                              cTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAg                                                                      >  1:270020/1‑151 (MQ=255)
                                                                           gaTAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTg                                                           >  1:68762/1‑150 (MQ=255)
                                                                             aaaccccTCCCCTACCCGCCCCCTTTTAACTGTTGTTTACAGGGGGCTTGTTTTTTATGACGGAAAAAAAAACCCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCAc                                                        <  2:90260/147‑1 (MQ=255)
                                                                                tCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGa                                                     >  2:41467/1‑151 (MQ=255)
                                                                                        caacccGCCCCTTTGTTACATTTGTGAAAAGGGGGCTTTTTTTTAAGGCCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAAt                                             <  2:46493/147‑1 (MQ=255)
                                                                                                 cccTTGTTACAGTGGTGTAAAGGGGGCTTGTTTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGacac                                    <  2:68762/151‑1 (MQ=255)
                                                                                                            gTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAg                         >  2:283217/1‑151 (MQ=255)
                                                                                                            gTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAg                         >  2:30160/1‑151 (MQ=255)
                                                                                                                                    ttATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACgg  <  2:169867/151‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                            
GTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTTAAACGG  >  NC_000913/1622676‑1622956

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: