Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	920,335	A→C	35.9%	P88P (CCT→CCG)	mpaA ←	murein peptide amidase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	920,335	0	A	C	35.9%	41.4 / 52.1	64	P88P (CCT→CCG)	mpaA	murein peptide amidase A
Reads supporting (aligned to +/- strand):  ref base A (16/25);  new base C (13/10);  total (29/35)
Fisher's exact test for biased strand distribution p-value = 2.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGGCCCGTAACCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGGAGAGCGTGACGACCGAAGAGTTTTCATC  >  minE/920267‑920388                                                                     |                                                      ttGGCCCGTAACCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAaagaa                                                       >  1:208983/1‑69 (MQ=255) ttGGCCCGTAACCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAaagaa                                                       >  1:637091/1‑69 (MQ=255) ttGGCCCGTAACCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAaagaa                                                       >  1:2326075/1‑69 (MQ=255) ttGGCCCGTAACCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAaagaa                                                       >  1:3403504/1‑69 (MQ=255) ttGGCCCGTAACCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAaagaa                                                       >  1:1340749/1‑69 (MQ=255) ttGGCCCGTAACCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAaagaa                                                       <  1:3395955/69‑1 (MQ=255) ttGGCCCGTAACCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAaagaa                                                       >  1:3231695/1‑69 (MQ=255) ttGGCCCGTAACCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAaagaa                                                       <  1:3108140/69‑1 (MQ=255) ttGGCCCGTAACCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAaagaa                                                       >  1:1788082/1‑69 (MQ=255) ttGGCCCGTAACCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAaagaa                                                       >  1:2207564/1‑69 (MQ=255) ttGGCCCGTAACCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAaagaa                                                       <  1:2331819/69‑1 (MQ=255) ttGGCCCGTAACCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAaagaa                                                       <  1:189931/69‑1 (MQ=255)       cGTAACCCCAACTGTCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAAAGAAGGTGTc                                                 >  1:559290/1‑69 (MQ=255)       cGTAACCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAAAGAAGGTGTc                                                 >  1:299942/1‑69 (MQ=255)       cGTAACCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAAAGAAGGTGTc                                                 >  1:2686259/1‑69 (MQ=255)       cGTAACCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAAAGAAGGTGTc                                                 >  1:2676900/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:2322887/1‑67 (MQ=255)        gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGgtgt                                                  >  1:1084515/1‑67 (MQ=255)        gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGgtgt                                                  >  1:2726594/1‑67 (MQ=255)        gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGgtgt                                                  >  1:2761532/1‑67 (MQ=255)        gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGgtgt                                                  >  1:2787902/1‑67 (MQ=255)        gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGgtgt                                                  >  1:2873168/1‑67 (MQ=255)        gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGgtgt                                                  >  1:2891522/1‑67 (MQ=255)        gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGgtgt                                                  >  1:409116/1‑67 (MQ=255)        gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGgtgt                                                  >  1:941457/1‑67 (MQ=255)        gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGgtgt                                                  >  1:1039746/1‑67 (MQ=255)        gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGgtgt                                                  >  1:1131508/1‑67 (MQ=255)        gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGgtgt                                                  >  1:1463846/1‑67 (MQ=255)        gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGGTGTCa                                                <  1:2311809/69‑1 (MQ=255)        gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGGTGTCa                                                <  1:907838/69‑1 (MQ=255)        gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGGTGTCa                                                <  1:8260/69‑1 (MQ=255)        gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGGTGTCa                                                <  1:1090285/69‑1 (MQ=255)        gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGGTGTCa                                                <  1:1193843/69‑1 (MQ=255)        gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGGTGTCa                                                <  1:3247688/69‑1 (MQ=255)        gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGGTGTCa                                                <  1:2198039/69‑1 (MQ=255)        gTAACCCTAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGGCGGCGCAAAGACGGTGTCa                                                <  1:2223533/69‑1 (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:984432/69‑1 (MQ=255)                                    cacacacaGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGGAGAGCGTgac                    <  1:2270645/69‑1 (MQ=255)                                    cacacacaGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGGAGAGCGTgac                    <  1:1092537/69‑1 (MQ=255)                                    cacacacaGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGGAGAGCGTgac                    <  1:2182803/69‑1 (MQ=255)                                        cacaGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTc                                           >  1:1381504/1‑42 (MQ=255)                                        cacaGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGgag                            >  1:19283/1‑57 (MQ=255)                                        cacaGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGgag                            >  1:957548/1‑57 (MQ=255)                                          caGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGgag                            <  1:1545411/55‑1 (MQ=255)                                          caGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGgag                            <  1:947623/55‑1 (MQ=255)                                          caGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGgag                            <  1:2322666/55‑1 (MQ=255)                                          caGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGgag                            <  1:897877/55‑1 (MQ=255)                                          caGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGgag                            <  1:873882/55‑1 (MQ=255)                                          caGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGgag                            <  1:2026814/55‑1 (MQ=255)                                          caGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGgag                            <  1:2764197/55‑1 (MQ=255)                                          caGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGgag                            <  1:647559/55‑1 (MQ=255)                                          caGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGgag                            <  1:2005543/55‑1 (MQ=255)                                          caGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGgag                            <  1:488169/55‑1 (MQ=255)                                          caGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGgag                            <  1:2416600/55‑1 (MQ=255)                                          caGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGgag                            <  1:3336908/55‑1 (MQ=255)                                          caGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGgag                            <  1:1949809/55‑1 (MQ=255)                                          caGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGgag                            <  1:3234850/55‑1 (MQ=255)                                          caGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGgag                            <  1:2656712/55‑1 (MQ=255)                                          caGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGgag                            <  1:3073732/55‑1 (MQ=255)                                          caGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGgag                            <  1:3000402/55‑1 (MQ=255)                                           aGCACCACATGATGACGGCGCAAAGACGGTGTCAACGTCCGCAGAGCACAGgag                            <  1:3205112/54‑1 (MQ=255)                                                      atgaCGGCGCAAAGAAGGTGTCAACGTCCGTAGAGCACAGGAGAGCGTGACGACCGAAGAGTTTtcatc  >  1:3316518/1‑69 (MQ=255)                                                      atgaCGGCGCAAAGAAGGTGTCAACGTCCGTAGAGCACAGGAGAGCGTGACGACCGAAGAGTTTtcatc  >  1:2377338/1‑69 (MQ=255)                                                      atgaCGGCGCAAAGAAGGTGTCAACGTCCGTAGAGCACAGGAGAGCGTGACGACCGAAGAGTTTtcatc  >  1:1442475/1‑69 (MQ=255)                                                      atgaCGGCGCAAAGAAGGTGTCAACGTCCGCAGAGca                                  <  1:2568126/37‑1 (MQ=255)                                                      atgaCGGCGCAAAGAAGGTGTCAACGTCCGCAGAGca                                  <  1:2644549/37‑1 (MQ=255)                                                      atgaCGGCGCAAAGAAGGTGTCAACGTCCGCAGAGca                                  <  1:2740209/37‑1 (MQ=255)                                                      atgaCGGCGCAAAGAAGGTGTCAACGTCCGCAGAGca                                  <  1:307818/37‑1 (MQ=255)                                                                     |                                                      TTGGCCCGTAACCCCAACTGGCAGCCGTCAGGATTCACACACAGCACCACATGATGACGGCGCAAAGAAGGTGTCAACGTCCGCAGAGCACAGGAGAGCGTGACGACCGAAGAGTTTTCATC  >  minE/920267‑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: ‑