Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,835,770	G→A	35.1%	L43L (CTG→TTG)	yjeJ ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,835,770	0	G	A	35.1%	63.3 / 49.5	88	L43L (CTG→TTG)	yjeJ	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (44/13);  new base A (16/15);  total (60/28)
Fisher's exact test for biased strand distribution p-value = 1.77e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.84e-01

CAGCAGCATCCAGTTTATGCTTTTGATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAATTTTCAGCGCCA  >  minE/2835707‑2835832                                                                |                                                               cagcagcaTCCAGTTTATGCTTTTGATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTcc                                                            >  1:285728/1‑68 (MQ=255) cagcagcaTCCAGTTTATGCTTTTGATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTcc                                                            >  1:1772146/1‑68 (MQ=255) cagcagcaTCCAGTTTATGCTTTTGATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTcc                                                            >  1:1060752/1‑68 (MQ=255) cagcagcaTCCAGTTTATGCTTTTGATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTcc                                                            >  1:1598183/1‑68 (MQ=255) cagcagcaTCCAGTTTATGCTTTTGATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTcc                                                            >  1:1921166/1‑68 (MQ=255) cagcagcaTCCAGTTTATGCTTTTGATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTcc                                                            >  1:718292/1‑68 (MQ=255) cagcagcaTCCAGTTTATGCTTTTGATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTcc                                                            >  1:1386356/1‑68 (MQ=255) cagcagcaTCCAGTTTATGCTTTTGATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTcc                                                            >  1:3286111/1‑68 (MQ=255) cagcagcaTCCAGTTTATGCTTTTGATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTcc                                                            >  1:2859777/1‑68 (MQ=255) cagcagcaTCCAGTTTATGCTTTTGATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTcc                                                            >  1:2863859/1‑68 (MQ=255) cagcagcaTCCAGTTTATGCTTTTGATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTcc                                                            >  1:85046/1‑68 (MQ=255) cagcagcaTCCAGTTTATGCTTTTGATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTcc                                                            >  1:876421/1‑68 (MQ=255) cagcagcaTCCAGTTTATGCTTTTGATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTcc                                                            >  1:109794/1‑68 (MQ=255) cagcagcaTCCAGTTTATGCTTTTGATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTcc                                                            >  1:2077789/1‑68 (MQ=255)           cAGTTTATGCTTTTGATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGtc                                                              >  1:2314903/1‑55 (MQ=255)           cAGTTTATGCTTTTGATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTCCATgacagaca                                                  >  1:2636693/1‑68 (MQ=255)           cAGTTTATGCTTTTGATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTCCATgacagaca                                                  >  1:1940320/1‑68 (MQ=255)           cAGTTTATGCTTTTGATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTCCATgacagaca                                                  >  1:2124448/1‑68 (MQ=255)           cAGTTTATGCTTTTGATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTCCATgacagaca                                                  >  1:1865111/1‑68 (MQ=255)           cAGTTTATGCTTTTGATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTCCATgacagaca                                                  >  1:1660779/1‑68 (MQ=255)           cAGTTTATGCTTTTGATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTCCATgacagaca                                                  >  1:581859/1‑68 (MQ=255)                tATGCTTTTGATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTCCATGACAGACATAAAg                                             >  1:2675024/1‑68 (MQ=255)                         gATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTCCATGACAGACATAAAGAACAACGAt                                    <  1:1793838/68‑1 (MQ=255)                         gATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTCCATGACAGACATAAAGAACAACGAt                                    <  1:3206593/68‑1 (MQ=255)                         gATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTCCATGACAGACATAAAGAACAACGAt                                    <  1:2815551/68‑1 (MQ=255)                            gCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCtt                                <  1:2785635/69‑1 (MQ=255)                             cAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTCCATGACAGACATAAAGAACAACGAt                                    <  1:1223755/64‑1 (MQ=255)                                    cTTTCCAGCGCAATCCGAAGGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCTTTATTGCGa                        <  1:1895385/69‑1 (MQ=255)                                    cTTTCCAGCGCAATCAGAAGGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCTTTATTGCGa                        <  1:1075207/69‑1 (MQ=255)                                    cTTTCCAGCGCAATCAGAAGGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCTTTATTGCGa                        <  1:2186196/69‑1 (MQ=255)                                    cTTTCCAGCGCAATCAGAAGGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCTTTATTGCGa                        <  1:3017882/69‑1 (MQ=255)                                    cTTTCCAGCGCAATCAGAAGGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCTTTATTGCGa                        <  1:2115893/69‑1 (MQ=255)                                    cTTTCCAGCGCAATCAGAAGGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCTTTATTGCGa                        <  1:3013931/69‑1 (MQ=255)                                    cTTTCCAGCGCAATCAGAAGGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCTTTATTGCGa                        <  1:3300830/69‑1 (MQ=255)                                    cTTTCCAGCGCAATCAGAAGGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCTTTATTGCGa                        <  1:3306691/69‑1 (MQ=255)                                    cTTTCCAGCGCAATCAGAAGGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCTTTATTGCGa                        <  1:184580/69‑1 (MQ=255)                                    cTTTCCAGCGCAATCAGAAGGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCTTTATTGCGa                        <  1:1250562/69‑1 (MQ=255)                                    cTTTCCAGCGCAATCAGAAGGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCTTTATTGCGa                        <  1:833351/69‑1 (MQ=255)                                    cTTTCCAGCGCAATCAGAAGGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCTTTATTGCGa                        <  1:1381421/69‑1 (MQ=255)                                    cTTTCCAGCGCAATCAGAAGGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCTTTATTGCGa                        <  1:186908/69‑1 (MQ=255)                                     tttCCAGCGCAATCAGAAGGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCTTTTTTGCGa                        <  1:435976/68‑1 (MQ=255)                                         cAGCGCAATCAGCAGGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGt                     <  1:1575438/67‑1 (MQ=255)                                         cAGCGCAATCAGCAGGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGt                     <  1:3029510/67‑1 (MQ=255)                                         cAGCGCAATCAGCAGGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGt                     <  1:1947381/67‑1 (MQ=255)                                                  cagcagGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAAtt            >  1:133303/1‑67 (MQ=255)                                                  cagcagGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAAtt            >  1:2872810/1‑67 (MQ=255)                                                  cagcagGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAAtt            >  1:2804332/1‑67 (MQ=255)                                                  cagcagGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAAtt            >  1:1432957/1‑67 (MQ=255)                                                  cagcagGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAAtt            >  1:937015/1‑67 (MQ=255)                                                  cagcagGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAAtt            >  1:2048714/1‑67 (MQ=255)                                                  cagcagGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAAtt            >  1:646030/1‑67 (MQ=255)                                                  cagcagGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAAtt            >  1:3209122/1‑67 (MQ=255)                                                  cagcagGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAAtt            >  1:2472444/1‑67 (MQ=255)                                                  cagcagGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAAtt            >  1:3234590/1‑67 (MQ=255)                                                  cagcagGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAAtt            >  1:1887928/1‑67 (MQ=255)                                                  cagcagGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAAtt            >  1:2133710/1‑67 (MQ=255)                                                  cagcagGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAAtt            >  1:1944136/1‑67 (MQ=255)                                                  cagcagGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAAtt            >  1:2068928/1‑67 (MQ=255)                                                  cagcacGTCGCGCAGTTCCATGACAGACATAAAGaa                                           >  1:2707153/1‑36 (MQ=255)                                                   agcagGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAATTTTc         >  1:430229/1‑69 (MQ=255)                                                     aagGTCGCGCAATTCCATGGCAGACATAAAGaacaac                                       >  1:243512/2‑37 (MQ=38)                                                     aagGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCtt                                >  1:2758838/2‑44 (MQ=255)                                                     aagGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCtt                                >  1:1122239/2‑44 (MQ=255)                                                     aagGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCtt                                >  1:1164265/2‑44 (MQ=255)                                                     aagGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCtt                                >  1:1254525/2‑44 (MQ=255)                                                     aagGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCtt                                >  1:1807422/2‑44 (MQ=255)                                                     aagGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCtt                                >  1:1930840/2‑44 (MQ=255)                                                     aagGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCtt                                >  1:1965304/2‑44 (MQ=255)                                                     aagGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCtt                                >  1:2145464/2‑44 (MQ=255)                                                     aagGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCtt                                >  1:2331850/2‑44 (MQ=255)                                                     aagGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCtt                                >  1:2932418/2‑44 (MQ=255)                                                     aagGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCtt                                >  1:3086633/2‑44 (MQ=255)                                                     aagGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCtt                                >  1:2543540/2‑44 (MQ=255)                                                     aagGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCtt                                >  1:2539159/2‑44 (MQ=255)                                                     aagGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCtt                                >  1:3273485/2‑44 (MQ=255)                                                     aagGTCGCGCAATTCCATGGCAGAAATAAAGAACAACGATTCtt                                >  1:681150/2‑44 (MQ=37)                                                      agGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAATTTTCAgc      >  1:1602043/1‑69 (MQ=255)                                                      agGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCTTTATTGCGa                        <  1:3210170/51‑1 (MQ=255)                                                       ggTCGCGCAGTTCCATGACAGACATAAAGAACAACGAt                                    <  1:561747/38‑1 (MQ=255)                                                         tCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAATTTTCAGCGcc   <  1:160801/69‑1 (MQ=255)                                                         tCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAATTTTCAGCGcc   <  1:1978802/69‑1 (MQ=255)                                                         tCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAATTTTCAGCGcc   <  1:110536/69‑1 (MQ=255)                                                         tCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAATTTTCAGCGcc   <  1:93916/69‑1 (MQ=255)                                                          cgcgCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAATTTTCAGCGCCa  >  1:1430813/1‑69 (MQ=255)                                                          cgcgCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAATTTTCAGCGCCa  >  1:406889/1‑69 (MQ=255)                                                          cgcgCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAATTTTCAGCGCCa  >  1:3228114/1‑69 (MQ=255)                                                          cgcgCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAATTTTCAGCGCCa  >  1:989661/1‑69 (MQ=255)                                                          cgcgCAGCTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAATTTTCAGCGCCa  >  1:973026/1‑69 (MQ=255)                                                                |                                                               CAGCAGCATCCAGTTTATGCTTTTGATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAATTTTCAGCGCCA  >  minE/2835707‑2835832

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