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

TGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTT  >  NC_000913/1622671‑1622950
                                                                                                                                                   |                                                                                                                                      
tGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGaaaa                                                                                                                                    >  2:55327/1‑150 (MQ=255)
              ataatTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTaa                                                                                                                      >  1:144484/1‑150 (MQ=255)
              ataatTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTaa                                                                                                                      >  1:380147/1‑150 (MQ=255)
                    tttttttctgtttttcgggCCGCTTTTTCCAGCCAAAAATGGGATGATAAATGAACGGAGATAACCCCTCCCCTACCCGCCCCTTTGTTACATTTGTGTACAGGGGGCTTGTTTTTTAGGACGGCAAAAAAAAACCGCCAGTAAACcggc                                                                                                                <  2:604968/147‑1 (MQ=255)
                       ttttttgtttttcgggCCGTATTTTACAGCCAAAAAGGGGATAACTAATGAAGGGAAATAACCCCCCCCCAACCCGCCCCTTTGTTACTGTGGTGTAAAGGGGGCTTTTTTTTTATGCCGCAAAAAAAAAACCGCCAGTAAACCGGCGGt                                                                                                             <  1:260428/150‑1 (MQ=255)
                               ttttctggCCGCTTTTTTCAGCCAAAAAGGGGATGACAAATGAACGGAAATAACCCCTCCCCTACCCGCCCCTTTGTTACAGTTGTGTAAAAGGGGCTTTTTTTTTATGCCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAAt                                                                                                         <  2:543073/146‑1 (MQ=255)
                               ttttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCtt                                                                                                     >  1:271503/1‑150 (MQ=255)
                               tgttctgcCCGCTTTTTCCACCCAAAAAGGGAATGATAAAGGAAGGGAAATACCCCCCCCCCAACCCGCCCCTTTTTTACATTTGTGAACAGGGGGCTTTTTTTTTAGGCGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCt                                                                                                      <  1:238075/147‑1 (MQ=255)
                                 ttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGc                                                                                                   >  1:155121/1‑150 (MQ=255)
                                 ttctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGc                                                                                                   >  1:539197/1‑150 (MQ=255)
                                  tctggCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCa                                                                                                  >  2:21634/1‑150 (MQ=255)
                                       ccGCGTTTTCCACCCAAAAAGGGGATAATTAATGAAGGGAAATACCCCCTCACCAACCCCCCCCTTTTTTACAGTGTTGAAAAGGGGGCGTTTTTTTTAGGCGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGAt                                                                                             <  1:129594/150‑1 (MQ=255)
                                                  aGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGttt                                                                                 >  2:187209/1‑151 (MQ=255)
                                                        aaaGGGGATGATTAAGGAACGGAAATACCCCCTCACCTACCCGCCCCTTTTTTACGGTTGTGAAAAGGGGGCTTGATTTTTATGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGttt                                                                                 <  1:28627/145‑1 (MQ=255)
                                                                         aagggAGATAACCCCCCCCCAACCCGGCCCTTTTTTACAGTTGTTTACAGGGGGCTTTTTTTTTAGGGCGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCAt                                                            <  2:256653/146‑1 (MQ=255)
                                                                                  accccccccccTACCCGCCCCTTTTTTACAGTTGTGTAAAGGGGGCTTGTTTTTTAGGACGGCAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACt                                                 <  2:485158/142‑1 (MQ=255)
                                                                                   aTCCCCCCCCCAACCCGCCCCTTTTTTACATTGTGGACAAGGGGGCTTTTTTTTTAGGAGGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTcc                                                       <  1:74584/145‑1 (MQ=255)
                                                                                          aCCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAAACGGCGGGGAATGCTTGCATGGGTAGATTTGTGTTTTTCTTTTACGCTAACAGGCATTTTCCTGCACTGATaa                                            >  1:212352/1‑148 (MQ=255)
                                                                                                         ttGTTACAGTTGTGTACAGGGGGCTTGTTTTTTATGACGGAAAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAg                            <  1:55321/150‑1 (MQ=255)
                                                                                                         ttGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGt                           >  2:159228/1‑151 (MQ=255)
                                                                                                              aCAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCAt                      >  1:347463/1‑151 (MQ=255)
                                                                                                                       tACAAGGAGCCTGATTTTTATGACGGCGAAAAAAAACCGGCAGGAAACCGGCGGTGAATGGTTGCATGGCTGGATTTAAGTTTTGCTTTTACGCCAACAGGCATTTTCCTGCACTGATAACGAATCGTTTACAAAGTAGCCTCAGTTTtc              >  1:184226/1‑150 (MQ=255)
                                                                                                                                   tGATTTTTATGAAGGAGAAAAAAAAAGGCCAGTAAACCGGGGGTGAATGTTTGGATGGATAGATTTGTGTTTTGGTTTTACGATAAAAGGGATTTTTCTGGAATGATAACGAATCGTTGACACAGTAGCATCAGTTTTTTTAATGAATGtt  >  1:66904/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                      
TGGCCGTACTGGGAATAATTTTATTTTCTGGTTTTCTGGCCGCGTATTTCAGCCACAAATGGGATGACTAATGAACGGAGATAATCCCTCACCTAACCGGCCCCTTGTTACAGTTGTGTACAAGGGGCCTGATTTTTATGACGGCGAAAAAAAACCGCCAGTAAACCGGCGGTGAATGCTTGCATGGATAGATTTGTGTTTTGCTTTTACGCTAACAGGCATTTTCCTGCACTGATAACGAATCGTTGACACAGTAGCATCAGTTTTCTCAATGAATGTT  >  NC_000913/1622671‑1622950

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