Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,419,579	G→A	15.6%	R128R (CGC→CGT)	ccmH ←	heme lyase, CcmH subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,419,579	0	G	A	15.6%	77.1 / 7.3	32	R128R (CGC→CGT)	ccmH	heme lyase, CcmH subunit
Reads supporting (aligned to +/- strand):  ref base G (16/11);  new base A (2/3);  total (18/14)
Fisher's exact test for biased strand distribution p-value = 6.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

AGATAAACAACATATCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTGGCAG  >  minE/1419511‑1419644                                                                     |                                                                  aGATAAACAACATATCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCAcgcg                                                                 >  1:2331838/1‑71 (MQ=255)    tAAACAACATATCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTa                                                               <  1:655706/70‑1 (MQ=255)     aaaCAACATATCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTAcc                                                             >  1:3880029/1‑71 (MQ=255)     aaaCAACATATCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTAcc                                                             >  1:3348701/1‑71 (MQ=255)      aacaacATATCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTAccc                                                            <  1:3955378/71‑1 (MQ=255)      aacaacATATCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTAccc                                                            <  1:3382915/71‑1 (MQ=255)        caacaTATCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACACGTAccc                                                            >  1:3924345/1‑69 (MQ=255)                     aCGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACACGTACCCGACGCCGCGAGCgg                                              >  1:2744386/1‑70 (MQ=255)                           aCCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGta                                         <  1:996099/69‑1 (MQ=255)                           aCCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGta                                         <  1:4056113/69‑1 (MQ=255)                           aCCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACACGTACCCGACGCCGCGAGCGGGCGta                                         <  1:1838826/69‑1 (MQ=255)                            ccTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTatat                                      <  1:1943858/71‑1 (MQ=255)                            ccTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTatat                                      >  1:3210294/1‑71 (MQ=255)                              ttCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCg                                               >  1:3499037/1‑60 (MQ=255)                                ccGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGAcc                                  >  1:3723260/1‑71 (MQ=255)                                ccGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACACGTACCCGACGCCGCGAGCGGGCGTATATGAc                                   <  1:1588011/70‑1 (MQ=255)                                     aCGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGta                                         >  1:868262/1‑59 (MQ=255)                                     aCGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCa                                <  1:3269573/68‑1 (MQ=255)                                                aGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAAc                              >  1:3916401/1‑59 (MQ=255)                                                 ggAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCgccaatgccaa                  >  1:3000844/1‑70 (MQ=255)                                                  gAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGta                                         <  1:1765176/46‑1 (MQ=255)                                                  gAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCgccaatgccaa                  >  1:670103/1‑69 (MQ=255)                                                  gAAACGCTTCCGGCACCACGCGTACCCGACGCCGAGAGCGGGCGTATATGACCCAACCgccaatgccaa                  >  1:1251247/1‑69 (MQ=255)                                                    aaCGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCgccaatgccaa                  <  1:1632276/67‑1 (MQ=255)                                                    aaCGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGc              <  1:1041361/71‑1 (MQ=255)                                                    aaCGCTTCCGGCACCACACGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGc              <  1:1027886/71‑1 (MQ=255)                                                     aCGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGa                                    >  1:661411/1‑48 (MQ=255)                                                     aCGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGcc             <  1:4133363/71‑1 (MQ=255)                                                       gCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCac           >  1:3965312/1‑71 (MQ=255)                                                       gCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCac           >  1:845141/1‑71 (MQ=255)                                                             ggCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTgg     >  1:3925887/1‑71 (MQ=255)                                                                 ccacGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTGgcag  >  1:2020014/1‑70 (MQ=255)                                                                   acGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCacta         >  1:2255851/1‑61 (MQ=255)                                                                     gcgTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTgg     <  1:744991/63‑1 (MQ=255)                                                                     |                                                                  AGATAAACAACATATCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTGGCAG  >  minE/1419511‑1419644

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