Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,525,136	A→G	47.7%	G79G (GGT→GGC)	yfcJ ←	predicted transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,525,136	0	A	G	47.7%	‑1.9 / 56.9	65	G79G (GGT→GGC)	yfcJ	predicted transporter
Reads supporting (aligned to +/- strand):  ref base A (29/5);  new base G (19/12);  total (48/17)
Fisher's exact test for biased strand distribution p-value = 4.66e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.94e-01

AGCAACGCGCCGCCAGCCAGACCACAAGCTAACATTCCCTGAAGCGCCGAACGTTTTGCACCATATTGATCGGCCAGTCGCCCGGCGTAACCACGCGTCAGCACCGTAGCCAGAAACTGA  >  minE/1525077‑1525196                                                            |                                                             agcaACGCGCCGCCAGCCAGACCACAAGCTAACATTCCCTGAAGCGCCGAACGTTTTGCACCATATTGa                                                     >  1:310256/1‑69 (MQ=255) agcaACGCGCCGCCAGCCAGACCACAAGCTAACATTCCCTGAAGCGCCGAACGTTTTGCACCATATTGa                                                     >  1:88634/1‑69 (MQ=255)               cgccagACCACAAGCTAACATCCCCTGAAGCGCCGAACGTTTTGCGCCATATTGATCGGCCAGTCGccc                                       >  1:3287439/2‑69 (MQ=255)               cgccagACCACAAGCTAACATCCCCTGAAGCGCCGAACGTTTTGCGCCATATTGATCGGCCAGTCGccc                                       >  1:1175392/2‑69 (MQ=255)               cgccagACCACAAGCTAACATCCCCTGAAGCGCCGAACGTTTTGCGCCATATTGATCGGCCAGTCGccc                                       >  1:2201914/2‑69 (MQ=255)               cgccagACCACAAGCTAACATCCCCTGAAGCGCCGAACGTTTTGCGCCATATTGATCGGCCAGTCGccc                                       >  1:1005381/2‑69 (MQ=255)               cgccagACCACAAGCTAACATCCCCTGAAGCGCCGAACGTTTTGCGCCATATTGATCGGCCAGTCGccc                                       >  1:3066933/2‑69 (MQ=255)               cgccagACCACAAGCTAACATCCCCTGAAGCGCCGAACGTTTTGCGCCATATTGATCGGCCAGTCGccc                                       >  1:839512/2‑69 (MQ=255)               cgccagACCACAAGCTAACATCCCCTGAAGCGCCGAACGTTTTGCGCCATATTGATCGGCCAGTCGccc                                       >  1:1611927/2‑69 (MQ=255)               cgccagACCACAAGCTAACATCCCCTGAAGCGCCGAACGTTTTGCGCCATATTGATCGGCCAGTCGccc                                       >  1:178129/2‑69 (MQ=255)               cgccagACCACAAGCTAACATCCCCTGAAGCGCCGAACGTTTTGCGCCATATTGATCGGCCAGTCGccc                                       >  1:1621707/2‑69 (MQ=255)               cgccagACCACAAGCTAACATCCCCTGAAGCGCCGAACGTTTTGCGCCATATTGATCGGCCAGTCGcc                                        >  1:1635898/2‑68 (MQ=255)               ccccagACCACAAGCTAACATCCCCTGAAGCGCCGAACGTTTTGCGCCATATTGATCGGCCAGTCGccc                                       >  1:2794662/3‑69 (MQ=255)                       cacaAGCTAACATTCCCTGAAGCGCCGAACGTTTTGCACCATATTGATCGGCCAGTCGCCCGGCGTAAc                               <  1:1455429/69‑1 (MQ=255)                       cacaAGCTAACATTCCCTGAAGCGCCGAACGTTTTGCACCATATTGATCGGCCAGTCGCCCGGCGTAAc                               <  1:549107/69‑1 (MQ=255)                           aGCTAACATTCCCTGAAGCGCCGAACGTTTTGCACCATATTGATCGGCCAGt                                            >  1:1875075/1‑52 (MQ=255)                           aGCTAACATTCCCTGAAGCGCCGAACGTTTTGCACCATATTGATCGGCCAGTCGCCCGGCGt                                  >  1:2230022/1‑62 (MQ=255)                           aGCTAACATTCCCTGAAGCGCCGAACGTTTTGCACCATATTGATCGGCCAGTCGCCCGGCGTa                                 >  1:2331324/1‑63 (MQ=255)                           aGCTAACATTCCCTGAAGCGCCGAACGTTTTGCACCATATTGATCGGCCAGTCGCCCGGCGTa                                 >  1:2823977/1‑63 (MQ=255)                           aGCTAACATTCCCTGAAGCGCCGAACGTTTTGCACCATATTGATCGGCCAGTCGCCCGGCGTa                                 >  1:3435307/1‑63 (MQ=255)                           aGCTAACATTCCCTGAAGCGCCGAACGTTTTGCACCATATTGATCGGCCAGTCGCCCGGCGTa                                 >  1:2776221/1‑63 (MQ=255)                           aGCTAACATTCCCTGAAGCGCCGAACGTTTTGCACCATATTGATCGGCCAGTCGCCCGGCGTa                                 >  1:311252/1‑63 (MQ=255)                           aGCTAACATTCCCTGAAGCGCCGAACGTTTTGCACCATATTGATCGGCCAGTCGCCCGGCGTa                                 >  1:493424/1‑63 (MQ=255)                           aGCTAACATTCCCTGAAGCGCCGAACGTTTTGCACCATATTGATCGGCCAGTCGCCCGGCGTa                                 >  1:2507103/1‑63 (MQ=255)                           aGCTAACATTCCCTGAAGCGCCGAACGTTTTGCACCATATTGATCGGCCAGTCGCCCGGCGTa                                 >  1:2335186/1‑63 (MQ=255)                           aGCTAACATTCCCTGAAGCGCCGAACGTTTTGCACCATATTGATCGGCCAGTCGCCCGGCGTa                                 >  1:631955/1‑63 (MQ=255)                           aGCTAACATTCCCTGAAGCGCCGAACGTTTTGCACCATATTGATCGGCCAGTCGCCCGGCGTa                                 >  1:1521287/1‑63 (MQ=255)                           aGCTAACATTCCCTGAAGCGCCGAACGTTTTGCACCATATTGATCGGCCAGTCGCCCGGCGTa                                 >  1:2250047/1‑63 (MQ=255)                           aGCTAACATTCCCTGAAGCGCCGAACGTTTTGCACCATATTGATCGGCCAGTCGCCCGGCGTAACCAcg                           >  1:910980/1‑69 (MQ=255)                           aGCTAACATTCCCTGAAGCGCCGAACGTTTTGCACCATATTGATCGGCCAGTCGCCCGGCGTAACCAcg                           >  1:124354/1‑69 (MQ=255)                           aGCTAACATTCCCTGAAGCGCCGAACGTTTTGCACCATATTGATCGGCCAGTCGCCCGGCGTAACCAc                            >  1:2777646/1‑68 (MQ=255)                                     cccTGAAGCGCCGAACGTTTTGCACCATATTGATCGGCCAGTCGCCCGGCGTAACCACGCGTCAGCAc                  >  1:2561084/1‑68 (MQ=255)                                     cccTGAAGCGCCGAACGTTTTGCACCATATTGATCGGCCAGTCGCCCGGCGTAACCACGCGTCAGCAc                  >  1:2764059/1‑68 (MQ=255)                                     cccTGAAGCGCCGAACGTTTTGCACCATATTGATCGGCCAGTCGCCCGGCGTAACCACGCGTCAGCAc                  >  1:559303/1‑68 (MQ=255)                                     cccTGAAGCGCCGAACGTTTTGCACCATATTGATCGGCCAGTCGCCCGGCGTAACCACGCGTCAGCAc                  >  1:3291333/1‑68 (MQ=255)                                     cccTGAAGCGCCGAACGTTTTGCACCATATTGATCGGCCAGTCGCCCGGCGTAACCACGCGTCAGCAc                  >  1:897934/1‑68 (MQ=255)                                     cccTGAAGCGCCGAACGTTTTGCACCATATTGATCGGCCAGTCGCCCGGCGTAACCACGCGTCAGCAc                  >  1:719364/1‑68 (MQ=255)                                     cccTGAAGCGCCGAACGTTTTGCACCATATTGATCGGCCAGTCGCCCGGCGTAACCACGCGTCAGCAc                  >  1:2547747/1‑68 (MQ=255)                                        tGAAGCGCCGAACGTTTTGCACCATATTGATCGGCCAGTCGCCCGGCGTAACCACGCGTCAGCACCGTa              <  1:702783/69‑1 (MQ=255)                                        tGAAGCGCCGAACGTTTTGCACCATATTGATCGGCCAGTCGCCCGGCGTAACCACGCGTCAGCACCGTa              <  1:328305/69‑1 (MQ=255)                                        tGAAGCGCCGAACGTTTTGCACCATATTGATCGGCCAGTCGCCCGGCGTAACCACGCGTCAGCACCGTa              <  1:1809117/69‑1 (MQ=255)                                         gAAGCGCCGAACGTTTTGCACCATATTGATCGGCCAGTCGCCCGGCGTAACCACGCGTCAGCACCGTa              >  1:3296360/1‑68 (MQ=255)                                         gAAGCGCCGAACGTTTTGCACCATATTGATCGGCCAGTCGCCCGGCGTAACCACGCGTCAGCACCGTa              >  1:2287443/1‑68 (MQ=255)                                         gAAGCGCCGAACGTTTTGCACCATATTGATCGGCCAGTCGCCCGGCGTAACCACGCGTCAGCACCGTa              >  1:2004137/1‑68 (MQ=255)                                         gAAGCGCCGAACGTTTTGCACCATATTGATCGGCCAGTCGCCCGGCGTAACCACGCGTCAGCACCGTa              >  1:2767400/1‑68 (MQ=255)                                                    cGTTTTGCGCCATATTGATCGGCCAGTCGCCCGGCGTAACCACGCGTCAGCACCGTAGCCAGAAACTGa  <  1:821370/69‑1 (MQ=255)                                                    cGTTTTGCGCCATATTGATCGGCCAGTCGCCCGGCGTAACCACGCGTCAGCACCGTAGCCAGAAACTGa  <  1:47380/69‑1 (MQ=255)                                                    cGTTTTGCGCCATATTGATCGGCCAGTCGCCCGGCGTAACCACGCGTCAGCACCGTAGCCAGAAACTGa  <  1:2685108/69‑1 (MQ=255)                                                    cGTTTTGCGCCATATTGATCGGCCAGTCGCCCGGCGTAACCACGCGTCAGCACCGTAGCCAGAAACTGa  <  1:3033347/69‑1 (MQ=255)                                                    cGTTTTGCGCCATATTGATCGGCCAGTCGCCCGGCGTAACCACGCGTCAGCACCGTAGCCAGAAACTGa  <  1:1068847/69‑1 (MQ=255)                                                    cGTTTTGCGCCATATTGATCGGCCAGTCGCCCGGCGTAACCACGCGTCAGCACCGTAGCCAGAAACTGa  <  1:1444288/69‑1 (MQ=255)                                                    cGTTTTGCGCCATATTGATCGGCCAGTCGCCCGGCGTAACCACGCGTCAGCACCGTAGCCAGAAACTGa  <  1:1504492/69‑1 (MQ=255)                                                    cGTTTTGCGCCATATTGATCGGCCAGTCGCCCGGCGTAACCACGCGTCAGCACCGTAGCCAGAAACTGa  <  1:1618393/69‑1 (MQ=255)                                                    cGTTTTGCGCCATATTGATCGGCCAGTCGCCCGGCGTAACCACGCGTCAGCACCGTAGCCAGAAACTGa  <  1:2299910/69‑1 (MQ=255)                                                    cGTTTTGCGCCATATTGATCGGCCAGTCGCCCGGCGTAACCACGCGTCAGCACCGTAGCCAGAAACTGa  <  1:2681211/69‑1 (MQ=255)                                                    cGTTTTGCGCCATATTGATCGGCCAGTCGCCCGGCGTAACCACGCGTCAGCACCGTAGCCAGAAACTGa  <  1:2986420/69‑1 (MQ=255)                                                     gTTTTGCGCCATATTGATCGGCCAGTCGCCCGGCGTAACCACGCGTCAGCACCGTAGCCAGAAACTGa  <  1:1525986/68‑1 (MQ=255)                                                       tttGCGCCATATTGATCGGCCAGTCGCCCGGCGTAACCACGCGTCAGCACCGTAGCCAGaaa      >  1:2125560/1‑62 (MQ=255)                                                       tttGCGCCATATTGATCGGCCAGTCGCCCGGCGTAACCACGCGTCAGCACCGTAGCCAGAAAc     >  1:2446294/1‑63 (MQ=255)                                                       tttGCGCCATATTGATCGGCCAGTCGCCCGGCGTAACCACGCGTCAGCACCGTAGCCAGAAACt    >  1:496288/1‑64 (MQ=255)                                                       tttGCGCCATATTGATCGGCCAGTCGCCCGGCGTAACCACGCGTCAGCACCGTAGCCAGAAACt    >  1:2482217/1‑64 (MQ=255)                                                       tttGCGCCATATTGATCGGCCAGTCGCCCGGCGTAACCACGCGTCAGCACCGTAGCCAGAAACt    >  1:1654766/1‑64 (MQ=255)                                                       tttGCGCCATATTGATCGGCCAGTCGCCCGGCGTAACCACGCGTCAGCACCGTAGCCAGAAACt    >  1:2886469/1‑64 (MQ=255)                                                       tttGCGCCATATTGATCGGCCAGTCGCCCGGCGTAACCACGCGTCAGCACCGTAGCCAGAAACt    >  1:3372236/1‑64 (MQ=255)                                                       tttGCGCCATATTGATCGGCCAGTCGCCCGGCGTAACCACGCGTCAGCACCGTAGCCAGAAACt    >  1:3206975/1‑64 (MQ=255)                                                            |                                                             AGCAACGCGCCGCCAGCCAGACCACAAGCTAACATTCCCTGAAGCGCCGAACGTTTTGCACCATATTGATCGGCCAGTCGCCCGGCGTAACCACGCGTCAGCACCGTAGCCAGAAACTGA  >  minE/1525077‑1525196

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