Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,086,717	C→T	80.8%	G164G (GGC→GGT)	yqjG →	predicted S‑transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,086,717	0	C	T	80.8%	79.7 / 11.6	57	G164G (GGC→GGT)	yqjG	predicted S‑transferase
Reads supporting (aligned to +/- strand):  ref base C (9/2);  new base T (16/30);  total (25/32)
Fisher's exact test for biased strand distribution p-value = 6.78e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.07e-01

CATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTCAATACCGCGTTTGATGCGCTGGGCGCGAAAGCGGGTGATTACTACCCTCCAGCCCTGCAAACAAAAATTGACG  >  minE/2086654‑2086766                                                                |                                                  cATCGTCAGCAACGAATCATCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa                                              <  1:1617730/69‑1 (MQ=39) cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTTGGTGCGaa                                              <  1:3262766/69‑1 (MQ=255) cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa                                              <  1:2822952/69‑1 (MQ=255) cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa                                              <  1:2073722/69‑1 (MQ=255) cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa                                              <  1:2139288/69‑1 (MQ=255) cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa                                              <  1:222357/69‑1 (MQ=255) cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa                                              <  1:2486261/69‑1 (MQ=255) cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa                                              <  1:2634181/69‑1 (MQ=255) cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa                                              <  1:2758427/69‑1 (MQ=255) cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa                                              <  1:27937/69‑1 (MQ=255) cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa                                              <  1:1032548/69‑1 (MQ=255) cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa                                              <  1:3107966/69‑1 (MQ=255) cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa                                              <  1:3281553/69‑1 (MQ=255) cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa                                              <  1:3332760/69‑1 (MQ=255) cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa                                              <  1:526156/69‑1 (MQ=255) cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa                                              <  1:542721/69‑1 (MQ=255) cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa                                              <  1:593182/69‑1 (MQ=255) cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa                                              <  1:767026/69‑1 (MQ=255) cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa                                              <  1:1927085/69‑1 (MQ=255) cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa                                              <  1:1056881/69‑1 (MQ=255) cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa                                              <  1:1309470/69‑1 (MQ=255) cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa                                              <  1:1606973/69‑1 (MQ=255) cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa                                              <  1:1005549/69‑1 (MQ=255) cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa                                              <  1:1646979/69‑1 (MQ=255) cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa                                              <  1:1664808/69‑1 (MQ=255) cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa                                              <  1:1836846/69‑1 (MQ=255) cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa                                              <  1:1876440/69‑1 (MQ=255) cATCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCCCGTTTGATGCGCTGGGTGCGaa                                              <  1:2772323/69‑1 (MQ=255)  aTCGTCAGCAACGAATCAGCAGAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa                                              <  1:2317907/68‑1 (MQ=255)                   agcagAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGAAAGCTGGAGATtac                                 >  1:153232/1‑64 (MQ=38)                   agcagAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGAAAGCTGGAGATtac                                 >  1:619579/1‑64 (MQ=38)                   agcagAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGAAAGCTGGAGATtac                                 >  1:71843/1‑64 (MQ=38)                   agcagAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGAAAGCTGGAGATtac                                 >  1:429786/1‑64 (MQ=38)                   agcagAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGAAAGCTGGAGATtac                                 >  1:1124460/1‑64 (MQ=38)                   agcagAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGAAAGCTGGAGATtac                                 >  1:3295588/1‑64 (MQ=38)                   agcagAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGAAAGCTGGAGATtac                                 >  1:1958383/1‑64 (MQ=38)                   agcagAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGAAAGCTGGAGATtac                                 >  1:752555/1‑64 (MQ=38)                   agcagAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGAAAGCTGGAGATtac                                 >  1:1469635/1‑64 (MQ=38)                   agcagAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGAAAGCTGGAGATtac                                 >  1:3092583/1‑64 (MQ=38)                   agcagAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGAAAGCTGGAGATtac                                 >  1:2228339/1‑64 (MQ=38)                   agcagAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGAAAGCTGGAGATtac                                 >  1:1878711/1‑64 (MQ=38)                   agcagAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGAAAGCTGGAGATtac                                 >  1:2093880/1‑64 (MQ=38)                   agcagAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGAAAGCTGGAGATtac                                 >  1:2573056/1‑64 (MQ=38)                   agcagAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGAAAGCTGGAGATtac                                 >  1:872820/1‑64 (MQ=38)                   agcagAAATCATACGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGAAAGCTGGAGATtac                                 >  1:2449748/1‑64 (MQ=38)                               acGCATGTTCAATACCGCGTTTGATGCGCTGGGTGCGaa                                              <  1:731982/38‑1 (MQ=37)                                cGCATGTTCAATACCGCGTTTGATGCGCTGGGCGCGAAAGCGGGTGATtact                                <  1:895850/52‑1 (MQ=255)                                cGCATGTTCAATACCGCGTTTGATGCGCTGGGCGCGAAAGCGGGTGATtact                                <  1:899736/52‑1 (MQ=255)                                cGCATGTTCAATACCGCGTTTGATGCGCTGGGCGCGAAAGCGGGTGATTACTACCCTCCAGCCCTGca                >  1:1072671/1‑68 (MQ=255)                                      ttCAATACCGCGTTTGATGCGCTGGGCGCGAAAGCGGGTGATTACTACCCTCCCGCCCTGCAAACaaaa         >  1:1102931/1‑69 (MQ=255)                                      ttCAATACCGCGTTTGATGCGCTGGGCGCGAAAGCGGGTGATTACTACCCTCCAGCCCTGCAAACaaaa         >  1:577242/1‑69 (MQ=255)                                      ttCAATACCGCGTTTGATGCGCTGGGCGCGAAAGCGGGTGATTACTACCCTCCAGCCCTGCAAACaaaa         >  1:1360394/1‑69 (MQ=255)                                      ttCAATACCGCGTTTGATGCGCTGGGCGCGAAAGCGGGTGATTACTACCCTCCAGCCCTGCAAACaaaa         >  1:2791111/1‑69 (MQ=255)                                      ttCAATACCGCGTTTGATGCGCTGGGCGCGAAAGCGGGTGATTACTACCCTCCAGCCCTGCAAACaaaa         >  1:2556926/1‑69 (MQ=255)                                      ttCAATACCGCGTTTGATGCGCTGGGCGCGAAAGCGGGTGATTACTACCCTCCAGCCCTGCAAACaaaa         >  1:1809984/1‑69 (MQ=255)                                             ccGCGTTTGATGCGCTGGGCGCGAAAGCGGGTGATTACTACCCTCCAGCCCTGCAAACAAAAATTGACg  >  1:705647/1‑69 (MQ=255)                                             ccGCGTTTGATGCGCTGGGCGCGAAAGCGGGTGATTACTACCCTCCAGCCCTGCAAACAAAAATTGACg  >  1:193676/1‑69 (MQ=255)                                                                |                                                  CATCGTCAGCAACGAATCAGCAGAAATCATCCGCATGTTCAATACCGCGTTTGATGCGCTGGGCGCGAAAGCGGGTGATTACTACCCTCCAGCCCTGCAAACAAAAATTGACG  >  minE/2086654‑2086766

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