Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	1,537,415	0	C	G	27.3%	29.9 / 10.7	22	A388G (GCG→GGG)	mdfA	Multidrug transporter MdfA
Reads supporting (aligned to +/- strand):  ref base C (16/0);  new base G (0/6);  total (16/6)
Fisher's exact test for biased strand distribution p-value = 1.34e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TCTCCAGCGAACAGAGCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCC‑T‑GGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCG  >  ACB122/1537283‑1537536                                                                                                                                       |                                                                                                                          tctcCAGCGAACAGAGCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCC‑T‑GGCGGCgggg                                                                                                                       >  7:684/1‑137 (MQ=255)   tcCAGCGAACAGAGCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCC‑T‑GGCGGCgggggg                                                                                                                     >  7:165755/1‑136 (MQ=255)   tcCAGCGAACAGAGCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCC‑T‑GGCGGCGGGGGt                                                                                                                     >  4:28385/1‑137 (MQ=255)             agagCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCC‑T‑GGCGGCGGGGGTGGCGGGGGCg                                                                                                           >  7:279604/1‑137 (MQ=255)                            tttCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCC‑T‑GGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCg                                                                                            >  3:254997/1‑137 (MQ=255)                            tttCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTCCCAGCTTGCTGCATTTTGCCC‑T‑GGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCg                                                                                            >  3:259118/1‑137 (MQ=255)                                   cATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCC‑T‑GGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGt                                                                                     >  5:272754/1‑137 (MQ=255)                                             atgatCACGATTGCCTTGCTGGGGTTTGGAGGATCGCTGCTGGCGATGATTGGCGCAGGTGCAAGCTTGCTGCATTTTGCCC‑T‑GGCGGCGGGGGTGTCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTggtggg                                                                           >  7:60636/1‑137 (MQ=255)                                                   aCGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCC‑T‑GGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGggcg                                                                     >  3:245146/1‑137 (MQ=255)                                                         gCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCC‑T‑GGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTgg                                                               >  7:180606/1‑137 (MQ=255)                                                                  gggTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCTAGCTTGCTGCATTTTGCCC‑T‑GGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAggggg                                                      >  1:77356/1‑136 (MQ=255)                                                                  gggTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCC‑T‑GGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGc                                                      >  8:253729/1‑137 (MQ=255)                                                                  gggTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCC‑T‑GGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGc                                                      >  6:217467/1‑137 (MQ=255)                                                                                gctgctGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCC‑T‑GGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAAt                                        >  5:211569/1‑137 (MQ=255)                                                                                       gCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCC‑T‑GGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTAtt                                 >  4:34694/1‑137 (MQ=255)                                                                                                          ccAGCTTGCTGCATTTTGCCC‑T‑GGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTcg              >  1:207191/1‑137 (MQ=255)                                                                                                                 gcggctttttCCCC‑G‑GGGGGGGGGGGGGGGGGGGGCGTGGGCGCGGTGGCCGGTGGGGGTTGGGGGGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTaa       <  8:62208/131‑1 (MQ=255)                                                                                                                  ggtctttttGCCCGG‑GGGGGGGGGGGGGGGGGGGGGGTGGGCGCGGGGGCCGGTGTGGTTTGTGGGGGGGGGGGGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTaa       <  3:135001/132‑1 (MQ=255)                                                                                                                  ggggaTTTTGCCCGG‑GGGGGGGGGGGGGGGGGGGGCGTGGGCGCTGTGGCGGCTGTGGTTTGGGGGGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTaa       <  5:166093/133‑1 (MQ=255)                                                                                                                   gcatttttGCCT‑G‑GGGGGGGGGGGGGGGGGGGGGGTGGGCGCTGTGGCCGGTGGGGTTGGGGGTGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAAc     <  7:251105/133‑1 (MQ=255)                                                                                                                     gaTTTTGCCCGG‑GGGGGGGGGGGGGGGGGGGGGGTGGGCGCGGGGGCCGCGGGGGTTTGGGGGGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAAcc    <  3:28385/136‑1 (MQ=255)                                                                                                                       ttttGCCC‑TGGGGGGGGGGGGGGGGGGGGGGGTTGGGGGTGTGGCGGGGGTGGTTTGTGGGGGTGGGGCGGCGTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCg  <  2:262285/137‑1 (MQ=255)                                                                                                                                       |                                                                                                                          TCTCCAGCGAACAGAGCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTGCCC‑T‑GGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCG  >  ACB122/1537283‑1537536

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: ‑