Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,505,836	A→G	22.6%	G501G (GGT→GGC)	yfbS ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,505,836	0	A	G	22.6%	53.8 / 10.8	44	G501G (GGT→GGC)	yfbS	predicted transporter
Reads supporting (aligned to +/- strand):  ref base A (23/11);  new base G (5/5);  total (28/16)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-02

AAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAGCCCTTTCACCGCCAGCGCGACACCTCCCGTTTTCTGT  >  minE/1505773‑1505890                                                                |                                                       aaaTAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGcc                                                   <  1:2917957/69‑1 (MQ=255) aaaTAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGcc                                                   <  1:1738058/69‑1 (MQ=255)  aaTAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGcc                                                   <  1:1571196/68‑1 (MQ=255)    tAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAAt                                                >  1:2951204/1‑69 (MQ=255)     aGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATg                                               <  1:1335371/69‑1 (MQ=255)             aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCa                                           >  1:2069350/1‑65 (MQ=255)             aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                                       >  1:1955936/1‑69 (MQ=255)             aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                                       >  1:2642875/1‑69 (MQ=255)             aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                                       >  1:2521224/1‑69 (MQ=255)             aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                                       >  1:2672126/1‑69 (MQ=255)             aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                                       >  1:2240270/1‑69 (MQ=255)             aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                                       >  1:323409/1‑69 (MQ=255)             aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                                       >  1:2050538/1‑69 (MQ=255)             aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                                       >  1:196487/1‑69 (MQ=255)             aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                                       >  1:421036/1‑69 (MQ=255)             aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                                       >  1:1815605/1‑69 (MQ=255)             aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                                       >  1:504616/1‑69 (MQ=255)             aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                                       >  1:648902/1‑69 (MQ=255)             aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                                       >  1:713820/1‑69 (MQ=255)             aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                                       >  1:781215/1‑69 (MQ=255)             aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAg                                       >  1:1030110/1‑69 (MQ=255)             aaCCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATAAg                                       >  1:880651/1‑69 (MQ=255)                     acaagacaaACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAGCCCTTTc                                <  1:2837577/68‑1 (MQ=255)                         gacaaACAAACAACCCAGCATCATATGTGGCCCGTAACCGCCGCCAATGTCCATCAGCCCTTTCACCGc                           <  1:2158003/69‑1 (MQ=255)                         gacaaACAAACAACCCAGCATCATATGTGGCCCGTAACCGCCGCCAATGTCCATCAGCCCTTTCACCGc                           <  1:1551177/69‑1 (MQ=255)                         gacaaACAAACAACCCAGCATCATATGTGGCCCGTAACCGCCGCCAATGTCCATCAGCCCTTTCACCGc                           <  1:1829722/69‑1 (MQ=255)                          acaaacaaacaaCCCAGCATCATATGTGGCCCGTAACCGCCGCCAATGTCCATCAGCCCTTTCACCGc                           <  1:468041/68‑1 (MQ=255)                           caaacaaaCAACCCAGCATCATATGTGGCCCGTAACCGCCGCCAATGTCCATCAGCCCTTTCACCGCCa                         >  1:2798308/1‑69 (MQ=255)                           caaacaaaCAACCCAGCATCATATGTGGCCCGTAACCGCCGCCAATGTCCATCAGCCCTTTCACCGCCa                         >  1:695049/1‑69 (MQ=255)                           caaacaaaCAACCCAGCATCATATGTGGCCCGTAACCGCCGCCAATGTCCATCAGCCCTTTCACCGCCa                         >  1:2944992/1‑69 (MQ=255)                           caaacaaaCAACCCAGCATCATATGTGGCCCGTAACCGCCGCCAATGTCCATCAGCCCTTTCACCGCCa                         >  1:2982651/1‑69 (MQ=255)                           caaacaaaCAACCCAGCATCATATGTGGCCCGTAACCGCCGCCAATGTCCATCAGCCCTTTCACCGCCa                         >  1:2139193/1‑69 (MQ=255)                               caaacaCCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAGCCCTTTCACCGCCAgcgc                     <  1:3138759/69‑1 (MQ=255)                                               atatGTGGCCCGTAACCACCGCCAATGTCCATCAGCCCTTTCACCGCCAGCGCGACACCTCCCGTTTTc     >  1:1268125/1‑69 (MQ=255)                                               atatGTGGCCCGTAACCACCGCCAATGTCCATCAGCCCTTTCACCGCCAGCGCGACACCTCCCGTTTTc     >  1:2547391/1‑69 (MQ=255)                                               atatGTGGCCCGTAACCACCGCCAATGTCCATCAGCCCTTTCACCGCCAGCGCGACACCTCCCGTTTTc     >  1:730763/1‑69 (MQ=255)                                               atatGTGGCCCGTAACCACCGCCAATGTCCATCAGCCCTTTCACCGCCAGCGCGACACCTCCCGTTTTc     >  1:2325473/1‑69 (MQ=255)                                                tatGTGGCCCGTAACCGCCGCCAATGTCCATCAGCCCTTTCACCGc                           <  1:3180439/46‑1 (MQ=255)                                                  tgtgGCCCGTAACCACCGCCAATGTCCATCAGCCCTTTCACCGCCAGCGCGACACCTCCCGTTTTCTGt  <  1:1292476/69‑1 (MQ=255)                                                  tgtgGCCCGTAACCACCGCCAATGTCCATCAGCCCTTTCACCGCCAGCGCGACACCTCCCGTTTTCTGt  <  1:2921975/69‑1 (MQ=255)                                                  tgtgGCCCGTAACCACCGCCAATGTCCATCAGCCCTTTCACCGCCAGCGCGACACCTCCCGTTTTCTGt  <  1:1285552/69‑1 (MQ=255)                                                  tgtgGCCCGTAACCACCGCCAATGTCCATCAGCCCTTTCACCGCCAGCGCGACACCTCCCGTTTTCTGt  >  1:858191/1‑69 (MQ=255)                                                  tgtgGCCCGTAACCACCGCCAATGTCCATCAGCCCTTTCACCGCCAGCGCGACACCTCCCGTTTTCTGt  <  1:2466504/69‑1 (MQ=255)                                                  tgtgGCCCGTAACCACCGCCAATGTCCATCAGACCTTTCACCGCCAGCGCGACACCTCCCGTTTTCTGt  <  1:1288718/69‑1 (MQ=255)                                                                |                                                       AAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCATCAGCCCTTTCACCGCCAGCGCGACACCTCCCGTTTTCTGT  >  minE/1505773‑1505890

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