Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,422,091	C→A	59.4%	V122V (GTG→GTT)	ccmF ←	heme lyase, CcmF subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,422,091	0	C	A	59.4%	21.0 / 61.5	64	V122V (GTG→GTT)	ccmF	heme lyase, CcmF subunit
Reads supporting (aligned to +/- strand):  ref base C (20/6);  new base A (22/16);  total (42/22)
Fisher's exact test for biased strand distribution p-value = 1.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

AGAGGTAAAGAGAATGAACAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGC  >  minE/1422025‑1422157                                                                   |                                                                   agagGTAAAGAGAATGAACAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCcac                                                                  <  1:2625728/69‑1 (MQ=255) agagGTAAAGAGAATGAACAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCcac                                                                  <  1:2261351/69‑1 (MQ=255)       aaaGAGAATGAACAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAata                                                              >  1:274328/1‑67 (MQ=255)       aaaGAGAATGAACAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAata                                                              >  1:2724758/1‑67 (MQ=255)       aaaGAGAATGAACAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAata                                                              >  1:2126540/1‑67 (MQ=255)        aaGAGAATGAACAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATcc                                                           <  1:991411/69‑1 (MQ=255)        aaGAGAATGAACAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATcc                                                           <  1:1439656/69‑1 (MQ=255)               tGAACAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGaa                                                    <  1:1097097/69‑1 (MQ=255)               tGAACAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGaa                                                    <  1:967848/69‑1 (MQ=255)                 aaCAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAAt                                                   >  1:2094663/1‑68 (MQ=255)                 aaCAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATa                                                  >  1:2456322/1‑69 (MQ=255)                 aaCAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATa                                                  >  1:2258919/1‑69 (MQ=255)                 aaCAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATa                                                  >  1:2498122/1‑69 (MQ=255)                 aaCAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATa                                                  >  1:3168609/1‑69 (MQ=255)                 aaCAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATa                                                  >  1:3115248/1‑69 (MQ=255)                 aaCAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATa                                                  >  1:2910830/1‑69 (MQ=255)                 aaCAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATa                                                  >  1:2004696/1‑69 (MQ=255)                 aaCAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATa                                                  >  1:2793764/1‑69 (MQ=255)                     gCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCt                                              <  1:2845568/69‑1 (MQ=255)                     gCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCt                                              <  1:274988/69‑1 (MQ=255)                     gCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCt                                              <  1:2763876/69‑1 (MQ=255)                     gCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCt                                              <  1:1440724/69‑1 (MQ=255)                     gCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCt                                              <  1:3244380/69‑1 (MQ=255)                     gCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCt                                              <  1:330218/69‑1 (MQ=255)                     gCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCt                                              <  1:3394754/69‑1 (MQ=255)                     gCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCt                                              <  1:615369/69‑1 (MQ=255)                     gCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCt                                              <  1:2280282/69‑1 (MQ=255)                     gCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCt                                              <  1:1328176/69‑1 (MQ=255)                     gCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCt                                              <  1:2016442/69‑1 (MQ=255)                     gCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGAACGCGGGCAACAATATCCAGCGGAATACGCt                                              <  1:1419468/69‑1 (MQ=255)                      cAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCt                                              <  1:197918/68‑1 (MQ=255)                       aaaaaGCCGACACTGACCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCt                                              <  1:849496/67‑1 (MQ=255)                               gACACTGACCATCCCCATTATCGCCAGTACACGGTCCACAATATCCAGCGGaa                                                    <  1:36728/53‑1 (MQ=255)                                acacTGACCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCt                                              <  1:2485649/58‑1 (MQ=255)                                    tGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAGTCGCCa                               <  1:440148/69‑1 (MQ=255)                                      aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc                              >  1:1792434/1‑68 (MQ=255)                                      aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc                              >  1:2338288/1‑68 (MQ=255)                                      aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc                              >  1:1745356/1‑68 (MQ=255)                                      aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc                              >  1:1717680/1‑68 (MQ=255)                                      aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc                              >  1:1500279/1‑68 (MQ=255)                                      aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc                              >  1:377029/1‑68 (MQ=255)                                      aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc                              >  1:1434402/1‑68 (MQ=255)                                      aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc                              >  1:617359/1‑68 (MQ=255)                                      aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc                              >  1:691688/1‑68 (MQ=255)                                      aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc                              >  1:807859/1‑68 (MQ=255)                                      aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc                              >  1:903492/1‑68 (MQ=255)                                      aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc                              >  1:1165456/1‑68 (MQ=255)                                      aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc                              >  1:1891646/1‑68 (MQ=255)                                      aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc                              >  1:29065/1‑68 (MQ=255)                                      aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc                              >  1:2898629/1‑68 (MQ=255)                                      aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc                              >  1:2863371/1‑68 (MQ=255)                                      aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc                              >  1:2240015/1‑68 (MQ=255)                                      aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc                              >  1:1288206/1‑68 (MQ=255)                                      aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc                              >  1:1240805/1‑68 (MQ=255)                                      aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc                              >  1:2760495/1‑68 (MQ=255)                                      aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAATCGCTAc                              >  1:1186313/1‑68 (MQ=255)                                      aCCATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCTGACTAAAAACCGCTAc                              >  1:36191/1‑68 (MQ=255)                                       ccATCCCCATTATCGCCAGTACGCGGGCAACAATATCCAGCGGAATACGCt                                              <  1:1897217/51‑1 (MQ=255)                                            cccATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAg                       >  1:2183644/1‑69 (MQ=255)                                                  aTCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAg                 >  1:2777035/1‑69 (MQ=255)                                                  aTCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAg                 >  1:2787557/1‑69 (MQ=255)                                                                 gCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGTTCCAGCCGCTCATCagcagc  >  1:1007775/1‑69 (MQ=255)                                                                 gCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCagcagc  >  1:3137915/1‑69 (MQ=255)                                                                 gCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCagcagc  >  1:2719014/1‑69 (MQ=255)                                                                 gCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCagcagc  >  1:392166/1‑69 (MQ=255)                                                                 gCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCagcagc  >  1:991978/1‑69 (MQ=255)                                                                   |                                                                   AGAGGTAAAGAGAATGAACAGCAAAAAGCCGACACTGACCATCCCCATTATCGCCAGTACACGGGCCACAATATCCAGCGGAATACGCTGACTAAAAATCGCCACCGCAAAGGTCCAGCCGCTCATCAGCAGC  >  minE/1422025‑1422157

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