Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,993,595	G→T	31.3%	A296A (GCG→GCT)	nupG →	nucleoside transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,993,595	0	G	T	31.3%	102.3 / 51.0	83	A296A (GCG→GCT)	nupG	nucleoside transporter
Reads supporting (aligned to +/- strand):  ref base G (31/26);  new base T (18/8);  total (49/34)
Fisher's exact test for biased strand distribution p-value = 2.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.01e-01

CTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTT  >  minE/1993529‑1993658                                                                   |                                                                cttAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCt                                                               >  1:2879858/1‑69 (MQ=255)         gCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCtt                                                       >  1:2653062/1‑69 (MQ=255)         gCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCtt                                                       >  1:3207673/1‑69 (MQ=255)         gCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCtt                                                       >  1:2755181/1‑69 (MQ=255)         gCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCtt                                                       >  1:2157944/1‑69 (MQ=255)         gCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCtt                                                       >  1:2877665/1‑69 (MQ=255)         gCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCtt                                                       >  1:573723/1‑69 (MQ=255)         gCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTCCGCTGTTTGCtt                                                       >  1:3172875/1‑69 (MQ=255)                 aTTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGAc                                               >  1:130811/1‑69 (MQ=255)                    aaGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCTCTGTTTGCTTACGGCGAccc                                             <  1:3113816/68‑1 (MQ=255)                    aaGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCTCTGTTTGCTTACGGCGAccc                                             <  1:3216943/68‑1 (MQ=255)                    aaGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCTCTGTTTGCTTACGGCGAccc                                             <  1:2948251/68‑1 (MQ=255)                    aaGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCTCTGTTTGCTTACGGCGAccc                                             <  1:2890695/68‑1 (MQ=255)                    aaGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCTCTGTTTGCTTACGGCGAccc                                             <  1:1499210/68‑1 (MQ=255)                    aaGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCTCTGTTTGCTTACGGCGAccc                                             <  1:2071815/68‑1 (MQ=255)                    aaGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCTCTGTTTGCTTACGGCGAccc                                             <  1:2108466/68‑1 (MQ=255)                    aaGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCTCTGTTTGCTTACGGCGAccc                                             <  1:2270089/68‑1 (MQ=255)                    aaGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGAccc                                             <  1:494219/68‑1 (MQ=255)                    aaGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGAccc                                             <  1:1538549/68‑1 (MQ=255)                    aaGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGAccc                                             <  1:2664432/68‑1 (MQ=255)                       aaCGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCTCTGTTTGCTTACGGCGACCCGa                                           >  1:2520474/1‑67 (MQ=255)                       aaCGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCTCTGTTTGCTTACGGCGACCCGAc                                          >  1:2960307/1‑68 (MQ=255)                       aaCGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCTCTGTTTGCTTACGGCGACCCGACt                                         >  1:1815024/1‑69 (MQ=255)                       aaCGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCTCTGTTTGCTTACGGCGACCCGACt                                         >  1:2614559/1‑69 (MQ=255)                       aaCGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCTCTGTTTGCTTACGGCGACCCGACt                                         >  1:1080585/1‑69 (MQ=255)                       aaCGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCTCTGTTTGCTTACGGCGACCCGACt                                         >  1:938823/1‑69 (MQ=255)                       aaCGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCTCTGTTTGCTTACGGCGACCCGACt                                         >  1:2372091/1‑69 (MQ=255)                       aaCGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCTCTGTTTGCTTACGGCGACCCGACt                                         >  1:2290887/1‑69 (MQ=255)                       aaCGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCTCTGTTTGCTTACGGCGACCCGACt                                         >  1:2221009/1‑69 (MQ=255)                       aaCGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCTCTGTTTGCTTACGGCGACCCGACt                                         >  1:865972/1‑69 (MQ=255)                       aaCGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCTCTGTTTGCTTACGGCGACCCGACt                                         >  1:757861/1‑69 (MQ=255)                       aaCGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCTCTGTTTGCTTACGGCGACCCGACt                                         >  1:733526/1‑69 (MQ=255)                       aaCGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCTCTGTTTGCTTACGGCGACCCGACt                                         >  1:1319240/1‑69 (MQ=255)                       aaCGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCTCTGTTTGCTTACGGCGACCCGACt                                         >  1:1460265/1‑69 (MQ=255)                       aaCGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCTCTGTTTGCTTACGGCGACCCGACt                                         >  1:2041819/1‑69 (MQ=255)                       aaCGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCTCTGTTTGCTTACGGCGACCCGACt                                         >  1:1942547/1‑69 (MQ=255)                       aaCGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCTCTGTTTGCTTACGGCGACCCGACt                                         >  1:1513240/1‑69 (MQ=255)                       aaCGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCTCTGTTTGCTTACGGCGACCCGACt                                         >  1:1524686/1‑69 (MQ=255)                             atgatgATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTc                                   >  1:884100/1‑69 (MQ=255)                             atgatgATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTc                                   >  1:438684/1‑69 (MQ=255)                             atgatgATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTc                                   >  1:954252/1‑69 (MQ=255)                             atgatgATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTc                                   >  1:1439722/1‑69 (MQ=255)                             atgatgATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTc                                   >  1:217643/1‑69 (MQ=255)                             atgatgATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTc                                   >  1:2081299/1‑69 (MQ=255)                             atgatgATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTc                                   >  1:1877866/1‑69 (MQ=255)                             atgatgATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTc                                   >  1:1809477/1‑69 (MQ=255)                                     cAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGgtactgt                           >  1:165994/1‑69 (MQ=255)                                        tATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGgtactgtact                        <  1:2507217/69‑1 (MQ=255)                                        tATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGgtactgtact                        <  1:119960/69‑1 (MQ=255)                                        tATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGgtactgtact                        <  1:2439695/69‑1 (MQ=255)                                        tATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGgtactgtact                        <  1:1475777/69‑1 (MQ=255)                                        tATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGgtactgtact                        <  1:3256463/69‑1 (MQ=255)                                        tATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGgtactgtact                        <  1:1654624/69‑1 (MQ=255)                                         aTTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGgtactgtact                        <  1:1012746/68‑1 (MQ=255)                                          ttGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTActgc                      >  1:3048756/1‑69 (MQ=255)                                            gtgGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGAc                                          >  1:2141037/1‑47 (MQ=255)                                            gtgGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTActgct                     >  1:1291946/1‑68 (MQ=255)                                            gtgGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTActgct                     >  1:3201991/1‑68 (MQ=255)                                            gtgGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTActgct                     >  1:399262/1‑68 (MQ=255)                                            gtgGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTActgct                     >  1:2852231/1‑68 (MQ=255)                                            gtgGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTActgct                     >  1:2953259/1‑68 (MQ=255)                                               gCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGt                  <  1:2477907/68‑1 (MQ=255)                                                    gATCCTGGGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTc            <  1:843790/69‑1 (MQ=255)                                                    gATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTc            >  1:3257743/1‑69 (MQ=255)                                                    gATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTc            <  1:2473192/69‑1 (MQ=255)                                                    gATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTc            <  1:804328/69‑1 (MQ=255)                                                    gATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTc            >  1:1293701/1‑69 (MQ=255)                                                    gATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTc            >  1:1298152/1‑69 (MQ=255)                                                    gATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTc            <  1:574300/69‑1 (MQ=255)                                                    gATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTc            >  1:1864013/1‑69 (MQ=255)                                                    gATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTc            <  1:3127375/69‑1 (MQ=255)                                                    gATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTc            >  1:3104058/1‑69 (MQ=255)                                                    gATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGGTCGGTACTGTACTGCTGGTACTGTc            >  1:3139626/1‑69 (MQ=255)                                                              tttGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGtt  <  1:1624486/69‑1 (MQ=255)                                                              tttGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGtt  <  1:2403296/69‑1 (MQ=255)                                                              tttGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGtt  <  1:1452366/69‑1 (MQ=255)                                                              tttGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGtt  <  1:2636825/69‑1 (MQ=255)                                                              tttGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGtt  <  1:1758932/69‑1 (MQ=255)                                                              tttGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGtt  <  1:2987555/69‑1 (MQ=255)                                                              tttGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATCATCGtt  <  1:439532/69‑1 (MQ=255)                                                              tttGCGCTGTTTGCTTACGGCGACCCGACTACGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGtt  <  1:1181131/69‑1 (MQ=255)                                                               ttGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGtt  <  1:317698/68‑1 (MQ=255)                                                               ttGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGtt  <  1:30686/68‑1 (MQ=255)                                                                   |                                                                CTTAAGCCGCTACGGTATTAAGAACGTAATGATGATCAGTATTGTGGCGTGGATCCTGCGTTTTGCGCTGTTTGCTTACGGCGACCCGACTCCGTTCGGTACTGTACTGCTGGTACTGTCGATGATCGTT  >  minE/1993529‑1993658

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