Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,993,837	T→C	46.9%	I377T (ATT→ACT)	nupG →	nucleoside transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,993,837	0	T	C	46.9%	1.6 / 54.4	64	I377T (ATT→ACT)	nupG	nucleoside transporter
Reads supporting (aligned to +/- strand):  ref base T (9/25);  new base C (15/15);  total (24/40)
Fisher's exact test for biased strand distribution p-value = 7.15e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCTTCGGCTGTATCCTCGGCGGCATCGTGAGCGGTAAAGTTGTTGAGATGTACACCCAAAACGGCATTACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTT  >  minE/1993770‑1993900                                                                    |                                                                ggCTTCGGCTGTATCCTCGGCGGCATCGTGAGCGGTAAAGTTGTTGAGATGTACACCCAAAACGGCAtt                                                                <  1:2800700/69‑1 (MQ=255)     tCGGCTGTATCCTCGGCGGCATCGTGAGCGGTAAAGTTGTTGAGATGTACACCCAAAACGGCATTAcc                                                             <  1:2212684/68‑1 (MQ=255)             aTCCTCGGCGGCATCGTGAGCGGTAAAGTGGTTGAGATGTACACCTAAAACGGCACTACCGACTGGCAg                                                    >  1:374088/1‑69 (MQ=255)             aTCCTCGGCGGCATCGTGAGCGGTAAAGTGGTTGAGATGTACACCCAAAACGGCACTAcc                                                             >  1:2666023/1‑60 (MQ=255)             aTCCTCGGCGGCATCGTGAGCGGTAAAGTGGTTGAGATGTACACCCAAAACGGCACTACCGACTGGCAg                                                    >  1:1426522/1‑69 (MQ=255)             aTCCTCGGCGGCATCGTGAGCGGTAAAGTGGTTGAGATGTACACCCAAAACGGCACTACCGACTGGCAg                                                    >  1:2006810/1‑69 (MQ=255)             aTCCTCGGCGGCATCGTGAGCGGTAAAGTGGTTGAGATGTACACCCAAAACGGCACTACCGACTGGCAg                                                    >  1:2189478/1‑69 (MQ=255)             aTCCTCGGCGGCATCGTGAGCGGTAAAGTGGTTGAGATGTACACCCAAAACGGCACTACCGACTGGCAg                                                    >  1:1772706/1‑69 (MQ=255)             aTCCTCGGCGGCATCGTGAGCGGTAAAGTGGTTGAGATGTACACCCAAAACGGCACTACCGACTGGCAg                                                    >  1:2549695/1‑69 (MQ=255)             aTCCTCGGCGGCATCGTGAGCGGTAAAGTGGTTGAGATGTACACCCAAAACGGCACTACCGACTGGCAg                                                    >  1:2601300/1‑69 (MQ=255)             aTCCTCGGCGGCATCGTGAGCGGTAAAGTGGTTGAGATGTACACCCAAAACGGCACTACCGACTGGCAg                                                    >  1:1436357/1‑69 (MQ=255)             aTCCTCGGCGGCATCGTGAGCGGTAAAGTGGTTGAGATGTACACCCAAAACGGCACTACCGACTGGCAg                                                    >  1:2773734/1‑69 (MQ=255)             aTCCTCGGCGGCATCGTGAGCGGTAAAGTGGTTGAGATGTACACCCAAAACGGCACTACCGACTGGCAg                                                    >  1:2922787/1‑69 (MQ=255)             aTCCTCGGCGGCATCGTGAGCGGTAAAGTGGTTGAGATGTACACCCAAAACGGCACTACCGACTGGCAg                                                    >  1:1173929/1‑69 (MQ=255)             aTCCTCGGCGGCATCGTGAGCGGTAAAGTGGTTGAGATGTACACCCAAAACGGCACTACCGACTGGCAg                                                    >  1:3386596/1‑69 (MQ=255)             aTCCTCGGCGGCATCGTGAGCGGTAAAGTGGTTGAGATGTACACCCAAAACGGCACTACCGACTGGCAg                                                    >  1:687349/1‑69 (MQ=255)             aTCCTCGGCGGCATCGTGAGCGGTAAAGTGGTTGAGATGTACACCCAAAACGGCACTACCGACTGGCAg                                                    >  1:964028/1‑69 (MQ=255)              tCCTCGGCGGCATCGTGAGCGGTAAAGTTGTTGAGATGTACACCCAAAACGGCATTACCGACTGGCAGa                                                   >  1:2353839/1‑69 (MQ=255)                  cggcggCATCGTGAGCGGTAAAGTTGTTGAGATGTACACCCAAAACGGCATTACCGACTGGCTGACCg                                                <  1:2339676/68‑1 (MQ=255)                            gTGAGCGGTAAAGTTGTTGAGATGTACACCCAAAACGGCATTACCGACTGGCAg                                                    >  1:642237/1‑54 (MQ=255)                            gTGAGCGGTAAAGTTGTTGAGATGTACACCCAAAACGGCATTACCGACTGGCAg                                                    >  1:1970069/1‑54 (MQ=255)                            gTGAGCGGTAAAGTTGTTGAGATGTACACCCAAAACGGCATTACCGACTGGCAg                                                    >  1:96325/1‑54 (MQ=255)                            gTGAGCGGTAAAGTTGTTGAGATGTACACCCAAAACGGCATTACCGACTGGCAg                                                    >  1:2039189/1‑54 (MQ=255)                            gTGAGCGGTAAAGTTGTTGAGATGTACACCCAAAACGGCATTACCGACTGGCAg                                                    >  1:1663847/1‑54 (MQ=255)                            gTGAGCGGTAAAGTTGTTGAGATGTACACCCAAAACGGCATTACCGACTGGCAg                                                    >  1:3297661/1‑54 (MQ=255)                            gTGAGCGGTAAAGTTGTTGAGATGTACACCCAAAACGGCATTACCGACTGGCAg                                                    >  1:2189329/1‑54 (MQ=255)                              gAGCGGTAAAGTTGTTGAGATGTACACCCAAAACGGCATTACCGACTGGCAGACCGTATGGTTGAtttt                                   >  1:2051656/1‑69 (MQ=255)                                 cGGTAAAGTGGTTGAGATGTACACCCAAAACGGCACTACCGACTGGCAGACCGTATGGCTGATTTTct                                 <  1:3008174/68‑2 (MQ=255)                                 cGGTAAAGTGGTTGAGATGTACACCCAAAACGGCACTACCGACTGGCAGACCGTATGGCTGATTTTCg                                 <  1:2810172/68‑1 (MQ=255)                                 cGGTAAAGTGGTTGAGATGTACACCCAAAACGGCACTACCGACTGGCAGACCGTATGGCTGATTTTCg                                 <  1:117217/68‑1 (MQ=255)                                 cGGTAAAGTGGTTGAGATGTACACCCAAAACGGCACTACCGACTGGCAGACCGTATGGCTGATTTTCg                                 <  1:1421420/68‑1 (MQ=255)                                 cGGTAAAGTGGTTGAGATGTACACCCAAAACGGCACTACCGACTGGCAGACCGTATGGCTGATTTTCg                                 <  1:2843570/68‑1 (MQ=255)                                 cGGTAAAGTGGTTGAGATGTACACCCAAAACGGCACTACCGACTGGCAGACCGTATGGCTGATTTTCg                                 <  1:3126880/68‑1 (MQ=255)                                 cGGTAAAGTGGTTGAGATGTACACCCAAAACGGCACTACCGACTGGCAGACCGTATGGCTGATTTTCg                                 <  1:1001325/68‑1 (MQ=255)                                 cGGTAAAGTGGTTGAGATGTACACCCAAAACGGCACTACCGACTGGCAGACCGTATGGCTGATTTTCg                                 <  1:2578786/68‑1 (MQ=255)                                 cGGTAAAGTGGTTGAGATGTACACCCAAAACGGCACTACCGACTGGCAGACCGTATGGCTGATTTTCg                                 <  1:1762259/68‑1 (MQ=255)                                 cGGTAAAGTGGTTGAGATGTACACCCAAAACGGCACTACCGACTGGCAGACCGTATGGCTGATTTTCg                                 <  1:2210093/68‑1 (MQ=255)                                 cGGTAAAGTGGTTGAGATGTACACCCAAAACGGCACTACCGACTGGCAGACCGTATGGCTGATTTTCg                                 <  1:1849119/68‑1 (MQ=255)                                 cGGTAAAGTGGTTGAGATGTACACCCAAAACGGCACTACCGACTGGCAGACCGTATGGCTGATTTTCg                                 <  1:841373/68‑1 (MQ=255)                                 cGGTAAAGTGGTTGAGATGTACACCCAAAACGGCACTACCGACTGGCAGACCGTATGGCTGATTTTCg                                 <  1:948939/68‑1 (MQ=255)                                  ggTAAAGTGGTTGAGATGTACACCCAAAACGGCACTACCGACTGGCAGACCGTATGGCTGATTTTCg                                 <  1:1930717/67‑1 (MQ=255)                                        gttgttGAGATGTACACCCAAAACGGCATTACCGACTGGCAGATCGTATGGTTGATTTTCGCTGGTTa                          <  1:70642/68‑1 (MQ=255)                                        gttgttGAGATGTACACCCAAAACGGCATTACCGACTGGCAGACCGTATGGTTGATTTTCTCTGGTTa                          <  1:1123483/68‑1 (MQ=255)                                        gttgttGAGATGTACACCCAAAACGGCATTACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTa                          <  1:49341/68‑1 (MQ=255)                                        gttgttGAGATGTACACCCAAAACGGCATTACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTa                          <  1:1943535/68‑1 (MQ=255)                                        gttgttGAGATGTACACCCAAAACGGCATTACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTa                          <  1:809679/68‑1 (MQ=255)                                        gttgttGAGATGTACACCCAAAACGGCATTACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTa                          <  1:1203594/68‑1 (MQ=255)                                        gttgttGAGATGTACACCCAAAACGGCATTACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTa                          <  1:2916939/68‑1 (MQ=255)                                        gttgttGAGATGTACACCCAAAACGGCATTACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTa                          <  1:1301331/68‑1 (MQ=255)                                        gttgttGAGATGTACACCCAAAACGGCATTACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTa                          <  1:2947070/68‑1 (MQ=255)                                        gttgttGAGATGTACACCCAAAACGGCATTACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTa                          <  1:1862185/68‑1 (MQ=255)                                        gttgttGAGATGTACACCCAAAACGGCATTACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTa                          <  1:1723797/68‑1 (MQ=255)                                        gttgttGAGATGTACACCCAAAACGGCATTACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTa                          <  1:2591891/68‑1 (MQ=255)                                        gttgttGAGATGTACACCCAAAACGGCATTACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTa                          <  1:1815639/68‑1 (MQ=255)                                        gttgttGAGATGTACACCCAAAACGGCATTACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTa                          <  1:2270046/68‑1 (MQ=255)                                                    tACACCCAAAACGGCATTACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTg             <  1:3114249/69‑1 (MQ=255)                                                    tACACCCAAAACGGCACTACCGACTGGCAGACCGTATGGCTGATTTTCg                                 <  1:3347994/49‑1 (MQ=39)                                                               cGGCATTACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGtt  <  1:1114555/69‑1 (MQ=255)                                                               cGGCATTACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGtt  <  1:2079528/69‑1 (MQ=255)                                                               cGGCATTACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGtt  <  1:2122383/69‑1 (MQ=255)                                                               cGGCATTACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGtt  <  1:2569538/69‑1 (MQ=255)                                                               cGGCATTACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGtt  <  1:591047/69‑1 (MQ=255)                                                               cGGCATTACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGtt  <  1:510875/69‑1 (MQ=255)                                                               cGGCATTACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGtt  <  1:1145404/69‑1 (MQ=255)                                                               cGGCATTACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGtt  <  1:3425593/69‑1 (MQ=255)                                                                    |                                                                GGCTTCGGCTGTATCCTCGGCGGCATCGTGAGCGGTAAAGTTGTTGAGATGTACACCCAAAACGGCATTACCGACTGGCAGACCGTATGGTTGATTTTCGCTGGTTACTCCGTGGTTCTGGCCTTCGCGTT  >  minE/1993770‑1993900

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