Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	124,043	A→G	39.1%	I641I (ATT→ATC)	gcd ←	glucose dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	124,043	0	A	G	39.1%	17.0 / 28.5	46	I641I (ATT→ATC)	gcd	glucose dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (20/8);  new base G (10/8);  total (30/16)
Fisher's exact test for biased strand distribution p-value = 3.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.13e-01

ATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAG  >  minE/123988‑124101                                                        |                                                           acggTGAGAGGAACGGGTTGAGCGTGACGCCATACGGTACACCGTACTTCGGCTGGATGCCGGATTCCg                                               <  1:3266283/67‑1 (MQ=255) acggTGAGAGGAACGGGTTGAGCGTGACGCCATACGGTACACCGTACTGCGGCTGGATGCCGGATTCCg                                               <  1:1556627/67‑1 (MQ=255) acggTGAGAGGAACGGGTTGAGCGTGACGCCATACGGTACACCGTACTGCGGCTGGATGCCGGATTCCg                                               <  1:2748375/67‑1 (MQ=255) acggTGAGAGGAACGGGTTGAGCGTGACGCCATACGGTACACCGTACTGCGGCTGGATGCCGGATTCCg                                               <  1:2016570/67‑1 (MQ=255) acggTGAGAGGAACGGGTTGAGCGTGACGCCATACGGTACACCGTACTGCGGCTGGATGCCGGATTCCg                                               <  1:2621210/67‑1 (MQ=255) acggTGAGAGGAACGGCTTGAGCGTGACGCCATACGGTACACCGTACTGCGGCTGGATGCCGGATTCCg                                               <  1:3061365/67‑1 (MQ=255)  cggTGAGAGGAACGGGTTGAGCGTGACGCCATACGGTACACCGTACTGCGGCTGGATGCCGGATTCCg                                               <  1:3429774/67‑1 (MQ=255)   ggTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTa                                             >  1:3210415/1‑69 (MQ=255)   ggTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTa                                             >  1:810315/1‑69 (MQ=255)      gagaGGAACTGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTAccc                                          >  1:1415204/1‑69 (MQ=255)      gagaGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTAccc                                          >  1:3289391/1‑69 (MQ=255)      gagaGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTAccc                                          >  1:1545603/1‑69 (MQ=255)      gagaGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTAccc                                          >  1:1916665/1‑69 (MQ=255)      gagaGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTAccc                                          >  1:219122/1‑69 (MQ=255)           gAACGGGTTGAGCGTGACGCCATACGGTACACCGTACTGCGGCTGGATGCCGGATTCCg                                               <  1:1562898/59‑1 (MQ=255)            aaCGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCt                                    >  1:93626/1‑69 (MQ=255)            aaCGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCt                                    >  1:2237151/1‑69 (MQ=255)            aaCGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCt                                    >  1:1990503/1‑69 (MQ=255)                ggTTGAGCGTGACACCATACGGTACACCGTACTGTGGGTGAATGCCGGATTCCGTACCCGTGCCTTTgg                                >  1:2764125/1‑69 (MQ=255)                ggTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTgg                                >  1:1341986/1‑69 (MQ=255)                           acacCATACGGTACACCGTACTGTGGCTGAATGcc                                                       >  1:490259/1‑35 (MQ=255)                           acacCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTAGCCGTGCCTTTGGCATCTTTCGgc                     >  1:3368928/1‑69 (MQ=255)                           acacCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGt                                        >  1:3217745/1‑50 (MQ=255)                           acacCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGgc                     >  1:2552602/1‑69 (MQ=255)                           acacCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGgc                     >  1:73530/1‑69 (MQ=255)                           acacCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGgc                     >  1:2980884/1‑69 (MQ=255)                           acacCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGgc                     >  1:755174/1‑69 (MQ=255)                           acacCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGgc                     >  1:1280068/1‑69 (MQ=255)                                            gTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGcc    <  1:464655/69‑1 (MQ=255)                                            gTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGcc    <  1:1168369/69‑1 (MQ=255)                                            gTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGcc    <  1:2692565/69‑1 (MQ=255)                                            gTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGcc    <  1:2551200/69‑1 (MQ=255)                                            gTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGcc    <  1:1872213/69‑1 (MQ=255)                                            gTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGcc    <  1:1054851/69‑1 (MQ=255)                                               cTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGcc    <  1:536388/66‑1 (MQ=255)                                               cTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAg  <  1:2327679/68‑1 (MQ=255)                                                    gCTGGATGCCGGATTCCGTACCCGTGGCTTTGGCATCt                           >  1:789524/1‑38 (MQ=38)                                                    gCTGGATGCCGGATTCCGTACCCGTGCCTTTGGCATCt                           >  1:299982/1‑38 (MQ=255)                                                    gCTGGATGCCGGATTCCGTACCCGTGCCTTTGGCATCt                           >  1:2865160/1‑38 (MQ=255)                                                    gCTGGATGCCGGATTCCGTACCCGTGCCTTTGGCATCt                           >  1:2853568/1‑38 (MQ=255)                                                    gCTGGATGCCGGATTCCGTACCCGTGCCTTTGGCATCt                           >  1:488611/1‑38 (MQ=255)                                                    gCTGGATGCCGGATTCCGTACCCGTGCCTTTGGCATCt                           >  1:647097/1‑38 (MQ=255)                                                    gCTGGATGCCGGATTCCGTACCCGTGCCTTTGGCATCt                           >  1:1305758/1‑38 (MQ=255)                                                    gCTGGATGCCGGATTCCGTACCCGTGCCTTTGGCATCt                           >  1:781612/1‑38 (MQ=255)                                                    gCTGGATGCCGGATTCCGTACCCGTGCCTTTGGCATCt                           >  1:118960/1‑38 (MQ=255)                                                    gCTGGATGCCGGATTCCGTACCCGTGCCTTTGGCATCt                           >  1:811292/1‑38 (MQ=255)                                                        |                                                           ATGGTGAGAGGAACGGGTTGAGCGTGACACCATACGGTACACCGTACTGTGGCTGAATGCCGGATTCCGTACCCGTGCCTTTGGCATCTTTCGGCTGCTCCATCGGGTTGCCAG  >  minE/123988‑124101

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