Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	920,323	A→G	38.9%	R92R (CGT→CGC)	mpaA ←	murein peptide amidase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	920,323	0	A	G	38.9%	34.3 / 51.7	77	R92R (CGT→CGC)	mpaA	murein peptide amidase A
Reads supporting (aligned to +/- strand):  ref base A (21/26);  new base G (13/17);  total (34/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.61e-01

CCATTAGCATTGGCCCGTAACCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGGAGAGCGTGACGACCGAAGAGTTTTCATC  >  minE/920258‑920388                                                                  |                                                                  ccATTAGCATTGGCCCGTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCgg                                                                <  1:2164305/69‑1 (MQ=255) ccATTAGCATTGGCCCGTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCgg                                                                <  1:2194438/69‑1 (MQ=255) ccATTAGCATTGGCCCGTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCgg                                                                <  1:476617/69‑1 (MQ=255) ccATTAGCATTGGCCCGTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCgg                                                                <  1:3432658/69‑1 (MQ=255) ccATTAGCATTGGCCCGTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCgg                                                                <  1:324947/69‑1 (MQ=255) ccATTAGCATTGGCCCGTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCgg                                                                <  1:2782723/69‑1 (MQ=255)  cATTAGCATTGGCCCGTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCgg                                                                <  1:3189570/68‑1 (MQ=255)  cATTAGCATTGGCCCGTAACCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCgg                                                                <  1:413081/68‑1 (MQ=255)          ttGGCCCGTAACCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAaagaa                                                       <  1:3395955/69‑1 (MQ=255)          ttGGCCCGTAACCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAaagaa                                                       >  1:208983/1‑69 (MQ=255)          ttGGCCCGTAACCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAaagaa                                                       <  1:3108140/69‑1 (MQ=255)          ttGGCCCGTAACCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAaagaa                                                       >  1:3231695/1‑69 (MQ=255)          ttGGCCCGTAACCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAaagaa                                                       <  1:189931/69‑1 (MQ=255)          ttGGCCCGTAACCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAaagaa                                                       >  1:2207564/1‑69 (MQ=255)          ttGGCCCGTAACCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAaagaa                                                       >  1:1788082/1‑69 (MQ=255)          ttGGCCCGTAACCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAaagaa                                                       >  1:3403504/1‑69 (MQ=255)          ttGGCCCGTAACCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAaagaa                                                       >  1:1340749/1‑69 (MQ=255)          ttGGCCCGTAACCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAaagaa                                                       >  1:2326075/1‑69 (MQ=255)          ttGGCCCGTAACCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAaagaa                                                       >  1:637091/1‑69 (MQ=255)          ttGGCCCGTAACCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAaagaa                                                       <  1:2331819/69‑1 (MQ=255)                cGTAACCCCAACTGTCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAAAGAAGGTGTc                                                 >  1:559290/1‑69 (MQ=255)                cGTAACCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAAAGAAGGTGTc                                                 >  1:2676900/1‑69 (MQ=255)                cGTAACCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAAAGAAGGTGTc                                                 >  1:2686259/1‑69 (MQ=255)                cGTAACCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAAAGAAGGTGTc                                                 >  1:299942/1‑69 (MQ=255)                cGTAACCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAAAGAAGGTGTc                                                 >  1:1827228/1‑69 (MQ=255)                cGTAACCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAAAGAAGGTGTc                                                 >  1:1825369/1‑69 (MQ=255)                cGTAACCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAAAGAAGGTGTc                                                 >  1:1179222/1‑69 (MQ=255)                 gTAACCCTAACTGGCAGCCGTGAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGgtgt                                                  >  1:2081149/1‑67 (MQ=255)                 gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGgtgt                                                  >  1:2891522/1‑67 (MQ=255)                 gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGgtgt                                                  >  1:2873168/1‑67 (MQ=255)                 gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGgtgt                                                  >  1:2787902/1‑67 (MQ=255)                 gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGgtgt                                                  >  1:2761532/1‑67 (MQ=255)                 gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGgtgt                                                  >  1:2726594/1‑67 (MQ=255)                 gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGgtgt                                                  >  1:409116/1‑67 (MQ=255)                 gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGgtgt                                                  >  1:941457/1‑67 (MQ=255)                 gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGgtgt                                                  >  1:1084515/1‑67 (MQ=255)                 gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGgtgt                                                  >  1:1039746/1‑67 (MQ=255)                 gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGgtgt                                                  >  1:1463846/1‑67 (MQ=255)                 gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGgtgt                                                  >  1:1131508/1‑67 (MQ=255)                 gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGgtgt                                                  >  1:2322887/1‑67 (MQ=255)                 gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGGTGTCa                                                <  1:8260/69‑1 (MQ=255)                 gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGGTGTCa                                                <  1:1090285/69‑1 (MQ=255)                 gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGGTGTCa                                                <  1:2311809/69‑1 (MQ=255)                 gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGGTGTCa                                                <  1:907838/69‑1 (MQ=255)                 gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGGTGTCa                                                <  1:1193843/69‑1 (MQ=255)                 gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGGTGTCa                                                <  1:2223533/69‑1 (MQ=255)                 gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGGTGTCa                                                <  1:2198039/69‑1 (MQ=255)                 gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGGTGTCa                                                <  1:3247688/69‑1 (MQ=255)                 gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCAAATGATGgcggcg                                                              >  1:377522/1‑55 (MQ=255)                    aCCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAAAGAAGGTGTCAACg                                             >  1:66984/1‑69 (MQ=255)                    aCCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAAAGAAGGTGTCAACg                                             >  1:3312060/1‑69 (MQ=255)                             ggCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGGTGTCa                                                <  1:2128245/57‑1 (MQ=255)                                          aTTCACACACAGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGGAGAGCGt                       >  1:2889312/1‑69 (MQ=255)                                             cacacacaGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGGAGAGCGTgac                    <  1:1092537/69‑1 (MQ=255)                                             cacacacaGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGGAGAGCGTgac                    <  1:984432/69‑1 (MQ=255)                                             cacacacaGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGGAGAGCGTgac                    <  1:2270645/69‑1 (MQ=255)                                             cacacacaGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGGAGAGCGTgac                    <  1:2182803/69‑1 (MQ=255)                                                 cacaGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTc                                           >  1:1381504/1‑42 (MQ=255)                                                 cacaGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGgag                            >  1:957548/1‑57 (MQ=255)                                                 cacaGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGgag                            >  1:19283/1‑57 (MQ=255)                                                   caGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGgag                            <  1:2416600/55‑1 (MQ=255)                                                   caGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGgag                            <  1:2026814/55‑1 (MQ=255)                                                   caGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGgag                            <  1:947623/55‑1 (MQ=255)                                                   caGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGgag                            <  1:897877/55‑1 (MQ=255)                                                   caGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGgag                            <  1:873882/55‑1 (MQ=255)                                                   caGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGgag                            <  1:2656712/55‑1 (MQ=255)                                                   caGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGgag                            <  1:647559/55‑1 (MQ=255)                                                   caGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGgag                            <  1:2322666/55‑1 (MQ=255)                                                   caGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGgag                            <  1:488169/55‑1 (MQ=255)                                                   caGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGgag                            <  1:2764197/55‑1 (MQ=255)                                                   caGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGgag                            <  1:3000402/55‑1 (MQ=255)                                                   caGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGgag                            <  1:1545411/55‑1 (MQ=255)                                                   caGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGgag                            <  1:3336908/55‑1 (MQ=255)                                                   caGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGgag                            <  1:3073732/55‑1 (MQ=255)                                                   caGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGgag                            <  1:3234850/55‑1 (MQ=255)                                                   caGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGgag                            <  1:1949809/55‑1 (MQ=255)                                                   caGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGgag                            <  1:2005543/55‑1 (MQ=255)                                                    aGCACCACATGATGACGGCGCAAAGACGGTGTCAACGTCCGCAGAGCACAGgag                            <  1:3205112/54‑1 (MQ=255)                                                               atgaCGGCGCAAAGAAGGTGTCAACGTCCGTAGAGCACAGGAGAGCGTGACGACCGAAGAGTTTtcatc  >  1:1442475/1‑69 (MQ=255)                                                               atgaCGGCGCAAAGAAGGTGTCAACGTCCGTAGAGCACAGGAGAGCGTGACGACCGAAGAGTTTtcatc  >  1:3316518/1‑69 (MQ=255)                                                               atgaCGGCGCAAAGAAGGTGTCAACGTCCGTAGAGCACAGGAGAGCGTGACGACCGAAGAGTTTtcatc  >  1:2377338/1‑69 (MQ=255)                                                               atgaCGGCGCAAAGAAGGTGTCAACGTCCGCAGAGca                                  <  1:2740209/37‑1 (MQ=255)                                                               atgaCGGCGCAAAGAAGGTGTCAACGTCCGCAGAGca                                  <  1:307818/37‑1 (MQ=255)                                                               atgaCGGCGCAAAGAAGGTGTCAACGTCCGCAGAGca                                  <  1:2644549/37‑1 (MQ=255)                                                               atgaCGGCGCAAAGAAGGTGTCAACGTCCGCAGAGca                                  <  1:2568126/37‑1 (MQ=255)                                                                  |                                                                  CCATTAGCATTGGCCCGTAACCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGGAGAGCGTGACGACCGAAGAGTTTTCATC  >  minE/920258‑920388

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