Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	AM260480	1,842,557	0	C	G	27.8%	28.2 / 12.4	18	R721R (CGG→CGC)	h16_B1632	conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (7/6);  new base G (3/2);  total (10/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 = 3.44e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGCGGACGCGTTCCGGCTTGGCGCCGAGCGCCGTGCGGATCCGCTCCAGCATTGCCGTCTGCTCGACAGTGGGCCTGGCGCCAAGCAGGCGCTCATCGACACGCAGCTTCTCCGCCAGCGCCTGCAGCCGCTGGCGGTCGGCCTGATACCTGAGCAAGGCCGGCAGCGTCGCGAGCGCCTGTTCGGGCGCGTACCGCGTGATCAGCTCCTGCCTGCCGCGGATCTCGCGCCAGGCGTCCGGCTCAAGTTC  >  AM260480/1842430‑1842679                                                                                                                                |                                                                                                                           cgcgGACGCGTTCCGGCTTGGCGCCGAGCGCCGTGCGGATCCGCTCCAGCATTGCCGTCTGCTCGACAGTGGGCCTGGCGCCAAGCAGGCGCTCATCGACACGCAGCTTCTCCGCCAGCGCCTGCAGCCGCTgg                                                                                                                      <  1:69641/134‑1 (MQ=255)    ggACGCGTTCCGGCTTGGCGCCGAGCGCCGTGCGGATCCGCTCCAGCATTGCCGTCTGCTCGACAGTGGGCCTGGCGCCAAGCAGGCGCTCATCGACACGCAGCTTCTCCGCCAGCGCCTGCAGc                                                                                                                            <  2:837112/125‑1 (MQ=255)    ggACGCGTTCCGGCTTGGCGCCGAGCGCCGTGCGGATCCGCTCCAGCATTGCCGTCTGCTCGACAGTGGGCCTGGCGCCAAGCAGGCGCTCATCGACACGCAGCTTCTCCGCCAGCGCCTGCAGc                                                                                                                            >  1:837112/1‑125 (MQ=255)                               ccGTGCGGATCCGCTCCAGCATTGCCGTCTGCTCGACAGTGGGCCTGGCGCCAAGCAGGCGCTCATCGACACGCAGCTTCTCCGCCAGCGGCTGCAGGCGCTGGCGGTCGGCCTGATACCTGAGCAAGGCCgg                                                                                         >  1:212907/1‑133 (MQ=255)                               ccGTGCGGATCCGCTCCAGCATTGCCGTCTGCTCGACAGTGGGCCTGGCGCCAAGCAGGCGCTCATCGACACGCAGCTTCTCCGCCAGCGCCTGCAGCCGCTgg                                                                                                                      >  1:775512/1‑104 (MQ=255)                               ccGTGCGGATCCGCTCCAGCATTGCCGTCTGCTCGACAGTGGGCCTGGCGCCAAGCAGGCGCTCATCGACACGCAGCTTCTCCGCCAGCGCCTGCAGCCGCTgg                                                                                                                      <  2:775512/104‑1 (MQ=255)                               ccGTGCGGATCCGCTCCAGCATTGCCGTCTGCTCGACAGTGGGCCTGGCGCCAAGCAGGCGCTCATCGACACGCAGCTTCTCCGCCAGCGCCTGCAGCCGCTGGc                                                                                                                     >  1:46967/1‑105 (MQ=255)                               ccGTGCGGATCCGCTCCAGCATTGCCGTCTGCTCGACAGTGGGCCTGGCGCCAAGCAGGCGCTCATCGACACGCAGCTTCTCCGCCAGCGCCTGCAGCCGCTGGc                                                                                                                     <  2:46967/105‑1 (MQ=255)                                                                cGACAGTGGGCCTGGCGCCAAGCAGGCGCTCATCGACACGCAGCTTCTCCGCCAGCGGCTGCAGGCGCTGGCGGTCGGCCTGATACCTGAGCAAGGCCGGCAGCGTCGCGAGCGCCTGtt                                                                     <  1:625724/120‑1 (MQ=255)                                                                cGACAGTGGGCCTGGCGCCAAGCAGGCGCTCATCGACACGCAGCTTCTCCGCCAGCGGCTGCAGGCGCTGGCGGTCGGCCTGATACCTGAGCAAGGCCGGCAGCGTCGCGAGCGCCTGtt                                                                     >  2:625724/1‑120 (MQ=255)                                                                                 ccAAGCAGGCGCTCATCGACACGCAGCTTCTCCGCCAGCGCCTGCAGCCGCTGGCGGTCGGCCTGATACCTGAGCAAGGCCGGCAGCGTCGCGAGCGCCTGTTCGGGCGCGTACCGCGTGATCAGCTcctgcct                                      >  1:1018736/1‑134 (MQ=255)                                                                                                           cTTCTCCGCCAGCGCCTGCAGCCGCTGGCGGTCGGCCTGATACCTGAGCAAGGCCGGCAGCGTCGCGAGCGcc                                                                         <  2:108874/73‑1 (MQ=255)                                                                                                           cTTCTCCGCCAGCGCCTGCAGCCGCTGGCGGTCGGCCTGATACCTGAGCAAGGCCGGCAGCGTCGCGAGCGcc                                                                         >  1:108874/1‑73 (MQ=255)                                                                                                           cTTCTCCGCCAGCGCCTGCAGCCGCTGGCGGTCGGCCTGATACCTGAGCAAGGCCGGCAGCGTCGCGAGCGCCTGTTCGGGCGCGTACcgcg                                                      >  1:667084/1‑92 (MQ=255)                                                                                                           cTTCTCCGCCAGCGCCTGCAGCCGCTGGCGGTCGGCCTGATACCTGAGCAAGGCCGGCAGCGTCGCGAGCGCCTGTTCGGGCGCGTACcgcg                                                      <  2:667084/92‑1 (MQ=255)                                                                                                             tctcCGCCAGCGCCTGCAGCCGCTGGCGGTCGGCCTGATACCTGAGCAAGGCCGGCAGCGTCGCGAGCGcc                                                                         <  2:351669/71‑1 (MQ=255)                                                                                                             tctcCGCCAGCGCCTGCAGCCGCTGGCGGTCGGCCTGATACCTGAGCAAGGCCGGCAGCGTCGCGAGCGcc                                                                         >  1:351669/1‑71 (MQ=255)                                                                                                                 cGCCAGCGGCTGCAGGCGCTGGCGGTCGGCCTGATACCTGAGCAAGGCCGGCAGCGTCGCGAGCGCCTGTTCGGGCGCGTACcgcg                                                      >  1:400049/1‑86 (MQ=255)                                                                                                                 cGCCAGCGGCTGCAGGCGCTGGCGGTCGGCCTGATACCTGAGCAAGGCCGGCAGCGTCGCGAGCGCCTGTTCGGGCGCGTACcgcg                                                      <  2:400049/86‑1 (MQ=255)                                                                                                                     aGCGCCTGCAGCCGCTGGCGGTCGGCCTGATACCTGAGCAAGGCCGGCAGCGTCGCGAGCGCCTGTTCGGGCGCGTACCGCGTGATCAGCTCCTGCCTGCCGCGGATCTCGCGCCAGGCGTCCGGCTCAAGTTc  >  1:77136/1‑134 (MQ=255)                                                                                                                                |                                                                                                                           CGCGGACGCGTTCCGGCTTGGCGCCGAGCGCCGTGCGGATCCGCTCCAGCATTGCCGTCTGCTCGACAGTGGGCCTGGCGCCAAGCAGGCGCTCATCGACACGCAGCTTCTCCGCCAGCGCCTGCAGCCGCTGGCGGTCGGCCTGATACCTGAGCAAGGCCGGCAGCGTCGCGAGCGCCTGTTCGGGCGCGTACCGCGTGATCAGCTCCTGCCTGCCGCGGATCTCGCGCCAGGCGTCCGGCTCAAGTTC  >  AM260480/1842430‑1842679

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