Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,506,571	G→A	54.5%	G256G (GGC→GGT)	yfbS ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,506,571	0	G	A	54.5%	12.4 / 83.0	99	G256G (GGC→GGT)	yfbS	predicted transporter
Reads supporting (aligned to +/- strand):  ref base G (29/16);  new base A (43/11);  total (72/27)
Fisher's exact test for biased strand distribution p-value = 1.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAA  >  minE/1506514‑1506635                                                          |                                                                 aaCGTTCACGATAACGCGACGAAACCGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                       >  1:1489817/1‑69 (MQ=255) aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGGTGACGTTAGCGCCATAACGCTc                                                       >  1:417315/1‑69 (MQ=255) aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATACCGCTc                                                       >  1:2186837/1‑69 (MQ=255) aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                       >  1:2561837/1‑69 (MQ=255) aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                       >  1:2082504/1‑69 (MQ=255) aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                       >  1:1984645/1‑69 (MQ=255) aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                       >  1:2186833/1‑69 (MQ=255) aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                       >  1:2234645/1‑69 (MQ=255) aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                       >  1:667626/1‑69 (MQ=255) aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                       >  1:1847352/1‑69 (MQ=255) aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                       >  1:1833306/1‑69 (MQ=255) aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                       >  1:1812344/1‑69 (MQ=255) aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                       >  1:2479736/1‑69 (MQ=255) aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                       >  1:1740732/1‑69 (MQ=255) aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                       >  1:2532766/1‑69 (MQ=255) aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                       >  1:885560/1‑69 (MQ=255) aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                       >  1:2629179/1‑69 (MQ=255) aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                       >  1:1490620/1‑69 (MQ=255) aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                       >  1:1389183/1‑69 (MQ=255) aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                       >  1:587839/1‑69 (MQ=255) aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                       >  1:425160/1‑69 (MQ=255) aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                       >  1:1306921/1‑69 (MQ=255) aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                       >  1:3100386/1‑69 (MQ=255) aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                       >  1:900909/1‑69 (MQ=255) aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                       >  1:3251158/1‑69 (MQ=255) aaCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTc                                                       >  1:101581/1‑69 (MQ=255)  aCGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGCTAGCGCCAt                                                              <  1:1783241/61‑1 (MQ=255)      tCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTa                                                  <  1:1106277/69‑1 (MQ=255)            tAACGCGACGAAAACGCC‑CCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTgag                                           <  1:2536807/69‑1 (MQ=255)             aaCGCGACGAAATCGCCGCCAGCGTTCAACACCGTTGACGTTAGCGCCATAACGCTCACGTAATTTgag                                           <  1:1167207/69‑1 (MQ=255)             aaCGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTgag                                           <  1:218324/69‑1 (MQ=255)              aCGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTgag                                           <  1:1360836/68‑1 (MQ=255)                                 aGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCaa                       <  1:1562906/69‑1 (MQ=255)                                 aGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCaa                       <  1:3325288/69‑1 (MQ=255)                                 aGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCaa                       <  1:1944691/69‑1 (MQ=255)                                 aGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCaa                       <  1:2450911/69‑1 (MQ=255)                                 aGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCaa                       <  1:2846184/69‑1 (MQ=255)                                 aGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCaa                       <  1:1353016/69‑1 (MQ=255)                                  gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACTTAATTTGAGATCGTCCAGTCGTTGACCAAt                      >  1:2404824/1‑69 (MQ=255)                                  gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTgaga                                          >  1:1975993/1‑49 (MQ=255)                                  gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTg                            >  1:94951/1‑63 (MQ=255)                                  gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGCCCAAt                      >  1:1492457/1‑69 (MQ=255)                                  gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCaa                       >  1:1265655/1‑68 (MQ=255)                                  gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                      >  1:3059620/1‑69 (MQ=255)                                  gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                      >  1:2589399/1‑69 (MQ=255)                                  gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                      >  1:3018104/1‑69 (MQ=255)                                  gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                      >  1:895347/1‑69 (MQ=255)                                  gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                      >  1:891384/1‑69 (MQ=255)                                  gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                      >  1:3224354/1‑69 (MQ=255)                                  gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                      >  1:411442/1‑69 (MQ=255)                                  gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                      >  1:345296/1‑69 (MQ=255)                                  gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                      >  1:3355321/1‑69 (MQ=255)                                  gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                      >  1:3429928/1‑69 (MQ=255)                                  gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                      >  1:818378/1‑69 (MQ=255)                                  gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                      >  1:814195/1‑69 (MQ=255)                                  gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                      >  1:770109/1‑69 (MQ=255)                                  gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                      >  1:1996005/1‑69 (MQ=255)                                  gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                      >  1:757971/1‑69 (MQ=255)                                  gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                      >  1:1902827/1‑69 (MQ=255)                                  gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                      >  1:22446/1‑69 (MQ=255)                                  gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                      >  1:687173/1‑69 (MQ=255)                                  gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                      >  1:2546606/1‑69 (MQ=255)                                  gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                      >  1:2403023/1‑69 (MQ=255)                                  gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                      >  1:624061/1‑69 (MQ=255)                                  gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                      >  1:246128/1‑69 (MQ=255)                                  gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                      >  1:1664365/1‑69 (MQ=255)                                  gCGTTCAACACCGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                      >  1:201513/1‑69 (MQ=255)                                  gCGTTCAACACCGATGACGTTAGCACCATAACGCTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                      >  1:2421059/1‑69 (MQ=255)                                  gCGTTCAACACAGATGACGTTAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAAt                      >  1:1746644/1‑69 (MQ=255)                                      tCAACACCGATGACGTTAGCACCATAACGTTCACGTaa                                                 >  1:51494/1‑38 (MQ=38)                                      tCAACACCGATGACGTTAGCACCATAACGTTCACGTaa                                                 >  1:757350/1‑38 (MQ=38)                                      tCAACACCGATGACGTTAGCACCATAACGTTCACGTaa                                                 >  1:578591/1‑38 (MQ=38)                                      tCAACACCGATGACGTTAGCACCATAACGTTCACGTaa                                                 >  1:1084989/1‑38 (MQ=38)                                      tCAACACCGATGACGTTAGCACCATAACGTTCACGTaa                                                 >  1:3325019/1‑38 (MQ=38)                                      tCAACACCGATGACGTTAGCACCATAACGTTCACGTaa                                                 >  1:1564483/1‑38 (MQ=38)                                      tCAACACCGATGACGTTAGCACCATAACGTTCACGTaa                                                 >  1:1899741/1‑38 (MQ=38)                                      tCAACACCGATGACGTTAGCACCATAACGTTCACGTaa                                                 >  1:2272811/1‑38 (MQ=38)                                      tCAACACCGATGACGTTAGCACCATAACGTTCACGTaa                                                 >  1:2387348/1‑38 (MQ=38)                                      tCAACACCGATGACGTTAGCACCATAACGTTCACGTaa                                                 >  1:2701024/1‑38 (MQ=38)                                      tCAACACCGATGACGTTAGCACCATAACGTTCACGTaa                                                 >  1:2984966/1‑38 (MQ=38)                                      tCAACACCGATGACGTTAGCACCATAACGTTCACGTaa                                                 >  1:327901/1‑38 (MQ=38)                                           acCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGa             <  1:134659/69‑1 (MQ=255)                                              gATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATcc          <  1:1590121/69‑1 (MQ=255)                                               aTGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATcc          <  1:626901/68‑1 (MQ=255)                                                 gACGTTTGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCgg        >  1:441444/1‑68 (MQ=255)                                                 gACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCaa                       <  1:1243018/53‑1 (MQ=255)                                                 gACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCggg       >  1:597353/1‑69 (MQ=255)                                                 gACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCggg       >  1:2054223/1‑69 (MQ=255)                                                      tAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGaa  <  1:2084769/69‑1 (MQ=255)                                                      tAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGaa  <  1:959691/69‑1 (MQ=255)                                                      tAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGaa  <  1:1138439/69‑1 (MQ=255)                                                      tAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGaa  <  1:3380156/69‑1 (MQ=255)                                                      tAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGaa  <  1:2052085/69‑1 (MQ=255)                                                      tAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGaa  <  1:2928920/69‑1 (MQ=255)                                                      tAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGaa  <  1:3046915/69‑1 (MQ=255)                                                      tAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGaa  <  1:236819/69‑1 (MQ=255)                                                      tAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGaa  <  1:2425349/69‑1 (MQ=255)                                                      tAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGaa  <  1:2924991/69‑1 (MQ=255)                                                      tAGCACCATAACGTTCACGTAATTTGAGATCGTCCAGTCGTTGACCAATCATTGGCGATCCGGGGCGaa  <  1:526215/69‑1 (MQ=255)                                                          |                                                                 AACGTTCACGATAACGCGACGAAAACGCCGCCAGCGTTCAACACCGATGACGTTAGCGCCATAACGCTCACGTAATTTGAGATCATCCAGCCGTTGACCAATCATTGGCGATCCGGGGCGAA  >  minE/1506514‑1506635

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