Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,969,046	C→A	52.4%	T47N (ACC→AAC) ‡	yggG →	predicted peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,969,046	0	C	A	52.4%	‑4.3 / 100.5	82	T47N (ACC→AAC) ‡	yggG	predicted peptidase
Reads supporting (aligned to +/- strand):  ref base C (25/14);  new base A (24/19);  total (49/33)
Fisher's exact test for biased strand distribution p-value = 5.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGATGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCG  >  minE/1968988‑1969108                                                           |                                                               ctctCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGATGCGCAGGTGAAAAcc                                                               >  1:2384250/1‑60 (MQ=255) ctctCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGATGCGCAGGTGAAAAcc                                                               >  1:3163300/1‑60 (MQ=255) ctctCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGATGCGCAGGTGAAAAcc                                                               >  1:3124636/1‑60 (MQ=255) ctctCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGATGCGCAGGTGAAAAcc                                                               >  1:2860894/1‑60 (MQ=255) ctctCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGATGCGCAGGTGAAAAcc                                                               >  1:2474302/1‑60 (MQ=255) ctctCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGATGCGCAGGTGAAAAcc                                                               >  1:448233/1‑60 (MQ=255) ctctCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGATGCGCAGGTGAAAAcc                                                               >  1:733338/1‑60 (MQ=255) ctctCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGATGCGCAGGTGAAAAcc                                                               >  1:147595/1‑60 (MQ=255)   ctcATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGATGCGCAGGTGAAAACCCTGAGCGAt                                                      <  1:3191133/67‑1 (MQ=255)   ctcATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGATGCGCAGGTGAAAACCCTGAGCGAt                                                      <  1:1732432/67‑1 (MQ=255)        caGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGATGCGCAGGTGAAAACCCTGAGCGATCAGGCa                                                <  1:41883/68‑1 (MQ=255)        caGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGATGCGCAGGTGAAAACCCTGAGCGATCAGGCa                                                <  1:417275/68‑1 (MQ=255)        caGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGATGCGCAGGTGAAAACCCTGAGCGATCAGGCa                                                <  1:2094042/68‑1 (MQ=255)        caGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGATGCGCAGGTGAAAACCCTGAGCGATCAGGCa                                                <  1:534440/68‑1 (MQ=255)                  aGCTTTTCAGGCTTACAGTTTGAGTGATGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGgaga                                      >  1:821692/1‑68 (MQ=255)                  aGCTTTTCAGGCTTACAGTTTGAGTGATGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGgaga                                      >  1:2770788/1‑68 (MQ=255)                   gCTTTTCAGGCTTACAGTTTGAGTGATGCGCAGGTGAAAACCCTGAGCg                                                        >  1:2090877/1‑49 (MQ=255)                   gCTTTTCAGGCTTACAGTTTGAGTGATGCGCAGGTGAAAACCCTGAGCg                                                        >  1:846724/1‑49 (MQ=255)                   gCTTTTCAGGCTTACAGTTTGAGTGATGCGCAGGTGAAAACCCTGAGCg                                                        >  1:3082810/1‑49 (MQ=255)                     ttttCAGGCTTACAGTTTGAGTGATGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATg                                    >  1:2358568/1‑67 (MQ=255)                      tttCAGGCTTACAGTTTGAGTGATGCGCAGGTGAAAAAACTGAGCGATCAGGCATgcc                                            <  1:2844454/58‑3 (MQ=38)                      tttCAGGCTTACAGTTTGAGTGATGCGCAGGTGAAAAAACTGAGCGATCAGGCATgcc                                            <  1:2406913/58‑3 (MQ=38)                      tttCAGGCTTACAGTTTGAGTGATGCGCAGGTGAAAAAACTGAGCGATCAGGCATgcc                                            <  1:588951/58‑3 (MQ=38)                      tttCAGGCTTACAGTTTGAGTGATGCGCAGGTGAAAAAACTGAGCGATCAGGCATgcc                                            <  1:2171939/58‑3 (MQ=38)                      tttCAGGCTTACAGTTTGAGTGATGCGCAGGTGAAAAAACTGAGCGATCAGGCATgcc                                            <  1:591392/58‑3 (MQ=38)                                  aGTTTGAGTGATGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGcgac                      <  1:1968658/68‑1 (MQ=255)                                  aGTTTGAGTGATGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGcgac                      <  1:1435413/68‑1 (MQ=255)                                  aGTTTGAGTGATGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGcgac                      <  1:1990900/68‑1 (MQ=255)                                      tGAGTGATGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGga                                        <  1:3237328/46‑1 (MQ=255)                                      tGAGTGATGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGga                                        <  1:1240331/46‑1 (MQ=255)                                      tGAGTGATGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGga                                        <  1:781799/46‑1 (MQ=255)                                      tGAGTGATGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGga                                        <  1:2489271/46‑1 (MQ=255)                                       gAGTGATGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGga                                        <  1:862186/45‑1 (MQ=255)                                          tGATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGAtt                  <  1:2991797/64‑1 (MQ=255)                                          tGATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGAtt                  <  1:3385010/64‑1 (MQ=255)                                          tGATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGAtt                  <  1:29211/64‑1 (MQ=255)                                          tGATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGAtt                  <  1:2897352/64‑1 (MQ=255)                                          tGATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGAtt                  <  1:2850079/64‑1 (MQ=255)                                          tGATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGAtt                  <  1:1104304/64‑1 (MQ=255)                                          tGATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGAtt                  <  1:3392692/64‑1 (MQ=255)                                          tGATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGAtt                  <  1:2718744/64‑1 (MQ=255)                                          tGATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGAtt                  <  1:450899/64‑1 (MQ=255)                                          tGATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGAtt                  <  1:1906844/64‑1 (MQ=255)                                          tGATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGAtt                  <  1:1888614/64‑1 (MQ=255)                                          tGATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGAtt                  <  1:1815304/64‑1 (MQ=255)                                          tGATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGAtt                  <  1:2131257/64‑1 (MQ=255)                                          tGATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGAtt                  <  1:1823782/64‑1 (MQ=255)                                           gATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca            >  1:1045587/1‑69 (MQ=255)                                           gATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca            >  1:266202/1‑69 (MQ=255)                                           gATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca            >  1:2003705/1‑69 (MQ=255)                                           gATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca            >  1:3141251/1‑69 (MQ=255)                                           gATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca            >  1:955619/1‑69 (MQ=255)                                           gATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca            >  1:594836/1‑69 (MQ=255)                                           gATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca            >  1:216827/1‑69 (MQ=255)                                           gATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca            >  1:944942/1‑69 (MQ=255)                                           gATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca            >  1:578339/1‑69 (MQ=255)                                           gATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca            >  1:542250/1‑69 (MQ=255)                                           gATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca            >  1:1366852/1‑69 (MQ=255)                                           gATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca            >  1:2597019/1‑69 (MQ=255)                                           gATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca            >  1:501843/1‑69 (MQ=255)                                            aTGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca            >  1:109581/1‑68 (MQ=255)                                            aTGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca            >  1:3205535/1‑68 (MQ=255)                                            aTGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca            >  1:83824/1‑68 (MQ=255)                                            aTGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca            >  1:3147771/1‑68 (MQ=255)                                            aTGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca            >  1:305351/1‑68 (MQ=255)                                            aTGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca            >  1:1490780/1‑68 (MQ=255)                                            aTGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca            >  1:1590745/1‑68 (MQ=255)                                            aTGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca            >  1:2153467/1‑68 (MQ=255)                                            aTGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca            >  1:267305/1‑68 (MQ=255)                                            aTGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca            >  1:2670166/1‑68 (MQ=255)                                            aTGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca            >  1:2627907/1‑68 (MQ=255)                                               cgcAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCgccagcc         >  1:2410956/1‑68 (MQ=255)                                               cgcAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCgccagcc         >  1:1372812/1‑68 (MQ=255)                                               cgcAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCgccagcc         >  1:1828444/1‑68 (MQ=255)                                               cgcAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCgccagcc         >  1:3185222/1‑68 (MQ=255)                                               cgcAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCgccagcc         >  1:2978034/1‑68 (MQ=255)                                               cgcAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCgcca            >  1:2272316/1‑65 (MQ=255)                                                     tGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCgccagcca        >  1:2402918/1‑63 (MQ=255)                                                     tGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAg    >  1:1267775/1‑67 (MQ=255)                                                     tGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGc   >  1:1343955/1‑68 (MQ=255)                                                     tGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCg  >  1:2458978/1‑69 (MQ=255)                                                     tGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCg  >  1:2632556/1‑69 (MQ=255)                                                     tGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCg  >  1:1868184/1‑69 (MQ=255)                                                     tGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCg  >  1:3393068/1‑69 (MQ=255)                                                     tGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCg  >  1:1620233/1‑69 (MQ=255)                                                     tGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCg  >  1:2784999/1‑69 (MQ=255)                                                     tGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCg  >  1:1535811/1‑69 (MQ=255)                                                     tGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCg  >  1:850942/1‑69 (MQ=255)                                                     tGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCg  >  1:3312044/1‑69 (MQ=255)                                                     tGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATAGCGCCAGCCAATAg    >  1:2695544/1‑67 (MQ=255)                                                           |                                                               CTCTCATCAGGAGCGGAAGCTTTTCAGGCTTACAGTTTGAGTGATGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCG  >  minE/1968988‑1969108

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