Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,031,911	T→A	42.6%	Q30L (CAG→CTG)	ydgJ ←	predicted oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,031,911	0	T	A	42.6%	18.9 / 54.1	54	Q30L (CAG→CTG)	ydgJ	predicted oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (10/21);  new base A (5/18);  total (15/39)
Fisher's exact test for biased strand distribution p-value = 5.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

CGTTGGCCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCGCATAACCATACCCAATCAA  >  minE/1031853‑1031975                                                           |                                                                 cgtTGGCCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGt                                                        <  1:280095/69‑1 (MQ=255) cgtTGGCCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGt                                                        <  1:3094988/69‑1 (MQ=255)       ccAGTCAGCTTTTACTTTTGTTTCATCGCTGCTTGAGATTACCGCCAGCTCCAGCCCgg                                                            <  1:953173/59‑1 (MQ=37)       ccAGTCAGCTTTTACTTTTGTTTCATCGCTGCTGGAGATTACCGCCAGCTCCAGCCCgg                                                            <  1:3278426/59‑1 (MQ=255)       ccAGTCAGCTTTTACTTTTGTTTCATCGCTGCTGGAGATTACCGCCAGCTCCAGCCCgg                                                            <  1:2974085/59‑1 (MQ=255)       ccAGTCAGCTTTTACTTTTGTTTCATCGCTGCTGGAGATTACCGCCAGCTCCAGCCCgg                                                            <  1:2999643/59‑1 (MQ=255)       ccAGTCAGCTTTTACTTTTGTTTCATCGCTGCTGGAGATTACCGCCAGCTCCAGCCCgg                                                            <  1:2843469/59‑1 (MQ=255)       ccAGTCAGCTTTTACTTTTGTTTCATCGCTGCTGGAGATTACCGCCAGCTCCAGCCCgg                                                            <  1:2606262/59‑1 (MQ=255)       ccAGTCAGCTTTTACTTTTGTTTCATCGCTGCTGGAGATTACCGCCAGCTCCAGCCCgg                                                            <  1:3108993/59‑1 (MQ=255)       ccAGTCAGCTTTTACTTTTGTTTCATCGCTGCTGGAGATTACCGCCAGCTCCAGCCCgg                                                            <  1:313282/59‑1 (MQ=255)       ccAGTCAGCTTTTACTTTTGTTTCATCGCTGCTGGAGATTACCGCCAGCTCCAGCCCgg                                                            <  1:653891/59‑1 (MQ=255)       ccAGTCAGCTTTTACTTTTGTTTCATCGCTGCTGGAGATTACCGCCAGCTCCAGCCCgg                                                            <  1:175938/59‑1 (MQ=255)       ccAGTCAGCTTTTACTTTTGTTTCATCGCTGCTGGAGATTACCGCCAGCTCCAGCCCgg                                                            <  1:1756959/59‑1 (MQ=255)       ccAGTCAGCTTTTACTTTTGTTTCATCGCTGCTGGAGATTACCGCCAGCTCCAGCCCgg                                                            <  1:1728097/59‑1 (MQ=255)       ccAGTCAGCTTTTACTTTTGTTTCATCGCTGCTGGAGATTACCGCCAGCTCCAGCCCgg                                                            <  1:1632427/59‑1 (MQ=255)       ccAGTCAGCTTTTACTTTTGTTTCATCGCTGCTGGAGATTACCGCCAGCTCCAGCCCgg                                                            <  1:1593192/59‑1 (MQ=255)       ccAGTCAGCTTTTACTTTTGTTTCATCGCTGCTGGAGATTACCGCCAGCTCCAGCCCgg                                                            <  1:1535815/59‑1 (MQ=255)       ccAGTCAGCTTTTACTTTTGTTTCATCGCTGCTGGAGATTACCGCCAGCTCCAGCCCgg                                                            <  1:1504242/59‑1 (MQ=255)       ccAGTCAGCTTTTACTTTTGTTTCATCGCTGCTGGAGATTACCGCCAGCTCCAGCCCgg                                                            <  1:3277636/59‑1 (MQ=255)       ccAGTCAGCTTTTACTTTTGTTTCATCGCTGCTGGAGATTACCGCCAGCTCCAGCCCgg                                                            <  1:1341078/59‑1 (MQ=255)                 tttACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAgggg                                         >  1:1764858/1‑68 (MQ=255)                 tttACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAgggg                                         >  1:1931781/1‑68 (MQ=255)                 tttACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAgggg                                         >  1:1535448/1‑68 (MQ=255)                 tttACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAgggg                                         >  1:1492257/1‑68 (MQ=255)                 tttACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAgggg                                         >  1:2751080/1‑68 (MQ=255)                 tttACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAgggg                                         >  1:1410591/1‑68 (MQ=255)                 tttACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAggg                                          >  1:3029916/1‑67 (MQ=255)                    aCTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGgcg                                       <  1:1296498/67‑1 (MQ=255)                                           gaTTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCTCATAAc               <  1:6441/68‑1 (MQ=255)                                           gaTTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCGCATAAc               <  1:371852/68‑1 (MQ=255)                                           gaTTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCGCATAAc               <  1:422345/68‑1 (MQ=255)                                           gaTTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCGCATAAc               <  1:360300/68‑1 (MQ=255)                                           gaTTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCGCATAAc               <  1:454177/68‑1 (MQ=255)                                           gaTTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCGCATAAc               <  1:478455/68‑1 (MQ=255)                                           gaTTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCGCATAAc               <  1:524076/68‑1 (MQ=255)                                           gaTTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCGCATAAc               <  1:640502/68‑1 (MQ=255)                                           gaTTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCGCATAAc               <  1:103791/68‑1 (MQ=255)                                           gaTTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCGCATAAc               <  1:964527/68‑1 (MQ=255)                                           gaTTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCGCATAAc               <  1:2991237/68‑1 (MQ=255)                                           gaTTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCGCATAAc               <  1:2626409/68‑1 (MQ=255)                                           gaTTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCGCATAAc               <  1:254090/68‑1 (MQ=255)                                           gaTTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCGCATAAc               <  1:2429922/68‑1 (MQ=255)                                           gaTTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCGCATAAc               <  1:2404588/68‑1 (MQ=255)                                           gaTTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCGCATAAc               <  1:2392501/68‑1 (MQ=255)                                           gaTTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCGCATAAc               <  1:2185010/68‑1 (MQ=255)                                           gaTTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCGCATAAc               <  1:1345561/68‑1 (MQ=255)                                               aCCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCGCATAACCATa           >  1:2929764/1‑68 (MQ=255)                                                      gctCCAGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCGCATAACCATACCCAATCa   >  1:3325767/3‑69 (MQ=255)                                                      gctCCAGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCGCATAACCATACCCAATCa   >  1:2900274/3‑69 (MQ=255)                                                      gctCCAGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCGCATAACCATACCCAATCa   >  1:2878876/3‑69 (MQ=255)                                                      gctCCAGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCGCATAACCATACCCAATCa   >  1:485922/3‑69 (MQ=255)                                                      gctCCAGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCGCATAACCATACCCAATCa   >  1:678399/3‑69 (MQ=255)                                                         ccTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCGCATAAc               >  1:585247/1‑54 (MQ=255)                                                         ccTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCGCATAACCATACCCAATCaa  >  1:3433361/1‑67 (MQ=255)                                                           |                                                                 CGTTGGCCAGTCGGCTTTTACTTTTGTTTCATCACTGCTGGAGATTACCGCCAGTTCCTGCCCGGGCGTGCCCGCAATCAGGGGCGCATGGAAGGTTTTGCTCGCATAACCATACCCAATCAA  >  minE/1031853‑1031975

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