Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	298,842	G→A	71.2%	S501S (TCG→TCA)	mdlA →	fused predicted multidrug transporter subunits and ATP‑binding components of ABC superfamily

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	298,842	0	G	A	71.2%	34.0 / 13.4	38	S501S (TCG→TCA)	mdlA	fused predicted multidrug transporter subunits and ATP‑binding components of ABC superfamily
Reads supporting (aligned to +/- strand):  ref base G (6/5);  new base A (5/22);  total (11/27)
Fisher's exact test for biased strand distribution p-value = 4.71e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGATCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAA  >  minE/298800‑298906                                           |                                                                 attaGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGc                                            >  1:1599996/1‑65 (MQ=255) attaGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGc                                            >  1:2914523/1‑65 (MQ=255) attaGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGc                                            >  1:2382338/1‑65 (MQ=255) attaGTCAACGCGGAAATCCTCATCCTTGATGATCCGCTTTCGGCGGTGGACGGACGCACTGAGc                                            >  1:343716/1‑65 (MQ=255)          cgcgGAAATCCTCATCCTTGATGATGCGCTTTCAGCGGTGgacggacg                                                    <  1:1301304/48‑1 (MQ=255)          cgcgGAAATCCTCATCCTTGATGATGCGCTTTCAGCGGTGgacggacg                                                    <  1:805203/48‑1 (MQ=255)          cgcgGAAATCCTCATCCTTGATGATGCGCTTTCAGCGGTGgacggacg                                                    <  1:76062/48‑1 (MQ=255)          cgcgGAAATCCTCATCCTTGATGATGCGCTTTCAGCGGTGgacggacg                                                    <  1:682824/48‑1 (MQ=255)          cgcgGAAATCCTCATCCTTGATGATGCGCTTTCAGCGGTGgacggacg                                                    <  1:577723/48‑1 (MQ=255)          cgcgGAAATCCTCATCCTTGATGATGCGCTTTCAGCGGTGgacggacg                                                    <  1:3189070/48‑1 (MQ=255)          cgcgGAAATCCTCATCCTTGATGATGCGCTTTCAGCGGTGgacggacg                                                    <  1:2947296/48‑1 (MQ=255)          cgcgGAAATCCTCATCCTTGATGATGCGCTTTCAGCGGTGgacggacg                                                    <  1:2702400/48‑1 (MQ=255)          cgcgGAAATCCTCATCCTTGATGATGCGCTTTCAGCGGTGgacggacg                                                    <  1:2652003/48‑1 (MQ=255)          cgcgGAAATCCTCATCCTTGATGATGCGCTTTCAGCGGTGgacggacg                                                    <  1:2325316/48‑1 (MQ=255)          cgcgGAAATCCTCATCCTTGATGATGCGCTTTCAGCGGTGgacggacg                                                    <  1:1864002/48‑1 (MQ=255)          cgcgGAAATCCTCATCCTTGATGATGCGCTTTCAGCGGTGgacggacg                                                    <  1:1596486/48‑1 (MQ=255)          cgcgGAAATCCTCATCCTTGATGATGCGCTTTCAGCGGTGgacggacg                                                    <  1:1457576/48‑1 (MQ=255)          cgcgGAAATCCTCATCCTTGATGATGCGCTTTCAGCGGTGgacggacg                                                    <  1:1496905/48‑1 (MQ=255)          cgcgGAAATCCTCATCCTTGATGATGCGCTTTCAGCGGTGgacggacg                                                    <  1:1365802/48‑1 (MQ=255)          cgcgGAAATCCTCATCCTTGATGATGCGCTTTCAGCGGTGgacggacg                                                    <  1:1859340/48‑1 (MQ=255)          cgcgGAAATCCTCATCCTTGATGATGCGCTTTCAGCGGTGgacggacg                                                    <  1:110799/48‑1 (MQ=255)          cgcgGAAATCCTCATCCTTGATGATGCGCTTTCAGCGGTGgacggacg                                                    <  1:1819373/48‑1 (MQ=255)          cgcgGAAATCCTCATCCTTGATGATGCGCTTTCAGCGGTGgacggacg                                                    <  1:1819518/48‑1 (MQ=255)          cgcgGAAATCCTCATCCTTGATGATGCGCTTTCAGCGGTGgacggacg                                                    <  1:1859011/48‑1 (MQ=255)          cgcgGAAATCCTCATCCTTGATGATGCGCTTTCAGCGGGGgacggacg                                                    <  1:97080/48‑1 (MQ=38)           gcgGAAATCCTCATCCTTGATGATGCGCTTTCAGCGGTGgacggacg                                                    <  1:855050/47‑1 (MQ=255)                 aTCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGc                                                   >  1:1728274/1‑42 (MQ=255)                 aTCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGc                                                   >  1:364708/1‑42 (MQ=255)                 aTCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGGTCCTGCATAACCTg                        <  1:1506413/69‑1 (MQ=255)                 aTCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGATCCTGCATAACCTg                        <  1:2849054/69‑1 (MQ=255)                 aTCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGATCCTGCATAACCTg                        <  1:632162/69‑1 (MQ=255)                      cATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGATCCTGCATAACCTGCGTCa                   <  1:343650/69‑1 (MQ=255)                                   gcgcTTTCGGCGGTGGACGGACGCACTGAGCACCAGATCCTGCATAACCTGCGTc                    <  1:1655810/55‑1 (MQ=255)                                       tttCAGCGGTGGACGGACGCACAGAGCACCAGATCCTGCATAACCTGCGTCAGTGGGGGCAGGgaagaa  >  1:3435300/1‑69 (MQ=255)                                       tttCAGCGGTGGACGGACGCACAGAGCACCAGATCCTGCATAACCTGCGTCAGTGGGGGCAGGgaagaa  >  1:452620/1‑69 (MQ=255)                                       tttCAGCGGTGGACGGACGCACAGAGCACCAGATCCTGCATAACCTGCGTCAGTGGGGGCAGGgaagaa  >  1:3212196/1‑69 (MQ=255)                                       tttCAGCGGTGGACGGACGCACAGAGCACCAGATCCTGCATAACCTGCGTCAGTGGGGGCAGGgaagaa  >  1:1806006/1‑69 (MQ=255)                                       tttCAGCGGTGGACGGACGCACAGAGCACCAGATCCTGCATAACCTGCGTCAGTGGGGGCAGGgaagaa  >  1:215418/1‑69 (MQ=255)                                           |                                                                 ATTAGTCAACGCGGAAATCCTCATCCTTGATGATGCGCTTTCGGCGGTGGACGGACGCACTGAGCACCAGATCCTGCATAACCTGCGTCAGTGGGGGCAGGGAAGAA  >  minE/298800‑298906

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