Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	339,753	C→T	37.7%	E32E (GAG→GAA)	copA ←	copper transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	339,753	0	C	T	37.7%	47.2 / 50.0	82	E32E (GAG→GAA)	copA	copper transporter
Reads supporting (aligned to +/- strand):  ref base C (7/44);  new base T (12/19);  total (19/63)
Fisher's exact test for biased strand distribution p-value = 1.45e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGCAGTGACCGCAGGACAG  >  minE/339693‑339815                                                             |                                                               cAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATc                                                        >  1:2665259/1‑69 (MQ=255) cAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATc                                                        >  1:3110731/1‑69 (MQ=255)    cTGTTCTGCACTGGGAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGTTCAACATCCg                                                      >  1:1652853/1‑68 (MQ=255)    cTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGTTc                                                              >  1:3298553/1‑60 (MQ=255)    cTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGTTCAACATCCgg                                                     >  1:169458/1‑69 (MQ=255)    cTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGTTCAACATCCgg                                                     >  1:2781151/1‑69 (MQ=255)    cTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGTTCAACATCCgg                                                     >  1:466510/1‑69 (MQ=255)    cTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGTTCAACATCCgg                                                     >  1:1625327/1‑69 (MQ=255)    cTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGTTCAACATCCgg                                                     >  1:3334172/1‑69 (MQ=255)    cTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGTTCAACATCCgg                                                     >  1:600462/1‑69 (MQ=255)    cTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGTTCAACATCCgg                                                     >  1:2222688/1‑69 (MQ=255)    cTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGTTCAACATCCgg                                                     >  1:1811000/1‑69 (MQ=255)    cTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGTTCAACATCCg                                                      >  1:6503/1‑68 (MQ=255)    cTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGTTCAACATCCg                                                      >  1:1032689/1‑68 (MQ=255)      gTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGAc                                                   >  1:3287250/1‑69 (MQ=255)      gTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGAc                                                   >  1:1793661/1‑69 (MQ=255)         cTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGc                                                 <  1:3380730/68‑1 (MQ=255)              cTGGCATTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGt                                              >  1:436049/1‑66 (MQ=255)              cTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTg                                               >  1:2883328/1‑65 (MQ=255)              cTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTg                                               >  1:1031308/1‑65 (MQ=255)                ggCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGTTCAACATCCGGGCGCTGTTCGAg                                         <  1:2320856/69‑1 (MQ=255)                ggCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGTTCAACATCCGGGCGCTGTTCGAg                                         <  1:972122/69‑1 (MQ=255)                ggCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGTTCAACATCCGGGCGCTGTTCGAg                                         <  1:2029869/69‑1 (MQ=255)                ggCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGTTCAACATCCGGGCGCTGTTCGAg                                         <  1:2371687/69‑1 (MQ=255)                ggCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGTTCAACATCCGGGCGCTGTTCGAg                                         <  1:608435/69‑1 (MQ=255)                ggCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGTTCAACATCCGGGCGCTGTTCGAg                                         <  1:2799445/69‑1 (MQ=255)                ggCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGTTCAACATCCGGGCGCTGTTCGAg                                         <  1:3159874/69‑1 (MQ=255)                ggCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGTTCAACATCCGGGCGCTGTTCGAg                                         <  1:323821/69‑1 (MQ=255)                ggCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGTTCAACATCCGGGCGCTGTTCGAg                                         <  1:583685/69‑1 (MQ=255)                ggCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGTTCAACATCCGGGCGCTGTTCGAg                                         <  1:1790424/69‑1 (MQ=255)                ggCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGTTCAACATCCGGGCGCTGTTCGAg                                         <  1:1667409/69‑1 (MQ=255)                ggCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGTTCAACATCCGGGCGCTGTTCGAg                                         <  1:1484756/69‑1 (MQ=255)                ggCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGTTCAACATCCGGGCGCTGTTCGAg                                         <  1:834346/69‑1 (MQ=255)                ggCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGTTCAACATCCGGGCGCTGTTCGAg                                         <  1:1815781/69‑1 (MQ=255)                ggCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGTTCAACATCCGGGCGCTGTTCGAg                                         <  1:1283800/69‑1 (MQ=255)                ggCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGTTCAACATCCGGGCGCTGTTCGAg                                         <  1:1182578/69‑1 (MQ=255)                ggCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGTTCAACATCCGGGCGCTGTTCGAg                                         <  1:1163203/69‑1 (MQ=255)                ggCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGTTCAACATCCGGGCGCTGTTCGAg                                         <  1:1128581/69‑1 (MQ=255)                ggCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGTTCAACATCCGGGCGCTGTTCGAg                                         <  1:1930759/69‑1 (MQ=255)                         ggTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGActttcttt                                <  1:1963729/69‑1 (MQ=255)                         ggTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGActttcttt                                <  1:588327/69‑1 (MQ=255)                         ggTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGActttcttt                                <  1:1151470/69‑1 (MQ=255)                         ggTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGActttcttt                                <  1:3221672/69‑1 (MQ=255)                         ggTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGActttcttt                                <  1:306809/69‑1 (MQ=255)                         ggTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGActttcttt                                <  1:2943828/69‑1 (MQ=255)                         ggTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGActttcttt                                <  1:1584144/69‑1 (MQ=255)                         ggTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGActttcttt                                <  1:2557850/69‑1 (MQ=255)                         ggTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGActttcttt                                <  1:245439/69‑1 (MQ=255)                           tAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCa                              <  1:1543447/69‑1 (MQ=255)                           tAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCAGGACGCTGTTCAAGACTTTCTTTCa                              <  1:1768198/69‑1 (MQ=255)                             aCGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGActttcttt                                <  1:253535/65‑1 (MQ=255)                               gTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCAcgcg                          <  1:1005588/69‑1 (MQ=255)                               gTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCAcgcg                          <  1:3366482/69‑1 (MQ=255)                               gTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCAcgcg                          <  1:703147/69‑1 (MQ=255)                               gTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCAcgcg                          <  1:2149268/69‑1 (MQ=255)                               gTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCAcgcg                          <  1:1863311/69‑1 (MQ=255)                               gTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCAcgcg                          <  1:1364639/69‑1 (MQ=255)                               gTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCAcgcg                          <  1:2013203/69‑1 (MQ=255)                               gTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCAcgcg                          <  1:1535979/69‑1 (MQ=255)                               gTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCAcgcg                          <  1:636661/69‑1 (MQ=255)                               gTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCAcgcg                          <  1:1624146/69‑1 (MQ=255)                                    cTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCAGGCGTTTaa                     <  1:1118273/69‑1 (MQ=255)                                    cTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTaa                     <  1:2812583/69‑1 (MQ=255)                                         gTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGc                  <  1:536690/67‑1 (MQ=255)                                         gTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGc                  <  1:2428594/67‑1 (MQ=255)                                         gTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGc                  <  1:3371408/67‑1 (MQ=255)                                         gTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGc                  <  1:2209867/67‑1 (MQ=255)                                         gTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGc                  <  1:2589174/67‑1 (MQ=255)                                         gTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGc                  <  1:1665392/67‑1 (MQ=255)                                         gTGATAGACACATCAGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGc                  <  1:1414364/67‑1 (MQ=255)                                          tGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGc                  <  1:2434558/66‑1 (MQ=255)                                                       cGCCTGCTCAACATCCGGACGCTGTTCAAGTCTTTCTTTCACGCGTTTAACGCAGTGACCGCAGGACAg  <  1:1756210/69‑1 (MQ=255)                                                       cGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGCAGTGACCGCAGGATAg  <  1:56281/69‑1 (MQ=255)                                                       cGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGCAGTGACCGCAGGACAg  <  1:3281898/69‑1 (MQ=255)                                                       cGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGCAGTGACCGCAGGACAg  <  1:2963232/69‑1 (MQ=255)                                                       cGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGCAGTGACCGCAGGACAg  <  1:2794639/69‑1 (MQ=255)                                                       cGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGCAGTGACCGCAGGACAg  <  1:2530559/69‑1 (MQ=255)                                                       cGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGCAGTGACCGCAGGACAg  <  1:2375397/69‑1 (MQ=255)                                                       cGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGCAGTGACCGCAGGACAg  <  1:2138617/69‑1 (MQ=255)                                                       cGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGCAGTGACCGCAGGACAg  <  1:769954/69‑1 (MQ=255)                                                       cGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGCAGTGACCGCAGGACAg  <  1:855899/69‑1 (MQ=255)                                                       cGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGCAGTGACCGCAGGACAg  <  1:1006219/69‑1 (MQ=255)                                                             |                                                               CAGCTGTTCTGCACTGGCAGTCCCGGTAACGTGCGCTTCAGTGATAGACACATCCGCCTGCTCAACATCCGGACGCTGTTCAAGACTTTCTTTCACGCGTTTAACGCAGTGACCGCAGGACAG  >  minE/339693‑339815

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