Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	712,419	C→T	35.5%	G339G (GGG→GGA)	pyrC ←	dihydro‑orotase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	712,419	0	C	T	35.5%	60.5 / 50.0	87	G339G (GGG→GGA)	pyrC	dihydro‑orotase
Reads supporting (aligned to +/- strand):  ref base C (23/33);  new base T (26/5);  total (49/38)
Fisher's exact test for biased strand distribution p-value = 1.14e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

GTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTG  >  minE/712353‑712478                                                                   |                                                            gTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTcccc                                                           >  1:2495100/1‑69 (MQ=255) gTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTcccc                                                           >  1:2507705/1‑69 (MQ=255) gTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTcccc                                                           >  1:154891/1‑69 (MQ=255)  tAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCg                                                          >  1:827215/1‑69 (MQ=255)  tAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCg                                                          >  1:3111248/1‑69 (MQ=255)  tAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCg                                                          >  1:1150256/1‑69 (MQ=255)  tAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCg                                                          >  1:1221007/1‑69 (MQ=255)  tAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCg                                                          >  1:2180597/1‑69 (MQ=255)  tAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCg                                                          >  1:63649/1‑69 (MQ=255)  tAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCg                                                          >  1:1637783/1‑69 (MQ=255)        aaCCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAgg                                                    <  1:1559871/69‑1 (MQ=255)        aaCCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAgg                                                    <  1:592064/69‑1 (MQ=255)        aaCCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAgg                                                    <  1:2286331/69‑1 (MQ=255)        aaCCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAgg                                                    <  1:2298448/69‑1 (MQ=255)        aaCCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAgg                                                    <  1:281394/69‑1 (MQ=255)        aaCCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAgg                                                    <  1:1553024/69‑1 (MQ=255)        aaCCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAgg                                                    <  1:1229498/69‑1 (MQ=255)        aaCCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAgg                                                    <  1:130345/69‑1 (MQ=255)        aaCCTAATCACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGAGAgg                                                    <  1:2858582/69‑1 (MQ=255)         aCCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTccc                                                            >  1:3006988/1‑60 (MQ=255)         aCCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGa                                                   <  1:811427/69‑1 (MQ=255)         aCCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGa                                                   >  1:3027665/1‑69 (MQ=255)         aCCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGa                                                   >  1:310700/1‑69 (MQ=255)         aCCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGa                                                   >  1:3031687/1‑69 (MQ=255)         aCCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGa                                                   >  1:1321780/1‑69 (MQ=255)         aCCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGa                                                   <  1:1893386/69‑1 (MQ=255)         aCCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGa                                                   >  1:1937625/1‑69 (MQ=255)             aaTGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATgg                                               <  1:3131372/69‑1 (MQ=255)             aaTGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATgg                                               <  1:3044277/69‑1 (MQ=255)             aaTGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATgg                                               <  1:986120/69‑1 (MQ=255)             aaTGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATgg                                               <  1:173448/69‑1 (MQ=255)                gACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAgg                                                    <  1:2559627/61‑1 (MQ=255)                gACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAgg                                                    <  1:2891030/61‑1 (MQ=255)                gACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAgg                                                    <  1:586526/61‑1 (MQ=255)                gACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAgg                                                    <  1:2993483/61‑1 (MQ=255)                gACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAgg                                                    <  1:1616658/61‑1 (MQ=255)                gACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAgg                                                    <  1:1258830/61‑1 (MQ=255)                gACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAgg                                                    <  1:939786/61‑1 (MQ=255)                     caGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCTCCGGCGAgg                                                    >  1:2442271/1‑56 (MQ=255)                     caGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCTCCGGCGAg                                                     >  1:660411/1‑55 (MQ=255)                     caGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCTCCGGCGAGGa                                                   >  1:1821835/1‑57 (MQ=255)                     caGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCTCCGGCGAGGa                                                   >  1:1142597/1‑57 (MQ=255)                     caGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCTCCGGCGAGGa                                                   >  1:3164147/1‑57 (MQ=255)                     caGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCTCCGGCGAGGa                                                   >  1:2774519/1‑57 (MQ=255)                     caGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCTCCGGCGAGGa                                                   >  1:2561701/1‑57 (MQ=255)                     caGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCTCCGGCGAGGa                                                   >  1:1702328/1‑57 (MQ=255)                     caGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCTCCGGCGAGGa                                                   >  1:362116/1‑57 (MQ=255)                     caGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCTCCGGCGAGGa                                                   >  1:474382/1‑57 (MQ=255)                     caGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCTCCGGCGAGGa                                                   >  1:478470/1‑57 (MQ=255)                     caGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCTCCGGCGAGGa                                                   >  1:573843/1‑57 (MQ=255)                          aGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAgg                                                    >  1:659520/1‑51 (MQ=255)                             tACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCTCCGGCGAGGAATGGCACCAGCGTGtcatca                               <  1:2059913/69‑1 (MQ=255)                             tACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCTCCGGCGAGGAATGGCACCAGCGTGtcatca                               <  1:1613270/69‑1 (MQ=255)                             tACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCTCCGGCGAGGAATGGCACCAGCGTGtcatca                               <  1:1279749/69‑1 (MQ=255)                             tACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCTCCGGCGAGGAATGGCACCAGCGTGtcatca                               <  1:1293159/69‑1 (MQ=255)                             tACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCTCCGGCGAGGAATGGCACCAGCGTGtcatca                               <  1:861115/69‑1 (MQ=255)                             tACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGtcatca                               >  1:2756585/1‑69 (MQ=255)                             tACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGtcatca                               >  1:2752508/1‑69 (MQ=255)                              aCGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGtcatca                               >  1:677678/1‑68 (MQ=255)                                      tATTGTTTAACGGACCAGGGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGc                       <  1:1953921/68‑1 (MQ=255)                                      tATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGc                       <  1:671832/68‑1 (MQ=255)                                      tATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGc                       <  1:1548370/68‑1 (MQ=255)                                      tATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGc                       <  1:2486546/68‑1 (MQ=255)                                      tATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGc                       <  1:1982603/68‑1 (MQ=255)                                      tATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGc                       <  1:2049629/68‑1 (MQ=255)                                               aCGGACCAGCGTACCGTTTCCCCGGCGAGGAATgg                                               <  1:3050137/35‑1 (MQ=255)                                               aCGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTc              >  1:1551386/1‑68 (MQ=255)                                               aCGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTc              >  1:2098487/1‑68 (MQ=255)                                               aCGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTc              >  1:2881276/1‑68 (MQ=255)                                               aCGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTc              >  1:915698/1‑68 (MQ=255)                                               aCGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTc              >  1:3246039/1‑68 (MQ=255)                                               aCGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTc              >  1:3346001/1‑68 (MQ=255)                                               aCGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTc              >  1:865704/1‑68 (MQ=255)                                                        cGTACCGTTTCTCCGGGGAGGAATGGCACCAGCGTCTCATCAGTCAGTGCGATGCTTTCAGCAACCTGt    >  1:1061122/1‑69 (MQ=255)                                                        cGTACCGTTTCTCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGt    >  1:77216/1‑69 (MQ=255)                                                        cGTACCGTTTCTCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGt    >  1:879039/1‑69 (MQ=255)                                                        cGTACCGTTTCTCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGt    >  1:995736/1‑69 (MQ=255)                                                        cGTACCGTTTCTCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGt    >  1:2881854/1‑69 (MQ=255)                                                        cGTACCGTTTCTCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGt    >  1:2073347/1‑69 (MQ=255)                                                        cGTACCGTTTCTCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGt    >  1:2030161/1‑69 (MQ=255)                                                        cGTACCGTTTCTCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGt    >  1:1952087/1‑69 (MQ=255)                                                        cGTACCGTTTCTCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGt    >  1:1795510/1‑69 (MQ=255)                                                        cGTACCGTTTCTCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGt    >  1:167968/1‑69 (MQ=255)                                                        cGTACCGTTTCTCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGt    >  1:1555636/1‑69 (MQ=255)                                                        cGTACCGTTTCTCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGt    >  1:1462405/1‑69 (MQ=255)                                                        cGTACCGTTTCTCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGAAGCTTTCAGCAACCTGt    >  1:2128131/1‑69 (MQ=255)                                                        cGTACCGTTTCTCCGGCGAGGAATGGCACAAGCGTGTCATCAGTCAGTGCGATGCTTTCa             >  1:904537/1‑60 (MQ=255)                                                        cGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGt    <  1:2425486/69‑1 (MQ=255)                                                        cGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGt    <  1:2193283/69‑1 (MQ=255)                                                        cGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGt    <  1:1120940/69‑1 (MQ=255)                                                          tACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTg  <  1:2985162/69‑1 (MQ=255)                                                                   |                                                            GTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTG  >  minE/712353‑712478

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑