Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,456,7931.T20.0% 32.0 / 10.0 15intergenic (+44/‑137)gapB/PP_2150glyceraldehyde‑3‑phosphate dehydrogenase/thiamine biosynthesis lipoprotein
Reads supporting (aligned to +/- strand):  ref base . (5/7);  new base T (2/1);  total (7/8)
Fisher's exact test for biased strand distribution p-value = 5.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.95e-01
Rejected as consensus: Frequency below/above cutoff threshold.

ACGAATTCGGTTATAGCTGCCAGGTGGTTCGTGTGATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGC  >  NC_002947/2456670‑2456933
                                                                                                                            |                                                                                                                                            
aCGAATTCGGTTATAGCTGCCAGGTGGTTCGTGTGATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAAtgtgtg                                                                                                                       >  8:68662/1‑147 (MQ=255)
               gCTGCCAGGTGGTTCGTGTGATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTa                                                                                                                                    <  3:140798/119‑1 (MQ=255)
               gCTGCCAGGTGGTTCGTGTGATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTa                                                                                                                                    >  4:140798/1‑119 (MQ=255)
                      ggtggtTCGTGTGATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTcgc                                                                                                 >  7:103639/1‑147 (MQ=255)
                       gtggtTCGTGTGATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTcgcg                                                                                                <  1:261329/147‑1 (MQ=255)
                        tggtTCGTGTGATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTcgcg                                                                                                <  8:103639/146‑1 (MQ=255)
                                        aGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTAc                                                                               <  3:133855/147‑1 (MQ=255)
                                          aCATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACa                                                                              >  3:229392/1‑147 (MQ=255)
                                                         aCCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGcc                                                                <  4:229392/146‑1 (MQ=255)
                                                           ccGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTg                                                            <  1:209967/147‑1 (MQ=255)
                                                                               aTTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCAGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGCGTAAACCCGCAACAAAAGCATTATCGCTGACCCAGTAGGACCGGTTTTACCcgc                                         >  8:63802/1‑146 (MQ=255)
                                                                                    cTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCCCTCCGACAGGATTATAGCTGGCCCGGTAGGAGCGGGTTTACCCGCGAagag                                   >  4:15406/1‑147 (MQ=39)
                                                                                             cAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAAACCCCAACAACAAGATTATAGCGGGCCCTGTAGGGGCGGGTTTACCCGCGAAAAGGCCAGCAc                           >  4:5295/1‑147 (MQ=255)
                                                                                                          ggAACCTTCGGGTTCCCGTTTTTTTTTATTCTTAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGaaa             <  3:15406/147‑1 (MQ=255)
                                                                                                                     gTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGTATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTgc  <  7:68662/147‑1 (MQ=255)
                                                                                                                            |                                                                                                                                            
ACGAATTCGGTTATAGCTGCCAGGTGGTTCGTGTGATGGAAGACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGC  >  NC_002947/2456670‑2456933

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