Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	811,467	C→T	53.3%	E517E (GAG→GAA)	ychM ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	811,467	0	C	T	53.3%	4.7 / 64.0	75	E517E (GAG→GAA)	ychM	predicted transporter
Reads supporting (aligned to +/- strand):  ref base C (16/19);  new base T (20/20);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 8.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.37e-01

ATGCCAGCGCGAGCCATAGTGCGCAGTGGCTGGAATTCCACGTTGCACACGCGCAGTTCACATCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACGGCA  >  minE/811402‑811533                                                                  |                                                                   aTGCCTGCGCGAGCCATAGTGCGCAGTGGCTGGAATTCTAGGTTGCACACGCGCAGTTCACATCCTTc                                                                  <  1:130949/68‑1 (MQ=255) aTGCCTGCGCGAGCCATAGTGCGCAGTGGCTGGAATTCCAGGTTGCACACGCGCAGTTCACATCCTTc                                                                  <  1:273654/68‑1 (MQ=255) aTGCCTGCGCGAGCCATAGTGCGCAGTGGCTGGAATTCCAGGTTGCACACGCGCAGTTCACATCCTTc                                                                  <  1:2778946/68‑1 (MQ=255) aTGCCTGCGCGAGCCATAGTGCGCAGTGGCTGGAATTCCAGGTTGCACACGCGCAGTTCACATCCTTc                                                                  <  1:2627626/68‑1 (MQ=255) aTGCCTGCGCGAGCCATAGTGCGCAGTGGCTGGAATTCCAGGTTGCACACGCGCAGTTCACATCCTTc                                                                  <  1:298004/68‑1 (MQ=255) aTGCCTGCGCGAGCCATAGTGCGCAGTGGCTGGAATTCCAGGTTGCACACGCGCAGTTCACATCCTTc                                                                  <  1:2971945/68‑1 (MQ=255)  tGCCTGCGCGAGCCATAGTGCGCAGTGGCTGGAATTCCAGGTTGCACACGCGCAGTTCACATCCTTc                                                                  <  1:1388054/67‑1 (MQ=255)  tGCCTGCGCGAGCCATAGTGCGCAGTGGCTGGAATTCCAGGTTGCACACGCGCAGTTCACATCCTTc                                                                  <  1:907977/67‑1 (MQ=255)       gcgcGAGCCATAGTGCGCAGTGGCTGGAATTCCACGTTGCACACGCGCAGTTCACATCCCTCGTGCAGa                                                           <  1:2292892/69‑1 (MQ=255)       gcgcGAGCCATAGTGCGCAGTGGCTGGAATTCCACGTTGCACACGCGCAGTTCACATCCCTCGGGCAGa                                                           <  1:2396517/69‑1 (MQ=255)       gcgcGAGCCATAGTGCGCAGTGGCTGGAATTCCACGTTGCACACGCGCAGTTCACATCCCTCGGGCAGa                                                           <  1:1216256/69‑1 (MQ=255)        cgcgAGCCATAGTGCGCAGTGGCTGGAATTCCAGGTTGCACACGCGCAGTTCACATCCTTCGGGCAGAc                                                          >  1:2979344/1‑69 (MQ=255)        cgcgAGCCATAGTGCGCAGTGGCTGGAATTCCAGGTTGCACACGCGCAGTTCACATCCTTCGGGCAGAc                                                          >  1:3332086/1‑69 (MQ=255)        cgcgAGCCATAGTGCGCAGTGGCTGGAATTCCAGGTTGCACACGCGCAGTTCACATCCTTCGGGCAGAc                                                          >  1:462102/1‑69 (MQ=255)        cgcgAGCCATAGTGCGCAGTGGCTGGAATTCCAGGTTGCACACGCGCAGTTCACATCCTTCGGGCAGAc                                                          >  1:1670358/1‑69 (MQ=255)        cgcgAGCCATAGTGCGCAGTGGCTGGAATTCCAGGTTGCACACGCGCAGTTCACATCCTTCGGGCAGAc                                                          >  1:878146/1‑69 (MQ=255)        cgcgAGCCATAGTGCGCAGTGGCTGGAATTCCAGGTTGCACACGCGCAGTTCACATCCTTCGGGCAGAc                                                          >  1:118726/1‑69 (MQ=255)        cgcgAGCCATAGTGCGCAGTGGCTGGAATTCCAGGTTGCACACGCGCAGTTCACATCCTTCGGGCAGAc                                                          >  1:850948/1‑69 (MQ=255)                  aGTGCGCAGTGGCTGGAATTCCAGGTTGCACACGCGCAGTTCACATCCTTCGGGCAGACGCTTCACaaa                                                <  1:2312572/69‑1 (MQ=255)                  aGTGCGCAGTGGCTGGAATTCCAGGTTGCACACGCGCAGTTCACATCCTTCGGGCAGACGCTTCACaaa                                                <  1:1680999/69‑1 (MQ=255)                  aGTGCGCAGTGGCTGGAATTCCAGGTTGCACACGCGCAGTTCACATCCTTCGGGCAGACGCTTCACaaa                                                <  1:2080926/69‑1 (MQ=255)                  aGTGCGCAGTGGCTGGAATTCCAGGTTGCACACGCGCAGTTCACATCCTTCGGGCAGACGCTTCACaaa                                                <  1:1970626/69‑1 (MQ=255)                  aGTGCGCAGTGGCTGGAATTCCAGGTTGCACACGCGCAGTTCACATCCTTCGGGCAGACGCTTCACaaa                                                <  1:2670040/69‑1 (MQ=255)                  aGTGCGCAGTGGCTGGAATTCCAGGTTGCACACGCGCAGTTCACATCCTTCGGGCAGACGCTTCACaaa                                                <  1:3198329/69‑1 (MQ=255)                  aGTGCGCAGTGGCTGGAATTCCAGGTTGCACACGCGCAGTTCACATCCTTCGGGCAGACGCTTCACaaa                                                <  1:2824585/69‑1 (MQ=255)                  aGTGCGCAGTGGCTGGAATTCCAGGTTGCACACGCGCAGTTCACATCCTTCGGGCAGACGCTTCACaaa                                                <  1:1528543/69‑1 (MQ=255)                  aGTGCGCAGTGGCTGGAATTCCAGGTTGCACACGCGCAGTTCACATCCTTCGGGCAGACGCTTCACaaa                                                <  1:1478372/69‑1 (MQ=255)                  aGTGCGCAGTGGCTGGAATTCCAGGTTGCACACGCGCAGTTCACATCCTTCGGGCAGACGCTTCACaaa                                                <  1:1439668/69‑1 (MQ=255)                  aGTGCGCAGTGGCTGGAATTCCAGGTTGCACACGCGCAGTTCACATCCTTCGGGCAGACGCTTCACaaa                                                <  1:1416325/69‑1 (MQ=255)                    tGCGCAGTGGCTGGAATTCCAGGTTGCACACGCGCAGTTCACATCCTTCGGGCAGACGCTTCACaaa                                                <  1:2625741/67‑1 (MQ=255)                     gcgcAGTGGCTGGAATTCCACGTTGCACACGCGCAGTTCACATCCCTCGGGCAGACGCTTCACAAAACg                                             >  1:626793/1‑69 (MQ=255)                     gcgcAGTGGCTGGAATTCCACGTTGCACACGCGCAGTTCACATCCCTCGGGCAGACGCTTCACAAAACg                                             >  1:2023617/1‑69 (MQ=255)                     gcgcAGTGGCTGGAATTCCACGTTGCACACGCGCAGTTCACATCCCTCGGGCAGACGCTTCACAAAACg                                             >  1:3341928/1‑69 (MQ=255)                     gcgcAGTGGCTGGAATTCCACGTTGCACACGCGCAGTTCACATCCCTCGGGCAGACGCTTCACAAAACg                                             >  1:1453585/1‑69 (MQ=255)                     gcgcAGTGGCTGGAATTCCACGTTGCACACGCGCAGTTCACATCCCTCGGGCAGACGCTTCACAAAACg                                             >  1:988150/1‑69 (MQ=255)                     gcgcAGTGGCTGGAATTCCACGTTGCACACGCGCAGTTCACATCCCTCGGGCAGACGCTTCACAAAACg                                             >  1:2428617/1‑69 (MQ=255)                           tGGCTGGAATTCCACGTTGCACACGCGCAGTTCACATCCCTCGGGCAGACGCTTCACAAAACGCTGGaa                                       >  1:3065935/1‑69 (MQ=255)                           tGGCTGGAATTCCACGTTGCACACGCGCAGTTCACATCCCTCGGGCAGACGCTTCACAAAACGCTGGaa                                       >  1:2947900/1‑69 (MQ=255)                           tGGCTGGAATTCCACGTTGCACACGCGCAGTTCACATCCCTCGGGCAGACGCTTCACAAAACGCTGGaa                                       >  1:2987467/1‑69 (MQ=255)                           tGGCTGGAATTCCACGTTGCACACGCGCAGTTCACATCCCTCGGGCAGACGCTTCACAAAACGCTGGaa                                       >  1:2368259/1‑69 (MQ=255)                           tGGCTGGAATTCCACGTTGCACACGCGCAGTTCACATCCCTCGGGCAGACGCTTCACAAAACGCTGGaa                                       >  1:870742/1‑69 (MQ=255)                           tGGCTGGAATTCCACGTTGCACACGCGCAGTTCACATCCCTCGGGCAGACGCTTCACAAAACGCTGGaa                                       >  1:1632060/1‑69 (MQ=255)                           tGGCTGGAATTCCACGTTGCACACGCGCAGTTCACATCCCTCGGGCAGACGCTTCACAAAACGCTGGaa                                       >  1:1517960/1‑69 (MQ=255)                           tGGCTGGAATTCCACGTTGCACACGCGCAGTTCACATCCCTCGGGCAGACGCTTCACAAAACGCTGGaa                                       >  1:1240633/1‑69 (MQ=255)                                      ccACGTTGCACACGCGCAGTTCACATCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGac                            >  1:275480/1‑69 (MQ=255)                                                           cacaTCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGaa       >  1:1167063/1‑69 (MQ=255)                                                            acaTCCTTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGaa       >  1:1431058/1‑68 (MQ=255)                                                            acaTCCTTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGaa       >  1:1304277/1‑68 (MQ=255)                                                            acaTCCTTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAAc      >  1:3302431/1‑69 (MQ=255)                                                            acaTCCTTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAAc      >  1:3085344/1‑69 (MQ=255)                                                            acaTCCTTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAAc      >  1:3391917/1‑69 (MQ=255)                                                            acaTCCTTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAAc      >  1:568852/1‑69 (MQ=255)                                                            acaTCCTTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAAc      >  1:590008/1‑69 (MQ=255)                                                            acaTCCTTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAAc      >  1:2607443/1‑69 (MQ=255)                                                            acaTCCTTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAAc      >  1:1100804/1‑69 (MQ=255)                                                            acaTCCTTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAAc      >  1:1342960/1‑69 (MQ=255)                                                            acaTCCTTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAAc      >  1:2497084/1‑69 (MQ=255)                                                            acaTCCTTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAAc      >  1:2231422/1‑69 (MQ=255)                                                            acaTCCTTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAAc      >  1:1939847/1‑69 (MQ=255)                                                              aTCCCTCGGGCATACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACgg    <  1:2263242/69‑1 (MQ=255)                                                              aTCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACgg    <  1:1936524/69‑1 (MQ=255)                                                              aTCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACgg    <  1:1336876/69‑1 (MQ=255)                                                              aTCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACgg    <  1:1695049/69‑1 (MQ=255)                                                              aTCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACgg    <  1:1711469/69‑1 (MQ=255)                                                              aTCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACgg    <  1:861038/69‑1 (MQ=255)                                                              aTCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACgg    <  1:755112/69‑1 (MQ=255)                                                              aTCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACgg    <  1:685745/69‑1 (MQ=255)                                                              aTCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACgg    <  1:668307/69‑1 (MQ=255)                                                              aTCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACgg    <  1:3042525/69‑1 (MQ=255)                                                              aTCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACgg    <  1:448209/69‑1 (MQ=255)                                                              aTCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACgg    <  1:2017248/69‑1 (MQ=255)                                                              aTCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACgg    <  1:244107/69‑1 (MQ=255)                                                              aTCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACgg    <  1:3268727/69‑1 (MQ=255)                                                               tCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACgg    <  1:757192/68‑1 (MQ=255)                                                               tCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACgg    <  1:994051/68‑1 (MQ=255)                                                                cccTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACGGc   >  1:3163635/1‑68 (MQ=255)                                                                cccTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACGGCa  >  1:1929861/1‑69 (MQ=255)                                                                cccTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACGGCa  >  1:1883532/1‑69 (MQ=255)                                                                cccTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACGGCa  >  1:1799793/1‑69 (MQ=255)                                                                cccTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACGGCa  >  1:3356275/1‑69 (MQ=255)                                                                cccTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACGGCa  >  1:3183650/1‑69 (MQ=255)                                                                  |                                                                   ATGCCAGCGCGAGCCATAGTGCGCAGTGGCTGGAATTCCACGTTGCACACGCGCAGTTCACATCCCTCGGGCAGACGCTTCACAAAACGCTGGAACGCATCAAGACCACCAGCATCAAGTACCGGAACGGCA  >  minE/811402‑811533

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑