Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,419,063	C→T	29.4%	A300A (GCG→GCA)	ccmH ←	heme lyase, CcmH subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,419,063	0	C	T	29.4%	78.4 / 33.2	78	A300A (GCG→GCA)	ccmH	heme lyase, CcmH subunit
Reads supporting (aligned to +/- strand):  ref base C (18/37);  new base T (11/12);  total (29/49)
Fisher's exact test for biased strand distribution p-value = 3.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

TCAACATCATCTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTC  >  minE/1419002‑1419125                                                              |                                                               tCAACATCATCTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATAc                                                          >  1:640825/1‑68 (MQ=255) tCAACATCATCTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACa                                                         >  1:2110155/1‑69 (MQ=255) tCAACATCATCTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACa                                                         >  1:932527/1‑69 (MQ=255) tCAACATCATCTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACa                                                         >  1:2448775/1‑69 (MQ=255) tCAACATCATCTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACa                                                         >  1:1844217/1‑69 (MQ=255) tCAACATCATCTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACa                                                         >  1:1456803/1‑69 (MQ=255) tCAACATCATCTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACa                                                         <  1:322843/69‑1 (MQ=255) tCAACATCATCTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACa                                                         >  1:737561/1‑69 (MQ=255)         atcTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAg                                                  >  1:1076911/1‑68 (MQ=255)         atcTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAg                                                  >  1:1007424/1‑68 (MQ=255)         atcTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAg                                                  >  1:2691490/1‑68 (MQ=255)         atcTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAg                                                  >  1:1627222/1‑68 (MQ=255)         atcTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAg                                                  >  1:1581006/1‑68 (MQ=255)         atcTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAg                                                  >  1:1110800/1‑68 (MQ=255)         atcTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAg                                                  >  1:894459/1‑68 (MQ=255)         atcTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAg                                                  >  1:763046/1‑68 (MQ=255)         atcTCCAACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAg                                                  >  1:437335/1‑68 (MQ=255)              ccACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACg                                            >  1:89383/1‑69 (MQ=255)              acACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACg                                            >  1:1861820/2‑69 (MQ=255)                aCGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCTCACGGa                                          <  1:2825431/69‑1 (MQ=255)                aCGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGa                                          <  1:2925769/69‑1 (MQ=255)                aCGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGa                                          <  1:829395/69‑1 (MQ=255)                aCGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGa                                          <  1:293472/69‑1 (MQ=255)                aCGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGa                                          <  1:2796787/69‑1 (MQ=255)                aCGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGa                                          <  1:3079513/69‑1 (MQ=255)                aCGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGa                                          <  1:2408556/69‑1 (MQ=255)                aCGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACAAGGa                                          <  1:3421798/69‑1 (MQ=255)                   cgGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGatat                                       <  1:1480012/69‑1 (MQ=255)                   cgGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGatat                                       <  1:2292541/69‑1 (MQ=255)                   cgGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGatat                                       <  1:2128819/69‑1 (MQ=255)                         cGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTa                                  <  1:2789326/68‑1 (MQ=255)                          ggCTTCGCCAAATCGCTGCTGCTCAAAGGCTTTAAACGCATa                                                           <  1:231903/42‑1 (MQ=255)                          ggCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATa                                                           <  1:50850/42‑1 (MQ=255)                          ggCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATa                                                           <  1:1318080/42‑1 (MQ=255)                          ggCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATa                                                           <  1:1853717/42‑1 (MQ=255)                          ggCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATa                                                           <  1:1883487/42‑1 (MQ=255)                          ggCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATa                                                           <  1:18945/42‑1 (MQ=255)                          ggCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATa                                                           <  1:1901474/42‑1 (MQ=255)                          ggCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATa                                                           <  1:1902319/42‑1 (MQ=255)                          ggCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATa                                                           <  1:1762768/42‑1 (MQ=255)                          ggCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATa                                                           <  1:254554/42‑1 (MQ=255)                          ggCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATa                                                           <  1:1906176/42‑1 (MQ=255)                          ggCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATa                                                           <  1:1258177/42‑1 (MQ=255)                          ggCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATa                                                           <  1:1990608/42‑1 (MQ=255)                          ggCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATa                                                           <  1:2312785/42‑1 (MQ=255)                          ggCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATa                                                           <  1:1593954/42‑1 (MQ=255)                          ggCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATa                                                           <  1:3284454/42‑1 (MQ=255)                          ggCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATa                                                           <  1:331606/42‑1 (MQ=255)                              tCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATa                                                           <  1:609534/38‑1 (MQ=255)                                 ccAAATCGTTGCTGCTCAAAGGCATTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGtt                         >  1:3082894/1‑69 (MQ=255)                                 ccAAATCGTTGCTGCTCAAAGGCATTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGtt                         >  1:3058070/1‑69 (MQ=255)                                 ccAAATCGTTGCTGCTCAAAGGCATTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGtt                         >  1:2992264/1‑69 (MQ=255)                                 ccAAATCGTTGCTGCTCAAAGGCATTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGtt                         >  1:214881/1‑69 (MQ=255)                                 ccAAATCGTTGCTGCTCAAAGGCATTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGtt                         >  1:2906538/1‑69 (MQ=255)                                 ccAAATCGTTGCTGCTCAAAGGCATTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGtt                         >  1:1314563/1‑69 (MQ=255)                                 ccAAATCGTTGCTGCTCAAAGGCATTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGtt                         >  1:2483798/1‑69 (MQ=255)                                 ccAAATCGTTGCTGCTCAAAGGCATTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGtt                         >  1:1181174/1‑69 (MQ=255)                                 ccAAATCGTTGCTGCTCAAAGGCATTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGtt                         >  1:2247670/1‑69 (MQ=255)                                 ccAAATCGTTGCTGCTCAAAGGCATTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGtt                         >  1:2192357/1‑69 (MQ=255)                                 acAAATCGTTGCTGCTCAAAGGCATTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGt                          >  1:3058074/2‑68 (MQ=255)                                              gctCAAAGGCATTAAACGCATACATGCTTAGCACACGGa                                          <  1:810974/39‑1 (MQ=255)                                                    aGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGc       <  1:2565719/68‑1 (MQ=255)                                                    aGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGc       <  1:2225789/68‑1 (MQ=255)                                                        aTTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGTTCTCACCAGCTGGCGCAGCAGTTc  <  1:1596765/69‑1 (MQ=255)                                                        aTTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGTTCTCACCAACTGGCGCAGCAGTTc  <  1:1649500/69‑1 (MQ=255)                                                        aTTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGTTCTCACCAACTGGCGCAGCAGTTc  <  1:979732/69‑1 (MQ=255)                                                        aTTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGTTCTCACCAACTGGCGCAGCAGTTc  <  1:523255/69‑1 (MQ=255)                                                        aTTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGTTCTCACCAACTGGCGCAGCAGTTc  <  1:512668/69‑1 (MQ=255)                                                        aTTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGTTCTCACCAACTGGCGCAGCAGTTc  <  1:1663608/69‑1 (MQ=255)                                                        aTTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGTTCTCACCAACTGGCGCAGCAGTTc  <  1:174010/69‑1 (MQ=255)                                                        aTTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGTTCTCACCAACTGGCGCAGCAGTTc  <  1:3035326/69‑1 (MQ=255)                                                        aTTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGTTCTCACCAACTGGCGCAGCAGTTc  <  1:1809885/69‑1 (MQ=255)                                                        aTTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGTTCTCACCAACTGGCGCAGCAGTTc  <  1:1992957/69‑1 (MQ=255)                                                        aTTAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGTTCTCACCAACTGGCGCAGCAGTTc  <  1:2186001/69‑1 (MQ=255)                                                        aTTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTc  <  1:3386092/69‑1 (MQ=255)                                                        aTTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTc  <  1:1748453/69‑1 (MQ=255)                                                        aTTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTc  <  1:2647015/69‑1 (MQ=255)                                                         ttAAATGCATACATGCTTAACACACGGATATTGCTATGGTCCGTTCTCACCAACTGGCGCAGCAGTTc  <  1:1947206/68‑1 (MQ=255)                                                              |                                                               TCAACATCATCTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTC  >  minE/1419002‑1419125

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