Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	1,969,065	T→C	45.2%	C53C (TGT→TGC)	yggG →	predicted peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	1,969,065	0	T	C	45.2%	10.3 / 69.9	84	C53C (TGT→TGC)	yggG	predicted peptidase
Reads supporting (aligned to +/- strand):  ref base T (34/12);  new base C (24/14);  total (58/26)
Fisher's exact test for biased strand distribution p-value = 3.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

AGCTTTTCAGGCTTACAGTTTGAGTGATGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAACTAT  >  minE/1969005‑1969133                                                             |                                                                     aGCTTTTCAGGCTTACAGTTTGAGTGATGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGgaga                                                               >  1:821692/1‑68 (MQ=255) aGCTTTTCAGGCTTACAGTTTGAGTGATGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGgaga                                                               >  1:2770788/1‑68 (MQ=255)    ttttCAGGCTTACAGTTTGAGTGATGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATg                                                             >  1:2358568/1‑67 (MQ=255)                 aGTTTGAGTGATGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGcgac                                               <  1:1968658/68‑1 (MQ=255)                 aGTTTGAGTGATGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGcgac                                               <  1:1990900/68‑1 (MQ=255)                 aGTTTGAGTGATGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGcgac                                               <  1:1435413/68‑1 (MQ=255)                     tGAGTGATGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGga                                                                 <  1:2489271/46‑1 (MQ=255)                     tGAGTGATGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGga                                                                 <  1:3237328/46‑1 (MQ=255)                     tGAGTGATGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGga                                                                 <  1:781799/46‑1 (MQ=255)                     tGAGTGATGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGga                                                                 <  1:1240331/46‑1 (MQ=255)                      gAGTGATGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGga                                                                 <  1:862186/45‑1 (MQ=255)                         tGATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGAtt                                           <  1:2850079/64‑1 (MQ=255)                         tGATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGAtt                                           <  1:2131257/64‑1 (MQ=255)                         tGATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGAtt                                           <  1:2897352/64‑1 (MQ=255)                         tGATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGAtt                                           <  1:29211/64‑1 (MQ=255)                         tGATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGAtt                                           <  1:2718744/64‑1 (MQ=255)                         tGATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGAtt                                           <  1:1815304/64‑1 (MQ=255)                         tGATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGAtt                                           <  1:1906844/64‑1 (MQ=255)                         tGATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGAtt                                           <  1:1888614/64‑1 (MQ=255)                         tGATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGAtt                                           <  1:2991797/64‑1 (MQ=255)                         tGATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGAtt                                           <  1:1823782/64‑1 (MQ=255)                         tGATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGAtt                                           <  1:3385010/64‑1 (MQ=255)                         tGATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGAtt                                           <  1:3392692/64‑1 (MQ=255)                         tGATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGAtt                                           <  1:450899/64‑1 (MQ=255)                         tGATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGAtt                                           <  1:1104304/64‑1 (MQ=255)                          gATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca                                     >  1:3141251/1‑69 (MQ=255)                          gATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca                                     >  1:266202/1‑69 (MQ=255)                          gATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca                                     >  1:501843/1‑69 (MQ=255)                          gATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca                                     >  1:542250/1‑69 (MQ=255)                          gATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca                                     >  1:578339/1‑69 (MQ=255)                          gATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca                                     >  1:594836/1‑69 (MQ=255)                          gATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca                                     >  1:944942/1‑69 (MQ=255)                          gATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca                                     >  1:1045587/1‑69 (MQ=255)                          gATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca                                     >  1:2597019/1‑69 (MQ=255)                          gATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca                                     >  1:955619/1‑69 (MQ=255)                          gATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca                                     >  1:1366852/1‑69 (MQ=255)                          gATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca                                     >  1:2003705/1‑69 (MQ=255)                          gATGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca                                     >  1:216827/1‑69 (MQ=255)                           aTGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca                                     >  1:3147771/1‑68 (MQ=255)                           aTGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca                                     >  1:109581/1‑68 (MQ=255)                           aTGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca                                     >  1:83824/1‑68 (MQ=255)                           aTGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca                                     >  1:1490780/1‑68 (MQ=255)                           aTGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca                                     >  1:1590745/1‑68 (MQ=255)                           aTGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca                                     >  1:3205535/1‑68 (MQ=255)                           aTGCGCAGGTGAAAAAACTGAGCGATCAGGCATGCCAGGAGATGGACAGCAAAGCGACGATTGCgcca                                     >  1:305351/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:3185222/1‑68 (MQ=255)                              cgcAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCgccagcc                                  >  1:1372812/1‑68 (MQ=255)                              cgcAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCgccagcc                                  >  1:2978034/1‑68 (MQ=255)                              cgcAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCgccagcc                                  >  1:2410956/1‑68 (MQ=255)                              cgcAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCgccagcc                                  >  1:1828444/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:1620233/1‑69 (MQ=255)                                    tGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCg                           >  1:3312044/1‑69 (MQ=255)                                    tGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCg                           >  1:3393068/1‑69 (MQ=255)                                    tGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCg                           >  1:1535811/1‑69 (MQ=255)                                    tGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCg                           >  1:1868184/1‑69 (MQ=255)                                    tGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCg                           >  1:2632556/1‑69 (MQ=255)                                    tGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCg                           >  1:2784999/1‑69 (MQ=255)                                    tGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCg                           >  1:2458978/1‑69 (MQ=255)                                    tGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCg                           >  1:850942/1‑69 (MQ=255)                                    tGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATAGCGCCAGCCAATAg                             >  1:2695544/1‑67 (MQ=255)                                                     cAGGCATGTCAGGAGATGGACAGCAAGGCGTCGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCt           >  1:547388/1‑68 (MQ=255)                                                     cAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCgccagcca                                 >  1:329042/1‑46 (MQ=255)                                                     cAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCt           >  1:544244/1‑68 (MQ=255)                                                     cAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCt           >  1:2220671/1‑68 (MQ=255)                                                     cAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCt           >  1:2726558/1‑68 (MQ=255)                                                        gCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGAc        <  1:2348821/68‑1 (MQ=255)                                                        gCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGAc        <  1:1523749/68‑1 (MQ=255)                                                        gCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGAc        <  1:2537264/68‑1 (MQ=255)                                                        gCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGAc        <  1:54892/68‑1 (MQ=255)                                                          aTGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACaa      >  1:2940602/1‑68 (MQ=255)                                                          aTGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACaa      >  1:1342550/1‑68 (MQ=255)                                                          aTGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAAc     >  1:2597823/1‑69 (MQ=255)                                                          aTGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAAc     >  1:472750/1‑69 (MQ=255)                                                          aTGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAAc     >  1:1214470/1‑69 (MQ=255)                                                          aTGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAAc     >  1:2067420/1‑69 (MQ=255)                                                          aTGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGACAATAGCGAATACGCTAAACGTCTGACAAc     >  1:1673884/1‑69 (MQ=255)                                                             tCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAACTAt  >  1:2721465/1‑69 (MQ=255)                                                             tCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAACTAt  >  1:2830290/1‑69 (MQ=255)                                                             tCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAACTAt  >  1:928289/1‑69 (MQ=255)                                                             tCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAACTAt  >  1:320238/1‑69 (MQ=255)                                                             tCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAACTAt  >  1:2413970/1‑69 (MQ=255)                                                             |                                                                     AGCTTTTCAGGCTTACAGTTTGAGTGATGCGCAGGTGAAAACCCTGAGCGATCAGGCATGTCAGGAGATGGACAGCAAGGCGACGATTGCGCCAGCCAATAGCGAATACGCTAAACGTCTGACAACTAT  >  minE/1969005‑1969133

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