Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,985,313	A→G	41.4%	V282V (GTA→GTG)	creD →	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,985,313	0	A	G	41.4%	9.9 / 27.3	41	V282V (GTA→GTG)	creD	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (6/18);  new base G (5/12);  total (11/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.93e-01

TGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATT  >  minE/2985250‑2985376                                                                |                                                                tGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAac                                                               <  1:606502/66‑1 (MQ=255) tGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAac                                                               <  1:3163471/66‑1 (MQ=255)                 ggCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCAGATCAATa                                             <  1:901478/68‑1 (MQ=255)                  gCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATa                                             <  1:3299579/67‑1 (MQ=255)                     aTGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCaa                                         >  1:1158773/1‑68 (MQ=255)                     aTGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCaa                                         >  1:1272382/1‑68 (MQ=255)                     aTGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCaa                                         >  1:38788/1‑68 (MQ=255)                     aTGATACCGGCTGGGAAAACTTCCCCGCGTTTAGCGTCGCAGTGACGACGCCAGCCGATCAATACCaa                                         >  1:3401258/1‑68 (MQ=255)                     aTGAAACCGGCTGGGAAAACTTCCCCGCGTTTCGCGTCGCAGTGACGACGCCAGCCGATCAATACCaa                                         >  1:2910273/1‑68 (MQ=255)                         tACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAac                                                               <  1:8885/42‑1 (MQ=255)                           ccGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGCCGCCAg                                                       >  1:740123/1‑48 (MQ=255)                           ccGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACt                                   >  1:2877790/1‑68 (MQ=255)                           ccGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTg                                  >  1:2334391/1‑69 (MQ=255)                           ccGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTg                                  >  1:548187/1‑69 (MQ=255)                           ccGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTg                                  >  1:1106598/1‑69 (MQ=255)                                 gggAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAAtt                                       <  1:1109056/58‑1 (MQ=255)                                 gggAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAAtt                                       <  1:3068771/58‑1 (MQ=255)                                 gggAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAAtt                                       <  1:812297/58‑1 (MQ=255)                                 gggAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAAtt                                       <  1:2776084/58‑1 (MQ=255)                                       aCTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTa                      <  1:75636/69‑1 (MQ=255)                                       aCTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTa                      <  1:532113/69‑1 (MQ=255)                                        cTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTa                      >  1:2382057/1‑68 (MQ=255)                                                cGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCa             <  1:2361010/69‑1 (MQ=255)                                                cGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCa             <  1:2620354/69‑1 (MQ=255)                                                        gTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTctgct      <  1:1731287/68‑1 (MQ=255)                                                        gTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTctgct      <  1:2000694/68‑1 (MQ=255)                                                        gTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTctgct      <  1:531148/68‑1 (MQ=255)                                                        gTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTctgct      <  1:3028181/68‑1 (MQ=255)                                                        gTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTctgct      <  1:1694492/68‑1 (MQ=255)                                                           gCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCTATTCTGCTGAtt  <  1:2768919/69‑1 (MQ=255)                                                           gCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCTATTCTGCTGAtt  <  1:2773086/69‑1 (MQ=255)                                                           gCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGAtt  <  1:1410077/69‑1 (MQ=255)                                                           gCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGAtt  <  1:1849894/69‑1 (MQ=255)                                                           gCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGAtt  <  1:663646/69‑1 (MQ=255)                                                           gCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGAtt  <  1:2814798/69‑1 (MQ=255)                                                           gCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGAtt  <  1:599354/69‑1 (MQ=255)                                                           gCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGAtt  <  1:2339340/69‑1 (MQ=255)                                                           gCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGAtt  <  1:533424/69‑1 (MQ=255)                                                           gCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGAtt  <  1:336371/69‑1 (MQ=255)                                                           gCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGAtt  <  1:2890506/69‑1 (MQ=255)                                                           gCAGTGACGACGCCAGCCGATCAATACCAATTAACGGACCGGGCGACTAAGTACGCCATTCTGCTGAtt  <  1:2880793/69‑1 (MQ=255)                                                                |                                                                TGAGCGTTTTGCTTCAGGCAATGATACCGGCTGGGAAAACTTCCCGGCGTTTAGCGTCGCAGTAACGACGCCAGCCGATCAATACCAATTAACTGACCGGGCGACTAAGTACGCCATTCTGCTGATT  >  minE/2985250‑2985376

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