Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	2,128,109	C→T	35.0%	A178A (GCC→GCT)	yhbV →	predicted protease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	2,128,109	0	C	T	35.0%	39.3 / 30.2	57	A178A (GCC→GCT)	yhbV	predicted protease
Reads supporting (aligned to +/- strand):  ref base C (8/29);  new base T (9/11);  total (17/40)
Fisher's exact test for biased strand distribution p-value = 7.80e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGTAACCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCTTTACCGCGCGTTCGGAAGACCGCCCGAAAGATGAGTGTGAA  >  minE/2128049‑2128157                                                             |                                                 tCGTAACCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCtt                                          <  1:677123/69‑1 (MQ=255) tCGTAACCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCtt                                          <  1:631438/69‑1 (MQ=255) tCGTAACCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCtt                                          <  1:1397317/69‑1 (MQ=255) tCGTAACCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCtt                                          <  1:3165385/69‑1 (MQ=255) tCGTAACCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCtt                                          <  1:2726345/69‑1 (MQ=255) tCGTAACCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCtt                                          <  1:2646008/69‑1 (MQ=255) tCGTAACCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCtt                                          <  1:2331572/69‑1 (MQ=255)  cGTAACCAGTTTGAAGTGGAAGTCCTGCGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCttt                                         >  1:1244151/1‑69 (MQ=255)  cGTAACCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCttt                                         >  1:890378/1‑69 (MQ=255)  cGTAACCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCttt                                         >  1:870142/1‑69 (MQ=255)  cGTAACCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCttt                                         >  1:1578626/1‑69 (MQ=255)  cGTAACCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCttt                                         >  1:2295240/1‑69 (MQ=255)      aCCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCtt                                          <  1:2655582/64‑1 (MQ=255)      aCCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCtt                                          <  1:493118/64‑1 (MQ=255)      aCCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCtt                                          <  1:3328465/64‑1 (MQ=255)      aCCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCtt                                          <  1:3126623/64‑1 (MQ=255)      aCCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCtt                                          <  1:2232703/64‑1 (MQ=255)      aCCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCtt                                          <  1:1487865/64‑1 (MQ=255)      aCCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCtt                                          <  1:1909585/64‑1 (MQ=255)      aCCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCtt                                          <  1:1415669/64‑1 (MQ=255)               aaGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCTTTACCGCGCGTTCg                             <  1:2283154/68‑1 (MQ=255)                aGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCTTTACCGCGCGTTCgg                            <  1:2408129/68‑1 (MQ=255)                aGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCTTTACCGCGCGTTCgg                            <  1:2418542/68‑1 (MQ=255)                aGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCTTTACCGCGCGTTCgg                            <  1:399640/68‑1 (MQ=255)                aGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCTTTACCGCGCGTTCgg                            <  1:3375974/68‑1 (MQ=255)                aGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCTTTACCGCGCGTTCgg                            <  1:2696484/68‑1 (MQ=255)                aGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCTTTACCGCGCGTTCgg                            <  1:2892300/68‑1 (MQ=255)                  tGGAAGTCTTGAGCTACGGTCATCTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGtt                               >  1:65518/1‑63 (MQ=255)                  tGGAAGTCTTGAGCTACGGTCATCTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGTTc                              >  1:30204/1‑64 (MQ=255)                  tGGAAGTCTTGAGCTACGGTCATCTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGTTCGGAAg                         >  1:437275/1‑69 (MQ=255)                  tGGAAGTCTTGAGCTACGGTCATCTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGTTCGGAAg                         >  1:3233521/1‑69 (MQ=255)                  tGGAAGTCTTGAGCTACGGTCATCTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGTTCGGAAg                         >  1:1936838/1‑69 (MQ=255)                  tGGAAGTCTTGAGCTACGGTCATCTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGTTCGGAAg                         >  1:872816/1‑69 (MQ=255)                  tGGAAGTCTTGAGCTACGGTCATCTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGTTCGGAAg                         >  1:401645/1‑69 (MQ=255)                  tGGAAGTCTTGAGCTACGGTCATCTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGTTCGGAAg                         >  1:963389/1‑69 (MQ=255)                  tGGAAGTCTTGAGCGACGGTCATCTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGTTCGGAAg                         >  1:3199274/1‑69 (MQ=255)                    gAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCTTTACCGCGCGTTCgg                            >  1:1922316/1‑64 (MQ=255)                          cTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCTTTACCGCGCGTTCGGAAg                         >  1:1400269/1‑61 (MQ=255)                          cTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCTTTACCGCGCGTTCGGAAGACCGccc                  >  1:1025287/1‑68 (MQ=255)                               caaCGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCTTTACCGCGCGTTCGGAAGACCGCCCGaaa              <  1:2445872/65‑1 (MQ=255)                                     tCATCTGCCGCTGGCCTACTCCGCCCGCTGCTTTACCGCGCGTTCg                             <  1:966459/46‑1 (MQ=255)                                         cTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGTTCGGaa                          <  1:671184/45‑1 (MQ=38)                                         cTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGTTCGGaa                          <  1:1333154/45‑1 (MQ=38)                                         cTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGTTCGGaa                          <  1:1366254/45‑1 (MQ=38)                                         cTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGTTCGGaa                          <  1:81154/45‑1 (MQ=38)                                         cTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGTTCGGaa                          <  1:1323364/45‑1 (MQ=38)                                         cTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGTTCGGaa                          <  1:3166482/45‑1 (MQ=38)                                         cTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGTTCGGaa                          <  1:2635930/45‑1 (MQ=38)                                         cTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGTTCGGaa                          <  1:2582118/45‑1 (MQ=38)                                         cTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGTTCGGaa                          <  1:2549651/45‑1 (MQ=38)                                         cTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGTTCGGaa                          <  1:2464271/45‑1 (MQ=38)                                         cTGCCGCTGGCCTACTCTGCTCGCTGCTTTACCGCGCGTTCGGaa                          <  1:2392731/45‑1 (MQ=38)                                         cTGCCGCTGGCCTACTCCGCCCGCTGCTTTACCGCGCGTTCGGAAGACCGCCCGAAAGATGAGTGTGaa  <  1:3354134/69‑1 (MQ=255)                                         cTGCCGCTGGCCTACTCCGCCCGCTGCTTTACCGCGCGTTCGGAAGACCGCCCGAAAGATGAGTGTGaa  <  1:3247407/69‑1 (MQ=255)                                         cTGCCGCTGGCCTACTCCGCCCGCTGCTTTACCGCGCGTTCGGAAGACCGCCCGAAAGATGAGTGTGaa  <  1:2986948/69‑1 (MQ=255)                                         cTGCCGCTGGCCTACTCCGCCCGCTGCTTTACCGCGCGTTCGGAAGACCGCCCGAAAGATGAGTGTGaa  <  1:1290604/69‑1 (MQ=255)                                         cTGCCGCTGGCCTACTCCGCCCGCTGCTTTACCGCGCGTTCGGAAGACCGCCCGAAAGATGAGTGTGaa  <  1:2274914/69‑1 (MQ=255)                                                             |                                                 TCGTAACCAGTTTGAAGTGGAAGTCCTGAGCTACGGTCATCTGCCGCTGGCCTACTCCGCCCGCTGCTTTACCGCGCGTTCGGAAGACCGCCCGAAAGATGAGTGTGAA  >  minE/2128049‑2128157

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