Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,619,653	A→G	50.7%	A362A (GCT→GCC)	zntA ←	zinc, cobalt and lead efflux system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,619,653	0	A	G	50.7%	‑1.8 / 65.2	73	A362A (GCT→GCC)	zntA	zinc, cobalt and lead efflux system
Reads supporting (aligned to +/- strand):  ref base A (18/18);  new base G (31/6);  total (49/24)
Fisher's exact test for biased strand distribution p-value = 2.79e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

AATCCACTCCTGCCAGCTGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGATACGGCTGAAACGGTCGATAAACCGCTCAAT  >  minE/2619596‑2619712                                                          |                                                            aaTCCACTCCTGCCAGCTGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCAt                                                  <  1:1210718/69‑1 (MQ=255) aaTCCACTCCTGCCAGCTGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCAt                                                  <  1:533125/69‑1 (MQ=255) aaTCCACTCCTGCCAGCTGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCAt                                                  <  1:2764050/69‑1 (MQ=255)         ccTGCCAGCTGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCgg                                          >  1:948058/1‑69 (MQ=255)         ccTGCCAGCTGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCgg                                          >  1:2915643/1‑69 (MQ=255)         ccTGCCAGCTGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCgg                                          >  1:2088101/1‑69 (MQ=255)         ccTGCCAGCTGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCgg                                          >  1:2815065/1‑69 (MQ=255)         ccTGCCAGCTGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCgg                                          >  1:2228865/1‑69 (MQ=255)         ccTGCCAGCTGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCgg                                          >  1:458685/1‑69 (MQ=255)         ccTGCCAGCTGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCgg                                          >  1:15577/1‑69 (MQ=255)         ccTGCCAGCTGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCgg                                          >  1:154882/1‑69 (MQ=255)         ccTGCCAGCTGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCgg                                          >  1:150843/1‑69 (MQ=255)         ccTGCCAGCTGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCgg                                          >  1:807664/1‑69 (MQ=255)         ccTGCCAGCTGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCgg                                          >  1:837761/1‑69 (MQ=255)         ccTGCCAGCTGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCgg                                          >  1:266928/1‑69 (MQ=255)         ccTGCCAGCTGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCgg                                          >  1:2641199/1‑69 (MQ=255)         ccTGCCAGCTGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCgg                                          >  1:2509760/1‑69 (MQ=255)                  tGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGa                                 >  1:846274/1‑69 (MQ=255)                  tGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGa                                 >  1:19884/1‑69 (MQ=255)                  tGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGa                                 >  1:1174669/1‑69 (MQ=255)                  tGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGa                                 >  1:940490/1‑69 (MQ=255)                  tGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGa                                 <  1:2211539/69‑1 (MQ=255)                  tGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGa                                 <  1:2207471/69‑1 (MQ=255)                                       gCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGATACGGCTGAAACGGTCGATaa            <  1:2815054/69‑1 (MQ=255)                                       gCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGATACGGCTGAAACGGTCGATaa            <  1:3213286/69‑1 (MQ=255)                                       gCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGATACGGCTGAAACGGTCGATaa            <  1:3202302/69‑1 (MQ=255)                                       gCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGATACGGCTGAAACGGTCGATaa            <  1:1607366/69‑1 (MQ=255)                                       gCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGATACGGCTGAAACGGTCGATaa            <  1:592156/69‑1 (MQ=255)                                       gCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGATACGGCTGAAACGGTCGATaa            <  1:135403/69‑1 (MQ=255)                                       gCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGATACGGCTGAAACGGTCGATaa            <  1:1277104/69‑1 (MQ=255)                                        cACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGATACGGCTGAAACGGTCGATaa            <  1:209127/68‑1 (MQ=255)                                           cAGCGTTACCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCg        <  1:3069154/69‑1 (MQ=255)                                           cAGCGTTACCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCg        <  1:1094997/69‑1 (MQ=255)                                           cAGCGTTACCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCg        <  1:913888/69‑1 (MQ=255)                                           cAGCGTTACCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCg        <  1:2481791/69‑1 (MQ=255)                                           cAGCGTTACCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCg        <  1:2477753/69‑1 (MQ=255)                                           cAGCGTTACCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCg        <  1:1403107/69‑1 (MQ=255)                                            aGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGa                                 <  1:153234/43‑1 (MQ=255)                                            aGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGa                                 <  1:2458816/43‑1 (MQ=255)                                            aGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGa                                 <  1:1995341/43‑1 (MQ=255)                                            aGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGa                                 <  1:1396610/43‑1 (MQ=255)                                            aGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGa                                 <  1:3066179/43‑1 (MQ=255)                                                ttaCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTc     >  1:2478615/3‑67 (MQ=255)                                                ttaCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTc     >  1:327510/3‑67 (MQ=255)                                                ttaCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTc     >  1:2181769/3‑67 (MQ=255)                                                ttaCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTc     >  1:1274166/3‑67 (MQ=255)                                                ttaCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTc     >  1:2321372/3‑67 (MQ=255)                                                ttaCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTc     >  1:2345198/3‑67 (MQ=255)                                                ttaCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTc     >  1:2763744/3‑67 (MQ=255)                                                ttaCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTc     >  1:2358850/3‑67 (MQ=255)                                                ttaCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTc     >  1:1068341/3‑67 (MQ=255)                                                 taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCaa   >  1:2386833/2‑68 (MQ=255)                                                 taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:1324010/2‑69 (MQ=255)                                                 taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:1379873/2‑69 (MQ=255)                                                 taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:814914/2‑69 (MQ=255)                                                 taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:1142794/2‑69 (MQ=255)                                                 taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:71371/2‑69 (MQ=255)                                                 taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:958432/2‑69 (MQ=255)                                                 taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:578750/2‑69 (MQ=255)                                                 taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:1747421/2‑69 (MQ=255)                                                 taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:534268/2‑69 (MQ=255)                                                 taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:468730/2‑69 (MQ=255)                                                 taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:1559124/2‑69 (MQ=255)                                                 taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:378010/2‑69 (MQ=255)                                                 taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:3102736/2‑69 (MQ=255)                                                 taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:2985705/2‑69 (MQ=255)                                                 taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:22133/2‑69 (MQ=255)                                                 taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:2904186/2‑69 (MQ=255)                                                 taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:2225678/2‑69 (MQ=255)                                                 taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:2562213/2‑69 (MQ=255)                                                 taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:239925/2‑69 (MQ=255)                                                 taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGCTCAAt  >  1:240914/2‑69 (MQ=255)                                                 taCCAGCAGGGCGACGGCCATAATCGCTGGCGTGTAAATACGGCTGAAACGGTCGATAAACCGATCAAt  >  1:981530/2‑69 (MQ=255)                                                          |                                                            AATCCACTCCTGCCAGCTGGCGGCAAACAGCAGCGGTGGCACCAGCGTCACCAGCAGAGCGACGGCCATAATCGCGGGCGTATAGATACGGCTGAAACGGTCGATAAACCGCTCAAT  >  minE/2619596‑2619712

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