Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,992,986	A→G	33.3%	E93E (GAA→GAG)	nupG →	nucleoside transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,992,986	0	A	G	33.3%	52.8 / 38.1	72	E93E (GAA→GAG)	nupG	nucleoside transporter
Reads supporting (aligned to +/- strand):  ref base A (24/24);  new base G (6/18);  total (30/42)
Fisher's exact test for biased strand distribution p-value = 4.78e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

TATGCCATTTGCCACACCATTGGCGCTATCACGCTGTTCATGGCGGCACAGGTCACGACACCGGAAGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCTTATATGCCAACGCTTGGGTTAATCAACACC  >  minE/1992921‑1993052                                                                  |                                                                   tatGCCATTTGCCACACCATTGGCGCTATCACGCTGTTCATGGCGGCACAGGTCACGACACCGGAAGCg                                                                 >  1:2749571/1‑69 (MQ=255)   tGCCATTTGCCACACCATTGGCGCTATCACGCTGTTCATGGCGGCACAGTTCACGACACCGGAAGCGAt                                                               >  1:2504305/1‑69 (MQ=255)   tGCCATTTGCCACACCATTGGCGCTATCACGCTGTTCATGGCGGCACAGGTCACGACACCGGAAGCGAt                                                               >  1:1079450/1‑69 (MQ=255)   tGCCATTTGCCACACCATTGGCGCTATCACGCTGTTCATGGCGGCACAGGTCACGACACCGGAAGCGAt                                                               >  1:2173140/1‑69 (MQ=255)   tGCCATTTGCCACACCATTGGCGCTATCACGCTGTTCATGGCGGCACAGGTCACGACACCGGAAGCGAt                                                               >  1:247826/1‑69 (MQ=255)   tGCCATTTGCCACACCATTGGCGCTATCACGCTGTTCATGGCGGCACAGGTCACGACACCGGAAGCGAt                                                               >  1:314084/1‑69 (MQ=255)   tGCCATTTGCCACACCATTGGCGCTATCACGCTGTTCATGGCGGCACAGGTCACGACACCGGAAGCGAt                                                               >  1:1691092/1‑69 (MQ=255)   tGCCATTTGCCACACCATTGGCCCTATCACGCAGTTCATGGCGGCACAGGTCACGACACCGTAAGCGAt                                                               >  1:1349554/1‑69 (MQ=255)          tGCCACACCATTGGCGCTATCACGCTGTTCATGGCGGCACAGGTCACGACACCGGAAGCGATGTTCCtt                                                        >  1:2095255/1‑69 (MQ=255)          tGCCACACCATTGGCGCTATCACGCTGTTCATGGCGGCACAGGTCACGACACCGGAAGCGATGTTCCtt                                                        >  1:2193784/1‑69 (MQ=255)          tGCCACACCATTGGCGCTATCACGCTGTTCATGGCGGCACAGGTCACGACACCGGAAGCGATGTTCCtt                                                        >  1:2835032/1‑69 (MQ=255)             cacaCCATTGGCGCTATCACGCTGTTCATGGCGGCACAGGTCACGACACCGGAAGCGATGTTCCTtgt                                                      <  1:131469/68‑1 (MQ=255)                      gtcgcTATCACGCTGTTCATGGCGGCACAGGTCACGACACCGGAAGCGATGTTCCTTGTGATATTGAtt                                            <  1:1008738/67‑1 (MQ=255)                      ggCGCTATCACGCTGTTCATGGCGGCACAGGTCACGACACCGGAAGCGATGTTCCTTGTGATATTGAtt                                            <  1:1901855/69‑1 (MQ=255)                            aTCACGCTGTTCATGGCGGCAGAGGTCACTACGCCGGAGGCGATGTTCCTTGTGatat                                                 <  1:3014382/58‑1 (MQ=255)                            aTCACGCTGTTCATGGCGGCAGAGGTCACTACGCCGGAGGCGATGTTCCTTGTGatat                                                 <  1:898578/58‑1 (MQ=255)                            aTCACGCTGTTCATGGCGGCAGAGGTCACTACGCCGGAGGCGATGTTCCTTGTGatat                                                 <  1:869498/58‑1 (MQ=255)                            aTCACGCTGTTCATGGCGGCAGAGGTCACTACGCCGGAGGCGATGTTCCTTGTGatat                                                 <  1:2671378/58‑1 (MQ=255)                            aTCACGCTGTTCATGGCGGCAGAGGTCACTACGCCGGAGGCGATGTTCCTTGTGatat                                                 <  1:273500/58‑1 (MQ=255)                            aTCACGCTGTTCATGGCGGCAGAGGTCACTACGCCGGAGGCGATGTTCCTTGTGatat                                                 <  1:672716/58‑1 (MQ=255)                            aTCACGCTGTTCATGGCGGCAGAGGTCACTACGCCGGAGGCGATGTTCCTTGTGatat                                                 <  1:3042181/58‑1 (MQ=255)                            aTCACGCTGTTCATGGCGGCAGAGGTCACTACGCCGGAGGCGATGTTCCTTGTGatat                                                 <  1:3078151/58‑1 (MQ=255)                            aTCACGCTGTTCATGGCGGCAGAGGTCACTACGCCGGAGGCGATGTTCCTTGTGatat                                                 <  1:3106324/58‑1 (MQ=255)                            aTCACGCTGTTCATGGCGGCAGAGGTCACTACGCCGGAGGCGATGTTCCTTGTGatat                                                 <  1:666582/58‑1 (MQ=255)                            aTCACGCTGTTCATGGCGGCAGAGGTCACTACGCCGGAGGCGATGTTCCTTGTGatat                                                 <  1:3317162/58‑1 (MQ=255)                            aTCACGCTGTTCATGGCGGCAGAGGTCACTACGCCGGAGGCGATGTTCCTTGTGatat                                                 <  1:3148631/58‑1 (MQ=255)                            aTCACGCTGTTCATGGCGGCAGAGGTCACTACGCCGGAGGCGATGTTCCTTGTGatat                                                 <  1:3243267/58‑1 (MQ=255)                            aTCACGCTGTTCATGGCGGCAGAGGTCACTACGCCGGAGGCGATGTTCCTTGTGatat                                                 <  1:3274723/58‑1 (MQ=255)                            aTCACGCTGTTCATGGCGGCAGAGGTCACTACGCCGGAGGCGATGTTCCTTGTGatat                                                 <  1:1461256/58‑1 (MQ=255)                            aTCACGCTGTTCATGGCGGCAGAGGTCACTACGCCGGAGGCGATGTTCCTTGTGatat                                                 <  1:2084169/58‑1 (MQ=255)                            aTCACGCTGTTCATGGCGGCAGAGGTCACTACGCCGGAGGCGATGTTCCTTGTGatat                                                 <  1:1970123/58‑1 (MQ=255)                            aTCACGCTGTTCATGGCGGCAGAGGTCACTACGCCGGAGGCGATGTTCCTTGTGatat                                                 <  1:182476/58‑1 (MQ=255)                                       cATGGCGGCACAGGTCACGACACCGGAAGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCTtat                             <  1:859169/67‑1 (MQ=255)                                       cATGGCGGCACAGGTCACGACACCGGAAGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCTtat                             <  1:1914254/67‑1 (MQ=255)                                       cATGGCGGCACAGGTCACGACACCGGAAGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCTtat                             <  1:1180710/67‑1 (MQ=255)                                       cATGGCGGCACAGGTCACGACACCGGAAGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCTtat                             <  1:2903605/67‑1 (MQ=255)                                       cATGGCGGCACAGGTCACGACACCGGAAGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCTTatat                           >  1:2712870/1‑69 (MQ=255)                                       cATGGCGGCACAGGTCACGACACCGGAAGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCTTatat                           >  1:2843502/1‑69 (MQ=255)                                           gcggcACAGGTCACGACACCGGAAGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCTta                              >  1:3376850/1‑62 (MQ=255)                                           gcggcACAGGTCACGACACCGGAAGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCTta                              >  1:2039458/1‑62 (MQ=255)                                           gcggcACAGGTCACGACACCGGAAGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCTTATATGCCa                       >  1:360843/1‑69 (MQ=255)                                                  aGGTCACTACGCCGGAGGCGATGTTCCTTGTGATATTGATTAACTCg                                      >  1:392784/1‑47 (MQ=255)                                                  aGGTCACTACGCCGGAGGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCTTATATGCCAACGCtt                 >  1:2328050/1‑68 (MQ=255)                                                  aGGTCACTACGCCGGAGGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCTTATATGCCAACGCtt                 >  1:349294/1‑68 (MQ=255)                                                  aGGTCACTACGCCGGAGGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCTTATATGCCAACGCtt                 >  1:2820826/1‑68 (MQ=255)                                                  aGGTCACTACGCCGGAGGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCTTATATGCCAACGCTTg                >  1:1835383/1‑69 (MQ=255)                                                  aGGTCACTACGCCGGAGGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCTTATATGCCAACGCTTg                >  1:2337267/1‑69 (MQ=255)                                                   ggTCACGACACCGGAAGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCt                                <  1:606347/52‑1 (MQ=255)                                                   ggTCACGACACCGGAAGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCt                                <  1:1305712/52‑1 (MQ=255)                                                   ggTCACGACACCGGAAGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCt                                <  1:351427/52‑1 (MQ=255)                                                   ggTCACGACACCGGAAGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCt                                <  1:774260/52‑1 (MQ=255)                                                   ggTCACGACACCGGAAGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCt                                <  1:3384343/52‑1 (MQ=255)                                                   ggTCACGACACCGGAAGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCt                                <  1:1145603/52‑1 (MQ=255)                                                   ggTCACGACACCGGAAGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCt                                <  1:1983896/52‑1 (MQ=255)                                                   ggTCACGACACCGGAAGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCt                                <  1:2725074/52‑1 (MQ=255)                                                   ggTCACGACACCGGAAGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCt                                <  1:2771771/52‑1 (MQ=255)                                                   ggTCACGACACCGGAAGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCt                                <  1:2910011/52‑1 (MQ=255)                                                    gTCACGACAGCGGAAGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCTTATATGCCAACg                    >  1:1340033/1‑63 (MQ=255)                                                    gTCACGACACCGGAAGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCTTATATGCCAACg                    >  1:1179107/1‑63 (MQ=255)                                                    gTCACGACACCGGAAGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCTTATATGCCAACg                    >  1:1754046/1‑63 (MQ=255)                                                    gTCACGACACCGGAAGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCTTATATGCCAACg                    >  1:1947414/1‑63 (MQ=255)                                                    gTCACGACACCGGAAGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCTTATATGCCAACg                    >  1:3015991/1‑63 (MQ=255)                                                     tCACGACACCGGAAGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCTTATATGCCAACGCTTggg              >  1:2944840/1‑68 (MQ=255)                                                     tCACGACACCGGAAGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCTTATATGCCAACGCTTggg              >  1:614999/1‑68 (MQ=255)                                                     tCACGACACCGGAAGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCTTATATGCCAACGCTTggg              >  1:349222/1‑68 (MQ=255)                                                          acacCGGAAGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCTTATATGCCAACGCTTGGGTTAATc        <  1:1311284/69‑1 (MQ=255)                                                                gAAGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCTTATATGCCAACGCTTGGGTTAATCAACAcc  <  1:1866738/69‑1 (MQ=255)                                                                gAAGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCTTATATGCCAACGCTTGGGTTAATCAACAcc  <  1:1582132/69‑1 (MQ=255)                                                                gAAGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCTTATATGCCAACGCTTGGGTTAATCAACAcc  <  1:2822215/69‑1 (MQ=255)                                                                gAAGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCTTATATGCCAACGCTTGGGTTAATCAACAcc  <  1:2989752/69‑1 (MQ=255)                                                                gAAGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCTTATATGCCAACGCTTGGGTTAATCAACAcc  <  1:1815406/69‑1 (MQ=255)                                                                gAAGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCTTATATGCCAACGCTTGGGTTAATCAACAcc  <  1:1709404/69‑1 (MQ=255)                                                                  |                                                                   TATGCCATTTGCCACACCATTGGCGCTATCACGCTGTTCATGGCGGCACAGGTCACGACACCGGAAGCGATGTTCCTTGTGATATTGATTAACTCGTTTGCTTATATGCCAACGCTTGGGTTAATCAACACC  >  minE/1992921‑1993052

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑