Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	190,756	G→A	49.5%	V19M (GTG→ATG)	yaeR →	predicted lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	190,756	0	G	A	49.5%	‑1.3 / 87.1	95	V19M (GTG→ATG)	yaeR	predicted lyase
Reads supporting (aligned to +/- strand):  ref base G (25/23);  new base A (36/11);  total (61/34)
Fisher's exact test for biased strand distribution p-value = 1.83e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGGAGACGTTATGCTGGGTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGGTGAGCAAAGCTTTCTACTGCGATATTCTTGGTTTCACGCTGCAAAGCG  >  minE/190692‑190804                                                                 |                                                 aGGAGACGTTATGCTGGGTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATg                                                >  1:2027876/1‑67 (MQ=255) aGGAGACGTTATGCTGGGTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGa                                               >  1:1046478/1‑68 (MQ=255) aGGAGACGTTATGCTGGGTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGa                                               >  1:984905/1‑68 (MQ=255) aGGAGACGTTATGCTGGGTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGa                                               >  1:916265/1‑68 (MQ=255) aGGAGACGTTATGCTGGGTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGa                                               >  1:596060/1‑68 (MQ=255) aGGAGACGTTATGCTGGGTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGa                                               >  1:37887/1‑68 (MQ=255) aGGAGACGTTATGCTGGGTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGa                                               >  1:2889294/1‑68 (MQ=255) aGGAGACGTTATGCTGGGTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGa                                               >  1:2806527/1‑68 (MQ=255) aGGAGACGTTATGCTGGGTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGa                                               >  1:2513533/1‑68 (MQ=255) aGGAGACGTTATGCTGGGTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGa                                               >  1:2305262/1‑68 (MQ=255) aGGAGACGTTATGCTGGGTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGa                                               >  1:2168767/1‑68 (MQ=255) aGGAGACGTTATGCTGGGTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGa                                               >  1:2045564/1‑68 (MQ=255) aGGAGACGTTATGCTGGGTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGa                                               >  1:1979035/1‑68 (MQ=255) aGGAGACGTTATGCTGGGTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGa                                               >  1:175196/1‑68 (MQ=255) aGGAGACGTTATGCTGGGTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGa                                               >  1:1698578/1‑68 (MQ=255) aGGAGACGTTATGCTGGGTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGa                                               >  1:1592675/1‑68 (MQ=255) aGGAGACGTTATGCTGGGTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGa                                               >  1:1268473/1‑68 (MQ=255) aGGAGACGTTATGCTGGGTTTAAAACAGGTTCACAATATTGCGATTATTGCGACGGATTATGCGATGa                                               >  1:2932554/1‑68 (MQ=255)        gTTATGCTGGGTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGGTGAGCAAAGCt                                       <  1:2347649/69‑1 (MQ=255)        gTTATGCTGGGTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGGTGAGCAAAGCt                                       <  1:2825278/69‑1 (MQ=255)           aTGCTGTGTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGGTGAGCAAAGCTTTc                                    <  1:950008/69‑1 (MQ=255)           aTGCTGGGTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGGTGAGCAAAGCTTTc                                    <  1:3403461/69‑1 (MQ=255)           aTGCTGGGTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGGTGAGCAAAGCTTTc                                    <  1:3213014/69‑1 (MQ=255)           aTGCTGGGTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGGTGAGCAAAGCTTTc                                    <  1:1198001/69‑1 (MQ=255)           aTGCTGGGTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGGTGAGCAAAGCTTTc                                    <  1:2967144/69‑1 (MQ=255)           aTGCTGGGTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGGTGAGCAAAGCTTTc                                    <  1:2068731/69‑1 (MQ=255)           aTGCTGGGTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGGTGAGCAAAGCTTTc                                    <  1:1616831/69‑1 (MQ=255)           aTGCTGGCTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGGTGAGCAAAGCTTTc                                    <  1:288697/69‑1 (MQ=255)              cTGGGTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGAGCAAAGCTTTCTAc                                 <  1:2261527/69‑1 (MQ=255)              cTGGGTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGAGCAAAGCTTTCTAc                                 <  1:1023247/69‑1 (MQ=255)              cTGGGTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGAGCAAAGCTTTCTAc                                 <  1:2253012/69‑1 (MQ=255)              cTGGGTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGAGCAAAGCTTTCTAc                                 <  1:2843705/69‑1 (MQ=255)              cTGGGTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGAGCAAAGCTTTCTAc                                 <  1:3131182/69‑1 (MQ=255)              cTGGGTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGAGCAAAGCTTTCTAc                                 <  1:587734/69‑1 (MQ=255)              cTGGGTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGAGCAAAGCTTTCTAc                                 <  1:838762/69‑1 (MQ=255)              cTGGGTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGAGCAAAGCTTTCTAc                                 <  1:859276/69‑1 (MQ=255)              cTGGGTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGAGCAAAGCTTTCTAc                                 <  1:1798451/69‑1 (MQ=255)               tGGGTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGAGCAAAGCTTTCTAc                                 <  1:2151992/68‑1 (MQ=255)                 ggTTTAAAACAGGTTCACCATATT‑CGATTATTGCGACGGATTATGCGGTGAGCAAAGCttt                                     >  1:1366101/1‑61 (MQ=255)                   tttAAAACAGGTTCACCATATTGCGATTATTGCGACTGATTATGCGGTGAGCAAAGCTTTCTACTGCGa                            <  1:1824096/69‑1 (MQ=255)                   tttAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGGTGAGCAAAGCTTTCTACTGCGa                            <  1:1771091/69‑1 (MQ=255)                   tttAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGGTGAGCAAAGCTTTCTACTGCGa                            <  1:1898874/69‑1 (MQ=255)                   tttAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGGTGAGCAAAGCTTTCTACTGCGa                            <  1:3089759/69‑1 (MQ=255)                   tttAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGGTGAGCAAAGCTTTCTACTGCGa                            <  1:1469345/69‑1 (MQ=255)                    ttAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGGTTAGCAAAGCTTTc                                    <  1:2944782/60‑1 (MQ=255)                    ttAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGAGCa                                            >  1:2033118/1‑52 (MQ=255)                    ttAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGAGCAAAGCTTTCTACTGCGat                           >  1:16087/1‑69 (MQ=255)                    ttAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGAGCAAAGCTTTCTACTGCGat                           >  1:1500661/1‑69 (MQ=255)                    ttAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGAGCAAAGCTTTCTACTGCGat                           >  1:2193531/1‑69 (MQ=255)                    ttAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGAGCAAAGCTTTCTACTGCGat                           >  1:2894197/1‑69 (MQ=255)                    ttAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGAGCAAAGCTTTCTACTGCGat                           >  1:2269001/1‑69 (MQ=255)                    ttAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGAGCAAAGCTTTCTACTGCGat                           >  1:311567/1‑69 (MQ=255)                    ttAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGAGCAAAGCTTTCTACTGCGat                           >  1:2163530/1‑69 (MQ=255)                    ttAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGAGCAAAGCTTTCTACTGCGat                           >  1:3256433/1‑69 (MQ=255)                    ttAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGAGCAAAGCTTTCTACTGCGat                           >  1:177524/1‑69 (MQ=255)                    ttAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGAGCAAAGCTTTCTACTGCGat                           >  1:431589/1‑69 (MQ=255)                    ttAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGAGCAAAGCTTTCTACTGCGat                           >  1:2658567/1‑69 (MQ=255)                    ttAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGAGCAAAGCTTTCTACTGCGat                           >  1:644027/1‑69 (MQ=255)                    ttAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGAGCAAAGCTTTCTACTGCGat                           >  1:2558775/1‑69 (MQ=255)                    ttAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGAGCAAAGCTTTCTACTGCGat                           >  1:2538195/1‑69 (MQ=255)                    ttAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGAGCAAAGCTTTCTACTGCGa                            >  1:1333258/1‑68 (MQ=255)                    ttAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGAGCAAAGCTTTCTACTGCGa                            >  1:454898/1‑68 (MQ=255)                    ttAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGAGCAAAGCTTTCAACTGCGat                           >  1:1926062/1‑69 (MQ=255)                     tAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGATGAGCAAAGCTTTCTAc                                 <  1:2421369/62‑1 (MQ=255)                           aGGTTCACCATATTGCGATTATTGCGACGGATTATGCGGTGAGCAAAGCTTTCTACTGCGATATTCTTg                    <  1:1590504/69‑1 (MQ=255)                           aGGTTCACCATATTGCGATTATTGCGACGGATTATGCGGTGAGCAAAGCTTTCTACTGCGATATTCTTg                    <  1:944596/69‑1 (MQ=255)                               tCACCATATTGCGATTATTGCGACGGATTATGCGGTGAGCAAAGCTTTCTACTGCGATATTCTTGGttt                <  1:746426/69‑1 (MQ=255)                               tCACCATATTGCGATTATTGCGACGGATTATGCGGTGAGCAAAGCTTTCTACTGCGATATTCTTGGttt                <  1:2805597/69‑1 (MQ=255)                                   cATATTGCGATTATTGCGACGGATTATGCGGTGAGCAAAGCTTTCTACTGCGATATTCTTGGTTTCAc             <  1:1608280/68‑1 (MQ=255)                                   cATATTGCGATTATTGCGACGGATTATGCGGTGAGCAAAGCTTTCTACTGCGATATTCTTGGTTTCAc             <  1:3356897/68‑1 (MQ=255)                                       ttGCGATTATTGCGACGGATTATGCGGTGAGCAAAGCTTTCTACTGCGATATTCTTGGTTTCACGCTg         >  1:373024/1‑68 (MQ=255)                                       ttGCGATTATTGCGACGGATTATGCGGTGAGCAAAGCTTTCTACTGCGATATTCTTGGTTTCACGCTg         >  1:2612257/1‑68 (MQ=255)                                       ttGCGATTATTGCGACGGATTATGCGGTGAGCAAAGCTTTCTACTGCGATATTCTTGGTTTCACGCTg         >  1:2729456/1‑68 (MQ=255)                                       ttGCGATTATTGCGACGGATTATGCGGTGAGCAAAGCTTTCTACTGCGATATTCTTGGTTTCACGCTg         >  1:2754613/1‑68 (MQ=255)                                       ttGCGATTATTGCGACGGATTATGCGGTGAGCAAAGCTTTCTACTGCGATATTCTTGGTTTCACGCTg         >  1:2965683/1‑68 (MQ=255)                                         gCGATTATTGCGACGGATTATGCGGTGAGCAAAGCTTTCTACTGCGATATTCTTGGTTTCACGCTGCaa      >  1:2505103/1‑69 (MQ=255)                                         gCGATTATTGCGACGGATTATGCGGTGAGCAAAGCTTTCTACTGCGATATTCTTGGTTTCACGCTGCaa      >  1:3227307/1‑69 (MQ=255)                                         gCGATTATTGCGACGGATTATGCGGTGAGCAAAGCTTTCTACTGCGATATTCTTGGTTTCACGCTGCaa      >  1:41627/1‑69 (MQ=255)                                         gCGATTATTGCGACGGATTATGCGGTGAGCAAAGCTTTCTACTGCGATATTCTTGGTTTCACGCTGCaa      >  1:2104076/1‑69 (MQ=255)                                         gCGATTATTGCGACGGATTATGCGGTGAGCAAAGCTTTCTACTGCGATATTCTTGGTTTCACGCTGCaa      <  1:573370/69‑1 (MQ=255)                                         gCGATTATTGCGACGGATTATGCGGTGAGCAAAGCTTTCTACTGCGATATTCTTGGTTTCACGCTGCaa      >  1:2723587/1‑69 (MQ=255)                                         gCGATTATTGCGACGGATTATGCGGTGAGCAAAGCTTTCTACTGCGATATTCTTGGTTTCACGCTGCaa      >  1:2376654/1‑69 (MQ=255)                                               attGCGACGGATTATGCGGTGAGCAAAGCTTTCTACTGCGATATTCTTGGTTTCACGCTGCAAAGCg  >  1:82512/1‑67 (MQ=255)                                               attGCGACGGATTATGCGGTGAGCAAAGCTTTCTACTGCGATATTCTTGGTTTCACGCTGCAAAGCg  >  1:754588/1‑67 (MQ=255)                                               attGCGACGGATTATGCGGTGAGCAAAGCTTTCTACTGCGATATTCTTGGTTTCACGCTGCAAAGCg  >  1:3006545/1‑67 (MQ=255)                                               attGCGACGGATTATGCGGTGAGCAAAGCTTTCTACTGCGATATTCTTGGTTTCACGCTGCAAAGCg  >  1:1864562/1‑67 (MQ=255)                                               attGCGACGGATTATGCGGTGAGCAAAGCTTTCTACTGCGATATTCTTGGTTTCACGCTGCAAAGCg  >  1:684879/1‑67 (MQ=255)                                               attGCGACGGATTATGCGGTGAGCAAAGCTTTCTACTGCGATATTCTTGGTTTCACGCTGCAAAGCg  >  1:501471/1‑67 (MQ=255)                                               attGCGACGGATTATGCGGTGAGCAAAGCTTTCTACTGCGATATTCTTGGTTTCACGCTGCAAAGCg  >  1:199640/1‑67 (MQ=255)                                               attGCGACGGATTATGCGGTGAGCAAAGCTTTCTACTGCGATATTCTTGGTTTCACGCTGCAAAGCg  >  1:3115167/1‑67 (MQ=255)                                               attGCGACGGATTATGCGGTGAGCAAAGCTTTCTACTGCGATATTCTTGGTTTCACGCTGCAAAGCg  >  1:2235826/1‑67 (MQ=255)                                               attGCGACGGATTATGCGGTGAGCAAAGCTTTCTACTGCGATATTCTTGGTTTCACGCTGCAAAGCg  >  1:1600434/1‑67 (MQ=255)                                               attGCGACGGATTATGCGGTGAGCAAAGCTTTCTACTGCGATATTCTTGGTTTCACGCTGCAAAGCg  >  1:2788942/1‑67 (MQ=255)                                               attGCGACGGATTATGCGGTGAGCAAAGCTTTCTACTGCGATATTCTTGGTTTCACGCTGCAAAGCg  >  1:2776557/1‑67 (MQ=255)                                               attGCGACGGATTATGCGGTGAGCAAAGCTTTCTACTGCGATATTCTTGGTTTCACGCTGCAAAGCg  >  1:2230768/1‑67 (MQ=255)                                                                 |                                                 AGGAGACGTTATGCTGGGTTTAAAACAGGTTCACCATATTGCGATTATTGCGACGGATTATGCGGTGAGCAAAGCTTTCTACTGCGATATTCTTGGTTTCACGCTGCAAAGCG  >  minE/190692‑190804

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