Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,361	T→C	24.0%	G631G (GGA→GGG) ‡	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,361	0	T	C	24.0%	80.6 / 25.1	50	G631G (GGA→GGG) ‡	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base T (21/17);  new base C (4/8);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 3.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.46e-01

GCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAG  >  minE/608296‑608428                                                                  |                                                                    gCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGaa                                                                <  1:2120626/71‑1 (MQ=255)     cGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgc                                                            <  1:5113031/71‑1 (MQ=255)       gAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgcgc                                                          <  1:3688099/71‑1 (MQ=255)         aGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgcgcgc                                                        >  1:980134/1‑71 (MQ=255)         aGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgcgcgc                                                        >  1:3105067/1‑71 (MQ=255)               agagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCac                                                  <  1:1399283/71‑1 (MQ=255)               agagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCac                                                  <  1:2490337/71‑1 (MQ=255)                  ggCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCacga                                                >  1:2161674/1‑70 (MQ=255)                   gCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCacgacg                                              <  1:3342709/71‑1 (MQ=255)                       cggCACCTTTGCGGTCACGACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCg                                          >  1:616796/3‑71 (MQ=255)                        ggCACCTTTGAGGTCACGACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCg                                          <  1:276826/69‑1 (MQ=255)                         gCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCacga                                                >  1:3379603/1‑63 (MQ=255)                         gCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCacga                                                >  1:3639095/1‑63 (MQ=255)                         gCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACgacgac                                            >  1:933964/1‑67 (MQ=255)                          cACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACgacgac                                            >  1:5403555/1‑66 (MQ=255)                          cACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCt                                       >  1:5464753/1‑71 (MQ=255)                          cACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCt                                       >  1:5213759/1‑71 (MQ=255)                            ccTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCttg                                     <  1:2552236/71‑3 (MQ=255)                               ttCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAg                                                     >  1:2966184/1‑52 (MQ=255)                               ttCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGgt                                  >  1:2646935/1‑71 (MQ=255)                               ttCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGgt                                  >  1:4138635/1‑71 (MQ=255)                                tgcggTCACGACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCgt                                  <  1:1831786/67‑1 (MQ=255)                                tCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGgt                                  <  1:2347125/70‑1 (MQ=255)                                 cTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCacga                                                >  1:2409864/1‑55 (MQ=255)                                    gTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGgtgt                                <  1:502969/68‑1 (MQ=255)                                    gTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTAc                             >  1:5106630/1‑71 (MQ=255)                                    gTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTAc                             >  1:3947601/1‑71 (MQ=255)                                    gTCACGACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTAc                             <  1:1035538/71‑1 (MQ=255)                                     tCACGACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCgta                                 >  1:4677054/1‑65 (MQ=255)                                       aCGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTggg                                    >  1:928333/1‑61 (MQ=255)                                       aCGACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTAccgg                          >  1:4037640/1‑69 (MQ=255)                                       aCGACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTAccg                           <  1:2079948/70‑2 (MQ=255)                                         gACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGtt                        <  1:3841163/71‑1 (MQ=255)                                         gACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGtt                        <  1:1515024/71‑1 (MQ=255)                                          aCCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGttt                       <  1:4754005/71‑1 (MQ=255)                                             gTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCt                    <  1:2869956/71‑1 (MQ=255)                                              tGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGtt                        >  1:5672208/1‑66 (MQ=255)                                               gCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTa                  >  1:4342200/1‑71 (MQ=255)                                                cATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGTTTTCtggc                 <  1:1104124/71‑4 (MQ=255)                                                 aTCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTgg                                     >  1:2202266/1‑50 (MQ=255)                                                    gggTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATa             <  1:1332168/71‑1 (MQ=255)                                                       ttAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCg                                          <  1:1314312/39‑1 (MQ=255)                                                         aGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACaa        >  1:4628740/1‑71 (MQ=255)                                                           cACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACg      >  1:5220220/1‑71 (MQ=255)                                                           cACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACg      <  1:3273572/71‑1 (MQ=255)                                                            aCCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGt     <  1:1479695/71‑1 (MQ=255)                                                             cccGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGt     >  1:2916497/1‑70 (MQ=255)                                                               cGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg  <  1:1196942/71‑1 (MQ=255)                                                               cGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg  <  1:4386550/71‑1 (MQ=255)                                                                 gTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGttt                       <  1:5369907/48‑1 (MQ=255)                                                                  |                                                                    GCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAG  >  minE/608296‑608428

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