Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,387,943	C→A	32.1%	L157L (CTG→CTT)	yigL ←	predicted hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,387,943	0	C	A	32.1%	102.0 / 52.8	84	L157L (CTG→CTT)	yigL	predicted hydrolase
Reads supporting (aligned to +/- strand):  ref base C (14/43);  new base A (9/18);  total (23/61)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAGGTAGAGAAACTGACGTTGACGCGATCGCCCCAACGAGCGTTAATCGCCTGCTCCAGCGGCAGCAGTTGTTCATGGGAATCGCAGGTGAAGAACACTTTGCTGACGCCTTCCGGCTCCAGTAATC  >  minE/2387877‑2388004                                                                   |                                                              tAAGGTAGAGAAACTGACGTTGACGCGATCGCCCCAACGAGCGTTAATCGCCTGCTCCAGCGGcagcag                                                             <  1:332329/69‑1 (MQ=255)      tAGAGAAACTGACGTTGACGCGATCGCCCCAACGAGCGTTAATCGCCTGCTCCAGCGGCAGCAgttgtt                                                        <  1:1315178/69‑1 (MQ=255)      tAGAGAAACTGACGTTGACGCGATCGCCCCAACGAGCGTTAATCGCCTGCTCCAGCGGCAGCAgttgtt                                                        <  1:59658/69‑1 (MQ=255)      tAGAGAAACTGACGTTGACGCGAGCGCCCCAACGAGCGTTAATCGCCTGCTCCAGCGGCAGCAgttgtt                                                        <  1:753500/69‑1 (MQ=255)       agagAAACTGACGTTGACGCGATCGCCCCAACGAGCGTTAATCGCCTGCTCCAGCGGCAGCAgttgtt                                                        <  1:962831/68‑1 (MQ=255)       agagAAACTGACGTTGACGCGATCGCCCCAACGAGCGTTAATCGCCTGCTCCAGCGGCAGCAGTTGTTc                                                       <  1:3096356/69‑1 (MQ=255)       agagAAACTGACGTTGACGCGATCGCCCCAACGAGCGTTAATCGCCTGCTCCAGCGGCAGCAGTTGTTc                                                       <  1:1138515/69‑1 (MQ=255)       agagAAACTGACGTTGACGCGATCGCCCCAACGAGCGTTAATCGCCTGCTCCAGCGGCAGCAGTTGTTc                                                       <  1:2975480/69‑1 (MQ=255)       agagAAACTGACGTTGACGCGATCGCCCCAACGAGCGTTAATCGCCTGCTCCAGCGGCAGCAGTTGTTc                                                       <  1:746949/69‑1 (MQ=255)       agagAAACTGACGTTGACGCGATCGCCCCAACGAGCGTTAATCGCCTGCTCCAGCGGCAGCAGTTGTTc                                                       <  1:275539/69‑1 (MQ=255)       agagAAACTGACGTTGACGCGATCGCCCCAACGAGCGTTAATCGCCTGCTCCAGCGGCAGCAGTTGTTc                                                       <  1:1578255/69‑1 (MQ=255)       agagAAACTGACGTTGACGCGATCGCCCCAACGAGCGTTAATCGCCTGCTCCAGCGGCAGCAGTTGTTc                                                       <  1:1601450/69‑1 (MQ=255)       agagAAACTGACGTTGACGCGATCGCCCCAACGAGCGTTAATCGCCTGCTCCAGCGGCAGCAGTTGTTc                                                       <  1:2798718/69‑1 (MQ=255)                gACGTTGACGCGATCGCCCCAACGAGCGTTAATCGCCTGCTCCAGCGGCAGCAGTTGTTCATGGGAAt                                               >  1:735479/1‑68 (MQ=255)                gACGTTGACGCGATCGCCCCAACGAGCGTTAATCGCCTGCTCCAGCGGCAGCAGTTGTTCATGGGAAt                                               >  1:2336042/1‑68 (MQ=255)                gACGTTGACGCGATCGCCCCAACGAGCGTTAATCGCCTGCTCCAGCGGCAGCAGTTGTTCATGGGAAt                                               >  1:2481170/1‑68 (MQ=255)                gACGTTGACGCGATCGCCCCAACGAGCGTTAATCGCCTGCTCCAGCGGCAGCAGTTGTTCATGGGAAt                                               >  1:82751/1‑68 (MQ=255)                gACGTTGACGCGATCGCCCCAACGAGCGTTAATCGCCTGCTCCAGCGGCAGCAGTTGTTCATGGGAAt                                               >  1:282254/1‑68 (MQ=255)                 acattGACGCGATCGCCCCAGCGAGCGTTAATCGCCTGCTCCAGTGGCAGAAGTTTTTCATGGGa                                                 <  1:2812961/62‑1 (MQ=255)                 acattGACGCGATCGCCCCAGCGAGCGTTAATCGCCTGCTCCAGTGGCAGAAGTTTTTCATGGGa                                                 <  1:547468/62‑1 (MQ=255)                 acattGACGCGATCGCCCCAGCGAGCGTTAATCGCCTGCTCCAGTGGCAGAAGTTTTTCATGGGa                                                 <  1:2593334/62‑1 (MQ=255)                 acattGACGCGATCGCCCCAGCGAGCGTTAATCGCCTGCTCCAGTGGCAGAAGTTTTTCATGGGa                                                 <  1:2574716/62‑1 (MQ=255)                 acattGACGCGATCGCCCCAGCGAGCGTTAATCGCCTGCTCCAGTGGCAGAAGTTTTTCATGGGa                                                 <  1:2890386/62‑1 (MQ=255)                 acattGACGCGATCGCCCCAGCGAGCGTTAATCGCCTGCTCCAGTGGCAGAAGTTTTTCATGGGa                                                 <  1:2896655/62‑1 (MQ=255)                 acattGACGCGATCGCCCCAGCGAGCGTTAATCGCCTGCTCCAGTGGCAGAAGTTTTTCATGGGa                                                 <  1:2921437/62‑1 (MQ=255)                 acattGACGCGATCGCCCCAGCGAGCGTTAATCGCCTGCTCCAGTGGCAGAAGTTTTTCATGGGa                                                 <  1:62592/62‑1 (MQ=255)                 acattGACGCGATCGCCCCAGCGAGCGTTAATCGCCTGCTCCAGTGGCAGAAGTTTTTCATGGGa                                                 <  1:6980/62‑1 (MQ=255)                 acattGACGCGATCGCCCCAGCGAGCGTTAATCGCCTGCTCCAGTGGCAGAAGTTTTTCATGGGa                                                 <  1:3088080/62‑1 (MQ=255)                 acattGACGCGATCGCCCCAGCGAGCGTTAATCGCCTGCTCCAGTGGCAGAAGTTTTTCATGGGa                                                 <  1:3254148/62‑1 (MQ=255)                 acattGACGCGATCGCCCCAGCGAGCGTTAATCGCCTGCTCCAGTGGCAGAAGTTTTTCATGGGa                                                 <  1:816452/62‑1 (MQ=255)                 acattGACGCGATCGCCCCAGCGAGCGTTAATCGCCTGCTCCAGTGGCAGAAGTTTTTCATGGGa                                                 <  1:1283673/62‑1 (MQ=255)                 acattGACGCGATCGCCCCAGCGAGCGTTAATCGCCTGCTCCAGTGGCAGAAGTTTTTCATGGGa                                                 <  1:1480747/62‑1 (MQ=255)                 acattGACGCGATCGCCCCAGCGAGCGTTAATCGCCTGCTCCAGTGGCAGAAGTTTTTCATGGGa                                                 <  1:927880/62‑1 (MQ=255)                 acattGACGCGATCGCCCCAGCGAGCGTTAATCGCCTGCTCCAGTGGCAGAAGTTTTTCATGGGa                                                 <  1:2378601/62‑1 (MQ=255)                 acattGACGCGATCGCCCCAGCGAGCGTTAATCGCCTGCTCCAGTGGCAGAAGTTTTTCATGGGa                                                 <  1:1732273/62‑1 (MQ=255)                 acattGACGCGATCGCCCCAGCGAGCGTTAATCGCCTGCTCCAGGGGCAGAAGTTTTTCATGGGa                                                 <  1:75239/62‑1 (MQ=255)                   gTTGACGCGATCGCCCCAACGAGCGTTAATCGCCTGCTCCAGCGTCAGCAGTTGTTCATGGGAATCGc                                            <  1:924151/68‑1 (MQ=255)                   gTTGACGCGATCGCCCCAACGAGCGTTAATCGCCTGCTCCAGCGGCAGCAGTTGTTCATGGGAATCGc                                            <  1:607117/68‑1 (MQ=255)                   gTTGACGCGATCGCCCCAACGAGCGTTAATCGCCTGCTCCAGCGGCAGCAGTTGTTCATGGGAATCGc                                            <  1:760810/68‑1 (MQ=255)                   gTTGACGCGATCGCCCCAACGAGCGTTAATCGCCTGCTCCAGCGGCAGCAGTTGTTCATGGGAATCGc                                            <  1:772042/68‑1 (MQ=255)                   gTTGACGCGATCGCCCCAACGAGCGTTAATCGCCTGCTCCAGCGGCAGCAGTTGTTCATGGGAATCGc                                            <  1:2077524/68‑1 (MQ=255)                   gTTGACGCGATCGCCCCAACGAGCGTTAATCGCCTGCTCCAGCGGCAGCAGTTGTTCATGGGAATCGc                                            <  1:1124255/68‑1 (MQ=255)                   gTTGACGCGATCGCCCCAACGAGCGTTAATCGCCTGCTCCAGCGGCAGCAGTTGTTCATGGGAATCGc                                            <  1:873174/68‑1 (MQ=255)                   gTTGACGCGATCGCCCCAACGAGCGTTAATCGCCTGCTCCAGCGGCAGCAGTTGTTCATGGGAATCGc                                            <  1:2138185/68‑1 (MQ=255)                   gTTGACGCGATCGCCCCAACGAGCGTTAATCGCCTGCTCCAGCGGCAGCAGTTGTTCATGGGAATCGc                                            <  1:3286555/68‑1 (MQ=255)                   gTTGACGCGATCGCCCCAACGAGCGTTAATCGCCTGCTCCAGCGGCAGCAGTTGTTCATGGGAATCGc                                            <  1:2318883/68‑1 (MQ=255)                   gTTGACGCGATCGCCCCAACGAGCGTTAATCGCCTGCTCCAG‑GGCAGCAGTTGTTCATGGCAATCGc                                            <  1:440306/67‑1 (MQ=255)                      gACGCGATCGCCCCAACGAGCGTTAATCGCGTGCTCCAGCGGCAGCAGTTGTTCATGGGAATCGCAgg                                         >  1:1401298/1‑68 (MQ=255)                      gACGCGATCGCCCCAACGAGCGTTAATCGCCTGCTCCAGCGGCAGCAGTTGTTCATGGGAATCGCAgg                                         >  1:952589/1‑68 (MQ=255)                      gACGCGATCGCCCCAACGAGCGTTAATCGCCTGCTCCAGCGGCAGCAGTTGTTCATGGGAATCGCAgg                                         >  1:3315420/1‑68 (MQ=255)                      gACGCGATCGCCCCAACGAGCGTTAATCGCCTGCTCCAGCGGCAGCAGTTGTTCATGGGAATCGCAgg                                         >  1:3155770/1‑68 (MQ=255)                      gACGCGATCGCCCCAACGAGCGTTAATCGCCTGCTCCAGCGGCAGCAGTTGTTCATGGGAATCGCAgg                                         >  1:2269768/1‑68 (MQ=255)                                    aaCGAGCGTTAATCGCCTGCTCCAGCGGCAGCAGTTGTTc                                                       <  1:3208885/40‑1 (MQ=255)                                     aCGAGCGTTAATCGCCTGCTCCAGCGGCAGCAGTTGTTCATGGGAATCGCAGGTGAAGAACACTTTGCt                         <  1:851286/69‑1 (MQ=255)                                     aCGAGCGTTAATCGCCTGCTCCAGCGGCAGCAGTTGTTCATGGGAATCGCAGGTGAAGAACACTTTGCt                         <  1:984880/69‑1 (MQ=255)                                     aCGAGCGTTAATCGCCTGCTCCAGCGGCAGCAGTTGTTCATGGGAATCGCAGGTGAAGAACACTTTGCt                         <  1:2260459/69‑1 (MQ=255)                                     aCGAGCGTTAATCGCCTGCTCCAGCGGCAGCAGTTGTTCATGGGAATCGCAGGTGAAGAACACTTTGCt                         <  1:2654812/69‑1 (MQ=255)                                     aCGAGCGTTAATCGCCTGCTCCAGCGGCAGCAGTTGTGCATGGGAATCGCAGGTGAAGAACACTTTGCt                         <  1:1101139/69‑1 (MQ=255)                                          cGTTAATCGCCTGCTCCAGCGGCAGCAGTTGTTCATGGGAATCGCAGGTGAAGAACACTTTGCTGACGc                    <  1:2250307/69‑1 (MQ=255)                                          cGTTAATCGCCTGCTCCAGCGGCAGCAGTTGTTCATGGGAATCGCAGGTGAAGAACACTTTGCTGACGc                    <  1:2359383/69‑1 (MQ=255)                                            ttAATCGCCTGCTCCAGCGGCAGCAGTTGTTCATGGGAATCGCAGGTGAAGAACACTTTGCTGACGc                    >  1:2883181/1‑67 (MQ=255)                                            ttAATCGCCTGCTCCAGCGGCAGCAGTTGTTCATGGGAATCGCAGGTGAAGAACACTTTGCTGACGCCt                  <  1:1250218/69‑1 (MQ=255)                                            ttAATCGCCTGCTCCAGCGGCAGCAGTTGTTCATGGGAATCGCAGGTGAAGAACACTTTGCTGACGCCt                  <  1:1515732/69‑1 (MQ=255)                                            ttAATCGCCTGCTCCAGCGGCAGCAGTTGTTCATGGGAATCGCAGGTGAAGAACACTTTGCTGACGCCt                  <  1:3318112/69‑1 (MQ=255)                                                 cGCCTGCTCCAGCGGCAGCAGTTGTTCATGGGAATCGCAGGTGAAGAACACTTTGCTGACGCCTTCCgg             >  1:2820948/1‑69 (MQ=255)                                                 cGCCTGCTCCAGCGGCAGCAGTTGTTCATGGGAATCGCAGGTGAAGAACACTTTGCTGACGCCTTCCgg             >  1:801598/1‑69 (MQ=255)                                                 cGCCTGCTCCAGCGGCAGCAGTTGTTCATGGGAATCGCAGGTGAAGAACACTTTGCTGACGCCTTCCgg             >  1:1890080/1‑69 (MQ=255)                                                    cTGCTCCATTGGCAGAAGTTTTTCATGGGAATCGCAGGTGAAGAACACTTTGCTGACGCCTTCCGGc            >  1:3366381/1‑67 (MQ=255)                                                    cTGCTCCATTGGCAGAAGTTTTTCATGGGAATCGCAGGTGAAGAACACTTTGCTGACGCCTTCCGGc            >  1:28854/1‑67 (MQ=255)                                                    cTGCTCCATTGGCAGAAGTTTTTCATGGGAATCGCAGGTGAAGAACACTTTGCTGACGCCTTCCGGc            >  1:1620164/1‑67 (MQ=255)                                                    cTGCTCCAGTGGCAGAAGTTTTTCATGGGAATCGCAGGTGAAGAACACTTTGCTGACGCCTTCCGGc            >  1:1957316/1‑67 (MQ=255)                                                    cTGCTCCAGTGGCAGAAGTTTTTCATGGGAATCGCAGGTGAAGAACACTTTGCTGACGCCTTCCGGc            >  1:969982/1‑67 (MQ=255)                                                    cTGCTCCAGTGGCAGAAGTTTTTCATGGGAATCGCAGGTGAAGAACACTTTGCTGACGCCTTCCGGc            >  1:1567358/1‑67 (MQ=255)                                                    cTGCTCCAGTGGCAGAAGTTTTTCATGGGAATCGCAGGTGAAGAACACTTTGCTGACGCCTTCCGGc            >  1:1703983/1‑67 (MQ=255)                                                    cTGCTCCAGTGGCAGAAGTTTTTCATGGGAATCGCAGGTGAAGAACACTTTGCTGACGCCTTCCGGc            >  1:3225325/1‑67 (MQ=255)                                                    cTGCTCCAGTGGCAGAAGTTTTTCATGGGAATCGCAGGTGAAGAACACTTTGCTGACGCCTTCCGGc            >  1:577209/1‑67 (MQ=255)                                                       cTCCAGCGGCAGCAGTTGTTCATGGGAATCGCAGGTGAAGAACACTTTGCt                         <  1:2867703/51‑1 (MQ=255)                                                         gcAGCGGCAGCAGTTGTTCATGGGAATCGCAGGTGAAGAACACTTTGCt                         <  1:1526044/48‑1 (MQ=255)                                                         ccAGCGGCAGCAGTTGTTCATGGGAATCGCAGGTGAAGAACACTTTGCt                         <  1:1784852/49‑1 (MQ=255)                                                         ccAGCGGCAGCAGTTGTTCATGGGAATCGCAGGTGAAGAACACTTTGCt                         <  1:3122685/49‑1 (MQ=255)                                                         ccAGCGGCAGCAGTTGTTCATGGGAATAGCAGGTGAAGAACACTTTGCt                         <  1:1425145/49‑1 (MQ=255)                                                            gCGGCAGCAGTTGTTCATGGGAATCGCAGGTGAAGAACACTTTGCTGACGCCTTCCGGCTCCAGTAATc  <  1:375515/69‑1 (MQ=255)                                                            gCGGCAGCAGTTGTTCATGGGAATCGCAGGTGAAGAACACTTTGCTGACGCCTTCCGGCTCCAGTAATc  <  1:2145318/69‑1 (MQ=255)                                                            gCGGCAGCAGTTGTTCATGGGAATCGCAGGTGAAGAACACTTTGCTGACGCCTTCCGGCTCCAGTAATc  <  1:2394000/69‑1 (MQ=255)                                                            gCGGCAGCAGTTGTTCATGGGAATCGCAGGTGAAGAACACTTTGCTGACGCCTTCCGGCTCCAGTAATc  <  1:2447697/69‑1 (MQ=255)                                                                   |                                                              TAAGGTAGAGAAACTGACGTTGACGCGATCGCCCCAACGAGCGTTAATCGCCTGCTCCAGCGGCAGCAGTTGTTCATGGGAATCGCAGGTGAAGAACACTTTGCTGACGCCTTCCGGCTCCAGTAATC  >  minE/2387877‑2388004

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