Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	997,596	C→G	42.1%	T69T (ACC→ACG)	ynfM →	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	997,596	0	C	G	42.1%	43.8 / 97.6	95	T69T (ACC→ACG)	ynfM	predicted transporter
Reads supporting (aligned to +/- strand):  ref base C (30/25);  new base G (35/5);  total (65/30)
Fisher's exact test for biased strand distribution p-value = 7.18e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

TGCACTTCTCTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTGTCCATTTCCACGGCGATGTTGGCTATTGGTTTGCTGTTTA  >  minE/997536‑997660                                                             |                                                                 tGCACTTCTCTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCcgcg                                                            >  1:1033373/1‑67 (MQ=255) tGCACTTCTCTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCcgcg                                                            >  1:501031/1‑67 (MQ=255) tGCACTTCTCTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCcgcg                                                            >  1:2623455/1‑67 (MQ=255) tGCACTTCTCTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCcgcg                                                            >  1:2905883/1‑67 (MQ=255) tGCACTTCTCTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCcgcg                                                            >  1:1722750/1‑67 (MQ=255)  gCACTTCTCTATTGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAAc                                                         >  1:1006319/1‑69 (MQ=255)  gCACTTCTCTATTGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAAc                                                         >  1:1968442/1‑69 (MQ=255)  gCACTTCTCTATTGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAAc                                                         >  1:2280636/1‑69 (MQ=255)  gCACTTCTCTATTGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAAc                                                         >  1:1785849/1‑69 (MQ=255)  gCACTTCTCTATTGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAAc                                                         >  1:1760076/1‑69 (MQ=255)  gCACTTCTCTATTGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAAc                                                         >  1:1756165/1‑69 (MQ=255)  gCACTTCTCTATTGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAAc                                                         >  1:3319515/1‑69 (MQ=255)  gCACTTCTCTATTGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAAc                                                         >  1:3241320/1‑69 (MQ=255)  gCACTTCTCTATTGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAAc                                                         >  1:3143182/1‑69 (MQ=255)  gCACTTCTCTATTGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAAc                                                         >  1:3204058/1‑69 (MQ=255)  gCACTTCTCTATTGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAAc                                                         >  1:1337833/1‑69 (MQ=255)  gCACTTCTCTATTGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAAc                                                         >  1:911876/1‑69 (MQ=255)  gCACTTCTCTATTGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAAc                                                         >  1:2996629/1‑69 (MQ=255)  gCACTTCTCTATTGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGACTTTGGCTTAACGCCCGCGaa                                                          >  1:2656387/1‑68 (MQ=255)       tctcTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACagtag                                                    <  1:1267277/69‑1 (MQ=255)       tctcTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACagtag                                                    <  1:733758/69‑1 (MQ=255)       tctcTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACagtag                                                    <  1:1909582/69‑1 (MQ=255)       tctcTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACagtag                                                    <  1:3194678/69‑1 (MQ=255)       tctcTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACagtag                                                    <  1:1678886/69‑1 (MQ=255)       tctcTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACagtag                                                    <  1:1672131/69‑1 (MQ=255)       tctcTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACagtag                                                    <  1:825562/69‑1 (MQ=255)       tctcTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACagtag                                                    <  1:3307881/69‑1 (MQ=255)       tctcTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACagtag                                                    <  1:1416768/69‑1 (MQ=255)        ctctATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACagtag                                                    <  1:103315/68‑1 (MQ=255)             ttGTTTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAACAGTAGTATTTc                                              >  1:2227199/1‑69 (MQ=255)             ttGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTc                                              >  1:2745964/1‑69 (MQ=255)             ttGTGTGCAACCTATCGTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAACAGTAGTATTTc                                              >  1:847178/1‑69 (MQ=255)             ttGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGTCTTAACGCCCGCGAACAGTAGTATTTc                                              >  1:2104548/1‑69 (MQ=255)             ttGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAACAGTAGTAttt                                               >  1:777747/1‑68 (MQ=255)             ttGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAACAGTAGTAttt                                               >  1:254129/1‑68 (MQ=255)             ttGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAACAGTAGTAttt                                               >  1:556764/1‑68 (MQ=255)             ttGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAACAGTAGTATTTc                                              >  1:2062395/1‑69 (MQ=255)             ttGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAACAGTAGTATTTc                                              >  1:914507/1‑69 (MQ=255)             ttGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAACAGTAGTATTTc                                              >  1:1207872/1‑69 (MQ=255)             ttGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAACAGTAGTATTTc                                              >  1:2993245/1‑69 (MQ=255)             ttGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAACAGTAGTATTTc                                              >  1:1585429/1‑69 (MQ=255)             ttGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAACAGTAGTATTTc                                              >  1:2692217/1‑69 (MQ=255)             ttGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAACAGTAGTATTTc                                              >  1:1675805/1‑69 (MQ=255)             ttGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAACAGTAGTATTTc                                              >  1:1840184/1‑69 (MQ=255)             ttGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAACAGTAGTATTTc                                              >  1:2562564/1‑69 (MQ=255)             ttGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAACAGTAGTATTTc                                              >  1:2481845/1‑69 (MQ=255)             ttGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAACAGTAGTATTTc                                              >  1:794286/1‑69 (MQ=255)             ttGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAACAGTAGTATTTc                                              >  1:2276344/1‑69 (MQ=255)             ttGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAACAGTAGTATTTc                                              >  1:1958118/1‑69 (MQ=255)             ttGTGTGCAACCTATCCTTCCGGTGCTTTCACAAGAGTTTGGCTTAACGCCCGCGAACAGTAGTATTTc                                              >  1:1930898/1‑69 (MQ=255)                 gtgCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACt                                           >  1:2256555/1‑68 (MQ=255)                 gtgCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTg                                          >  1:665632/1‑69 (MQ=255)                 gtgCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTg                                          >  1:3159921/1‑69 (MQ=255)                 gtgCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTg                                          >  1:923915/1‑69 (MQ=255)                 gtgCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTg                                          >  1:3138947/1‑69 (MQ=255)                 gtgCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTg                                          >  1:3125851/1‑69 (MQ=255)                 gtgCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTg                                          >  1:882835/1‑69 (MQ=255)                 gtgCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTg                                          >  1:2865941/1‑69 (MQ=255)                 gtgCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTg                                          >  1:1548838/1‑69 (MQ=255)                 gtgCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTg                                          >  1:2157255/1‑69 (MQ=255)                 gtgCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTg                                          >  1:2159005/1‑69 (MQ=255)                 gtgCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTg                                          >  1:2527035/1‑69 (MQ=255)                    cAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTGTcc                                       >  1:2446567/1‑69 (MQ=255)                    cAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTGTcc                                       >  1:1522378/1‑69 (MQ=255)                    cAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTGTcc                                       >  1:2146625/1‑69 (MQ=255)                    cAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTGTcc                                       >  1:1612416/1‑69 (MQ=255)                    cAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTGTcc                                       >  1:1817517/1‑69 (MQ=255)                    cAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTGTcc                                       >  1:1096765/1‑69 (MQ=255)                    cAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTGTcc                                       >  1:2334981/1‑69 (MQ=255)                    cAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTGTc                                        <  1:1022141/68‑1 (MQ=255)                    cAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTGTc                                        <  1:655152/68‑1 (MQ=255)                    cAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTGTc                                        <  1:2783062/68‑1 (MQ=255)                    cAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTGTc                                        <  1:2115582/68‑1 (MQ=255)                    cAGCCTATCCTTACGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTGTc                                        <  1:1852525/68‑1 (MQ=255)                             cTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTGTCCATTTCCACg                              >  1:1246199/1‑69 (MQ=255)                             cTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTGTCCATTTCCACg                              >  1:2490550/1‑69 (MQ=255)                             cTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTGTCCATTTCCACg                              >  1:1900442/1‑69 (MQ=255)                                   gTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTGTCCATTTCCACGGCGAt                         <  1:642634/68‑1 (MQ=255)                                      cTTTCGCAGGAGTTTGGCTTAACGCCCGCGAACAGTAGTATTTCACTGTCCATTTCCACGGCGATGTTg                     >  1:2488653/1‑69 (MQ=255)                                      cTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTGTCCATTTCCACGGCGATGTTg                     >  1:916539/1‑69 (MQ=255)                                      cTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTGTCCATTTCCACGGCGATGTTg                     >  1:1390630/1‑69 (MQ=255)                                                    tGGCTTAACGCCCGCGAACAGTAGTATTTCACTGTCCATTTCAACGGc                            <  1:731160/48‑1 (MQ=255)                                                    tGGCTTAACGCCCGCGAACAGTAGTATTTCACTGTCCATTTCAACGGc                            <  1:950662/48‑1 (MQ=255)                                                    tGGCTTAACGCCCGCGAACAGTAGTATTTCACTGTCCATTTCAACGGc                            <  1:1660730/48‑1 (MQ=255)                                                    tGGCTTAACGCCCGCGAACAGTAGTATTTCACTGTCCATTTCAACGGc                            <  1:2613661/48‑1 (MQ=255)                                                    tGGCTTAACGCCCGCGAACAGTAGTATTTCACTGTCCATTTCAACGGc                            <  1:2722676/48‑1 (MQ=255)                                                        ttAACCCCCGCGAACAGTAGTATTTCACTGTCCATTTCCACGGCGATGTTGGCTATTGGTTTGCTGtt    <  1:449263/68‑1 (MQ=255)                                                        ttAACCCCCGCGAACAGTAGTATTTCACTGTCCATTTCCACGGCGATGTTGGCTATTGGTTTGCTGtt    <  1:2638632/68‑1 (MQ=255)                                                        ttAACCCCCGCGAACAGTAGTATTTCACTGTCCATTTCCACGGCGATGTTGGCTATTGGTTTGCTGtt    <  1:1188096/68‑1 (MQ=255)                                                        ttAACCCCCGCGAACAGTAGTATTTCACTGTCCATTTCCACGGCGATGTTGGCTATTGGTTTGCTGtt    <  1:1181894/68‑1 (MQ=255)                                                        ttAACCCCCGCGAACAGTAGTATTTCACTGTCCATTTCCACGGCGATGTTGGCTATTGGTTTGCTGtt    <  1:3082446/68‑1 (MQ=255)                                                        ttAACCCCCGCGAACAGTAGTATTTCACTGTCCATTTCCACGGCGATGTTGGCTATTGGTTTGCTGtt    <  1:3150166/68‑1 (MQ=255)                                                         tAACCCCCGCGAACAGTAGTATTTCACTGTCCATTTCCACGGCGATGTTGGCTATTTGTTTGCTGtt    <  1:3167040/67‑1 (MQ=255)                                                         tAACCCCCGCGAACAGTAGTATTTCACTGTCCATTTCCACGGCGATGTTGGCTATTGGTTTGCTGTTTa  <  1:2133254/69‑1 (MQ=255)                                                         tAACCCCCGCGAACAGTAGTATTTCACTGTCCATTTCCACGGCGATGTTGGCTATTGGTTTGCTGTTTa  <  1:1061106/69‑1 (MQ=255)                                                             |                                                                 TGCACTTCTCTATTGTGTGCAGCCTATCCTTCCGGTGCTTTCGCAGGAGTTTGGCTTAACCCCCGCGAACAGTAGTATTTCACTGTCCATTTCCACGGCGATGTTGGCTATTGGTTTGCTGTTTA  >  minE/997536‑997660

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