Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,617,914	0	C	G	26.3%	27.2 / 12.6	19	intergenic (‑127/+46)	PP_1417/PP_1418	tricarboxylate transport protein TctB/tricarboxylate transport protein TctC
Reads supporting (aligned to +/- strand):  ref base C (13/1);  new base G (5/0);  total (18/1)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.91e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GATGCGATCTGTACGAAGCATCGTATCGGCTGTGCCAGGGGCCGCTTCGCGCCCGATCCCGACACAAGGCCGCTCCCACAGGAGTGGTGGTGAGTCCGACAGGCTTACTGGATCAGGCCGAATTCCTTGGCCAGCGCCTTGTAGTCAGCCACTTGCTTCTTCACATAGCCGTCCAGCTCTTCGCCGGTC  >  NC_002947/1617856‑1618044                                                           |                                                                                                                                   gATGCGATCTGTACGAAGCATCGTATCGGCTGTGCCAGGGGCCGCTTCGCGCCCGATCCCGACACAAGGCCGCTCCCACAGGAGTGGTGGTGAGTCCGACAGGCTTACTGGATCAGGCCGAATTCCTTGGCCAgcg                                                       >  2:274075/1‑136 (MQ=255)                    aTCGTATCGGCTGTGCCAGGGGCCGCTTCGCGCCCGATCCCGACACAAGGCCGCTCCCACAGGAGTGGTGGTGAGTCCGACAGGCTTACTGGATCAGGCCGAATTCCTTGGCCAGCGCCTTGTAGTCAGCCACTTg                                    >  1:26244/1‑136 (MQ=255)                         aTCGGCTGTGCCAGGGGCCGCTTCGCGCCCGATCCCGACACAAGGCCGCTCCCACAGGAGTGGTGGTGAGTCCGACAGGCTTACTGGATCAGGCCGAATTCCTTGGCCAGCGCCTTGTAGTCAGCCACTTGcttct                               >  1:414397/1‑136 (MQ=255)                         aTCGGCTGTGCCAGGGGCCGCTTCGCGCCCGATCCCGACACAAGGCCGCTCCAACAGGAGTGGTGGTGAGTCCGACAGGCTTACTGGATCAGGCCGAATTACTTGGCCAGCGCCTTGTAGTCAGCCACTTGcttct                               >  1:92991/1‑136 (MQ=255)                                gtgCCAGGGGCCGCTTCGCGCCCGATCCCGACACAAGGCCGCTCCCACAGGAGTGGTGGTGAGTCCGACAGGCTTACTGGATCAGGCCGAATTCCTTGGCCAGCGCCTTGTAGTCAGCCACTTGCTTCTTCACATa                        >  1:161784/1‑136 (MQ=255)                                gtgCCAGGGGCCGCTTCGCGCCCGATCCCGACACAAGGCCGCTCCCACAGGAGTGGTGGTGAGTCCGACAGGCTTACTGGATCAGGCCGAATTCCTTGGCCAGCGCCTTGTAGTCAGCCACTTGCTTCTTCACATa                        >  1:412181/1‑136 (MQ=255)                                          cagCTTCGCGCCCGATCGCGACGCAAGGCCGCTCCCACAGGAGTGGTGGTGAGTCCGACAGGCTTACTGGATCAGGCCGAATTCCTTGGCCAGCGCCTTGTAGTCAGCCACTTGCTTCTTCACATAGCCGTCCAGc              >  1:283075/3‑136 (MQ=255)                                              ttCGCGCCCGATCGCGACGCAAGGCCGCTCCCACAGGAGTGGTGGTGAGTCCGACAGGCTTACTGGATCAGGCCGAATTCCTTGGCCAGCGCCTTGTAGTCAGCCACTTGCTTCTTCACATAGCCGTCCAGCTCtt          >  1:132142/1‑136 (MQ=255)                                              ttCGCGCCCGATCCCGACACAAGGCCGCTCCCACAGGAGTGGTGGTGAGTCCGACAGGCTTACTGGATCAGGCCGAATTCCTTGGCCAGCGCCTTGTAGTCAGCCACTTGCTTCTTCACATAGCCGTCCAGCTCtt          >  1:235425/1‑136 (MQ=255)                                              ttCGCGCCCGATCCCGACACAAGGCCGCTCCCACAGGAGTGGTGGTGAGTCCGACAGGCTTACTGGATCAGGCCGAATTCCTTGGCCAGCGCCTTGTAGTCAGCCACTTGCTTCTTCACATAGCCGTCCAGCTCtt          >  1:440796/1‑136 (MQ=255)                                              ttCGCGCCCGATCCCGACACAAGGCCGCTCCCACAGGAGTGGTGGTGAGTCCGACAGGCTTACTGGATCAGGCCGAATTCCTTGGCCAGCGCCTTGTAGTCAGCCACTTGCTTCTTCACATAGCCGTCCAGCTCtt          <  2:26244/136‑1 (MQ=255)                                               tCGCGCCCGATCGCGACGCAAGGCCGCTCCCACTGGAGTGGTGGTGAGTCCGACAGGCTTACTGGATCAGGCCGTATTGCTTGGCCAGCGCGTTGTAGTCGGCCACTTGCTTCTTCACATAGCCGTCCAGCTCTTc         >  2:107609/1‑136 (MQ=38)                                               tCGCGCCCGATCGCGACGCAAGGCCGCTCCCACAGGAGTGGTGGTGAGTCCGACAGGCTTACTGGATCAGGCCGAATTCCTTGGCCAGCGCCTTGTAGTCAGCCACTTGCTTCTTCACATAGCCGTCCAGCTCTTc         >  1:217932/1‑136 (MQ=255)                                               tCGCGCCCGATCGCGACGCAAGGCCGCTCCCACAGGAGTGGTGGTGAGTCCGACAGGCTTACTGGATCAGGCCGAATTCCTTGGCCAGCGCCTTGTAGTCAGCCACTTGCTTCTTCACATAGCCGTCCAGCTCTTc         >  1:11186/1‑136 (MQ=255)                                               tCGCGCCCGATCCCGACACAAGGCCGCTCCCACAGGAGTGGTGGTGAGTCCGACAGGCTTACTGGATCAGGCCGAATTCCTTGGCCAGCGCCTTGTAGTCAGCCACTTGCTTCTTCACATAGCCGTCCAGCTCTTc         >  1:74740/1‑136 (MQ=255)                                               tCGCGCCCGATCCCGACACAAGGCCGCTCCCACAGGAGTGGTGGTGAGTCCGACAGGCTTACTGGATCAGGCCGAATTCCATGGCCAGCGCCTAGTAGTCAGCCACTTCCTTCTTCACATAGCCGTCCAGCTCtta         >  1:373371/1‑135 (MQ=255)                                                cgcgCCCGATCCAGACACAAGGCCGCTCCCACAGGAGTGGTGGTGAGTCCGACAGGCTTACTGGATCAGGCCGAATTCCTTGGCCAGCGCCTTGTAGTCAGCCACTTGCTTCTTCACATAGCCGTCCAGCTCTTCg        >  2:398870/1‑136 (MQ=255)                                                      cGATCCCGACACAAGGCCGCTCCCACAGGAGTGGTGGTGAGTCCGACAGGCTTACTGGATCAGGCCGAATTCCTTGGCCAGCGCCTTGTAGTCAGCCACTTGCTTCTTCACATAGCCGTCCAGCTCTTCGCCGGTc  >  2:80925/1‑136 (MQ=255)                                                      cGATCCCGACACAAGGCCGCTCCCACAGGAGTGGTGGTGAGTCCGACAGGCGTACAGGATCAGGCCGAATTCCTTGGCCAGCGCCTTGTAGTCAGCCACTTGCTTCATCACATAGCCGTCCAGCTCTTCGCCCGTc  >  2:290526/1‑136 (MQ=255)                                                           |                                                                                                                                   GATGCGATCTGTACGAAGCATCGTATCGGCTGTGCCAGGGGCCGCTTCGCGCCCGATCCCGACACAAGGCCGCTCCCACAGGAGTGGTGGTGAGTCCGACAGGCTTACTGGATCAGGCCGAATTCCTTGGCCAGCGCCTTGTAGTCAGCCACTTGCTTCTTCACATAGCCGTCCAGCTCTTCGCCGGTC  >  NC_002947/1617856‑1618044

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