Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,835,778	A→G	40.8%	V40A (GTC→GCC)	yjeJ ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,835,778	0	A	G	40.8%	27.0 / 54.2	76	V40A (GTC→GCC)	yjeJ	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (27/18);  new base G (16/15);  total (43/33)
Fisher's exact test for biased strand distribution p-value = 4.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGTTTATGCTTTTGATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAATTTTCAGCGCCAGTGCAATGAAATT  >  minE/2835717‑2835845                                                              |                                                                    cAGTTTATGCTTTTGATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTCCATgacagaca                                                               >  1:2636693/1‑68 (MQ=255) cAGTTTATGCTTTTGATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTCCATgacagaca                                                               >  1:1865111/1‑68 (MQ=255) cAGTTTATGCTTTTGATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTCCATgacagaca                                                               >  1:1660779/1‑68 (MQ=255) cAGTTTATGCTTTTGATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTCCATgacagaca                                                               >  1:1940320/1‑68 (MQ=255) cAGTTTATGCTTTTGATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTCCATgacagaca                                                               >  1:2124448/1‑68 (MQ=255) cAGTTTATGCTTTTGATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTCCATgacagaca                                                               >  1:581859/1‑68 (MQ=255)      tATGCTTTTGATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTCCATGACAGACATAAAg                                                          >  1:2675024/1‑68 (MQ=255)               gATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTCCATGACAGACATAAAGAACAACGAt                                                 <  1:2815551/68‑1 (MQ=255)               gATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTCCATGACAGACATAAAGAACAACGAt                                                 <  1:3206593/68‑1 (MQ=255)               gATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTCCATGACAGACATAAAGAACAACGAt                                                 <  1:1793838/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:2115893/69‑1 (MQ=255)                          cTTTCCAGCGCAATCAGAAGGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCTTTATTGCGa                                     <  1:2186196/69‑1 (MQ=255)                          cTTTCCAGCGCAATCAGAAGGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCTTTATTGCGa                                     <  1:3013931/69‑1 (MQ=255)                          cTTTCCAGCGCAATCAGAAGGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCTTTATTGCGa                                     <  1:3017882/69‑1 (MQ=255)                          cTTTCCAGCGCAATCAGAAGGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCTTTATTGCGa                                     <  1:186908/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:1381421/69‑1 (MQ=255)                          cTTTCCAGCGCAATCAGAAGGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCTTTATTGCGa                                     <  1:1250562/69‑1 (MQ=255)                          cTTTCCAGCGCAATCAGAAGGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCTTTATTGCGa                                     <  1:833351/69‑1 (MQ=255)                          cTTTCCAGCGCAATCAGAAGGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCTTTATTGCGa                                     <  1:1075207/69‑1 (MQ=255)                           tttCCAGCGCAATCAGAAGGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCTTTTTTGCGa                                     <  1:435976/68‑1 (MQ=255)                               cAGCGCAATCAGCAGGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGt                                  <  1:1575438/67‑1 (MQ=255)                               cAGCGCAATCAGCAGGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGt                                  <  1:1947381/67‑1 (MQ=255)                               cAGCGCAATCAGCAGGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGt                                  <  1:3029510/67‑1 (MQ=255)                                        cagcagGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAAtt                         >  1:1432957/1‑67 (MQ=255)                                        cagcagGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAAtt                         >  1:2804332/1‑67 (MQ=255)                                        cagcagGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAAtt                         >  1:133303/1‑67 (MQ=255)                                        cagcagGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAAtt                         >  1:2472444/1‑67 (MQ=255)                                        cagcagGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAAtt                         >  1:2872810/1‑67 (MQ=255)                                        cagcagGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAAtt                         >  1:937015/1‑67 (MQ=255)                                        cagcagGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAAtt                         >  1:3234590/1‑67 (MQ=255)                                        cagcagGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAAtt                         >  1:2133710/1‑67 (MQ=255)                                        cagcagGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAAtt                         >  1:1887928/1‑67 (MQ=255)                                        cagcagGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAAtt                         >  1:2068928/1‑67 (MQ=255)                                        cagcagGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAAtt                         >  1:2048714/1‑67 (MQ=255)                                        cagcagGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAAtt                         >  1:646030/1‑67 (MQ=255)                                        cagcagGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAAtt                         >  1:1944136/1‑67 (MQ=255)                                        cagcagGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAAtt                         >  1:3209122/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:1164265/2‑44 (MQ=255)                                           aagGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCtt                                             >  1:2932418/2‑44 (MQ=255)                                           aagGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCtt                                             >  1:3086633/2‑44 (MQ=255)                                           aagGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCtt                                             >  1:1122239/2‑44 (MQ=255)                                           aagGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCtt                                             >  1:1807422/2‑44 (MQ=255)                                           aagGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCtt                                             >  1:3273485/2‑44 (MQ=255)                                           aagGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCtt                                             >  1:2758838/2‑44 (MQ=255)                                           aagGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCtt                                             >  1:2543540/2‑44 (MQ=255)                                           aagGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCtt                                             >  1:2539159/2‑44 (MQ=255)                                           aagGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCtt                                             >  1:2331850/2‑44 (MQ=255)                                           aagGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCtt                                             >  1:2145464/2‑44 (MQ=255)                                           aagGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCtt                                             >  1:1254525/2‑44 (MQ=255)                                           aagGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCtt                                             >  1:1965304/2‑44 (MQ=255)                                           aagGTCGCGCAATTCCATGGCAGACATAAAGAACAACGATTCtt                                             >  1:1930840/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:110536/69‑1 (MQ=255)                                               tCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAATTTTCAGCGcc                <  1:160801/69‑1 (MQ=255)                                               tCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAATTTTCAGCGcc                <  1:1978802/69‑1 (MQ=255)                                               tCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAATTTTCAGCGcc                <  1:93916/69‑1 (MQ=255)                                                cgcgCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAATTTTCAGCGCCa               >  1:406889/1‑69 (MQ=255)                                                cgcgCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAATTTTCAGCGCCa               >  1:989661/1‑69 (MQ=255)                                                cgcgCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAATTTTCAGCGCCa               >  1:3228114/1‑69 (MQ=255)                                                cgcgCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAATTTTCAGCGCCa               >  1:1430813/1‑69 (MQ=255)                                                cgcgCAGCTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAATTTTCAGCGCCa               >  1:973026/1‑69 (MQ=255)                                                        tCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAATTTTCAGCGCCAGTGCaa         >  1:2463961/1‑67 (MQ=255)                                                         ccATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAATTTTCAGCGCCAGTGCAATGa      <  1:2138796/69‑1 (MQ=255)                                                         ccATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAATTTTCAGCGCCAGTGCAATGa      >  1:898608/1‑69 (MQ=255)                                                         ccATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAATTTTCAGCGCCAGTGCAATGa      >  1:1122299/1‑69 (MQ=255)                                                         caaTGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAATTTTCAGCGCCAGTGCAATGa      >  1:1560777/3‑69 (MQ=255)                                                           aTGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAATTTTCAGCGCCAGTGCAATGaa     <  1:2681508/68‑1 (MQ=255)                                                           aTGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAATTTTCAGCGCCAGTGCAATGaa     <  1:2008196/68‑1 (MQ=255)                                                           aTGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAATTTTCAGCGCCAGTGCAATGaa     <  1:1908061/68‑1 (MQ=255)                                                           aTGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAATTTTCAGCGCCAGTGCAATGaa     <  1:1141061/68‑1 (MQ=255)                                                              acagacaTAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAATTTTCAGCGCCAGTGCAATGAAatt  >  1:3305453/1‑68 (MQ=255)                                                              acagacaTAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAATTTTCAGCGCCAGTGCAATGAAatt  >  1:3322256/1‑68 (MQ=255)                                                              acagacaTAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAATTTTCAGCGCCAGTGCAATGAAatt  >  1:2871785/1‑68 (MQ=255)                                                              acagacaTAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAATTTTCAGCGCCAGTGCAATGAAatt  >  1:812419/1‑68 (MQ=255)                                                              |                                                                    CAGTTTATGCTTTTGATGCAGACGACTTTCCAGCGCAATCAGCAGGTCGCGCAGTTCCATGACAGACATAAAGAACAACGATTCTTTATTGCGAGGTTCTTTAATTTTCAGCGCCAGTGCAATGAAATT  >  minE/2835717‑2835845

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