Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,441,897	C→T	26.5%	L462L (TTG→TTA)	ilvA ←	threonine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,441,897	0	C	T	26.5%	55.5 / 21.1	34	L462L (TTG→TTA)	ilvA	threonine deaminase
Reads supporting (aligned to +/- strand):  ref base C (15/10);  new base T (4/5);  total (19/15)
Fisher's exact test for biased strand distribution p-value = 4.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

GTCGCCAAGTTCGAACGCCGCCAGTACGCGCCCGTAGTCGGTGCCATGGCTGCGATAGTGGAACAAAGAAATGTTCCAGTACGTACCCAGCGTGTTGAGGAAGCGCAGCAGCGCGCCCG  >  minE/2441834‑2441952                                                                |                                                        gTCGCCAAGTTCGAACGCCGCCAGTACGCGCCCGTAGTCGGTGCCATGGCTGCGATAGTGGAACAAAGaaa                                                  >  1:1193685/1‑71 (MQ=255) gTCGCCAAGTTCGAACGCCGCCAGTACGCGCCCGTAGTCGGTGCCATGGCTGCGATAGTGGAACAAAGaaa                                                  >  1:1745182/1‑71 (MQ=255)            cGAACGCCGCCAGTACGCGCCCGTAGTCGGTGCCATGGCTGCGATAGTGGAACAAAGaaa                                                  >  1:974836/1‑60 (MQ=255)            cGAACGCCGCCAGTACGCGCCCGTAGTCGGTGCCATGGCTGCGATAGTGGAACAAAGaaa                                                  >  1:3588095/1‑60 (MQ=255)            cGAACGCCGCCAGTACGCGCCCGTAGTCGGTGCCATGGCTGCGATAGTGGAACAAAGAAATGTTCCAg                                          >  1:2735608/1‑68 (MQ=255)                 gccgccAGTACGCGCCCGTAGTCGGTGCCATGGCTGCGATAGTGGAATAAAGAAATGTTCCAGTACGTAcc                                  <  1:139757/71‑1 (MQ=255)                      cAGTACGCGCCCGTAGTCGGTGCCATGGCTGCGATAGTGGAACAAAGAAATGTTCCAGTACGTACCCAGCg                             <  1:2456262/71‑1 (MQ=255)                         tACGCGCCCGTAGTCGGTGCCATGGCTGCGATAGTGGAACAAAGAAATGTTcc                                            <  1:1091237/53‑1 (MQ=255)                         tACGCGCCCGTAGTCGGTGCCATGGCTGCGATAGTGGAACAAAGAAATGTTcc                                            <  1:376537/53‑1 (MQ=255)                         tACGCGCCCGTAGTCGGTGCCATGGCTGCGATAGTGGAACAAAGAAATGTTcc                                            <  1:1882454/53‑1 (MQ=255)                          aCGCGCCCGTAGTCGGTGCCATGGCTGCGATAGTGGAACAAAGAAATGTTcc                                            <  1:1381268/52‑1 (MQ=255)                           cgcgCCCGTAGTCGGTGCCATGGCTGCGATAGTGGAATAAAGAAATGTTCCAGTACGTACCCAGCGTGTTg                        <  1:2798218/71‑1 (MQ=255)                           cgcgCCCGTAGTCGGTGCCATGGCTGCGATAGTGGAATAAAGAAATGTTCCAGTACGTACCCAGCGTGTTg                        <  1:1530270/71‑1 (MQ=255)                           cgcgCCCGTAGTCGGTGCCATGGCTGCGATAGTGGAATAAAGAAATGTTCCAGTACGTACCCAGCGTGTTg                        <  1:928822/71‑1 (MQ=255)                           cgcgCCCGTAGTCGGTGCCATGGCTGCGATAGTGGAACAAAGAAATGTTCCa                                           >  1:1881477/1‑52 (MQ=255)                           cgcgCCCGTAGTCGGTGCCATGGCTGCGATAGTGGAACAAAGAAATGTTCCAg                                          >  1:2314804/1‑53 (MQ=255)                           cgcgCCCGTAGTCGGTGCCATGGCTGCGATAGTGGAACAAAGAAATGTTCCAGTACGTACCCAGCGTGtt                         <  1:2670498/70‑1 (MQ=255)                                  gTAGTCGGTGCCATGGCTGCGATAGTGGAACAAAGAAATGTTCCAGTACGTCCCCAGCGTGTTGAg                      >  1:409445/1‑66 (MQ=255)                                  gTAGTCGGTGCCATGGCTGCGATAGTGGAACAAAGAAATGTTCCAGTACGTACCCAGCGTGTTGAg                      >  1:487418/1‑66 (MQ=255)                                  gTAGTCGGTGCCATGGCTGCGATAGTGGAACAAAGAAATGTTCCAGTACGTACCCAGCGTGTTGAg                      >  1:3241330/1‑66 (MQ=255)                                     gTCGGTGCCATGGCTGCGATAGTGGAACAAAGAAATGTTCCAGTACGTACCCAGCGTGTTGAGGAAgcg                >  1:3498044/1‑69 (MQ=255)                                      tCGGTGCCATGGCTGCGATAGTGGAACAAAGAAATGTTCCAGTACGTACCCAGCGTGTTGAGGAAg                  <  1:541741/66‑1 (MQ=255)                                      tCGGTGCCATGGCTCCGATAGTGGAACAAAGAAATGTTCCAGTACGTACCCAGCGTGTTGAGGAAg                  <  1:3544911/66‑1 (MQ=255)                                        ggTGCCATGGCTGCGATAGTGGAATAAAGAAATGTTCCAGTACGTACCCAGCGTGTTGAGGAAGCgcag             >  1:788514/1‑69 (MQ=255)                                        ggTGCCATGGCTGCGATAGTGGAACAAAGAAATGTTCCAGTACGTACCCAGCGTGtt                         <  1:3136550/57‑1 (MQ=255)                                        ggTGCCATGGCTGCGATAGTGGAACAAAGAAATGTTCCAGTACGTACCCAGCGTGTTGAGGAAGCgcag             >  1:3383237/1‑69 (MQ=255)                                        ggTGCCATGGCTGCGATAGTGGAACAAAGAAATGTTCCAGTACGTACCCAGCGTGTTGAGGAAGCgcag             >  1:3154596/1‑69 (MQ=255)                                            ccATGGCTGCGATAGTGGAACAAAGAAATGTTCCAGTACGTACCCAGCGTGTTGa                       >  1:2338997/1‑55 (MQ=255)                                               tGGCTGCGATAGTGGAATAAAGAAATGTTCCAGTACGTACCCa                                <  1:3256851/43‑1 (MQ=255)                                                ggCTGCGATAGTGGAATAAAGAAATGTTCCAGTACGTACCCAGCGTGTTGAGGAAGCGCAGCAgcgc       >  1:1953699/1‑67 (MQ=255)                                                ggCTGCGATAGTGGAATAAAGAAATGTTCCAGTACGTACCCAGCGTGTTGAGGAAGCGCAGCAGCGCGcc    >  1:3416958/1‑70 (MQ=255)                                                ggCTGCGATAGTGGAACAAAGAAATGTTCCAGTACGTACCCAGCgt                            >  1:1065292/1‑46 (MQ=255)                                                ggCTGAGATAGTTGAACAAAGAAATGTTCCAGTACGTACCCAGCGTGTTGAGGAAg                  <  1:3420295/56‑1 (MQ=255)                                                 gCTGCGATAGTGGAATAAAGAAATGTTCCAGTACGTACCCAGCGTGTTGAGGAAGCGCAGCAGCGCGCCCg  >  1:893099/1‑71 (MQ=255)                                                                |                                                        GTCGCCAAGTTCGAACGCCGCCAGTACGCGCCCGTAGTCGGTGCCATGGCTGCGATAGTGGAACAAAGAAATGTTCCAGTACGTACCCAGCGTGTTGAGGAAGCGCAGCAGCGCGCCCG  >  minE/2441834‑2441952

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