Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	1,537,421	0	T	G	30.0%	40.2 / 14.6	30	V390G (GTG→GGG)	mdfA	Multidrug transporter MdfA
Reads supporting (aligned to +/- strand):  ref base T (19/2);  new base G (0/9);  total (19/11)
Fisher's exact test for biased strand distribution p-value = 3.84e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CGAACAGAGCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTG‑CC‑C‑TGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCG  >  ACB122/1537290‑1537554                                                                                                                                       |                                                                                                                                      ctaaCAGTGAAAGGGGCTGGTTTGGGCCTTGGTGGGTATTATCCCTATTGCCTTGCTGGGGTTTGGCGGAGCGCGGGTGGGTATGATTGTCGCGGGTGCCTGTTTGCTGCTTTTTG‑CC‑C‑TGGGGGGGGGGGTGgggg                                                                                                                                  <  6:36019/135‑1 (MQ=255)      agagCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTG‑CC‑C‑TGGCGGCGGGGGTGGCGGGGGCg                                                                                                                             >  5:344910/1‑137 (MQ=255)           aaGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTG‑CC‑C‑TGGCGGCGGGGGTGGCGGGGGCGTTGgc                                                                                                                        >  8:109022/1‑137 (MQ=255)           aaGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTG‑CC‑C‑TGGCGGCGGGGGTGGCGGGGGCGTTGgc                                                                                                                        >  2:393784/1‑137 (MQ=255)               ggCTGGTTTCGGCCATGCGGGGTATGATCACTATTGCCTTGCTGGGGTTTGCCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTG‑CC‑C‑TGGCGGCGGGGGGGGCGGGGGCGTTGGCGCtg                                                                                                                    <  4:162640/137‑1 (MQ=255)                gCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTG‑CC‑C‑TGGCGGCGGGGGTGGCGGGGGCGTTGGCGCtgt                                                                                                                   >  2:54614/1‑137 (MQ=255)                gCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTG‑CC‑C‑TGGCGGCGGGGGTGGCGGGGGCGTTGGCGCtgt                                                                                                                   >  2:61023/1‑137 (MQ=255)                gCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTG‑CC‑C‑TGGCGGCGGGGGTGGCGGGGGCGTTGGCGCtgt                                                                                                                   >  5:193419/1‑137 (MQ=255)                    gTTGCGGCCAGGGGGGGGTTGATCACGATTGCCTGGGGGGGGTTTGGCGGCGCGCGGCTGGCTTTGTTGGGGGCCGGTGCCAGGTTGGTGCTTTTTC‑CC‑C‑TGGGGGGGGGGGTGGGGGGGGCGTTGGCGCTGTGGcc                                                                                                               <  8:105139/137‑1 (MQ=255)                     tttCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTG‑CC‑C‑TGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCg                                                                                                              >  1:222786/1‑137 (MQ=255)                     tttCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTG‑CC‑C‑TGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCg                                                                                                              >  3:130555/1‑137 (MQ=255)                     tttCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTG‑CC‑C‑TGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCg                                                                                                              >  3:345966/1‑137 (MQ=255)                                        gatCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTG‑CC‑C‑TGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTggtgggtg                                                                                           >  3:411057/1‑137 (MQ=255)                                            aCGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTG‑CC‑C‑TGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGggcg                                                                                       >  1:149893/1‑137 (MQ=255)                                            aCGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTG‑CC‑C‑TGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGggcg                                                                                       >  6:72806/1‑137 (MQ=255)                                            aCGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTG‑CC‑C‑TGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGggcg                                                                                       >  4:314706/1‑137 (MQ=255)                                            aCGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTG‑CC‑C‑TGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGggcg                                                                                       >  1:121737/1‑137 (MQ=255)                                             cGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGAGCCGGTGCCAGCTTGCTGCATTTTG‑CC‑C‑TGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGggcgg                                                                                      >  2:366960/1‑137 (MQ=255)                                                  gCCTTGCTGGGGTTTGGCGGAGCCCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTG‑CC‑C‑TGGCGGAGGGGGTGGCGGGGGCGTTGGCTCTGTGGCCGCTGTGGATTGTGGTGGGTGGGCGGCCTgg                                                                                 >  8:374771/1‑137 (MQ=255)                                                  gCCTTGCTGGGGTTTGGCGCAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTG‑CC‑C‑TGGCGGCGTGGGTGGCGCGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTgg                                                                                 >  7:55817/1‑137 (MQ=255)                                                                                          gcgcCGGTGCCAGCTTGCTGCATTTTG‑CC‑C‑TGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGggccg                                         >  6:311763/1‑137 (MQ=255)                                                                                                 tGCCAGCTTGCTGCATTTTG‑CC‑C‑TGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGggccgggggggg                                  >  7:209111/1‑131 (MQ=255)                                                                                                           tgggtttttG‑CC‑CGGGGGGGGGGGGGGGGGGGGGGGGTGGGCGCTGGGGCCGCGGTGGTTTGTGGTGGGTGGGCGGCCGGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTaa                         <  1:61023/132‑1 (MQ=255)                                                                                                           ggtgaTTTTGCCC‑C‑GGGGGGGGGGGGGGGGGGGGGCGTGGGCGCGGGGGCCGGTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTaa                         <  7:109022/133‑1 (MQ=255)                                                                                                           ggtctttttG‑CCGC‑GGGGGGGGGGGGGCGGGGGGGTGTTGGCGCGGGGGCGGCTGGGGTTGGTGGTGGGGGGGCGGCGTGGTGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTaa                         <  5:69344/132‑1 (MQ=255)                                                                                                           gctgtttttG‑CC‑CGGGGGGGGGGGGGGGGGGGGGGCGTGGGCGCGGGGGCCGCTGTGGTTTGGGGGGGGTGGGCGGCCGGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTaa                         <  2:228015/132‑1 (MQ=255)                                                                                                           gccgaTTTTG‑CC‑CGGGGGGGGGGGGGGGGGGGGGGGGTGGGCGCGGTGGCCGCTGTGGTTTGTGGGGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTaa                         <  3:314706/133‑1 (MQ=255)                                                                                                              taTTTTG‑CC‑CGGGGGGGGGGGGGGGGGGGGGGGGTGGGCGCGGGGGCCGGGGGGGTTTGTGGGGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAAcc                      <  5:242092/136‑1 (MQ=255)                                                                                                               ttttgG‑CC‑G‑TGGGGGGGGGGGGGGGGGGGGCGGTGGCGGTGGGGCGGCGGGGGTTTGTGGGGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCg                    <  8:128316/136‑1 (MQ=255)                                                                                                                       cccgtGGGGGGGGGGTGGGCGGGGGGTTTGGCGGGGGGGCGGGTGGGGTTTGTGGTGGTGGCGGGGCCTGGGGACGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCt               <  2:173359/133‑1 (MQ=19)                                                                                                                                 gggggggggggggggCGGTGGCGCGGGGGCCGCGGTGGTTTTGGGTGGGGGGGCGGCCGGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGgcc     <  6:188252/136‑1 (MQ=37)                                                                                                                                    ggggggCGGGGGGGTGGGCGCTGTGGCGGCTTGGGTTTGTGGTGGGTGGGCGGCGTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGgccgcg  <  7:103203/137‑1 (MQ=38)                                                                                                                                       |                                                                                                                                      CGAACAGAGCAAGGGGCTGGTTTCGGCCATGCTGGGGATGATCACGATTGCCTTGCTGGGGTTTGGCGGAGCGCTGCTGGCGATGATTGGCGCCGGTGCCAGCTTGCTGCATTTTG‑CC‑C‑TGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCG  >  ACB122/1537290‑1537554

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