Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,034,921	A→C	34.0%	A453A (GCA→GCC)	tolC →	transport channel

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,034,921	0	A	C	34.0%	78.1 / 74.2	97	A453A (GCA→GCC)	tolC	transport channel
Reads supporting (aligned to +/- strand):  ref base A (22/42);  new base C (26/7);  total (48/49)
Fisher's exact test for biased strand distribution p-value = 3.87e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGT  >  minE/2034853‑2034983                                                                     |                                                               gAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCa                                                                <  1:1030225/69‑1 (MQ=255) gAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCa                                                                <  1:2909386/69‑1 (MQ=255) gAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCa                                                                <  1:2108120/69‑1 (MQ=255)      ggATCTGCTGGCACTGAACATTGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCa                                                           <  1:1996943/69‑1 (MQ=255)      ggATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCa                                                           <  1:2768266/69‑1 (MQ=255)      ggATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCa                                                           <  1:9574/69‑1 (MQ=255)      ggATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCa                                                           <  1:1973980/69‑1 (MQ=255)      ggATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCa                                                           <  1:1371432/69‑1 (MQ=255)      ggATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCa                                                           <  1:1252647/69‑1 (MQ=255)         tCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAAc                                                        >  1:1415617/1‑69 (MQ=255)         tCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAAc                                                        >  1:2386659/1‑69 (MQ=255)         tCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAAc                                                        >  1:3151965/1‑69 (MQ=255)         tCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAAc                                                        >  1:550461/1‑69 (MQ=255)         tCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAAc                                                        >  1:465655/1‑69 (MQ=255)         tCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAAc                                                        >  1:2736073/1‑69 (MQ=255)         tCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAAc                                                        >  1:586779/1‑69 (MQ=255)         tCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAAc                                                        >  1:663706/1‑69 (MQ=255)         tCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTCCACCGCAAAc                                                        >  1:3418365/1‑69 (MQ=255)               ggCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCgg                                                   <  1:2304642/68‑1 (MQ=255)               ggCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCgg                                                   <  1:340694/68‑1 (MQ=255)                     gAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGaa                                            >  1:331037/1‑69 (MQ=255)                     gAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGaa                                            >  1:3031351/1‑69 (MQ=255)                        cAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCCCCGCAAACGCCGGAACAGAATGc                                         >  1:2615713/1‑69 (MQ=255)                        cAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCCCCGCAAACGCCGGAACAGAATGc                                         >  1:2653587/1‑69 (MQ=255)                        cAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCCCCGCAAACGCCGGAACAGAATGc                                         >  1:3380923/1‑69 (MQ=255)                        cAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCCCCGCAAACGCCGGAACAGAATGc                                         >  1:2515735/1‑69 (MQ=255)                        cAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCCCCGCAAACGCCGGAACAGAATGc                                         >  1:419832/1‑69 (MQ=255)                        cAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCCCCGCAAACGCCGGAACAGAATGc                                         >  1:1024167/1‑69 (MQ=255)                        cAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCCCCGCAAACGCCGGAACAGAATGc                                         >  1:2338703/1‑69 (MQ=255)                        cAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCCCCGCAAACGCCGGAACAGAATGc                                         >  1:3027906/1‑69 (MQ=255)                        cAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCCCCGCAAACGCCGGAACAGAATGc                                         >  1:469646/1‑69 (MQ=255)                        cAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCCCCGCAAACGCCGGAACAGAATGc                                         >  1:2288930/1‑69 (MQ=255)                        cAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCCCCGCAAACGCCGGAACAGAATGc                                         >  1:1366365/1‑69 (MQ=255)                        cAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCCCCGCAAACGCCGGAACAGAATGc                                         >  1:98651/1‑69 (MQ=255)                        cAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCCCCGCAAACGCCGGAACAGAATGc                                         >  1:2050175/1‑69 (MQ=255)                        cAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCCCCGCAAACGCCGGAACAGAATGc                                         >  1:1263584/1‑69 (MQ=255)                        cAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCCCCGCAAACGCCGGAACAGAATGc                                         >  1:144910/1‑69 (MQ=255)                        cAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCCCCGCAAACGCCGGAACAGAATGc                                         >  1:961713/1‑69 (MQ=255)                        cAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCCCCGCAAACGCCGGAACAGAATGc                                         >  1:1787474/1‑69 (MQ=255)                           tGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCCCCGCAAACGCCGGAACa                                               >  1:3053260/1‑60 (MQ=255)                           tGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCCCCGCAAACGCCGGAACAGAATGCTa                                       >  1:2940464/1‑68 (MQ=255)                           tGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCCCCGCAAACGCCGGAACAGAATGCTa                                       >  1:1905289/1‑68 (MQ=255)                           tGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCCCCGCAAACGCCGGAACAGAATGCTa                                       >  1:82826/1‑68 (MQ=255)                           tGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCCCCGCAAACGCCGGAACAGAATGCTAt                                      >  1:676253/1‑69 (MQ=255)                           tGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCCCCGCAAACGCCGGAACAGAATGCTAt                                      >  1:1535361/1‑69 (MQ=255)                           tGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCCCCGCAAACGCCGGAACAGAATGCTAt                                      >  1:3313490/1‑69 (MQ=255)                           tGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCCCCGCAAACGCCGGAACAGAATGCTAt                                      >  1:680124/1‑69 (MQ=255)                           tGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCCCCGCAAACGCCGGAACAGAATGCTAt                                      >  1:563282/1‑69 (MQ=255)                              gcTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCgg                                                   <  1:1758765/53‑1 (MQ=255)                                    cAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATg                              >  1:708371/1‑68 (MQ=255)                                    cAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATg                              >  1:349579/1‑68 (MQ=255)                                    cAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATg                              >  1:2892441/1‑68 (MQ=255)                                    cAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATg                              >  1:2201176/1‑68 (MQ=255)                                    cAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATg                              >  1:523073/1‑68 (MQ=255)                                    cAAACCGGTTTCCACTAAT‑CGGGAAACGTTGCACCGCAAACGCCGGAACAGGATGCTATTGCTGATGGt                            <  1:2302589/69‑1 (MQ=255)                                     aaaCCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGt                            <  1:52170/69‑1 (MQ=255)                                     aaaCCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGt                            <  1:566878/69‑1 (MQ=255)                                     aaaCCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGt                            <  1:320751/69‑1 (MQ=255)                                     aaaCCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGt                            <  1:137527/69‑1 (MQ=255)                                     aaaCCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGt                            <  1:3001028/69‑1 (MQ=255)                                     aaaCCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGt                            <  1:2644089/69‑1 (MQ=255)                                     aaaCCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGt                            <  1:2760258/69‑1 (MQ=255)                                     aaaCCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGt                            <  1:1484859/69‑1 (MQ=255)                                              tCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAAGAGAATGCTATTGCTGATGGt                            <  1:1554613/60‑1 (MQ=255)                                              tCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGt                            <  1:404562/60‑1 (MQ=255)                                              tCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGt                            <  1:239645/60‑1 (MQ=255)                                              tCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGt                            <  1:3216788/60‑1 (MQ=255)                                              tCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGt                            <  1:3258445/60‑1 (MQ=255)                                              tCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGt                            <  1:3267301/60‑1 (MQ=255)                                              tCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGt                            <  1:2868455/60‑1 (MQ=255)                                              tCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGt                            <  1:536475/60‑1 (MQ=255)                                              tCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGt                            <  1:3408997/60‑1 (MQ=255)                                              tCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGt                            <  1:2239427/60‑1 (MQ=255)                                              tCCACTAATCCGGAAAACGTTACACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGt                            <  1:2533403/60‑1 (MQ=255)                                                   tAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAg              <  1:1272861/69‑1 (MQ=255)                                                   tAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAg              <  1:2320087/69‑1 (MQ=255)                                                   tAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAg              <  1:2989968/69‑1 (MQ=255)                                                   tAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAg              <  1:1917315/69‑1 (MQ=255)                                                   tAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAg              <  1:1726208/69‑1 (MQ=255)                                                   tAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAg              <  1:394860/69‑1 (MQ=255)                                                   tAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAg              <  1:2396628/69‑1 (MQ=255)                                                   tAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAg              <  1:2586287/69‑1 (MQ=255)                                                   tAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAg              <  1:3053579/69‑1 (MQ=255)                                                   tAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAg              <  1:1144973/69‑1 (MQ=255)                                                   tAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAg              <  1:2902996/69‑1 (MQ=255)                                                    aaTCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAg              <  1:783211/68‑1 (MQ=255)                                                    aaTCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGc             >  1:1896333/1‑69 (MQ=255)                                                       ccGGAAAACGTTGCCCCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCg          <  1:1437844/69‑1 (MQ=255)                                                       ccGGAAAACGTTGCCCCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCg          <  1:1098200/69‑1 (MQ=255)                                                       ccGGAAAACGTTGCCCCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCg          <  1:1589349/69‑1 (MQ=255)                                                       ccGGAAAACGTTGCCCCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCg          <  1:1798556/69‑1 (MQ=255)                                                       ccGGAAAACGTTGCCCCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCg          <  1:616892/69‑1 (MQ=255)                                                       ccGGAAAACGTTGCCCCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCg          <  1:3435600/69‑1 (MQ=255)                                                       ccGGAAAACGTTGCCCCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATCCGCCTGATAGCCCg          <  1:2225641/69‑1 (MQ=255)                                                        cGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGt                            >  1:1700393/1‑50 (MQ=255)                                                        cGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGt                            >  1:1881497/1‑50 (MQ=255)                                                        cGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGt                            >  1:2722169/1‑50 (MQ=255)                                                               cGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGt  >  1:3325356/1‑69 (MQ=255)                                                               cGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGt  >  1:951070/1‑69 (MQ=255)                                                                  tGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGt                            <  1:3162043/40‑1 (MQ=255)                                                                     |                                                               GAGCAGGATCTGCTGGCACTGAACAATGCGCTGAGCAAACCGGTTTCCACTAATCCGGAAAACGTTGCACCGCAAACGCCGGAACAGAATGCTATTGCTGATGGTTATGCGCCTGATAGCCCGGCACCAGT  >  minE/2034853‑2034983

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