Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,416,365	C→A	43.8%	G454G (GGG→GGT)	yifK ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,416,365	0	C	A	43.8%	28.0 / 58.7	63	G454G (GGG→GGT)	yifK	predicted transporter
Reads supporting (aligned to +/- strand):  ref base C (22/13);  new base A (3/24);  total (25/38)
Fisher's exact test for biased strand distribution p-value = 5.79e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAATAAATTGCGCGTTTGGTTACGGTTTGTGCGTTTTGCTGCTTATTCCTCCAGTTTATGCGCTTTCCCGTGGCGATTAAGGCCAAAAACTTTATAAATCGCCGTCACCGCCAGCATAAAGATGA  >  minE/2416299‑2416423                                                                   |                                                           aaataaatTGCGCGTTTGGTTACGGTTTGTGCGTTTTGCTGCTTATTCCTCCAGTTTATGCGCTTTcc                                                           <  1:798961/68‑1 (MQ=255) aaataaatTGCGCGTTTGGTTACGGTTTGTGCGTTTTGCTGCTTATTCCTCCAGTTTATGCGCTTTcc                                                           <  1:2525542/68‑1 (MQ=255) aaataaatTGCGCGTTTGGTTACGGTTTGTGCGTTTTGCTGCTTATTCCTCCAGTTTATGCGCTTTcc                                                           <  1:3158987/68‑1 (MQ=255)      aatTGCGCGTTTGGTTACGGTTTGTGCGTTTTGCTGCTTATTCTTCCAGTTTATGCGCTTTACCGTGGc                                                     <  1:2726645/69‑1 (MQ=255)      aatTGCGCGTTTGGTTACGGTTTGTGCGTTTTGCTGCTTATTCCTCCAGTTTATGCGCTTTCCCGTGGc                                                     <  1:335283/69‑1 (MQ=255)      aatTGCGCGTTTGGTTACGGTTTGTGCGTTTTGCTGCTTATTCCTCCAGTTTATGCGCTTTCCCGTGGc                                                     <  1:2925158/69‑1 (MQ=255)      aatTGCGCGTTTGGTTACGGTTTGTGCGTTTTGCTGCTTATTCCTCCAGTTTATGCGCTTTACCGTGGc                                                     <  1:229317/69‑1 (MQ=255)      aatTGCGCGTTTGGTTACGGTTTGTGCGTTTTGCTGCTTATTCCTCCAGTTTATGCGCTTTACCGTGGc                                                     <  1:2723660/69‑1 (MQ=255)      aatTGCGCGTTTGGTTACGGTTTGTGCGTTTTGCTGCTTATTCCTCCAGTTTATGCGCTTTACCGTGGc                                                     <  1:2331398/69‑1 (MQ=255)      aatTGCGCGTTTGGTTACGGTTTGTGCGTTTTGCTGCTTATTCCTCCAGTTTATGCGCTTTACCGTGGc                                                     <  1:2525934/69‑1 (MQ=255)      aatTGCGCGTTTGGTTACGGTTTGTGCGTTTTGCTGCTTATTCCTCCAGTTTATGCGCTTTACCGTGGc                                                     <  1:2682571/69‑1 (MQ=255)      aatTGCGCGTTTGGTTACGGTTTGTGCGTTTTGCTGCTTATTCCTCCAGTTTATGCGCTTTACCGTGGc                                                     <  1:270988/69‑1 (MQ=255)      aatTGCGCGTTTGGTTACGGTTTGTGCGTTTTGCTGCTTATTCCTCCAGTTTATGCGCTTTACCGTGGc                                                     <  1:2754460/69‑1 (MQ=255)      aatTGCGCGTTTGGTTACGGTTTGTGCGTTTTGCTGCTTATTCCTCCAGTTTATGCGCTTTACCGTGGc                                                     <  1:2973606/69‑1 (MQ=255)      aatTGCGCGTTTGGTTACGGTTTGTGCGTTTTGCTGCTTATTCCTCCAGTTTATGCGCTTTACCGTGGc                                                     <  1:3015496/69‑1 (MQ=255)      aatTGCGCGTTTGGTTACGGTTTGTGCGTTTTGCTGCTTATTCCTCCAGTTTATGCGCTTTACCGTGGc                                                     <  1:3210296/69‑1 (MQ=255)      aatTGCGCGTTTGGTTACGGTTTGTGCGTTTTGCTGCTTATTCCTCCAGTTTATGCGCTTTACCGTGGc                                                     <  1:3264055/69‑1 (MQ=255)      aatTGCGCGTTTGGTTACGGTTTGTGCGTTTTGCTGCTTATTCCTCCAGTTTATGCGCTTTACCGTGGc                                                     <  1:3338482/69‑1 (MQ=255)      aatTGCGCGTTTGGTTACGGTTTGTGCGTTTTGCTGCTTATTCCTCCAGTTTATGCGCTTTACCGTGGc                                                     <  1:411179/69‑1 (MQ=255)      aatTGCGCGTTTGGTTACGGTTTGTGCGTTTTGCTGCTTATTCCTCCAGTTTATGCGCTTTACCGTGGc                                                     <  1:568657/69‑1 (MQ=255)      aatTGCGCGTTTGGTTACGGTTTGTGCGTTTTGCTGCTTATTCCTCCAGTTTATGCGCTTTACCGTGGc                                                     <  1:689976/69‑1 (MQ=255)      aatTGCGCGTTTGGTTACGGTTTGTGCGTTTTGCTGCTTATTCCTCCAGTTTATGCGCTTTACCGTGGc                                                     <  1:796343/69‑1 (MQ=255)      aatTGCGCGTTTGGTTACGGTTTGTGCGTTTTGCTGCTTATTCCTCCAGTTTATGCGCTTTACCGTGGc                                                     <  1:1099341/69‑1 (MQ=255)      aatTGCGCGTTTGGTTACGGTTTGTGCGTTTTGCTGCTTATTCCTCCAGTTTATGCGCTTTACCGTGGc                                                     <  1:1204966/69‑1 (MQ=255)      aatTGCGCGTTTGGTTACGGTTTGTGCGTTTTGCTGCTTATTCCTCCAGTTTATGCGCTTTACCGTGGc                                                     <  1:1971597/69‑1 (MQ=255)      aatTGCGCGTTTGGTTACGGTTTGTGCGTTTTGCTGCTTATTCCTCCAGTTTATGCGCTTTACCGTGGc                                                     <  1:1264920/69‑1 (MQ=255)      aatTGCGCGTTTGGTTACGGTTTGTGCGTTTTGCTGCTTATTCCTCCAGTTTATGCGCTTTACCGTGGc                                                     <  1:2119788/69‑1 (MQ=255)      aatTGCGCGTTTGGTTACGGTTTGTGCGTTTTGCTGCTTATTCCTCCAGTTTATGCGCTTTACCGTGGc                                                     <  1:1698456/69‑1 (MQ=255)      aatTGCGCGTTTGGTTACGGTTTGTGCGTTTTGCTGCTTATTCCTCCAGTTTATGCGCTTTACCGTGGc                                                     <  1:2196382/69‑1 (MQ=255)               tttGGTTACGGTTTGTGCGTTTTGCTGCTTATTCCTCCAGTTTATGCGCTTTTCCGTGGc                                                     <  1:2762242/60‑1 (MQ=255)                 tGGTTACGGTTTGTGCGTTTTGCTGCTTATTCCTCCAGTTTATGCGCTTTCCCGTGGCGATTAAGGcc                                           >  1:700501/1‑68 (MQ=255)                      aCGGTTTGTGCGTTTTGCTGCTTATTCCTCCAGTTTATGCGCTTTCCCGTGGCGATTAAGGCCAAAAAc                                     >  1:789970/1‑69 (MQ=255)                      aCGGTTTGTGCGTTTTGCTGCTTATTCCTCCAGTTTATGCGCTTTCCCGTGGCGATTAAGGCCAAAAAc                                     >  1:2781531/1‑69 (MQ=255)                      aCGGTTTGTGCGTTTTGCTGCTTATTCCTCCAGTTTATGCGCTTTCCCGTGGCGATTAAGGCCAAAAAc                                     >  1:491992/1‑69 (MQ=255)                                      gctgctTATTCCTCCAGTTTATGCGCTTTACCGTGGCGATTAAGCCCAAAAACTTTATAAATCGCCGTc                     >  1:710841/1‑69 (MQ=255)                                      gctgctTATTCCTCCAGTTTATGCGCTTTACCGTGGCGATTAAGCCCAAAAACTTTATAAATCGCCGTc                     >  1:3212445/1‑69 (MQ=255)                                      gctgctTATTCCTCCAGTTTATGCGCTTTACCGTGGCGATTAAGCCCAAAAACTTTATAAATCGCCGTc                     >  1:1206179/1‑69 (MQ=255)                                       ctgctTATTCCTCCAGTTTATGCGCTTTCCCGTGGCGATTAAGGCCAAAAACTTTATAAATCGCCGTc                     <  1:1435891/68‑1 (MQ=255)                                       ctgctTATTCCTCCAGTTTATGCGCTTTCCCGTGGCGATTAAGGCCAAAAACTTTATAAATCGCCGTc                     <  1:1198571/68‑1 (MQ=255)                                       ctgctTATTCCTCCAGTTTATGCGCTTTCCCGTGGCGATTAAGGCCAAAAACTTTATAAATCGCCGTc                     <  1:1974471/68‑1 (MQ=255)                                            tATTCCTCCAGTTTATGCGCTTTCCCGTGGCGATTAAGGCCAAAAACTTTATAAATCGCCGTCACCGcc               >  1:2698809/1‑69 (MQ=255)                                            tATTCCTCCAGTTTATGCGCTTTCCCGTGGCGATTAAGGCCAAAAACTTTATAAATCGCCGTCACCGcc               >  1:2610232/1‑69 (MQ=255)                                            tATTCCTCCAGTTTATGCGCTTTCCCGTGGCGATTAAGGCCAAAAACTTTATAAATCGCCGTCACCGcc               >  1:2229016/1‑69 (MQ=255)                                              ttCCTCCAGTTTATGCGCTTTCCCGTGGCGATTAAGGCCAAAAACTTTATAAATCGCCGt                      <  1:2638377/60‑1 (MQ=255)                                              ttCCTCCAGTTTATGCGCTTTCCCGTGGCGATTAAGGCCAAAAACTTTATAAATCGCCGt                      <  1:1138212/60‑1 (MQ=255)                                              ttCCTCCAGTTTATGCGCTTTCCCGTGGCGATTAAGGCCAAAAACTTTATAAATCGCCGt                      <  1:3059573/60‑1 (MQ=255)                                              ttCCTCCAGTTTATGCGCTTTCCCGTGGCGATTAAGGCCAAAAACTTTATAAATCGCCGt                      <  1:3158632/60‑1 (MQ=255)                                              ttCCTCCAGTTTATGCGCTTTCCCGTGGCGATTAAGGCCAAAAACTTTATAAATCGCCGt                      <  1:3131350/60‑1 (MQ=255)                                                   ccAGTTTATGCGCTTTCCCGTGGCGATTAAGGCCAAAAACTTTATAAATCGCCGTCACCTCCAGCa           >  1:1768096/1‑66 (MQ=255)                                                   ccAGTTTATGCGCTTTCCCGTGGCGATTAAGGCCAAAAACTTTATAAATCGCCGTCACCGCCAGCa           >  1:1602467/1‑66 (MQ=255)                                                   ccAGTTTATGCGCTTTCCCGTGGCGATTAAGGCCAAAAACTTTATAAATCGCCGTCACCGCCAGCa           >  1:1957209/1‑66 (MQ=255)                                                   ccAGTTTATGCGCTTTCCCGTGGAGATTAAGGCCAAAAACTTTATAAATCGCCGTCACCGCCAGCa           >  1:2533535/1‑66 (MQ=255)                                                       tttATGCGCTTTCCCGTGGCGATTAAGGCCAAAAACTTTATAAATCGCCGTCACCGCCAGCATAAAgat    <  1:2727074/69‑1 (MQ=255)                                                       tttATGCGCTTTCCCGTGGCGATTAAGGCCAAAAACTTTATAAATCGCCGTCACCGCCAGCATAAAgat    <  1:2731028/69‑1 (MQ=255)                                                       tttATGCGCTTTCCCGTGGCGATTAAGGCCAAAAACTTTATAAATCGCCGTCACCGCCAGCATAAAgat    <  1:68898/69‑1 (MQ=255)                                                         tATGCGCTTTCCCGTGGCGATTAAGGCCAAAAACTTTATAAATCGCCGTCACCGCCAGCATAAAgatga  >  1:2200497/1‑69 (MQ=255)                                                         tATGCGCTTTCCCGTGGCGATTAAGGCCAAAAACTTTATAAATCGCCGTCACCGCCAGCATAAAgatga  >  1:3131017/1‑69 (MQ=255)                                                         tATGCGCTTTCCCGTGGCGATTAAGGCCAAAAACTTTATAAATCGCCGTCACCGCCAGCATAAAgatga  >  1:2126537/1‑69 (MQ=255)                                                         tATGCGCTTTCCCGTGGCGATTAAGGCCAAAAACTTTATAAATCGCCGTCACCGCCAGCATAAAgatga  >  1:2596052/1‑69 (MQ=255)                                                         tATGCGCTTTCCCGTGGCGATTAAGGCCAAAAACTTTATAAATCGCCGTCACCGCCAGCATAAAgatga  >  1:1022623/1‑69 (MQ=255)                                                         tATGCGCTTTCCCGTGGCGATTAAGGCCAAAAACTTTATAAATCGCCGTCACCGCCAGCATAAAgatga  >  1:2748191/1‑69 (MQ=255)                                                         tATGCGCTTTCCCGTGGCGATTAAGGCCAAAAACTTTATAAATCGCCGTCACCGCCAGCATAAAgatga  >  1:1229125/1‑69 (MQ=255)                                                         tATGCGCTTTCCCGTGGCGATTAAGGCCAAAAACTTTATAAATCGCCGTCACCGCCAGCATAAAgatga  >  1:2892703/1‑69 (MQ=255)                                                         tATGCGCTTTCCCGTGGCGATTAAGGCCAAAAACTTTATAAATCGCCGTCACCGCCAGCATAAAgatga  >  1:2898776/1‑69 (MQ=255)                                                         tATGCGCTTTCCCGTGGCGATTAAGGCCAAAAACTTTATAAATCGCCGTCACCGCCAGCATAAAgatga  >  1:1939206/1‑69 (MQ=255)                                                         tATGCGCTTTCCCGTGGCGATTAAGGCCAAAAACTTTATAAATCGCCGTCACCGCCAGCATAAAgatg   >  1:1226079/1‑68 (MQ=255)                                                                   |                                                           AAATAAATTGCGCGTTTGGTTACGGTTTGTGCGTTTTGCTGCTTATTCCTCCAGTTTATGCGCTTTCCCGTGGCGATTAAGGCCAAAAACTTTATAAATCGCCGTCACCGCCAGCATAAAGATGA  >  minE/2416299‑2416423

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