Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	712,701	T→A	46.6%	R245R (CGA→CGT)	pyrC ←	dihydro‑orotase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	712,701	0	T	A	46.6%	20.3 / 98.3	90	R245R (CGA→CGT)	pyrC	dihydro‑orotase
Reads supporting (aligned to +/- strand):  ref base T (19/29);  new base A (38/4);  total (57/33)
Fisher's exact test for biased strand distribution p-value = 6.36e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGATGGGTA  >  minE/712644‑712768                                                          |                                                                    gCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAAcc                                                          <  1:1293933/69‑1 (MQ=255)             tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                                             >  1:1346573/1‑54 (MQ=255)             tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                                             >  1:44741/1‑54 (MQ=255)             tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                                             >  1:458084/1‑54 (MQ=255)             tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                                             >  1:1725645/1‑54 (MQ=255)             tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                                             >  1:496170/1‑54 (MQ=255)             tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                                             >  1:3286216/1‑54 (MQ=255)             tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                                             >  1:642007/1‑54 (MQ=255)             tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                                             >  1:1864551/1‑54 (MQ=255)             tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                                             >  1:656208/1‑54 (MQ=255)             tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                                             >  1:2854359/1‑54 (MQ=255)             tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                                             >  1:2213709/1‑54 (MQ=255)             tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                                             >  1:1450826/1‑54 (MQ=255)             tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                                             >  1:737517/1‑54 (MQ=255)             tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                                             >  1:79009/1‑54 (MQ=255)             tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                                             >  1:2267128/1‑54 (MQ=255)             tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                                             >  1:802877/1‑54 (MQ=255)             tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                                             >  1:1215188/1‑54 (MQ=255)             tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                                             >  1:2281040/1‑54 (MQ=255)             tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                                             >  1:2291837/1‑54 (MQ=255)             tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                                             >  1:1201002/1‑54 (MQ=255)             tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                                             >  1:2316588/1‑54 (MQ=255)             tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                                             >  1:1055158/1‑54 (MQ=255)             tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                                             >  1:2334025/1‑54 (MQ=255)             tttGCGATGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTaaa                                                             >  1:2405437/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:1052542/1‑69 (MQ=255)               tGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGtt                                            >  1:1452077/1‑69 (MQ=255)                   aTGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTAAAACCGCTGGCGACCAGTTCAc                                         >  1:530965/1‑68 (MQ=255)                   aTGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTAAAACCGCTGGCGACCAGTTCAc                                         >  1:905550/1‑68 (MQ=255)                   aTGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTAAAACCGCTGGCGACCAGTTCAc                                         >  1:2289606/1‑68 (MQ=255)                   aTGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTAAAACCGCTGGCGACCAGTTCAc                                         >  1:869374/1‑68 (MQ=255)                   aTGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTAAAACCGCTGGCGACCAGTTCAc                                         >  1:2713108/1‑68 (MQ=255)                   aTGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTAAAACCGCTGGCGACCAGTTCAc                                         >  1:2234050/1‑68 (MQ=255)                   aTGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTAAAACCGCTGGCGACCAGTTCAc                                         >  1:280798/1‑68 (MQ=255)                   aTGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTAAAACCGCTGGCGACCAGTTCAc                                         >  1:685578/1‑68 (MQ=255)                   aTGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTAAAACCGCTGGCGACCAGTTCAc                                         >  1:2006064/1‑68 (MQ=255)                   aTGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTAAAACCGCTGGCGACCAGTTCAc                                         >  1:656067/1‑68 (MQ=255)                   aTGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATACACGATTAAAACCGCTGGCGACCAGTTCAc                                         >  1:1809503/1‑68 (MQ=255)                   aTGACGTGCATGTGGCGCAGAATCCGTTCCGAGGAATAAACGATTAAAACCGCTGGCGACCAGTTCAc                                         >  1:1204339/1‑68 (MQ=255)                    tGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGc                                       <  1:79687/69‑1 (MQ=255)                    tGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGc                                       <  1:3158617/69‑1 (MQ=255)                    tGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGc                                       <  1:2478831/69‑1 (MQ=255)                    tGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGc                                       <  1:977858/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)                     gACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCa                                      <  1:411320/69‑1 (MQ=255)                     gACGTGCATGTGGCGCAGAATCAGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCa                                      <  1:2769053/69‑1 (MQ=255)                           cATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCt                                >  1:549426/1‑69 (MQ=255)                           cATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCt                                >  1:1218336/1‑69 (MQ=255)                               tgGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTg                            >  1:255133/1‑69 (MQ=255)                                          ccGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTAc                 <  1:322583/69‑1 (MQ=255)                                          ccGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTAc                 <  1:1714152/69‑1 (MQ=255)                                          ccGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTAc                 <  1:1858038/69‑1 (MQ=255)                                          ccGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTAc                 <  1:159870/69‑1 (MQ=255)                                          ccGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTAc                 <  1:42050/69‑1 (MQ=255)                                          ccGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTAc                 <  1:422348/69‑1 (MQ=255)                                          ccGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTAc                 <  1:1763099/69‑1 (MQ=255)                                          ccGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTAc                 <  1:1731607/69‑1 (MQ=255)                                          ccGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTAc                 <  1:314666/69‑1 (MQ=255)                                          ccGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTAc                 <  1:290564/69‑1 (MQ=255)                                          ccGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTAc                 <  1:2147780/69‑1 (MQ=255)                                          ccGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTAc                 <  1:1271121/69‑1 (MQ=255)                                          ccGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTAc                 <  1:25071/69‑1 (MQ=255)                                          ccGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTAc                 <  1:2354943/69‑1 (MQ=255)                                          ccGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATAAc                 <  1:16034/69‑1 (MQ=255)                                           cGTTCCGAGGAATACACGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGGATATTGCg                <  1:1486053/69‑1 (MQ=255)                                           cGTTCCGAGGAATACACGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGGATATTGCg                <  1:308720/69‑1 (MQ=255)                                           cGTTCCGAGGAATACACGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGGATATTGCg                <  1:632589/69‑1 (MQ=255)                                           cGTTCCGAGGAATACACGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGGATATTGCg                <  1:1377869/69‑1 (MQ=255)                                           cGTACCGAGGAATACTCGATTAAAAGCGCTGGCGACCAGTTCACGc                                       >  1:2068505/1‑46 (MQ=255)                                           cGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGc                                       >  1:38493/1‑46 (MQ=255)                                           cGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGc                                       >  1:1698515/1‑46 (MQ=255)                                           cGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGc                                       >  1:349852/1‑46 (MQ=255)                                           cGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGc                                       >  1:338086/1‑46 (MQ=255)                                           cGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGc                                       >  1:1994100/1‑46 (MQ=255)                                           cGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGc                                       >  1:2175936/1‑46 (MQ=255)                                           cGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGc                                       >  1:2715603/1‑46 (MQ=255)                                           cGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGc                                       >  1:2675238/1‑46 (MQ=255)                                           cGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGc                                       >  1:2575084/1‑46 (MQ=255)                                           cGTACCGAGGAATACTCGATTAAAACAGCTGGCGACCAGTTCACGc                                       >  1:3422339/1‑46 (MQ=255)                                           cGTAACGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGc                                       >  1:1694852/1‑46 (MQ=255)                                                         cTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAg          <  1:462920/61‑1 (MQ=255)                                                         cTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAg          <  1:1244313/61‑1 (MQ=255)                                                         cTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAg          <  1:2267446/61‑1 (MQ=255)                                                         cTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGATGGGTa  >  1:2446053/1‑69 (MQ=255)                                                          |                                                                    GCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGATGGGTA  >  minE/712644‑712768

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