Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,615,872	C→T	18.3%	G139G (GGC→GGT)	yhhS →	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,615,872	0	C	T	18.3%	101.4 / 12.3	65	G139G (GGC→GGT)	yhhS	predicted transporter
Reads supporting (aligned to +/- strand):  ref base C (6/47);  new base T (8/4);  total (14/51)
Fisher's exact test for biased strand distribution p-value = 1.97e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGGCGCGTCATCCTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTCACT  >  minE/2615812‑2615932                                                             |                                                             ggggCGCGTCATCCTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGt                                                      <  1:420052/69‑1 (MQ=255) ggggCGCGTCATCCTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGt                                                      <  1:2418081/69‑1 (MQ=255) ggggCGCGTCATCCTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGt                                                      <  1:1818937/69‑1 (MQ=255)    gcgcGTCATCCTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTg                                                     <  1:2556552/67‑1 (MQ=255)    gcgcGTCATCCTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTg                                                     <  1:1496594/67‑1 (MQ=255)         tCATCCTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTGGTTGGCt                                              >  1:2760799/1‑69 (MQ=255)         tCATCCTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTGGTTGGCt                                              >  1:1156380/1‑69 (MQ=255)         tCATCCTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTGGTTGGCt                                              >  1:2095460/1‑69 (MQ=255)         tCATCCTTGGGATTGGGCAAAGTTTAGCCGGAACGGGATCGACCCTATGGGGCGTTg                                                          >  1:949918/1‑57 (MQ=255)           aTCCTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCg                                            >  1:2750475/1‑69 (MQ=255)           aTCCTTGGGATTGGACAAAGTTTTGCCGGAACGGGATCGACCCTGTGGGGTGTTGGCGTGGTTGGTTCg                                            <  1:1103167/69‑1 (MQ=255)           aTCCTTGGGATTGGACAAAGTTTTGCCGGAACGGGATCGACCCTGTGGGGTGTTGGCGTGGTTGGTTCg                                            <  1:439600/69‑1 (MQ=255)           aTCCTTGGGATTGGACAAAGTTTTGCCGGAACGGGATCGACCCTGTGGGGTGTTGGCGTGGTTGGTTCg                                            <  1:1065998/69‑1 (MQ=255)              cTTGGGATTGGGCAAAGTTTTGCCGGAACTGGATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTg                                         <  1:780791/69‑1 (MQ=255)              cTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTg                                         <  1:1278747/69‑1 (MQ=255)              cTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTg                                         <  1:2867951/69‑1 (MQ=255)              cTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTg                                         <  1:2115302/69‑1 (MQ=255)              cTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTg                                         <  1:3013725/69‑1 (MQ=255)              cTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTg                                         <  1:3019057/69‑1 (MQ=255)              cTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTg                                         <  1:3094471/69‑1 (MQ=255)              cTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTg                                         <  1:3364955/69‑1 (MQ=255)              cTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTg                                         <  1:602091/69‑1 (MQ=255)              cTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTg                                         <  1:73154/69‑1 (MQ=255)                            aaGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGt                           >  1:82979/1‑69 (MQ=255)                             aGTTTTGCCGGAACGGGATCGACCCTGTGGGGTGTTGGCGTGGTTggtt                                              >  1:712405/1‑47 (MQ=255)                             aGTTTTGCCGGAACGGGATCGACCCTGTGGGGTGTTGGCGTGGTTGGTTCGCTACATATCGGGCGGGTg                          >  1:685470/1‑69 (MQ=255)                             aGTTTTGCCGGAACGGGATCGACCCTGTGGGGTGTTGGCGTGGTTGGTTCGCTACATATCGGGCGGGTg                          >  1:1801836/1‑69 (MQ=255)                             aGTTTTGCCGGAACGGGATCGACCCTGTGGGGTGTTGGCGTGGTTGGTTCGCTACATATCGGGCGGGTg                          >  1:282370/1‑69 (MQ=255)                             aGTTTTGCCGGAACGGGATCGACCCTGTGGGGTGTTGGCGTGGTTGGTTCGCTACATATCGGGCGGGTg                          >  1:2633798/1‑69 (MQ=255)                             aGTTTTGCCGGAACGGGATCGACCCTGTGGGGTGTTGGCGTGGTTGGTTCGCTACATATCGGGCGGGTg                          >  1:2463124/1‑69 (MQ=255)                             aGTTTTGCCGGAACGGGATCGACCCTGTGGGGTGTTGGCGTGGTTGGTTCGCTACATATCGGGCGGGTg                          >  1:257224/1‑69 (MQ=255)                             aGTTTTGCCGGAACGGGATCGACCCTGTGGGGTGTTCGCGt                                                      >  1:1015621/1‑41 (MQ=37)                                       gAACGGGATCGACCCTGTGGGGTGTTGGCGTGGTTGGTTCGCTACATATCGGGCGGGTGATTTCGTGGa                <  1:2316723/69‑1 (MQ=255)                                            ggATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGTAACGGc           <  1:1180184/69‑1 (MQ=255)                                            ggATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGc           <  1:2516071/69‑1 (MQ=255)                                            ggATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGc           <  1:592903/69‑1 (MQ=255)                                            ggATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGc           <  1:548287/69‑1 (MQ=255)                                            ggATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGc           <  1:112993/69‑1 (MQ=255)                                            ggATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGc           <  1:451382/69‑1 (MQ=255)                                            ggATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGc           <  1:1706325/69‑1 (MQ=255)                                            ggATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGc           <  1:1862796/69‑1 (MQ=255)                                            ggATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGc           <  1:2228810/69‑1 (MQ=255)                                            ggATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGc           <  1:3011885/69‑1 (MQ=255)                                               tCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTg                                         <  1:683493/36‑1 (MQ=255)                                                 gACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGt      <  1:3234041/69‑1 (MQ=255)                                                 gACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGt      <  1:3361016/69‑1 (MQ=255)                                                 gACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGt      <  1:382207/69‑1 (MQ=255)                                                 gACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGt      <  1:386341/69‑1 (MQ=255)                                                 gACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGt      <  1:1106214/69‑1 (MQ=255)                                                  aCCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGt      <  1:1828730/68‑1 (MQ=255)                                                   cccTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTc     <  1:1536501/68‑1 (MQ=255)                                                   cccTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTc     <  1:1698719/68‑1 (MQ=255)                                                   cccTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTc     <  1:2771600/68‑1 (MQ=255)                                                   cccTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTc     <  1:238828/68‑1 (MQ=255)                                                     cTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTCACt  <  1:127259/69‑1 (MQ=255)                                                     cTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTCACt  <  1:2028655/69‑1 (MQ=255)                                                     cTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTCACt  <  1:2772372/69‑1 (MQ=255)                                                     cTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTCACt  <  1:2439204/69‑1 (MQ=255)                                                        tGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTc     <  1:1430471/63‑1 (MQ=255)                                                        tGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTc     <  1:1222961/63‑1 (MQ=255)                                                        tGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTc     <  1:1853432/63‑1 (MQ=255)                                                        tGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTc     <  1:438612/63‑1 (MQ=255)                                                        tGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTc     <  1:1873764/63‑1 (MQ=255)                                                        tGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTc     <  1:2626004/63‑1 (MQ=255)                                                        tGGGGCGTCGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTc     <  1:2046949/63‑1 (MQ=255)                                                             |                                                             GGGGCGCGTCATCCTTGGGATTGGGCAAAGTTTTGCCGGAACGGGATCGACCCTATGGGGCGTTGGCGTGGTTGGCTCGCTGCATATCGGGCGGGTGATTTCGTGGAACGGCATTGTCACT  >  minE/2615812‑2615932

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