Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,328	C→G	26.1%	Q642H (CAG→CAC) ‡	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,328	0	C	G	26.1%	107.8 / 35.8	65	Q642H (CAG→CAC) ‡	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base C (30/18);  new base G (11/6);  total (41/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGT  >  minE/608259‑608396                                                                      |                                                                     ggAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCt                                                                     >  1:415944/1‑71 (MQ=255)       aCCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCAc                                                               >  1:5691076/1‑71 (MQ=255)        ccATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCAc                                                               >  1:3762863/1‑70 (MQ=255)        ccATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCAc                                                               >  1:3321072/1‑70 (MQ=255)        ccATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTgg                                                                   <  1:1633976/66‑1 (MQ=255)        ccATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCAc                                                               >  1:526394/1‑70 (MQ=255)         cATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGa                                                             >  1:196632/1‑71 (MQ=255)         cATCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGa                                                             >  1:489970/1‑71 (MQ=255)          aTCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGAc                                                            <  1:1258203/71‑1 (MQ=255)           tCTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGATCACGAc                                                            >  1:432341/1‑70 (MQ=255)            cTTTGGCGTAGGTGAACGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGAcca                                                          >  1:5093733/1‑70 (MQ=255)              ttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGt                                                         >  1:3836571/1‑70 (MQ=255)              ttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg                                                        >  1:1192396/1‑71 (MQ=255)              ttAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCCTg                                                        >  1:4201025/1‑71 (MQ=255)               tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGAcc                                                           <  1:635675/67‑1 (MQ=255)               tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGt                                                         <  1:2513253/69‑1 (MQ=255)               tAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTg                                                        >  1:3385579/1‑70 (MQ=255)                  cGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATc                                                    <  1:3595044/71‑1 (MQ=255)                   gTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCg                                                   >  1:5471197/1‑71 (MQ=255)                    tAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCg                                                   <  1:2509113/70‑1 (MQ=255)                      aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGAc                                                            >  1:4945308/1‑59 (MQ=255)                      aGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGAc                                                            >  1:3140718/1‑59 (MQ=255)                      aGCAAACGGCAGTTTAGAAACGGAAGTCAGAGAGGCTACTCCACCTTTCTGGTAACGACCGTGCATCGGGt                                                >  1:4693906/1‑71 (MQ=255)                          aaCGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGc                                            >  1:2058665/1‑71 (MQ=255)                           aCGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGcg                                           >  1:321985/1‑70 (MQ=255)                           aCGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGcg                                           <  1:945225/71‑2 (MQ=255)                           aCGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGcg                                           >  1:3138578/1‑70 (MQ=255)                           aCGGCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGcg                                           >  1:2844125/1‑70 (MQ=255)                           aCGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGt                                                                  >  1:2937546/1‑48 (MQ=255)                           aCGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGc                                            >  1:1711072/1‑70 (MQ=255)                             ggCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGc                                                       <  1:3407023/57‑1 (MQ=255)                             ggCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGc                                                       <  1:2386775/57‑1 (MQ=255)                              gCAGTTTCGCCACCGACGTCAATGAGGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGTTAGCGccc                                        <  1:1432672/71‑1 (MQ=255)                                aGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCgg                                      <  1:3073683/71‑1 (MQ=255)                                aGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCgg                                      <  1:1762056/71‑1 (MQ=255)                                      gCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGaa                                <  1:2120626/71‑1 (MQ=255)                                       cAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCAc                                                               >  1:1451238/1‑39 (MQ=255)                                          cGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgc                            <  1:5113031/71‑1 (MQ=255)                                           ggAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCAccc                                        >  1:2368787/1‑58 (MQ=255)                                           ggAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCAccc                                        >  1:1922208/1‑58 (MQ=255)                                            gAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgcgc                          <  1:3688099/71‑1 (MQ=255)                                              aGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgcgcgc                        >  1:980134/1‑71 (MQ=255)                                              aGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGgcgcgc                        >  1:3105067/1‑71 (MQ=255)                                                tCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGt                                                <  1:5463821/45‑1 (MQ=255)                                                    tgagGCCACGGCACCTTTGCGGTCACGACCATGCATCGGGtt                                               <  1:900613/41‑1 (MQ=255)                                                    agagGCTACTGCACCTTTCTGGTCACGACCGTGCAt                                                     >  1:2995951/1‑36 (MQ=255)                                                    agagGCTACTGCACCTTTCTGGTCACGACCGTGCAt                                                     >  1:3128837/1‑36 (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)                                                         cTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAg                                             <  1:1676881/39‑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:5213759/1‑71 (MQ=255)                                                               cACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCt       >  1:5464753/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)                                                                     tCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGgt  <  1:2347125/70‑1 (MQ=255)                                                                      cTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCacga                >  1:2409864/1‑55 (MQ=255)                                                                      |                                                                     GGAGATACCATCTTTAGCGTAAGCAAACGGCAGTTTAGCAACGGAAGTCAGAGAGGCTACTGCACCTTTCTGGTCACGACCGTGCATCGGGTTAGCACCCGGTCCGAACGGCGCGCCAGCACGACGACCGTCTGGGGT  >  minE/608259‑608396

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