Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	414,381	T→C	17.2%	E235E (GAA→GAG)	ybeX ←	predicteed ion transport

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	414,381	0	T	C	17.2%	145.6 / 18.4	92	E235E (GAA→GAG)	ybeX	predicteed ion transport
Reads supporting (aligned to +/- strand):  ref base T (39/37);  new base C (8/8);  total (47/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.46e-01

GTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGG  >  minE/414319‑414448                                                               |                                                                    gTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACttcttc                                                                >  1:1182000/1‑68 (MQ=255) gTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCa                                                               >  1:2604201/1‑69 (MQ=255) gTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCa                                                               >  1:2892593/1‑69 (MQ=255) gTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCa                                                               >  1:3092110/1‑69 (MQ=255) gTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCa                                                               >  1:3209205/1‑69 (MQ=255)      gCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCAACCTCTTCATCg                                                            >  1:1733360/1‑67 (MQ=255)      gCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCAACCTCTTCATCGc                                                           >  1:1846429/1‑68 (MQ=255)      gCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCAACCTCTTCATCGc                                                           >  1:2286061/1‑68 (MQ=255)      gCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCAACCTCTTCATCGc                                                           >  1:1262177/1‑68 (MQ=255)      gCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCAACCTCTTCATCGc                                                           >  1:1250008/1‑68 (MQ=255)      gCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCAACCTCTTCATCGc                                                           >  1:2127010/1‑68 (MQ=255)      gCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCAACCTCTTCATCGc                                                           >  1:75130/1‑68 (MQ=255)      gCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCAACCTCTTCATCGc                                                           >  1:1133087/1‑68 (MQ=255)       ccACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCt                                                          <  1:2237091/68‑1 (MQ=255)       ccACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCt                                                          <  1:2766059/68‑1 (MQ=255)       ccACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCt                                                          <  1:2684765/68‑1 (MQ=255)       ccACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCCCt                                                          <  1:271083/68‑1 (MQ=255)            cgcCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTg                                                    <  1:2378141/69‑1 (MQ=255)            cgcCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTg                                                    <  1:258661/69‑1 (MQ=255)            cgcCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTg                                                    <  1:773864/69‑1 (MQ=255)            cgcCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTg                                                    <  1:3139829/69‑1 (MQ=255)            cgcCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTg                                                    <  1:3200136/69‑1 (MQ=255)            cgcCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTg                                                    <  1:416853/69‑1 (MQ=255)            cgcCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTg                                                    <  1:580092/69‑1 (MQ=255)            cgcCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTg                                                    <  1:712451/69‑1 (MQ=255)              ccGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTggg                                                  <  1:2811650/69‑1 (MQ=255)              ccGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTggg                                                  <  1:2968679/69‑1 (MQ=255)              ccGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTggg                                                  <  1:913243/69‑1 (MQ=255)                   aGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCg                                             <  1:59540/69‑1 (MQ=255)                   aGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCg                                             <  1:391294/69‑1 (MQ=255)                   aGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCg                                             <  1:2278427/69‑1 (MQ=255)                   aGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCg                                             <  1:3336512/69‑1 (MQ=255)                   aGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCg                                             <  1:2171991/69‑1 (MQ=255)                   aGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCg                                             <  1:1720980/69‑1 (MQ=255)                   aGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCg                                             <  1:3246797/69‑1 (MQ=255)                   aGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCg                                             <  1:2620741/69‑1 (MQ=255)                   aGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCg                                             <  1:673345/69‑1 (MQ=255)                   aGATGCCCAAATGCCTGCATCACCAGAACACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCg                                             >  1:3039698/1‑69 (MQ=255)                    gttGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCg                                             <  1:1685743/66‑1 (MQ=255)                       gCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCg                                             >  1:1681711/1‑65 (MQ=255)                       gCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACg                                         >  1:1487577/1‑69 (MQ=255)                       gCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACg                                         >  1:1598747/1‑69 (MQ=255)                       gCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACg                                         >  1:2608929/1‑69 (MQ=255)                       gCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACg                                         >  1:349649/1‑69 (MQ=255)                       gCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACg                                         >  1:2815690/1‑69 (MQ=255)                       gCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACg                                         >  1:2787481/1‑69 (MQ=255)                       gCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACg                                         >  1:171286/1‑69 (MQ=255)                       gCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACg                                         >  1:525342/1‑69 (MQ=255)                       gCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACg                                         >  1:2072573/1‑69 (MQ=255)                       gCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACg                                         >  1:1937220/1‑69 (MQ=255)                       gCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACg                                         >  1:728470/1‑69 (MQ=255)                       gCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACg                                         >  1:1004109/1‑69 (MQ=255)                       gCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACg                                         >  1:1838353/1‑69 (MQ=255)                              tGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCa                                                               >  1:58468/1‑40 (MQ=255)                              tGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTaa                                                        >  1:2490812/1‑47 (MQ=255)                              tGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGt                                   >  1:62528/1‑68 (MQ=255)                              tGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGt                                   >  1:1183764/1‑68 (MQ=255)                              tGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGt                                   >  1:1820407/1‑68 (MQ=255)                              tGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGt                                   >  1:1593737/1‑68 (MQ=255)                              tGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGt                                   >  1:2619585/1‑68 (MQ=255)                              tGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGt                                   >  1:2657448/1‑68 (MQ=255)                                       cACCAGACCACCGATAGTGTCAACCTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTc                         <  1:1383111/69‑1 (MQ=255)                                       cACCAGACCACCGATAGTGTCAACCTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTc                         <  1:2911834/69‑1 (MQ=255)                                       cACCAGACCACCGATAGTGTCAACCTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTc                         <  1:1161460/69‑1 (MQ=255)                                       cACCAGACCACCGATAGTGTCAACCTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTc                         <  1:2940881/69‑1 (MQ=255)                                       cACCAGACCACCGATAGTGTCAACCTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTc                         <  1:2207558/69‑1 (MQ=255)                                       cACCAGACCACCGATAGTGTCAACCTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTc                         <  1:798127/69‑1 (MQ=255)                                       cACCAGACCACCGATAGTGTCAACCTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTc                         <  1:464254/69‑1 (MQ=255)                                       cACCAGACCACCGATAGTGTCAACCTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTc                         <  1:1174254/69‑1 (MQ=255)                                        aCCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTc                         >  1:405108/1‑68 (MQ=255)                                                 ccGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGcc               <  1:3231823/69‑1 (MQ=255)                                                 ccGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGcc               <  1:1309015/69‑1 (MQ=255)                                                 ccGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGcc               <  1:1159167/69‑1 (MQ=255)                                                 ccGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGcc               <  1:1357497/69‑1 (MQ=255)                                                 ccGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGcc               <  1:861356/69‑1 (MQ=255)                                                  cGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGcc               <  1:1159873/68‑1 (MQ=255)                                                         gtCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGcgcg       <  1:3389150/69‑1 (MQ=255)                                                         gtCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGcgcg       <  1:1626215/69‑1 (MQ=255)                                                         gtCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGcgcg       <  1:1877387/69‑1 (MQ=255)                                                         gtCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGcgcg       <  1:1100447/69‑1 (MQ=255)                                                         gtCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGcgcg       <  1:1968868/69‑1 (MQ=255)                                                         gtCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGcgcg       <  1:2012215/69‑1 (MQ=255)                                                          tCGACTTCTTCATCGCTAAAGTGGGTTCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGt      >  1:1513247/1‑69 (MQ=255)                                                          tCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTgcgc        >  1:1626994/1‑67 (MQ=255)                                                          tCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGcgcg       >  1:796808/1‑68 (MQ=255)                                                          tCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGt      >  1:1312257/1‑69 (MQ=255)                                                          tCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGt      >  1:3123427/1‑69 (MQ=255)                                                          tCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGt      >  1:1322180/1‑69 (MQ=255)                                                          tCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGt      >  1:2931989/1‑69 (MQ=255)                                                          tCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGt      >  1:1856024/1‑69 (MQ=255)                                                            gACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTAc    >  1:654281/1‑69 (MQ=255)                                                            gACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTAc    >  1:2178800/1‑69 (MQ=255)                                                            gACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTAc    >  1:3074892/1‑69 (MQ=255)                                                               ttcttcATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACgg  >  1:2134113/1‑68 (MQ=255)                                                               ttcttcATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACgg  >  1:2605884/1‑68 (MQ=255)                                                               |                                                                    GTTTCGCCACGCGCCGGAAGATGCCCAAATGCCTGCATCACCAGACCACCGATAGTGTCGACTTCTTCATCGCTAAAGTGGGTGCCGAACGCTTCGTTGAAGTCTTCAATGGAAGCCAGTGCGCGTACGG  >  minE/414319‑414448

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