Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,993,886	T→G	27.3%	V393V (GTT→GTG)	nupG →	nucleoside transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,993,886	0	T	G	27.3%	112.8 / 60.1	99	V393V (GTT→GTG)	nupG	nucleoside transporter
Reads supporting (aligned to +/- strand):  ref base T (31/41);  new base G (9/18);  total (40/59)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.17e-01

TACACCCAAAACGGCATTACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACA  >  minE/1993821‑1993948                                                                  |                                                               tACACCCAAAACGGCATTACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTg                                                             <  1:3114249/69‑1 (MQ=255)            cGGCATTACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGtt                                                  <  1:2079528/69‑1 (MQ=255)            cGGCATTACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGtt                                                  <  1:3425593/69‑1 (MQ=255)            cGGCATTACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGtt                                                  <  1:2122383/69‑1 (MQ=255)            cGGCATTACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGtt                                                  <  1:1114555/69‑1 (MQ=255)            cGGCATTACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGtt                                                  <  1:1145404/69‑1 (MQ=255)            cGGCATTACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGtt                                                  <  1:591047/69‑1 (MQ=255)            cGGCATTACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGtt                                                  <  1:510875/69‑1 (MQ=255)            cGGCATTACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGtt                                                  <  1:2569538/69‑1 (MQ=255)                  tACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGc                                            <  1:2561100/69‑1 (MQ=255)                  tACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGc                                            <  1:1964163/69‑1 (MQ=255)                  tACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGc                                            <  1:163546/69‑1 (MQ=255)                  tACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGc                                            <  1:1572934/69‑1 (MQ=255)                  tACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGc                                            <  1:1505302/69‑1 (MQ=255)                  tACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGc                                            <  1:1329544/69‑1 (MQ=255)                  tACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGc                                            <  1:125088/69‑1 (MQ=255)                  tACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGc                                            <  1:3278749/69‑1 (MQ=255)                  tACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGc                                            <  1:3435609/69‑1 (MQ=255)                  tACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGc                                            <  1:969422/69‑1 (MQ=255)                  tACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGc                                            <  1:1110745/69‑1 (MQ=255)                   aCCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTc                                                       >  1:3319539/1‑57 (MQ=255)                   aCCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGc                                            >  1:968521/1‑68 (MQ=255)                   aCCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGc                                            >  1:203572/1‑68 (MQ=255)                   aCCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGc                                            >  1:636220/1‑68 (MQ=255)                   aCCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGc                                            >  1:2678430/1‑68 (MQ=255)                   aCCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCg                                           >  1:2220297/1‑69 (MQ=255)                   aCCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCg                                           >  1:316801/1‑69 (MQ=255)                   aCCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCg                                           >  1:505835/1‑69 (MQ=255)                   aCCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCg                                           >  1:2339064/1‑69 (MQ=255)                   aCCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCg                                           >  1:550052/1‑69 (MQ=255)                   aCCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCg                                           >  1:2370824/1‑69 (MQ=255)                   aCCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCg                                           >  1:3351821/1‑69 (MQ=255)                   aCCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCg                                           >  1:1726885/1‑69 (MQ=255)                   aCCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCg                                           >  1:1223400/1‑69 (MQ=255)                   aCCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCAAGGCg                                           >  1:2616601/1‑69 (MQ=255)                   aCCGACAGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCg                                           >  1:437989/1‑69 (MQ=255)                     cGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGa                                          <  1:659117/68‑1 (MQ=255)                     cGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGa                                          <  1:3241198/68‑1 (MQ=255)                     cGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGa                                          <  1:247968/68‑1 (MQ=255)                     cGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGa                                          <  1:585759/68‑1 (MQ=255)                     cGACTGGCAGACCGAATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGa                                          <  1:1042241/68‑1 (MQ=255)                       aCTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGa                                          <  1:1518537/66‑1 (MQ=255)                           gCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCa                                    >  1:945548/1‑68 (MQ=255)                           gCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCa                                    >  1:1471767/1‑68 (MQ=255)                                ccGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGtt                                                  <  1:2713962/49‑1 (MQ=255)                                ccGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGtt                                                  <  1:1236049/49‑1 (MQ=255)                                  gTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATaaa                            <  1:555342/69‑1 (MQ=255)                                  gTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATaaa                            <  1:1683367/69‑1 (MQ=255)                                  gTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATaaa                            <  1:1022601/69‑1 (MQ=255)                                     tGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATaaa                            >  1:293958/1‑66 (MQ=255)                                     tGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATaaa                            >  1:97790/1‑66 (MQ=255)                                     tGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATaaa                            >  1:1285135/1‑66 (MQ=255)                                     tGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATaaa                            >  1:1046906/1‑66 (MQ=255)                                     tGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATaaa                            >  1:2711786/1‑66 (MQ=255)                                     tGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGAGATGTTCAAATATaaa                            >  1:398600/1‑66 (MQ=255)                                       gTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGt                       >  1:3110273/1‑69 (MQ=255)                                       gTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGt                       >  1:1317112/1‑69 (MQ=255)                                       gTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGt                       >  1:1138224/1‑69 (MQ=255)                                       gTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGt                       >  1:287113/1‑69 (MQ=255)                                           aTTTTCGCAGGTTACTCCGTGGTGCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCgt                   >  1:1362355/1‑69 (MQ=255)                                           aTTTTCGCAGGTTACTCCGTGGTGCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCgt                   >  1:2544149/1‑69 (MQ=255)                                           aTTTTCGCAGGTTACTCCGTGGTGCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCgt                   >  1:955196/1‑69 (MQ=255)                                           aTTTTCGCAGGTTACTCCGTGGTGCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCgt                   >  1:2978857/1‑69 (MQ=255)                                           aTTTTCGCAGGTTACTCCGTGGTGCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCgt                   >  1:2363428/1‑69 (MQ=255)                                           aTTTTCGCAGGTTACTCCGTGGTGCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCgt                   >  1:1566930/1‑69 (MQ=255)                                           aTTTTCGCAGGTTACTCCGTGGTGCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCgt                   >  1:3253032/1‑69 (MQ=255)                                           aTTTTCGCAGGTTACTCCGTGGTGCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCgt                   >  1:1802241/1‑69 (MQ=255)                                           aTTTTCGCAGGTTACTCCGTGGTGCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCgt                   >  1:182629/1‑69 (MQ=255)                                             tttCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCgtgt                 >  1:1525866/1‑69 (MQ=255)                                             tttCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCgtgt                 >  1:2172897/1‑69 (MQ=255)                                                           ccGTGGTTCTGGCCTTCGCTTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGca      >  1:660876/1‑66 (MQ=255)                                                              tGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGcacaca  <  1:504346/67‑1 (MQ=255)                                                              tGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGcacaca  <  1:943482/67‑1 (MQ=255)                                                              tGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGcacaca  <  1:95427/67‑1 (MQ=255)                                                              tGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGcacaca  <  1:2458877/67‑1 (MQ=255)                                                              tGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGcacaca  <  1:2971501/67‑1 (MQ=255)                                                              tGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGcacaca  <  1:114724/67‑1 (MQ=255)                                                              tGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGcacaca  <  1:2628125/67‑1 (MQ=255)                                                              tGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGcacaca  <  1:1469241/67‑1 (MQ=255)                                                              tGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGcacaca  <  1:234257/67‑1 (MQ=255)                                                              tGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGcacaca  <  1:2065492/67‑1 (MQ=255)                                                              tGGTGCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCgtg                  <  1:990031/51‑1 (MQ=255)                                                              tGGTGCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCgtg                  <  1:3262794/51‑1 (MQ=255)                                                              tGGTGCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCgtg                  <  1:108071/51‑1 (MQ=255)                                                              tGGTGCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCgtg                  <  1:1377617/51‑1 (MQ=255)                                                              tGGTGCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCgtg                  <  1:1432669/51‑1 (MQ=255)                                                              tGGTGCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCgtg                  <  1:185438/51‑1 (MQ=255)                                                              tGGTGCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCgtg                  <  1:1926619/51‑1 (MQ=255)                                                              tGGTGCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCgtg                  <  1:777249/51‑1 (MQ=255)                                                              tGGTGCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCgtg                  <  1:2010488/51‑1 (MQ=255)                                                              tGGTGCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCgtg                  <  1:2019036/51‑1 (MQ=255)                                                              tGGTGCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCgtg                  <  1:2134526/51‑1 (MQ=255)                                                              tGGTGCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCgtg                  <  1:2269514/51‑1 (MQ=255)                                                              tGGTGCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCgtg                  <  1:1020630/51‑1 (MQ=255)                                                              tGGTGCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCgtg                  <  1:2681399/51‑1 (MQ=255)                                                              tGGTGCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCgtg                  <  1:3047770/51‑1 (MQ=255)                                                              tGGTGCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCgtg                  <  1:3434929/51‑1 (MQ=255)                                                               ggTGCTGGCCTTCGCGTTCATGTCGATGTTCAAATATAAACACGTTCgtg                  <  1:2357556/50‑1 (MQ=255)                                                               ggTGCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCgtg                  <  1:28553/50‑1 (MQ=255)                                                                  |                                                               TACACCCAAAACGGCATTACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTTCATGGCGATGTTCAAATATAAACACGTTCGTGTCCCGACAGGCACACA  >  minE/1993821‑1993948

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