Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,373	G→T	24.9%	G627G (GGC→GGA)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,373	0	G	T	24.9%	105.7 / 22.3	52	G627G (GGC→GGA)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base G (22/17);  new base T (5/8);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 3.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.91e-01

AGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAA  >  minE/608304‑608443                                                                      |                                                                       aGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgcgcgc                                                                       >  1:980134/1‑71 (MQ=255) aGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgcgcgc                                                                       >  1:3105067/1‑71 (MQ=255)       agagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCac                                                                 <  1:2490337/71‑1 (MQ=255)       agagGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCac                                                                 <  1:1399283/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:4138635/1‑71 (MQ=255)                       ttCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGgt                                                 >  1:2646935/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:3947601/1‑71 (MQ=255)                            gTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTAc                                            >  1:5106630/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:1515024/71‑1 (MQ=255)                                 gACCATGCATCGGGTTAGCGCCCGGCGCGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGtt                                       <  1:3841163/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:3273572/71‑1 (MQ=255)                                                   cACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACg                     >  1:5220220/1‑71 (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)                                                           ccGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTAcc                               >  1:321984/1‑53 (MQ=255)                                                            cGAACGGTGTTCCGGCGCGACGACCGTCCGGCGTATTACCGGTTTTCTGGCCGTACACCACGTTAGAAGt             >  1:3633108/1‑70 (MQ=255)                                                            cGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTg            <  1:899285/71‑1 (MQ=255)                                                              aaCGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg                 <  1:5010659/64‑1 (MQ=255)                                                              aaCGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAg                 <  1:2862864/64‑1 (MQ=255)                                                              aaCGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGAt          >  1:3351370/1‑71 (MQ=255)                                                                cGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATgg        >  1:3361145/1‑71 (MQ=255)                                                                  gcgcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGt       >  1:1669343/1‑70 (MQ=255)                                                                  gcgcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGt       >  1:2058523/1‑70 (MQ=255)                                                                  gcgcgcCAGCACGACGACCGTCTGGGGTGTTAACCGTTTTCTTACCATACACAACGTTAGAAGTGATGGt       >  1:2981590/1‑70 (MQ=255)                                                                    gcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATaca                          <  1:2340756/49‑1 (MQ=255)                                                                    gcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATaca                          <  1:3983295/49‑1 (MQ=255)                                                                    gcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGt                    <  1:210484/55‑1 (MQ=255)                                                                    gcgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGt                    <  1:4115035/55‑1 (MQ=255)                                                                     cgcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAg    >  1:2490241/1‑70 (MQ=255)                                                                      gcCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGaa  >  1:781163/1‑71 (MQ=255)                                                                      |                                                                       AGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGTGTTACCCGTTTTCTTACCATACACAACGTTAGAAGTGATGGTCAGAA  >  minE/608304‑608443

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