Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,422,361	A→G	21.3%	D32D (GAT→GAC)	ccmF ←	heme lyase, CcmF subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,422,361	0	A	G	21.3%	33.4 / 4.2	28	D32D (GAT→GAC)	ccmF	heme lyase, CcmF subunit
Reads supporting (aligned to +/- strand):  ref base A (12/10);  new base G (2/4);  total (14/14)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.39e-01

GCCACAGACATAAACAGCAGCCAGGCAAACAAGCGGGAAGACGCCATCATGCGCGCATCTCCGCGCGCCACGCCCCATAGCGGATACACGGACAGCAGCAGCGCAATTCCCAGCGCC  >  minE/1422305‑1422421                                                         |                                                             gccaCAGACATAAACAGCAGCCAGGCAAACAAGCGGGAAGACGCCATCATGCGCGCATCTCCgcgcgc                                                   <  1:206035/68‑1 (MQ=255)       gACATAAACAGCAGCCAGGCAAACAAGCGGGAAGACGCCATCATGCGCGCATCTCCGCGCGCCACGccc                                            <  1:1421328/69‑1 (MQ=255)             aaCAGCAGCCAGGCAAACAAGCGGGAAGACGCCATCATGCGCGCATCTCcgcg                                                      >  1:749759/1‑53 (MQ=255)             aaCAGCAGCCAGGCAAACAAGCGGGAAGACGCCATCATGCGCGCATCTCCGCGCGCCACGCCCCATAg                                       >  1:126875/1‑68 (MQ=255)             aaCAGCAGCCAGGCAAACAAGCGGGAAGACGCCATCATGCGCGCATCTCCGCGCGCCACGCCCCATAg                                       >  1:321199/1‑68 (MQ=255)             aaCAGCAGCCAGGCAAACAAGCGGGAAGACGCCATCATGCGCGCATCTCCGCGCGCCACGCCCCATAg                                       >  1:3115182/1‑68 (MQ=255)                  cagcCAGGCAAACAAGCGGGAAGACGCCATCATGCGCGCATCTCCGCGCGCCACGCCCCATAGCGGATa                                 <  1:1550100/69‑1 (MQ=255)                  cagcCAGGCAAACAAGCGGGAAGACGCCATCATGCGCGCATCTCCGCGCGCCACGCCCCATAGCGGATa                                 <  1:2345791/69‑1 (MQ=255)                      cAGGCAAACAAGCGGGAAGACGCCATCATGCGCGCATCTCCGCGCGCCACGCCCCATAGCGGATACACg                             <  1:2225404/69‑1 (MQ=255)                        ggCAAACAAGCGGGAAGACGCCATCATGCGCGCATCTCCGCGCGCCACGCCCCATAGCGGATACACGGa                           <  1:421598/69‑1 (MQ=255)                         gCAAACAAGCGGGAAGACGCCATCATGCGCGCGTCCCCGCGCGCCACGCCCCATAGCg                                     <  1:1861709/58‑1 (MQ=39)                         gCAAACAAGCGGGAAGACGCCATCATGCGCGCGTCCCCGCGCGCCACGCCCCATAGCg                                     <  1:2027359/58‑1 (MQ=39)                         gCAAACAAGCGGGAAGACGCCATCATGCGCGCGTCCCCGCGCGCCACGCCCCATAGCg                                     <  1:1820002/58‑1 (MQ=39)                         gCAAACAAGCGCGAAGACGCCATCATGCGCGCGTCCCCGCGCGCCACGCCCCATAGCg                                     <  1:165512/58‑1 (MQ=38)                               aaGCGGGAAGACGCCATCATGCGCGCATCTCCGCGCGCCACGCCCCATAGCGGATACACGGAcagcagc                    <  1:2155982/69‑1 (MQ=255)                               aaGCGGGAAGACGCCATCATGCGCGCATCTCCGCGCGCCACGCCCCATAGCGGATACACGGAcagcagc                    <  1:1420157/69‑1 (MQ=255)                                      aaGACGCCATCATGCGCGCGTCCCCGCGCGCCACGCCCCATAGCGGATACACGGACAGCAGCAGCGCaa             >  1:3389474/1‑69 (MQ=255)                                      aaGACGCCATCATGCGCGCGTCCCCGCGCGCCACGCCCCATAGCGGATACACGGACAGCAGCAGCGCa              >  1:1260721/1‑68 (MQ=255)                                       aGACGCCATCATGCGCGCATCTCCGCGCGCCACGCCCCATAGCGGATACACGGAcagcagc                    <  1:1626331/61‑1 (MQ=255)                                              atcatGCGCGCATCTCCGCGCGCCCCGCCCCATAGCGGATACACGGACCGCAGCAGCGCAATTCCCAg      >  1:306136/1‑68 (MQ=255)                                              atcatGCGCGCATCTCCGCGCGCCACGCCCCATAGCGGATACACGGACAGCAGCAGCGCAATTCCCAgc     >  1:1459740/1‑69 (MQ=255)                                              atcatGCGCGCATCTCCGCGCGCCACGCCCCATAGCGGATACACGGACAGCAGCAGCGCAATTCCCAgc     >  1:2311218/1‑69 (MQ=255)                                              atcatGCGCGCATCTCCGCGCGCCACGCCCCATAGCGGATACACGGACAGCAGCAGCGCAATTCCCAgc     >  1:695387/1‑69 (MQ=255)                                               tcatGCGCGCATCTCCGCGCGCCACGCCCCATAGCGGATACACGGACAGCAGCAGCGCAATTCCCAgc     <  1:858475/68‑1 (MQ=255)                                                 atGCGCGCATCTCCGCGCGCCACGCCCCATAGCGGATACACGGACAGCAGCAGCGCAATTCCCAGCGcc  >  1:2568943/1‑69 (MQ=255)                                                 atGCGCGCATCTCCGCGCGCCACGCCCCATAGCGGATACACGGACAGCAGCAGCGCAATTCCCAGCGcc  >  1:657135/1‑69 (MQ=255)                                                 atGCGCGCATCTCCGCGCGCCACGCCCCATAGCGGATACACGGACAGCAGCAGCGCAATTCCCAGCGcc  >  1:705946/1‑69 (MQ=255)                                                 atGCGCGCATCTCCGCGCGCCACGCCCCATAGCGGATACACGGACAGCAGCAGCGCAATTCCCAGCGcc  >  1:1327333/1‑69 (MQ=255)                                                         |                                                             GCCACAGACATAAACAGCAGCCAGGCAAACAAGCGGGAAGACGCCATCATGCGCGCATCTCCGCGCGCCACGCCCCATAGCGGATACACGGACAGCAGCAGCGCAATTCCCAGCGCC  >  minE/1422305‑1422421

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