Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	383,184	C→T	28.6%	G198G (GGC→GGT)	entB →	isochorismatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	383,184	0	C	T	28.6%	76.4 / 29.4	73	G198G (GGC→GGT)	entB	isochorismatase
Reads supporting (aligned to +/- strand):  ref base C (20/32);  new base T (9/12);  total (29/44)
Fisher's exact test for biased strand distribution p-value = 7.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCTGGCCGATTTCAGCCGTGACGAGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCGGGTGGTGATGACTGAAGAATTACTGCCAGCACCTATCCCCGCCAGCAAAGCGGCGCTGCGTG  >  minE/383119‑383248                                                                  |                                                                 gcgcTGGCCGATTTCAGCCGTGACGAGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCgg                                                               >  1:976623/1‑69 (MQ=255) gcgcTGGCCGATTTCAGCCGTGACGAGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCgg                                                               >  1:741483/1‑69 (MQ=255) gcgcTGGCCGATTTCAGCCGTGACGAGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCgg                                                               >  1:2590472/1‑69 (MQ=255) gcgcTGGCCGATTTCAGCCGTGACGAGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCgg                                                               >  1:2313032/1‑69 (MQ=255) gcgcTGGCCGATTTCAGCCGTGACGAGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCgg                                                               >  1:2102381/1‑69 (MQ=255) gcgcTGGCCGATTTCAGCCGTGACGAGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCgg                                                               >  1:1923575/1‑69 (MQ=255)      ggCCGATTTCAGCCGTGACGAGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCGggtggt                                                          <  1:2545391/69‑1 (MQ=255)      ggCCGATTTCAGCCGTGACGAGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCGggtggt                                                          <  1:2766513/69‑1 (MQ=255)      ggCCGATTTCAGCCGTGACGAGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCGggtggt                                                          <  1:1376609/69‑1 (MQ=255)      ggCCGATTTCAGCCGTGACGAGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCGggtggt                                                          <  1:1226014/69‑1 (MQ=255)      ggCCGATTTCAGCCGTGACGAGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCGggtggt                                                          <  1:1161172/69‑1 (MQ=255)      ggCCGATTTCAGACGTGACGAGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCGggtggt                                                          <  1:2958660/69‑1 (MQ=255)           aTTTCAGCCGTGACGAGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCGggtggt                                                          <  1:2880441/64‑1 (MQ=255)              tCAGCCGTGACGAGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCGGGTGGtgat                                                       >  1:1170904/1‑64 (MQ=255)               cAGCCGTGACGAGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCGggtggt                                                          <  1:2936996/60‑1 (MQ=255)               cAGCCGTGACGAGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCGggtggt                                                          <  1:554891/60‑1 (MQ=255)                  ccGTGACGAGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCGGGTGGTGATGACTgaag                                               <  1:1505636/68‑1 (MQ=255)                  ccGTGACGAGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCGGGTGGTGATGACTgaag                                               <  1:1466331/68‑1 (MQ=255)                  ccGTGACGAGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCGGGTGGTGATGACTgaag                                               <  1:534253/68‑1 (MQ=255)                  ccGTGACGAGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCGGGTGGTGATGACTgaag                                               <  1:241062/68‑1 (MQ=255)                  ccGTGACGAGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCGGGTGGTGATGACTgaag                                               <  1:1443666/68‑1 (MQ=255)                  ccGTGACGAGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCGGGTGGTGATGACTgaag                                               <  1:1056975/68‑1 (MQ=255)                  ccGTGACGAGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCGGGTGGTGATGACTgaag                                               <  1:837224/68‑1 (MQ=255)                  ccGTGACGAGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCGGGTGGTGATGACTgaag                                               <  1:3269726/68‑1 (MQ=255)                       aCGAGCATTTGATGTCTCTGAAATATGTGGCCGGACGTTCTGGCCGGGTGGTGATGACTGAAGAATTAc                                         <  1:1626493/69‑1 (MQ=255)                       aCGAGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCGGGTGGTGATGACTGAAGAATTAc                                         <  1:120679/69‑1 (MQ=255)                       aCGAGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCGGGTGGTGATGACTGAAGAATTAc                                         <  1:2968321/69‑1 (MQ=255)                       aCGAGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCGGGTGGTGATGACTGAAGAATTAc                                         >  1:189243/1‑69 (MQ=255)                       aCGAGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCGGGTGGTGATGACTGAAGAATTAc                                         <  1:1287238/69‑1 (MQ=255)                        cGAGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCGGGTGGTGATGACTGAAGAATTAc                                         <  1:2704832/68‑1 (MQ=255)                          aGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGTCGGGTGGTGATGACCGAAGAATtgctgc                                      <  1:2925502/69‑1 (MQ=255)                          aGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGTCGGGTGGTGATGACCGAAGAATtgctgc                                      <  1:2914172/69‑1 (MQ=255)                          aGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGTCGGGTGGTGATGACCGAAGAATtgctgc                                      <  1:3342233/69‑1 (MQ=255)                          aGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGTCGGGTGGTGATGACCGAAGAATtgctgc                                      <  1:3364581/69‑1 (MQ=255)                          aGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGTCGGGTGGTGATGACCGAAGAATtgctgc                                      <  1:473388/69‑1 (MQ=255)                          aGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGTCGGGTGGTGATGACCGAAGAATtgctgc                                      <  1:2815483/69‑1 (MQ=255)                          aGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGTCGGGTGGTGATGACCGAAGAATtgctgc                                      <  1:543689/69‑1 (MQ=255)                          aGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGTCGGGTGGTGATGACCGAAGAATtgctgc                                      <  1:563857/69‑1 (MQ=255)                          aGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGTCGGGTGGTGATGACCGAAGAATtgctgc                                      <  1:2718425/69‑1 (MQ=255)                          aGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGTCGGGTGGTGATGACCGAAGAATtgctgc                                      <  1:1068202/69‑1 (MQ=255)                          aGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGTCGGGTGGTGATGACCGAAGAATtgctgc                                      <  1:2611146/69‑1 (MQ=255)                                tGATGTCGCTGAAATATGTGGCCGGACGTTCTGGcc                                                                 >  1:473898/1‑36 (MQ=255)                                tGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCGGGTGGTGATGACTGAAGAATTACTGCCAGCa                                 >  1:3406174/1‑68 (MQ=255)                                tGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCGGGTGGTGATGACTGAAGAATTACTGCCAGCa                                 >  1:899095/1‑68 (MQ=255)                                tGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCGGGTGGTGATGACTGAAGAATTACTGCCAGCa                                 >  1:1264764/1‑68 (MQ=255)                                tGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCGGGTGGTGATGACTGAAGAATTACTGCCAGCa                                 >  1:223125/1‑68 (MQ=255)                                tGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCGGGTGGTGATGACTGAAGAATTACTGCCAGCa                                 >  1:2122714/1‑68 (MQ=255)                                tGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCGGGTGGTGATGACTGAAGAATTACTGCCAGCa                                 >  1:3216993/1‑68 (MQ=255)                                tGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCGGGTGGTGATGACTGAAGAATTACTGCCAGCa                                 >  1:3189349/1‑68 (MQ=255)                                tGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCGGGTGGTGATGACTGAAGAATTACTGCCAGCa                                 >  1:1654657/1‑68 (MQ=255)                                     tCGCTGAAATATGTGGCCGGACGTTCTGGCCGGGTGGTGATGACTGAAGAATTACTGCCAg                                   <  1:1270687/61‑1 (MQ=255)                                     tCGCTGAAATATGTGGCCGGACGTTCTGGCCGGGTGGTGATGACTGAAGAATTACTGCCAg                                   <  1:3331056/61‑1 (MQ=255)                                     tCGCTGAAATATGTGGCCGGACGTTCTGGCCGGGTGGTGATGACTGAAGAATTACTGCCAg                                   <  1:1589192/61‑1 (MQ=255)                                     tCGCTGAAATATGTGGCCGGACGTTCTGGCCGGGTGGTGATGACTGAAGAATTACTGCCAg                                   <  1:209485/61‑1 (MQ=255)                                     tCGCTGAAATATGTGGCCGGACGTTCTGGCCGGGTGGTGATGACTGAAGAATTACTGCCAg                                   <  1:2463157/61‑1 (MQ=255)                                                     ccGGACGTTCTGGTCGGGTGGTGATGACCGAAGAATTGCTGCCAGCACCTATCCCCGCCAGCAAAgcgg           >  1:1514233/1‑69 (MQ=255)                                                     ccGGACGTTCTGGTCGGGTGGTGATGACCGAAGAATTGCTGCCAGCACCTATCCCCGCCAGCAAAgcgg           >  1:1745117/1‑69 (MQ=255)                                                     ccGGACGTTCTGGTCGGGTGGTGATGACCGAAGAATTGCTGCCAGCACCTATCCCCGCCAGCAAAgcgg           >  1:1114048/1‑69 (MQ=255)                                                     ccGGACGTTCTGGTCGGGTGGTGATGACCGAAGAATTGCTGCCAGCACCTATCCCCGCCAGCAAAgcgg           >  1:1955091/1‑69 (MQ=255)                                                     ccGGACGTTCTGGTCGGGTGGTGATGACCGAAGAATTGCTGCCAGCACCTATCCCCGCCAGCAAAgcgg           >  1:3362777/1‑69 (MQ=255)                                                     ccGGACGTTCTGGTCGGGTGGTGATGACCGAAGAATTGCTGCCAGCACCTATCCCCGCCAGCAAAgcgg           >  1:3141362/1‑69 (MQ=255)                                                     ccGGACGTTCTGGTCGGGTGGTGATGACCGAAGAATTGCTGCCAGCACCTATCCCCGCCAGCAAAgcgg           >  1:3005018/1‑69 (MQ=255)                                                     ccGGACGTTCTGGTCGGGTGGTGATGACCGAAGAATTGCTGCCAGCACCTATCCCCGCCAGCAAAgcgg           >  1:2771450/1‑69 (MQ=255)                                                     ccGGACGTTCTGGTCGGGTGGTGATGACCGAAGAATTGCTGCCAGCACCTATCCCCGCCAGCAAAgcgg           >  1:1168521/1‑69 (MQ=255)                                                      cGGACGTTCTGGCCGGGTGGTGATGACTGAAGAATTACTGCCAGCACCTATCCCCGCCAGCAAAgcggc          <  1:1943673/69‑1 (MQ=255)                                                      cGGACGTTCTGGCCGGGTGGTGATGACTGAAGAATTACTGCCAGCACCTATCCCCGCCAGCAAAgcggc          <  1:448800/69‑1 (MQ=255)                                                      cGGACGTTCTGGCCGGGTGGTGATGACTGAAGAATTACTGCCAGCACCTATCCCCGCCAGCAAAgcggc          <  1:2061625/69‑1 (MQ=255)                                                      cGGACGTTCTGGCCGGGTGGTGATGACTGAAGAATTACTGCCAGCACCTATCCCCGCCAGCAAAgcggc          <  1:778011/69‑1 (MQ=255)                                                        gACGTTCTGGTCGGGTGGTGATGACCGAAGAATtgctgc                                      <  1:2412802/39‑1 (MQ=37)                                                            ttCTGGCCGGGTGGTGATGACTGAAGAATTACTGCCAGCACCTATCCCCGCCAGCAAAGCGGCGCTGc     >  1:2410460/1‑68 (MQ=255)                                                            ttCTGGCCGGGTGGTGATGACTGAAGAATTACTGCCAGCACCTATCCCCGCCAGCAAAGCGGCGCTGc     >  1:1376368/1‑68 (MQ=255)                                                            ttCTGGCCGGGTGGTGATGACTGAAGAATTACTGCCAGCACCTATCCCCGCCAGCAAAGCGGCGCTGc     >  1:951206/1‑68 (MQ=255)                                                            ttCTGGCAGGGTGGTGATGACTGAAGAATTACTGCCAGCACCTATCCCCGCCAGCAAAGCGGCGCTGc     >  1:2025358/1‑68 (MQ=255)                                                              cTGGCCGGGTGGTGATGACTGAAGAATTACTGCCAGCACCTATCCCCGCCAGCAAAGCGGCGCTGCGTg  <  1:2017540/69‑1 (MQ=255)                                                                  |                                                                 GCGCTGGCCGATTTCAGCCGTGACGAGCATTTGATGTCGCTGAAATATGTGGCCGGACGTTCTGGCCGGGTGGTGATGACTGAAGAATTACTGCCAGCACCTATCCCCGCCAGCAAAGCGGCGCTGCGTG  >  minE/383119‑383248

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑