Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,628,431	0	C	.	25.0%	40.0 / 10.1	16	coding (1044/1758 nt)	PP_4944	carbamoyltransferase
Reads supporting (aligned to +/- strand):  ref base C (6/6);  new base . (2/2);  total (8/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.87e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGGCGATCACGTCTTCGTTGGAGTACGACGGGCCGAGGTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAGGCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCAGCGCACAGCCGCCGGCGAAGGCCAGCTTGCCGGTTTCGCGCAGGGTGTC  >  NC_002947/5628322‑5628562                                                                                                              |                                                                                                                                    cggcgATCACGTCTTCGTTGGAGTACGACGGGCCGAGGTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAGGCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGAc                                                                                                           <  1:309618/136‑1 (MQ=255)     gATCACGTCTTCGTTGGAGTACGACGGGCCGAGGTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAGGCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAAc                                                                                                       <  2:486304/136‑1 (MQ=255)                        tACGACGGGCCGAGGTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAGGCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGAc                                                                                    >  2:63643/1‑136 (MQ=255)                                       tAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAGGCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGAtcttct                                                                     <  1:236607/136‑1 (MQ=255)                                                           gACCGGAACGCCACGGGCGTGGGAAACATAGGCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCAgc                                                 >  2:6552/1‑136 (MQ=255)                                                            aCCGGAACGCCACGGGCGTGGGAAACATAGGCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCAgcg                                                >  1:355864/1‑136 (MQ=255)                                                               ggAACGCCACGGGCGTGGGAAACATAGGCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGAGGATCTTCTGGTTGAGCTTGACGTTCAGCGcac                                             >  2:223260/1‑136 (MQ=255)                                                                                   aaCATAGGCAGCGGCGCCAACGGCAGT‑CCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCAGCGCACAGCCGCCGGCGAAGGCCAGCt                        >  1:49864/1‑136 (MQ=255)                                                                                   aaCATAGGCAGCGGCGCCAACGGCAGT‑CCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCAGCGCACAGCCGCCGGCGAAGGCCAGCt                        >  1:52541/1‑136 (MQ=255)                                                                                        aGGCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCAGCGCACAGCCGCCGGCGAAGGCCAGCTTGcc                    >  2:121772/1‑136 (MQ=255)                                                                                               ggcgCCAACGGCAGTCCCGGAGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCAGCGCACAGCCGCAGGAGAAGGCCAGCTTGCCGGTTTcg             <  1:545608/136‑1 (MQ=255)                                                                                                 cgcCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCAGCGCACAGCCGCCGGCGAAGGCCAGCTTGCCGGTTTcgcg           >  1:409411/1‑136 (MQ=255)                                                                                                         gCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCAGCGCACAGCCGCCGGCGAAGGCCAGCTTGCCGGTTTCGCGCAGGgtgt   <  2:545477/136‑1 (MQ=255)                                                                                                         gCAGT‑CCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCAGCGCACAGCCGCCGGCGAAGGCCAGCTTGCCGGTTTCGCGCAGGGTGTc  <  2:49864/136‑1 (MQ=255)                                                                                                         gCAGT‑CCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCAGCGCACAGCCGCCGGCGAAGGCCAGCTTGCCGGTTTCGCGCAGGGTGTc  <  2:52541/136‑1 (MQ=255)                                                                                                          cAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCAGCGCACAGCCGCCGGCGAAGGCCAGCTTGCCGGTTTCGCGCAGGGTGTc  <  1:121772/136‑1 (MQ=255)                                                                                                              |                                                                                                                                    CGGCGATCACGTCTTCGTTGGAGTACGACGGGCCGAGGTAGACGTGCTCCATCTTCTCGACCGGAACGCCACGGGCGTGGGAAACATAGGCAGCGGCGCCAACGGCAGTCCCGGCGTCGCCGGAGGCCGGCTGGACGAACAGCTCTTTCACATCCGGACGGGCGATGATCTTCTGGTTGAGCTTGACGTTCAGCGCACAGCCGCCGGCGAAGGCCAGCTTGCCGGTTTCGCGCAGGGTGTC  >  NC_002947/5628322‑5628562

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