Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	240,243	C→T	11.7%	A236A (GCG→GCA)	araJ ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	240,243	0	C	T	11.7%	177.9 / 8.3	92	A236A (GCG→GCA)	araJ	predicted transporter
Reads supporting (aligned to +/- strand):  ref base C (47/34);  new base T (5/6);  total (52/40)
Fisher's exact test for biased strand distribution p-value = 5.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

TAGCATATTTCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTATGGCTTTACGTAGCTGAACCAGGCAAA  >  minE/240189‑240305                                                       |                                                               tAGCATATTTCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGaaaa                                                  >  1:2488473/1‑69 (MQ=255) tAGCATATTTCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGaaaa                                                  >  1:2458212/1‑69 (MQ=255) tAGCATATTTCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGaaaa                                                  >  1:1157579/1‑69 (MQ=255) tAGCATATTTCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGaaaa                                                  >  1:1006341/1‑69 (MQ=255) tAGCATATTTCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGaaaa                                                  >  1:698136/1‑69 (MQ=255) tAGCATATTTCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGaaaa                                                  >  1:921077/1‑69 (MQ=255) tAGCATATTTCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGaaaa                                                  >  1:2607924/1‑69 (MQ=255) tAGCATATTTCCCAGCACCATCCCTAGCACAACTAACATCATAATAAAGGTCATCGCCGTTTCCGaaaa                                                  >  1:617143/1‑69 (MQ=255)        tttCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGt                                                           >  1:2182530/1‑53 (MQ=255)        tttCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACGGGa                                            >  1:2080691/1‑68 (MQ=255)        tttCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCgg                                             >  1:2252229/1‑67 (MQ=255)        tttCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGa                                            >  1:2502342/1‑68 (MQ=255)        tttCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGa                                            >  1:2685907/1‑68 (MQ=255)        tttCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGa                                            >  1:2935576/1‑68 (MQ=255)        tttCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGa                                            >  1:1489269/1‑68 (MQ=255)        tttCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGa                                            >  1:530439/1‑68 (MQ=255)        tttCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGa                                            >  1:2619823/1‑68 (MQ=255)        tttCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGa                                            >  1:1146279/1‑68 (MQ=255)        tttCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGa                                            >  1:619077/1‑68 (MQ=255)        tttCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGa                                            >  1:628333/1‑68 (MQ=255)        ttccccAGCACCATCCCTAGCCCCACTAACATCATAATAAAGGTCATTGCCGTTTCCGaaaa                                                  >  1:651088/4‑62 (MQ=255)        ttccccAGCACCATCCCTAGCCCCACTAACATCATAATAAAGGTCATTGCCGTTTCCGAAAAACCGGa                                            >  1:1698240/4‑68 (MQ=255)        ttccccAGCACCATCCCTAGCCCCACTAACATCATAATAAAGGTCATTGCCGTTTCCGAAAAACCGGa                                            >  1:510347/4‑68 (MQ=255)        ttccccAGCACCATCCCTAGCCCCACTAACATCATAATAAAGGTCATTGCCGTTTCCGAAAAACCGGa                                            >  1:1724990/4‑68 (MQ=255)        ttccccAGCACCATCCCTAGCCCCACTAACATCATAATAAAGGTCATTGCCGTTTCCGAAAAACCGGa                                            >  1:1525415/4‑68 (MQ=255)          ccccAGCACCATCCCTAGCCCCACTAACATCATAATAAAGGTCATTGCCGTTTCCGAAAAACCGGAAAt                                         <  1:1126642/68‑1 (MQ=255)          ccccAGCACCATCCCTAGCCCCACTAACATCATAATAAAGGTCATTGCCGTTTCCGAAAAACCGGAAAt                                         <  1:686590/68‑1 (MQ=255)          ccccAGCACCATCCCTAGCCCCACTAACATCATAATAAAGGTCATTGCCGTTTCCGAAAAACCGGAAAt                                         <  1:779293/68‑1 (MQ=255)          ccccAGCACCATCCCTAGCCCCACTAACATCATAATAAAGGTCATTGCCGTTTCCGAAAAACCGGAAAt                                         <  1:881905/68‑1 (MQ=255)          ccccAGCACAATCCCTAGCCCCACTAACATCATAATAAAGGTCATTGCCGTTTCCGAAAAACCGGAAAt                                         <  1:3065101/68‑1 (MQ=255)            ccAGCACCATCCCTAGCCCCACTAACATCATAATAAAGGTCATTGCCGTTTCCGAAAAACCGGAAAt                                         <  1:2942522/67‑1 (MQ=255)            ccAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATa                                        <  1:2424924/68‑1 (MQ=255)            ccAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATa                                        <  1:1879808/68‑1 (MQ=255)            ccAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATa                                        <  1:1307242/68‑1 (MQ=255)             cAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATa                                        <  1:1781244/67‑1 (MQ=255)              aGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATaaa                                      <  1:1496303/68‑1 (MQ=255)              aGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATaaa                                      <  1:2417387/68‑1 (MQ=255)              aGCACCATCCCTAGCCCAACTAACATCATAAAAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATaaa                                      <  1:148944/68‑1 (MQ=255)                cACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAAcat                                   <  1:2718540/69‑1 (MQ=255)                    aTCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATg                               <  1:2749849/69‑1 (MQ=255)                    aTCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATg                               <  1:2938120/69‑1 (MQ=255)                    aTCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATg                               <  1:2893170/69‑1 (MQ=255)                    aTCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATg                               <  1:2764669/69‑1 (MQ=255)                    aTCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATg                               <  1:3230674/69‑1 (MQ=255)                    aTCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATg                               <  1:2273244/69‑1 (MQ=255)                    aTCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATg                               <  1:51405/69‑1 (MQ=255)                    aTCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATg                               <  1:2261442/69‑1 (MQ=255)                    aTCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATg                               <  1:1587379/69‑1 (MQ=255)                    aTCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATg                               <  1:700509/69‑1 (MQ=255)                    aTCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATg                               <  1:1438787/69‑1 (MQ=255)                    aTCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATg                               <  1:1310924/69‑1 (MQ=255)                     tCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATg                               <  1:3335219/68‑1 (MQ=255)                     tCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATg                               <  1:1937052/68‑1 (MQ=255)                     tCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATg                               <  1:951868/68‑1 (MQ=255)                     tCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATg                               <  1:1068730/68‑1 (MQ=255)                       ccTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGt                              >  1:969767/1‑67 (MQ=255)                       ccTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTa                             >  1:2849537/1‑68 (MQ=255)                       ccTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTa                             >  1:1223642/1‑68 (MQ=255)                       ccTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTa                             >  1:2634043/1‑68 (MQ=255)                       ccTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTa                             >  1:1843854/1‑68 (MQ=255)                       ccTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTa                             >  1:1851211/1‑68 (MQ=255)                       ccTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTAt                            >  1:2270848/1‑69 (MQ=255)                       ccTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTAt                            >  1:742400/1‑69 (MQ=255)                       ccTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTAt                            >  1:1984740/1‑69 (MQ=255)                       ccTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTAt                            >  1:400652/1‑69 (MQ=255)                       ccTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTAt                            >  1:2390502/1‑69 (MQ=255)                       ccTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTAt                            >  1:2677270/1‑69 (MQ=255)                       ccTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTAt                            >  1:2953455/1‑69 (MQ=255)                       ccTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTAt                            >  1:2492293/1‑69 (MQ=255)                       ccTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTAt                            >  1:2494838/1‑69 (MQ=255)                          aGCCCAACTAACATCATAATAAAGGTCATCGCCGttt                                                         >  1:325695/1‑37 (MQ=255)                            cccAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTATGGCtt                       >  1:11908/1‑69 (MQ=255)                            cccAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTATGGCtt                       >  1:13471/1‑69 (MQ=255)                            cccAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTATGGCtt                       >  1:1707455/1‑69 (MQ=255)                            cccAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACAGGAAATAAACATCATGTATGGCtt                       >  1:1366423/1‑69 (MQ=255)                                  tAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATg                               >  1:145364/1‑55 (MQ=255)                                  tAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATg                               >  1:3345858/1‑55 (MQ=255)                                  tAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATg                               >  1:1451033/1‑55 (MQ=255)                                  tAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATg                               >  1:1545609/1‑55 (MQ=255)                                  tAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATg                               >  1:2994310/1‑55 (MQ=255)                                  tAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATg                               >  1:2273992/1‑55 (MQ=255)                                  tAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTATGGCTTTACGTa                 <  1:546682/69‑1 (MQ=255)                                  tAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATCTATGGCTTTACGTa                 <  1:1261412/69‑1 (MQ=255)                                    aCATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAAtt                                        <  1:345032/44‑2 (MQ=255)                                        catAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTATGGCTTTACGTAGCTGaa           <  1:536209/69‑1 (MQ=255)                                               aaGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTATGGCTTTACGTAGCTGAACCAGGCa    <  1:3425838/69‑1 (MQ=255)                                               aaGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTATGGCTTTACGTAGCTGAACCAGGCa    <  1:2871866/69‑1 (MQ=255)                                               aaGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTATGGCTTTACGTAGCTGAACCAGGCa    <  1:3400447/69‑1 (MQ=255)                                               aaGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTATGGCTTTACGTAGCTGAACCAGGCa    <  1:3350245/69‑1 (MQ=255)                                                 ggTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTATGGCTTTACGTAGCTGAACCAGGCaaa  >  1:953080/1‑69 (MQ=255)                                                 ggTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTATGGCTTTACGTAGCTGAACCAGGCa    <  1:2602679/67‑1 (MQ=255)                                                      tCGCCGTTTCCGAAAAACCGGAAATAAACATCATg                               <  1:1795903/35‑1 (MQ=255)                                                       |                                                               TAGCATATTTCCCAGCACCATCCCTAGCCCAACTAACATCATAATAAAGGTCATCGCCGTTTCCGAAAAACCGGAAATAAACATCATGTATGGCTTTACGTAGCTGAACCAGGCAAA  >  minE/240189‑240305

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