Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	712,689	A→T	52.6%	G249G (GGT→GGA)	pyrC ←	dihydro‑orotase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	712,689	0	A	T	52.6%	11.0 / 112.1	95	G249G (GGT→GGA)	pyrC	dihydro‑orotase
Reads supporting (aligned to +/- strand):  ref base A (18/27);  new base T (46/4);  total (64/31)
Fisher's exact test for biased strand distribution p-value = 4.91e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

TGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACG  >  minE/712625‑712754                                                                 |                                                                  tGAAACAGCCCGCGCAGCCGCAGCTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTTCCg                                                                >  1:2357481/1‑68 (MQ=255) tGAAACAGCCCGCGCAGCCGCAGCTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTTCCg                                                                >  1:2357314/1‑68 (MQ=255) tGAAACAGCCCGCGCAGCCGCAGCTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTTCCg                                                                >  1:1589907/1‑68 (MQ=255) tGAAACAGCCCGCGCAGCCGCAGCTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTTCCg                                                                >  1:1541637/1‑68 (MQ=255) tGAAACAGCCCGCGCAGCCGCAGCTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTTCCg                                                                >  1:1535451/1‑68 (MQ=255) tGAAACAGCCCGCGCAGCCGCAGCTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTTCCg                                                                >  1:819668/1‑68 (MQ=255) tGAAACAGCCCGCGCAGCCGCAGCTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTTCCg                                                                >  1:190826/1‑68 (MQ=255) tGAAACAGCCCGCGCAGCCGCAGCTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATACGTTCCg                                                                >  1:1309276/1‑68 (MQ=255)    agcagcCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAgg                                                             <  1:2344561/68‑1 (MQ=255)                    gCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAAcc                                            <  1:1293933/69‑1 (MQ=255)                                tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                               >  1:458084/1‑54 (MQ=255)                                tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                               >  1:2316588/1‑54 (MQ=255)                                tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                               >  1:44741/1‑54 (MQ=255)                                tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                               >  1:79009/1‑54 (MQ=255)                                tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                               >  1:1864551/1‑54 (MQ=255)                                tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                               >  1:496170/1‑54 (MQ=255)                                tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                               >  1:2334025/1‑54 (MQ=255)                                tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                               >  1:642007/1‑54 (MQ=255)                                tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                               >  1:2291837/1‑54 (MQ=255)                                tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                               >  1:3286216/1‑54 (MQ=255)                                tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                               >  1:2213709/1‑54 (MQ=255)                                tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                               >  1:1725645/1‑54 (MQ=255)                                tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                               >  1:2267128/1‑54 (MQ=255)                                tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                               >  1:2281040/1‑54 (MQ=255)                                tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                               >  1:737517/1‑54 (MQ=255)                                tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                               >  1:1055158/1‑54 (MQ=255)                                tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                               >  1:1201002/1‑54 (MQ=255)                                tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                               >  1:2854359/1‑54 (MQ=255)                                tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                               >  1:1215188/1‑54 (MQ=255)                                tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                               >  1:1346573/1‑54 (MQ=255)                                tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                               >  1:802877/1‑54 (MQ=255)                                tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                               >  1:1450826/1‑54 (MQ=255)                                tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                               >  1:2405437/1‑54 (MQ=255)                                tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                               >  1:656208/1‑54 (MQ=255)                                tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTa                                                 >  1:1005318/1‑52 (MQ=255)                                tttGCGATGACGTGCATCTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                               >  1:2736996/1‑54 (MQ=255)                                  tGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGtt                              >  1:2409423/1‑69 (MQ=255)                                  tGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGtt                              >  1:1452077/1‑69 (MQ=255)                                  tGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGtt                              >  1:1052542/1‑69 (MQ=255)                                      aTGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTAAAACCGCTGGCGACCAGTTCAc                           >  1:280798/1‑68 (MQ=255)                                      aTGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTAAAACCGCTGGCGACCAGTTCAc                           >  1:2713108/1‑68 (MQ=255)                                      aTGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTAAAACCGCTGGCGACCAGTTCAc                           >  1:2289606/1‑68 (MQ=255)                                      aTGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTAAAACCGCTGGCGACCAGTTCAc                           >  1:2234050/1‑68 (MQ=255)                                      aTGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTAAAACCGCTGGCGACCAGTTCAc                           >  1:2006064/1‑68 (MQ=255)                                      aTGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTAAAACCGCTGGCGACCAGTTCAc                           >  1:1809503/1‑68 (MQ=255)                                      aTGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTAAAACCGCTGGCGACCAGTTCAc                           >  1:530965/1‑68 (MQ=255)                                      aTGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTAAAACCGCTGGCGACCAGTTCAc                           >  1:656067/1‑68 (MQ=255)                                      aTGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTAAAACCGCTGGCGACCAGTTCAc                           >  1:685578/1‑68 (MQ=255)                                      aTGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTAAAACCGCTGGCGACCAGTTCAc                           >  1:869374/1‑68 (MQ=255)                                      aTGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTAAAACCGCTGGCGACCAGTTCAc                           >  1:905550/1‑68 (MQ=255)                                      aTGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATAAACGATTAAAACCGCTGGCGACCAGTTCAc                           >  1:1204339/1‑68 (MQ=255)                                       tGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGc                         <  1:3158617/69‑1 (MQ=255)                                       tGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGc                         <  1:79687/69‑1 (MQ=255)                                       tGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGc                         <  1:2478831/69‑1 (MQ=255)                                       tGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGc                         <  1:977858/69‑1 (MQ=255)                                        gACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCa                        <  1:411320/69‑1 (MQ=255)                                        gACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCa                        <  1:1037082/69‑1 (MQ=255)                                        gACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCa                        <  1:1171346/69‑1 (MQ=255)                                        gACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCa                        <  1:1207432/69‑1 (MQ=255)                                        gACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCa                        <  1:2411646/69‑1 (MQ=255)                                        gACGTGCATGTGGCGCAGAATCAGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCa                        <  1:2769053/69‑1 (MQ=255)                                              cATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCt                  >  1:1218336/1‑69 (MQ=255)                                              cATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCt                  >  1:549426/1‑69 (MQ=255)                                                  tgGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTg              >  1:255133/1‑69 (MQ=255)                                                             ccGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTAc   <  1:422348/69‑1 (MQ=255)                                                             ccGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTAc   <  1:1271121/69‑1 (MQ=255)                                                             ccGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTAc   <  1:159870/69‑1 (MQ=255)                                                             ccGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTAc   <  1:1714152/69‑1 (MQ=255)                                                             ccGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTAc   <  1:1731607/69‑1 (MQ=255)                                                             ccGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTAc   <  1:1763099/69‑1 (MQ=255)                                                             ccGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTAc   <  1:1858038/69‑1 (MQ=255)                                                             ccGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTAc   <  1:42050/69‑1 (MQ=255)                                                             ccGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTAc   <  1:2147780/69‑1 (MQ=255)                                                             ccGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTAc   <  1:25071/69‑1 (MQ=255)                                                             ccGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTAc   <  1:290564/69‑1 (MQ=255)                                                             ccGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTAc   <  1:322583/69‑1 (MQ=255)                                                             ccGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTAc   <  1:2354943/69‑1 (MQ=255)                                                             ccGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTAc   <  1:314666/69‑1 (MQ=255)                                                             ccGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATAAc   <  1:16034/69‑1 (MQ=255)                                                              cGTTCCGAGGAATACACGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGGATATTGCg  <  1:1377869/69‑1 (MQ=255)                                                              cGTTCCGAGGAATACACGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGGATATTGCg  <  1:1486053/69‑1 (MQ=255)                                                              cGTTCCGAGGAATACACGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGGATATTGCg  <  1:632589/69‑1 (MQ=255)                                                              cGTTCCGAGGAATACACGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGGATATTGCg  <  1:308720/69‑1 (MQ=255)                                                              cGTACCGAGGAATACTCGATTAAAAGCGCTGGCGACCAGTTCACGc                         >  1:2068505/1‑46 (MQ=255)                                                              cGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGc                         >  1:2575084/1‑46 (MQ=255)                                                              cGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGc                         >  1:38493/1‑46 (MQ=255)                                                              cGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGc                         >  1:2715603/1‑46 (MQ=255)                                                              cGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGc                         >  1:1698515/1‑46 (MQ=255)                                                              cGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGc                         >  1:2175936/1‑46 (MQ=255)                                                              cGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGc                         >  1:338086/1‑46 (MQ=255)                                                              cGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGc                         >  1:1994100/1‑46 (MQ=255)                                                              cGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGc                         >  1:2675238/1‑46 (MQ=255)                                                              cGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGc                         >  1:349852/1‑46 (MQ=255)                                                              cGTACCGAGGAATACTCGATTAAAACAGCTGGCGACCAGTTCACGc                         >  1:3422339/1‑46 (MQ=255)                                                              cGTAACGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGc                         >  1:1694852/1‑46 (MQ=255)                                                                 |                                                                  TGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACG  >  minE/712625‑712754

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