Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,094	C→G	15.0%	M720I (ATG→ATC)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,094	0	C	G	15.0%	155.5 / 17.2	60	M720I (ATG→ATC)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base C (29/22);  new base G (2/7);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 7.56e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.09e-01

AGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTT  >  minE/608026‑608164                                                                     |                                                                       aGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCAt                                                                      >  1:1828542/1‑71 (MQ=255)         gAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTc                                                                >  1:2861031/1‑69 (MQ=255)         gAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCg                                                               <  1:2574956/70‑1 (MQ=255)          aaaCGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCCTCGag                                                             <  1:589931/71‑1 (MQ=255)           aaCGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCg                                                               >  1:1936877/1‑68 (MQ=255)            aCGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGag                                                             >  1:566372/1‑69 (MQ=255)            aCGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGag                                                             >  1:529629/1‑69 (MQ=255)              gTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGagca                                                           >  1:3640/1‑69 (MQ=255)                aCTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCAt                                                       >  1:3378394/1‑71 (MQ=255)                 cggCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCAtt                                                      <  1:1265563/69‑1 (MQ=12)                 cggCGTAGCCAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCAtt                                                      <  1:3645971/69‑1 (MQ=12)                   gCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCAttt                                                     >  1:2362056/1‑70 (MQ=255)                    cGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCa                                                   >  1:2963865/1‑71 (MQ=255)                        gCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGt                                               <  1:3744660/71‑1 (MQ=255)                          cAGAGACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCga                                                <  1:2729242/68‑2 (MQ=9)                           agagACACGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCg                                                 <  1:263141/66‑1 (MQ=9)                           aGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtca                                            <  1:385533/71‑1 (MQ=255)                           aGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtca                                            >  1:396873/1‑71 (MQ=255)                            gagacacGGATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATtcat                                           <  1:803100/68‑1 (MQ=9)                            gATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCAttt                                                     <  1:303918/61‑1 (MQ=255)                             aTACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcat                                           >  1:3380607/1‑70 (MQ=255)                               acacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc                                        >  1:1149639/1‑71 (MQ=255)                                cacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc                                        <  1:251882/70‑1 (MQ=255)                                 acGGATGGTCAGCTGCGGAAATTTTTCCGGGTTTTCCATCGCGTCGagcagc                                                         >  1:25087/1‑52 (MQ=255)                                   ggATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGttg                                    <  1:1812077/71‑2 (MQ=255)                                   ggATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGtt                                     >  1:2612834/1‑70 (MQ=255)                                        gTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCAttt                                                     <  1:731714/49‑1 (MQ=255)                                         tCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCAc                                                  <  1:160098/51‑1 (MQ=255)                                         tCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTc                              <  1:928150/71‑1 (MQ=255)                                          cAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGtt                               >  1:254837/1‑69 (MQ=255)                                           aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                            >  1:1010417/1‑71 (MQ=255)                                           aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                            <  1:3241491/71‑1 (MQ=255)                                           aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                            >  1:3155460/1‑71 (MQ=255)                                            gCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcatc                                          >  1:2916952/1‑56 (MQ=255)                                            gCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCa                             <  1:1912623/69‑1 (MQ=255)                                             cTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGt                          <  1:2538351/71‑1 (MQ=255)                                             cTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGt                          >  1:40621/1‑71 (MQ=255)                                             cTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGt                          <  1:3583132/71‑1 (MQ=255)                                               taggATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTgt                        >  1:1934025/3‑70 (MQ=255)                                               gCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGc                        >  1:3690568/1‑71 (MQ=255)                                                 ggATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCAc                                                  >  1:2967272/1‑43 (MQ=255)                                                  gATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGa                     <  1:1929460/71‑1 (MQ=255)                                                  gATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGa                     >  1:616967/1‑71 (MQ=12)                                                     atTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGt                                      <  1:1328372/51‑1 (MQ=255)                                                     atTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGacca                  >  1:564393/1‑71 (MQ=255)                                                      tttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccac                 <  1:2149654/71‑1 (MQ=255)                                                       ttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccacc                >  1:855116/1‑71 (MQ=255)                                                       ttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccacc                >  1:534188/1‑71 (MQ=255)                                                       ttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccacc                >  1:3251582/1‑71 (MQ=255)                                                        tttCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTGaccgcc                <  1:3278536/70‑1 (MQ=255)                                                            cGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCga           <  1:1278827/71‑1 (MQ=255)                                                            cGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCga           <  1:868636/71‑1 (MQ=255)                                                             gggTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcat                                           <  1:558239/38‑1 (MQ=255)                                                             gggTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgat          <  1:1627687/71‑1 (MQ=255)                                                               gTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatg         >  1:1553947/1‑70 (MQ=255)                                                                ttttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                            <  1:237535/50‑1 (MQ=255)                                                                ttttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAg                            <  1:1953432/50‑1 (MQ=255)                                                                ttttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatgga       >  1:1659397/1‑71 (MQ=255)                                                                  ttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatggat      >  1:1977977/1‑70 (MQ=255)                                                                   tCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGc    >  1:3689004/1‑71 (MQ=255)                                                                    ccATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTGAGGTGCTGACCACCTTCGATGGATGc    <  1:3419300/70‑1 (MQ=255)                                                                    ccATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGc    <  1:1213134/70‑1 (MQ=255)                                                                     cATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAgg                           <  1:904983/46‑1 (MQ=255)                                                                     cATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCtt  >  1:3171065/1‑71 (MQ=255)                                                                     |                                                                       AGCGAGTTGAAACGTACTGCGTAGCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTT  >  minE/608026‑608164

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