Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,780,0960TG53.3% 1.4 / 10.9 15intergenic (+55/+51)PP_1587/dapCNhaP‑type Na+(K+)/H+ antiporter/N‑succinyl‑L,L‑diaminopimelate aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (3/4);  new base G (8/0);  total (11/4)
Fisher's exact test for biased strand distribution p-value = 2.56e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-02
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGG  >  NC_002947/1779954‑1780243
                                                                                                                                              |                                                                                                                                                   
tCGTCAGCACCTGGTGGGGGAGGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGaaa                                                                                                                                                <  6:41298/148‑1 (MQ=255)
   tCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGGACAGCGCTGCGCGTGCAATCa                                                                                                                                             >  4:22193/1‑148 (MQ=255)
          cTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGt                                                                                                                                      <  1:70834/148‑1 (MQ=255)
              tggGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTAc                                                                                                                                    <  7:100403/146‑1 (MQ=255)
              tggGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTACCACAGGTACTGCGCTGCGCGTGAAATCAGCGGTGTGccg                                                                                                                                  >  7:114235/1‑147 (MQ=255)
                     tGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTAACCCACGTACGGCGCTGCGCGTGCAATGACGGCTGTACCTGAGGGAg                                                                                                                           >  3:119356/1‑148 (MQ=255)
                                  cAAGTACTGGTGGAGCTGTATATGAGGGCGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCCAGGCCCTTCGCGGGCGCGCCCGCTCCCACCGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCcgc                                                                                                              <  4:22673/148‑1 (MQ=38)
                                    aGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCGTGAAATCAGCGCTGTACCGGTGGGAGCGGGCGTGCCcgc                                                                                                              >  5:64679/1‑146 (MQ=255)
                                         tGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCCCTCCCACAAGGAATTCCCAGCGCGTGAAAACAACGCAGTACCAGGGGGAGCGGGCCTGCCCCCGAAcggg                                                                                                       >  2:51217/1‑148 (MQ=255)
                                           ggTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCCCAGGTACTGCGCTGCGCTTGTAATCAGCGCTGTACCTGGGGGAGCGGGCGTGCCCGCGGAGGGGtt                                                                                                     >  3:166134/1‑148 (MQ=255)
                                                  tGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCCCTCCCACAGGGAAAGCCCTGCGCTTGAAAACACCCCAGTACCGGTGGGCGCGGGCGTGCCCGCGCACGGCGTGATCacc                                                                                              >  8:34619/1‑146 (MQ=255)
                                                                gCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTAAAGCGCTGCGCGTGAAAACACCGCTGTACCTGTGGGCGCGGGCGTGCCCGCGAACGGGTTGATCACCCGTTTTGCacgcac                                                                                >  2:29313/1‑143 (MQ=37)
                                                                  gAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGCAGGGGTTGCTCAGCGGTTTTGCaggcaggc                                                                              >  3:62684/1‑148 (MQ=255)
                                                                        gggTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCCCAGGTACTGGGCGGGGCGTGTAATCAGCGCGGTACCTGTGGGAGCGGGCGTGCCCGCGCAGGGGGTGGTCAGCGGTTTTGGAGGAAGGCGCGCAt                                                                        >  8:115188/1‑148 (MQ=37)
                                                                                                              gcgcgcCCGCTCCCACAGGTACTGCGCTGCGCGTGAAATCAGCGCTGTACCGGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCaa                                  >  4:113668/1‑148 (MQ=255)
                                                                                                                                         tgcgctTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCAc       >  2:154264/1‑148 (MQ=255)
                                                                                                                                         tgcgctTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGAATTCGGCCAGCGGAGCAACCAGTGCCATGCGCACGCGCCCGGCAc       >  7:145972/1‑148 (MQ=255)
                                                                                                                                              ttGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCggg  >  5:186306/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
TCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGG  >  NC_002947/1779954‑1780243

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