Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	232,066	T→C	76.4%	E232E (GAA→GAG)	hemB ←	porphobilinogen synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	232,066	0	T	C	76.4%	62.7 / 13.4	51	E232E (GAA→GAG)	hemB	porphobilinogen synthase
Reads supporting (aligned to +/- strand):  ref base T (10/2);  new base C (19/19);  total (29/22)
Fisher's exact test for biased strand distribution p-value = 5.12e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

ACGCTCCAGCAGGTTTAACCATCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGTGATTCACGAATCGCCTCACGACGGTTCATTGGGTTCATCTGATAGCTTTTGCGGTCGCCTTTTAATG  >  minE/232008‑232130                                                           |                                                                 aCGCTCCAGCAGGTTTAACCATCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGTGATTCACGAATCg                                                        >  1:1142996/1‑69 (MQ=255)  cGCTCCAGCAGGTTTAACCATCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGTGATTCACGAATCGc                                                       >  1:1351547/1‑69 (MQ=255)  cGCTCCAGCAGGTTTAACCATCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGTGATTCACGAATCGc                                                       >  1:1717182/1‑69 (MQ=255)           agGTTTAACCATCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGTGATTCACGAATCGCCTCacgacg                                              >  1:3323372/1‑69 (MQ=255)           agGTTTAACCATCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGTGATTCACGAATCGCCTCacgacg                                              >  1:3355300/1‑69 (MQ=255)           agGTTTAACCATCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGTGATTCACGAATCGCCTCacgacg                                              >  1:652755/1‑69 (MQ=255)            ggTTTAACCATCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGCGACTCACGAATCGCCTCACGACgg                                             <  1:1017871/69‑1 (MQ=255)            ggTTTAACCATCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGCGACTCACGAATCGCCTCACGACgg                                             <  1:2627934/69‑1 (MQ=255)            ggTTTAACCATCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGCGACTCACGAATCGCCTCACGACgg                                             <  1:3074981/69‑1 (MQ=255)            ggTTTAACCATCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGCGACTCACGAATCGCCTCACGACgg                                             <  1:2602824/69‑1 (MQ=255)            ggTTTAACCATCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGCGACTCACGAATCGCCTCACGACgg                                             <  1:3329758/69‑1 (MQ=255)            ggTTTAACCATCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGCGACTCACGAATCGCCTCACGACgg                                             <  1:2387447/69‑1 (MQ=255)            ggTTTAACCATCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGCGACTCACGAATCGCCTCACGACgg                                             <  1:2380686/69‑1 (MQ=255)            ggTTTAACCATCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGCGACTCACGAATCGCCTCACGACgg                                             <  1:2289698/69‑1 (MQ=255)            ggTTTAACCATCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGCGACTCACGAATCGCCTCACGACgg                                             <  1:527754/69‑1 (MQ=255)            ggTTTAACCATCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGCGACTCACGAATCGCCTCACGACgg                                             <  1:1744341/69‑1 (MQ=255)            ggTTTAACCATCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGCGACTCACGAATCGCCTCACGACgg                                             <  1:762741/69‑1 (MQ=255)            ggTTTAACCATCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGCGACTCACGAATCGCCTCACGACgg                                             <  1:1353227/69‑1 (MQ=255)            ggTTTAACCATCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGCGACTCACGAATCGCCTCACGACgg                                             <  1:1222283/69‑1 (MQ=255)               ttAACCATCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGTGATTCACGAATCGCCTCACGACGGTTc                                          >  1:1922965/1‑69 (MQ=255)               ttAACCATCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGTGATTCACGAATCGCCTCACGACGGTTc                                          >  1:3134851/1‑69 (MQ=255)               ttAACCATCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGTGATTCACGAATCGCCTCACGACGGTTc                                          <  1:2970454/69‑1 (MQ=255)                     aTCAGGGAGTCTGCGCCCTGGGCTTCATCCAGCAGCGACTCACGAATCGCCTCACGACGGtt                                           >  1:440670/1‑62 (MQ=255)                     aTCAGGGAGTCTGCGCCCTGGGCTTAATCCAGCAGCGACTCACGAATCGCCTCACGACGGTTCATTggg                                    >  1:1471868/1‑69 (MQ=255)                     aTCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGCGACTCACGAATCGCCTCacga                                                >  1:1590496/1‑57 (MQ=255)                     aTCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGCGACTCACGAATCGCCTCACGACGGTTCATTggg                                    >  1:645650/1‑69 (MQ=255)                     aTCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGCGACTCACGAATCGCCTCACGACGGTTCATTggg                                    >  1:51140/1‑69 (MQ=255)                     aTCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGCGACTCACGAATCGCCTCACGACGGTTCATTggg                                    >  1:392395/1‑69 (MQ=255)                     aTCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGCGACTCACGAATCGCCTCACGACGGTTCATTggg                                    >  1:743601/1‑69 (MQ=255)                     aTCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGCGACTCACGAATCGCCTCACGACGGTTCATTggg                                    >  1:867235/1‑69 (MQ=255)                     aTCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGCGACTCACGAATCGCCTCACGACGGTTCATTggg                                    >  1:303427/1‑69 (MQ=255)                     aTCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGCGACTCACGAATCGCCTCACGACGGTTCATTggg                                    >  1:2729816/1‑69 (MQ=255)                     aTCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGCGACTCACGAATCGCCTCACGACGGTTCATTggg                                    >  1:2477301/1‑69 (MQ=255)                     aTCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGCGACTCACGAATCGCCTCACGACGGTTCATTgg                                     >  1:2224798/1‑68 (MQ=255)                     aTCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGCGACTCACGAATCGCCTCACGACGGTTCATTgg                                     >  1:779404/1‑68 (MQ=255)                     aTCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGCGACTCACGAATCGCCTCACGACGGTTCATTgg                                     >  1:2079654/1‑68 (MQ=255)                     aTCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGCGACTCACGAATCGCCTCACGACGGTTCATTgg                                     >  1:2171868/1‑68 (MQ=255)                     aTCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGCGACTCACGAATCGCCTCACGACGGTTCATTgg                                     >  1:2355526/1‑68 (MQ=255)                     aTCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGCGACTCACGAATCGCCTCACGACGGTTCATTgg                                     >  1:2566090/1‑68 (MQ=255)                     aTCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGCGACTCACGAATAGCCTCACGACGGTTCATTgg                                     >  1:2355320/1‑68 (MQ=255)                     aTCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGCGACTAACGAATCGCCTCACGACGGTTCATTggg                                    >  1:1056050/1‑69 (MQ=255)                                tGCGCCCTGGGCTTCATCCAGCAGCGACTCACGAATCGCCTCACGACgg                                             <  1:2619624/49‑1 (MQ=38)                                         ggCTTCATCCAGCAGCGACTCACGAATCGCCTCACGACGGTTc                                          <  1:476373/43‑1 (MQ=38)                                         ggCTTCATCCAGCAGCGACTCACGAATCGCCTCACGACGGTTc                                          <  1:2739946/43‑1 (MQ=38)                                         ggCTTCATCCAGCAGCGACTCACGAATCGCCTCACGACGGTTc                                          <  1:2601864/43‑1 (MQ=38)                                         ggCTTCATCCAGCAGCGACTCACGAATCGCCTCACGACGGTTc                                          <  1:598964/43‑1 (MQ=38)                                          gCTTCATCCAGCAGCGACTCACGAATCGCCTCACGACGGTTc                                          <  1:712442/42‑1 (MQ=38)                                                tCCAGCAGTGATTCACGAATCGCCTCACGACGGTTCATTGGGTTCATCTGATAGCTTTTGCGGTCGc           <  1:1630571/67‑1 (MQ=255)                                                tCCAGCAGTGAATCACGAATCGCCTCACGACGGTTCATTGGGTTCATCTGATAGCTTTTGCGGTCGc           <  1:1789526/67‑1 (MQ=255)                                                       gTGATTCACGAATCGCCTCACGACGGTTCATTGGGTTCATCTGATAGCTTTTGCGGTCGCCTTTTa     >  1:1758402/1‑66 (MQ=255)                                                       gTGATTCACGAATCGCCTCACGACGGTTCATTGGGTTCATCTGATAGCTTTTGCGGTCGCCTTTTAATg  >  1:2626791/1‑69 (MQ=255)                                                           |                                                                 ACGCTCCAGCAGGTTTAACCATCAGGCAGTCTGCGCCCTGGGCTTCATCCAGCAGTGATTCACGAATCGCCTCACGACGGTTCATTGGGTTCATCTGATAGCTTTTGCGGTCGCCTTTTAATG  >  minE/232008‑232130

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