Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	338,643	A→G	13.1%	R402R (CGT→CGC)	copA ←	copper transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	338,643	0	A	G	13.1%	175.8 / 12.2	98	R402R (CGT→CGC)	copA	copper transporter
Reads supporting (aligned to +/- strand):  ref base A (43/42);  new base G (9/4);  total (52/46)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

CTGGGCCTGGCGCACCATGCGAATGATTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACTGTCCCGGCATGGACGCTATCGCCTTCGCCTTTTTGCT  >  minE/338580‑338710                                                                |                                                                    ctGGGCCTGGCGCACCATACGAATAATCCGCGACAGCGTAGTATGGCTGCCAACCGCACTGGCGCGa                                                                  >  1:1052026/1‑67 (MQ=255) ctGGGCCTGGCGCACCATACGAATAATCCGCGACAGCGTAGTATGGCTGCCAACCGCACTGGCGCGa                                                                  >  1:2736133/1‑67 (MQ=255) ctGGGCCTGGCGCACCATACGAATAATCCGCGACAGCGTAGTATGGCTGCCAACCGCACTGGCGCGa                                                                  >  1:2671278/1‑67 (MQ=255) ctGGGCCTGGCGCACCATACGAATAATCCGCGACAGCGTAGTATGGCTGCCAACCGCACTGGCGCGa                                                                  >  1:3122859/1‑67 (MQ=255) ctGGGCCTGGCGCACCATACGAATAATCCGCGACAGCGTAGTATGGCTGCCAACCGCACTGGCGCGa                                                                  >  1:323534/1‑67 (MQ=255) ctGGGCCTGGCGCACCATACGAATAATCCGCGACAGCGTAGTATGGCTGCCAACCGCACTGGCGCGa                                                                  >  1:178595/1‑67 (MQ=255) ctGGGCCTGGCGCACCATACGAATAATCCGCGACAGCGTAGTATGGCTGCCAACCGCACTGGCGCGa                                                                  >  1:393086/1‑67 (MQ=255) ctGGGCCTGGCGCACCATACGAATAATCCGCGACAGCGTAGTATGGCTGCCAACCGCACTGGCGCGa                                                                  >  1:1504016/1‑67 (MQ=255) ctGGGACTGGCGCACCATACGAATAATCCGCGACAGCGTAGTATGGCTGCCAACCGCACTGGCGCGa                                                                  >  1:1461490/1‑67 (MQ=255)            gcACCATGCGAATGATTCGTGACAGCGTAGTATGGGTGCCAACCGCACTGGCACGAAACAGCACACTGc                                                     >  1:1586526/1‑69 (MQ=255)            gcACCATGCGAATGATTCGTGACAGCGTAGTATGGCTGCCACCCGCACTGGCACGAAACAGCACACTGc                                                     >  1:1454747/1‑69 (MQ=255)            gcACCATGCGAATGATTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACt                                                       >  1:1687182/1‑67 (MQ=255)            gcACCATGCGAATGATTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGc                                                     >  1:814090/1‑69 (MQ=255)            gcACCATGCGAATGATTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGc                                                     >  1:3277134/1‑69 (MQ=255)            gcACCATGCGAATGATTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGc                                                     <  1:3308340/69‑1 (MQ=255)            gcACCATGCGAATGATTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGc                                                     <  1:2009386/69‑1 (MQ=255)            gcACCATGCGAATGATTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGc                                                     >  1:2365174/1‑69 (MQ=255)            gcACCATGCGAATGATTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGc                                                     <  1:2867180/69‑1 (MQ=255)            gcACCATGCGAATGATTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGc                                                     >  1:2423610/1‑69 (MQ=255)            gcACCATGCGAATGATTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGc                                                     >  1:1338481/1‑69 (MQ=255)            gcACCATGCGAATGATTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGc                                                     >  1:494562/1‑69 (MQ=255)            gcACCATGCGAATGATTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAAAAGCACACTGc                                                     >  1:770194/1‑69 (MQ=255)                    cGAATGATTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGt                                             >  1:2441634/1‑69 (MQ=255)                    cGAATGATTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGt                                             >  1:466175/1‑69 (MQ=255)                    cGAATGATTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGt                                             >  1:346324/1‑69 (MQ=255)                    cGAATGATTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGt                                             >  1:3427776/1‑69 (MQ=255)                    cGAATGATTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGt                                             >  1:523060/1‑69 (MQ=255)                    cGAATGATTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGt                                             >  1:2821399/1‑69 (MQ=255)                    cGAATGATTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGt                                             >  1:1537746/1‑69 (MQ=255)                    cGAATGATTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGt                                             >  1:3359030/1‑69 (MQ=255)                    cGAATGATTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGt                                             >  1:2062798/1‑69 (MQ=255)                    cGAATGATTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGt                                             >  1:715090/1‑69 (MQ=255)                    cGAATGATTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGt                                             >  1:793665/1‑69 (MQ=255)                          aTTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACt                                       >  1:2849048/1‑69 (MQ=255)                          aTTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACt                                       >  1:938424/1‑69 (MQ=255)                          aTTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACt                                       >  1:3359232/1‑69 (MQ=255)                          aTTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACt                                       >  1:1399607/1‑69 (MQ=255)                          aTTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACt                                       >  1:1938935/1‑69 (MQ=255)                          aTTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACt                                       >  1:339355/1‑69 (MQ=255)                          aTTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACt                                       >  1:372609/1‑69 (MQ=255)                          aTTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACt                                       >  1:464238/1‑69 (MQ=255)                          aTTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACt                                       >  1:1562279/1‑69 (MQ=255)                          aTTCGTGACAGCGTAGTATGGATGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACt                                       >  1:2821754/1‑69 (MQ=255)                                     cGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACTGTCCCGGCacg                            <  1:2407331/69‑3 (MQ=255)                                     cGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACTGTCCCGGCATg                            <  1:3145601/69‑1 (MQ=255)                                     cGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACTGTCCCGGCATg                            <  1:623897/69‑1 (MQ=255)                                     cGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACTGTCCCGGCATg                            <  1:2748659/69‑1 (MQ=255)                                     cGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACTGTCCCGGCATg                            <  1:1054429/69‑1 (MQ=255)                                     cGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACTGTCCCGGCATg                            <  1:1056646/69‑1 (MQ=255)                                     cGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACTGTCCCGGCATg                            <  1:1119969/69‑1 (MQ=255)                                     cGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACTGTCCCGGCATg                            <  1:1268409/69‑1 (MQ=255)                                     cGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACTGTCCCGGCATg                            <  1:1464644/69‑1 (MQ=255)                                     cGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACTGTCCCGGCATg                            <  1:1661929/69‑1 (MQ=255)                                     cGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACTGTCCCGGCATg                            <  1:2100660/69‑1 (MQ=255)                                     cGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACTGTCCCGGCATg                            <  1:1841982/69‑1 (MQ=255)                                     cGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACTGTCCCGGCATg                            <  1:1686069/69‑1 (MQ=255)                                     cGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACTGTCCCGGCATg                            <  1:1779582/69‑1 (MQ=255)                                     cGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACTGTCCAGGCATg                            <  1:441056/69‑1 (MQ=255)                                      gtagtaTGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACTGTCCCGGCATg                            <  1:2065378/68‑1 (MQ=255)                                      gtagtaTGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACTGTCCCGGCATg                            <  1:2064827/68‑1 (MQ=255)                                      gtagtaTGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACTGTCCCGGCATg                            <  1:1903902/68‑1 (MQ=255)                                              gCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCtc                                              <  1:2879051/42‑2 (MQ=255)                                              gCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTg                                              <  1:1722120/42‑1 (MQ=255)                                              gCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTg                                              <  1:502474/42‑1 (MQ=255)                                              gCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTg                                              <  1:521029/42‑1 (MQ=255)                                              gCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTg                                              <  1:702428/42‑1 (MQ=255)                                              gCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTg                                              <  1:2692636/42‑1 (MQ=255)                                              gCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTg                                              <  1:3276074/42‑1 (MQ=255)                                              gCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTg                                              <  1:2935569/42‑1 (MQ=255)                                              gCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTg                                              <  1:3041450/42‑1 (MQ=255)                                              gCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTg                                              <  1:2150170/42‑1 (MQ=255)                                              gCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTg                                              <  1:3045715/42‑1 (MQ=255)                                              gCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTg                                              <  1:3108724/42‑1 (MQ=255)                                              gCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTg                                              <  1:2233712/42‑1 (MQ=255)                                              gCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACTGTCCCGGCATGGACGCTAtc                   >  1:2613143/1‑69 (MQ=255)                                              gCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACTGTCCCGGCATGGACGCTAtc                   >  1:68984/1‑69 (MQ=255)                                              gCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACTGTCCCGGCATGGACGCTAtc                   >  1:499794/1‑69 (MQ=255)                                              gCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACTGTCCCGGCATGGACGCTAtc                   >  1:1699290/1‑69 (MQ=255)                                              gCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACTGTCCCGGCATGGACGCTAtc                   >  1:1926524/1‑69 (MQ=255)                                              gCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACTGTCCCGGCATGGACGCTAtc                   >  1:1991519/1‑69 (MQ=255)                                              gCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGGCCTg                                              <  1:1385374/42‑1 (MQ=255)                                                tGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACTGTCCCGGCATGGACGCTAtcgc                 <  1:125212/69‑1 (MQ=255)                                                tGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACTGTCCCGGCATGGACGCTAtcgc                 <  1:371442/69‑1 (MQ=255)                                                tGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACTGTCCCGGCATGGACGCTAtcgc                 <  1:2798052/69‑1 (MQ=255)                                                tGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACTGTCCCGGCATGGACGCTAtcgc                 <  1:2713814/69‑1 (MQ=255)                                                tGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACTGTCCCGGCATGGACGCTAtcgc                 <  1:3405358/69‑1 (MQ=255)                                                tGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACTGTCCCGGCATGGACGCTAtcgc                 <  1:2413280/69‑1 (MQ=255)                                                tGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACTGTCCCGGCATGGACGCTAtcgc                 <  1:2425458/69‑1 (MQ=255)                                                     aCCGCACTGGCGCGAAACAGCACGCTACCGTCCTGCACCACTGTCCCGGCATgg                           <  1:989735/54‑1 (MQ=255)                                                     aCCGCACTGGCGCGAAACAGCACGCTACCGTCCTGCACCACTGTCCCGGCATgg                           <  1:600313/54‑1 (MQ=255)                                                     aCCGCACTGGCGCGAAACAGCACGCTACCGTCCTGCACCACTGTCCCGGCATgg                           <  1:420051/54‑1 (MQ=255)                                                     aCCGCACTGGCGCGAAACAGCACGCTACCGTCCTGCACCACTGTCCCGGCATgg                           <  1:2040368/54‑1 (MQ=255)                                                       cGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACTGTCCCGGCATGGACGCTAtcgccttcgcc          >  1:2907941/1‑69 (MQ=255)                                                       cGCACTAGCACGAAACAGCACACTGCCGTCCTGTACCACTGTCCCGGCATGGACGCTAtcgccttcgcc          >  1:1802/1‑69 (MQ=255)                                                             ggCACGAAACAGCACACTGCCGTCCTGTACCACTGTCCCGGCATGGACGCTATCGCCTTCGCCTTTTtg    >  1:253635/1‑69 (MQ=255)                                                             ggCACGAAACAGCACACTGCCGTCCTGTACCACTGTCCCGGCATGGACGCTATCGCCTTCGCCTTTTtg    >  1:361380/1‑69 (MQ=255)                                                             ggCACGAAACAGCACACTGCCGTCCTGTACCACTGTCCCGGCATGGACGCTATCGCCTTCGCCTTTTtg    >  1:3263155/1‑69 (MQ=255)                                                             ggCACGAAACAGCACACTGCCGTCCTGTAC‑‑CTGTCCCGGCATGGACGCTATCGCCTTCGCCTTTTtgct  >  1:895146/1‑69 (MQ=255)                                                                |                                                                    CTGGGCCTGGCGCACCATGCGAATGATTCGTGACAGCGTAGTATGGCTGCCAACCGCACTGGCACGAAACAGCACACTGCCGTCCTGTACCACTGTCCCGGCATGGACGCTATCGCCTTCGCCTTTTTGCT  >  minE/338580‑338710

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