Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,628,386	T→C	67.8%	P220P (CCA→CCG)	ypfG ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,628,386	0	T	C	67.8%	58.0 / 32.8	68	P220P (CCA→CCG)	ypfG	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (14/8);  new base C (28/18);  total (42/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCACAGCTAATCATCATCAGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGT  >  minE/1628321‑1628444                                                                  |                                                           tCACAGCTAATCATCATCAGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATg                                                           <  1:975939/67‑1 (MQ=255) tCACAGCTAATCATCATCAGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATg                                                           <  1:966513/67‑1 (MQ=255) tCACAGCTAATCATCATCAGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATg                                                           <  1:1030199/67‑1 (MQ=255) tCACAGCTAATCATCATCAGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATg                                                           <  1:399788/67‑1 (MQ=255) tCACAGCTAATCATCATCAGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATg                                                           <  1:395860/67‑1 (MQ=255) tCACAGCTAATCATCATCAGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATg                                                           <  1:3158880/67‑1 (MQ=255) tCACAGCTAATCATCATCAGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATg                                                           <  1:258446/67‑1 (MQ=255)                  caGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAg                                        >  1:985373/1‑69 (MQ=255)                  caGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAg                                        >  1:956431/1‑69 (MQ=255)                  caGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAg                                        >  1:2270898/1‑69 (MQ=255)                                     agCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGag                                                <  1:3088599/42‑1 (MQ=255)                                      gcgcAGTGACATTCACATCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTc                    >  1:391050/1‑69 (MQ=255)                                             gACATTCACCTCATGACGCAATGGATCAAGCGAGCAGCGCAGACCa                                    >  1:2944150/1‑46 (MQ=255)                                              aCATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATc             >  1:2888668/1‑68 (MQ=255)                                                aTTCACCTCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              <  1:1362702/65‑1 (MQ=255)                                                aTTCACCTCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              <  1:3137395/65‑1 (MQ=255)                                                aTTCACCTCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              <  1:3343613/65‑1 (MQ=255)                                                aTTCACCTCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              <  1:2771362/65‑1 (MQ=255)                                                aTTCACCTCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              <  1:351279/65‑1 (MQ=255)                                                aTTCACCTCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              <  1:377655/65‑1 (MQ=255)                                                aTTCACCTCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              <  1:425384/65‑1 (MQ=255)                                                aTTCACCTCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              <  1:2215407/65‑1 (MQ=255)                                                aTTCACCTCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              <  1:486458/65‑1 (MQ=255)                                                aTTCACCTCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              <  1:2096289/65‑1 (MQ=255)                                                aTTCACCTCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              <  1:1948788/65‑1 (MQ=255)                                                aTTCACCTCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              <  1:1442692/65‑1 (MQ=255)                                                aTTCACCTCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              <  1:782431/65‑1 (MQ=255)                                                aTTCACCTCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              <  1:819451/65‑1 (MQ=255)                                                aTTCACCTCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              <  1:1681432/65‑1 (MQ=255)                                                aTTCACCTCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              <  1:1023689/65‑1 (MQ=255)                                                 ttCACCTCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              <  1:2162789/64‑1 (MQ=255)                                                  tCACCTCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              <  1:2414318/63‑1 (MQ=255)                                                   cACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCaa       >  1:2244380/1‑69 (MQ=255)                                                   cACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCaa       >  1:393782/1‑69 (MQ=255)                                                       tCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              >  1:1010909/1‑58 (MQ=255)                                                       tCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              >  1:1358577/1‑58 (MQ=255)                                                       tCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              >  1:3096793/1‑58 (MQ=255)                                                       tCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              >  1:1015971/1‑58 (MQ=255)                                                       tCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              >  1:3377761/1‑58 (MQ=255)                                                       tCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              >  1:81655/1‑58 (MQ=255)                                                       tCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              >  1:1164926/1‑58 (MQ=255)                                                       tCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              >  1:1069727/1‑58 (MQ=255)                                                       tCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              >  1:604797/1‑58 (MQ=255)                                                       tCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              >  1:2312823/1‑58 (MQ=255)                                                       tCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              >  1:1517481/1‑58 (MQ=255)                                                       tCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              >  1:154363/1‑58 (MQ=255)                                                       tCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              >  1:1636845/1‑58 (MQ=255)                                                       tCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              >  1:1727356/1‑58 (MQ=255)                                                       tCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              >  1:1944028/1‑58 (MQ=255)                                                       tCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              >  1:2061405/1‑58 (MQ=255)                                                       tCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              >  1:2179851/1‑58 (MQ=255)                                                       tCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              >  1:226834/1‑58 (MQ=255)                                                       tCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              >  1:1301207/1‑58 (MQ=255)                                                       tCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              >  1:3068855/1‑58 (MQ=255)                                                       tCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              >  1:2527519/1‑58 (MQ=255)                                                       tCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              >  1:1248681/1‑58 (MQ=255)                                                       tCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              >  1:2667735/1‑58 (MQ=255)                                                       tCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              >  1:2807311/1‑58 (MQ=255)                                                       tCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              >  1:2808435/1‑58 (MQ=255)                                                       tCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              >  1:2858975/1‑58 (MQ=255)                                                       tCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGCCAGTTTCCATAAt              >  1:304776/1‑58 (MQ=255)                                                       tCACGACGTAACGGATCAAGCGAACAGCGCAGACCATTCATCCGACAGTTTCCATAAt              >  1:1171009/1‑58 (MQ=255)                                                        cACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCg   >  1:2344313/1‑68 (MQ=255)                                                        cACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCg   >  1:248232/1‑68 (MQ=255)                                                        cACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCg   >  1:190348/1‑68 (MQ=255)                                                        cACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGt  >  1:589045/1‑69 (MQ=255)                                                        cACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGt  >  1:3149325/1‑69 (MQ=255)                                                        cACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGt  >  1:1023876/1‑69 (MQ=255)                                                                  |                                                           TCACAGCTAATCATCATCAGCGCTTTATCATCAGTCAGCGCAGTGACATTCACCTCACGACGCAATGGATCAAGCGAGCAGCGCAGACCATTCATCCGCCAGTTTCCATAATCCAGCAAATCGT  >  minE/1628321‑1628444

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