Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,115	A→C	7.6%	R713R (CGT→CGG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,115	0	A	C	7.6%	202.8 / 12.9	79	R713R (CGT→CGG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand):  ref base A (36/37);  new base C (3/3);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

GCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATC  >  minE/608049‑608183                                                                   |                                                                     gCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACgg                                                                   <  1:3759950/70‑1 (MQ=255) gCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACgg                                                                   <  1:4059636/70‑1 (MQ=255) gCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACgg                                                                   <  1:2782407/70‑1 (MQ=255) gCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACgg                                                                   <  1:1554197/70‑1 (MQ=255) gCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGt                                                                  >  1:2263107/1‑71 (MQ=255)  ccAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGtt                                                                 >  1:678018/1‑71 (MQ=255)       tACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCAc                                                                     >  1:2117928/1‑62 (MQ=255)       tACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcatca                                                            >  1:2134546/1‑71 (MQ=255)         cacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACg                                                                    <  1:3307487/61‑1 (MQ=255)         cacGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc                                                           <  1:954018/70‑1 (MQ=255)           cGGATGGTCAGCT‑‑GGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGt                                                                  >  1:3782096/1‑59 (MQ=38)            ggATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTa                                                       <  1:2867888/71‑1 (MQ=255)            ggATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTa                                                       <  1:2862875/71‑1 (MQ=255)             gATTGTCAGGTTAGGATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGttga                                                      <  1:2025127/71‑3 (MQ=255)             gATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTaa                                                      <  1:3732417/71‑1 (MQ=255)                 gTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcatc                                                             >  1:2968005/1‑60 (MQ=255)                 gTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGtt                                                  <  1:1895377/71‑1 (MQ=255)                  tCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTc                                                 <  1:2209672/71‑1 (MQ=255)                   cAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGtt                                                  >  1:341219/1‑69 (MQ=255)                    aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACgg                                                                   >  1:873134/1‑51 (MQ=255)                    aGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCa                                                <  1:428259/70‑1 (MQ=255)                        taggATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTgt                                           >  1:2042885/3‑70 (MQ=255)                         cGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcat                                                              <  1:3068170/51‑1 (MQ=255)                         cGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACg                                                    <  1:1592085/61‑1 (MQ=255)                         cGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCt                                          <  1:2993827/71‑1 (MQ=255)                         aggATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTgtt                                          >  1:1240154/2‑69 (MQ=255)                          ggATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGATTCATTACGTTGACGTTGAGGTGTTg                                         <  1:808871/71‑1 (MQ=12)                           gATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGtt                                                                 <  1:3237633/46‑1 (MQ=255)                           gATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGa                                        >  1:87146/1‑71 (MQ=255)                           gATATTTTTCCGGGTGCTCGATGGCATCCAGCAGCATTTCCCGAtt                                                                 <  1:2344523/46‑1 (MQ=12)                            atatTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGtt                                                  <  1:2751272/60‑1 (MQ=255)                            atatTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGac                                       >  1:2458871/1‑71 (MQ=255)                               tttttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTtcat                                                              <  1:2823193/45‑1 (MQ=255)                                 tttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCa                                                >  1:1845417/1‑57 (MQ=255)                                 tttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCt                                          >  1:3168975/1‑63 (MQ=255)                                 tttCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGaccacc                                   >  1:747842/1‑70 (MQ=255)                                    ccGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTc                                <  1:1257222/70‑1 (MQ=255)                                    ccGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTc                                <  1:2698739/70‑1 (MQ=255)                                     cGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTg                                            <  1:2694289/57‑1 (MQ=255)                                     cGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGacc                                      <  1:1602284/63‑1 (MQ=255)                                       ggTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatg                            <  1:2205812/71‑1 (MQ=255)                                       ggTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatg                            <  1:757415/71‑1 (MQ=255)                                       ggTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatg                            <  1:662534/71‑1 (MQ=255)                                        gTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGtt                                                  <  1:3973057/48‑1 (MQ=255)                                        gTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGtt                                                  <  1:178005/48‑1 (MQ=255)                                        gTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatg                            >  1:1728911/1‑70 (MQ=255)                                         ttttCCATCGCGTCGAGCATCATTTCACGGTTCATCACGTTAACGtt                                                  <  1:3039242/47‑1 (MQ=255)                                         ttttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTc                                                 <  1:2731694/48‑1 (MQ=255)                                          tttCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc                                                           <  1:2931425/37‑1 (MQ=255)                                          tttCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc                                                           <  1:2926582/37‑1 (MQ=255)                                          tttCCATCGCGTCGAGCAGCATTTCACGGTTCATCAc                                                           <  1:1306678/37‑1 (MQ=255)                                           ttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatggatg                        >  1:319676/1‑71 (MQ=255)                                           ttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatggatg                        >  1:1585647/1‑71 (MQ=255)                                           ttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatggat                         >  1:3293199/1‑70 (MQ=255)                                           ttCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatg                            <  1:1791464/67‑1 (MQ=255)                                            tCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGc                       >  1:994847/1‑71 (MQ=255)                                             ccATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCt                      >  1:1545759/1‑71 (MQ=255)                                                tCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTg                                            <  1:637657/46‑1 (MQ=255)                                                  gcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCga                              >  1:1281253/1‑58 (MQ=255)                                                  gcgTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtg                 >  1:789870/1‑71 (MQ=255)                                                    gTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggt               >  1:3071801/1‑71 (MQ=255)                                                     tCGATCAGCATTTCACGGTTCATCACGTTAACGTTCa                                                <  1:1402645/37‑1 (MQ=255)                                                       gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatgga                          >  1:2311773/1‑57 (MQ=255)                                                       gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtg                 >  1:1126238/1‑66 (MQ=255)                                                       gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg             >  1:3085286/1‑70 (MQ=255)                                                       gAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGtggtgg             >  1:1333066/1‑70 (MQ=255)                                                           agcaTTTCCCGATTCATTACGTTGACGTTGAGGTGTTGACCGCCTTCGACATCCGCTTCGTGGTGGaaa          >  1:591115/1‑68 (MQ=12)                                                            gcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatgga                          >  1:65986/1‑52 (MQ=255)                                                            gcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatgga                          >  1:856523/1‑52 (MQ=255)                                                            gcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatgga                          >  1:58018/1‑52 (MQ=255)                                                            gcaTTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTaa       <  1:91753/71‑1 (MQ=255)                                                               tttCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTc                                <  1:2349467/43‑1 (MQ=255)                                                                 tCACGGTTCATCACGTTAACGTTCAGGTGCTGacca                                     >  1:511709/1‑36 (MQ=255)                                                                 tCACGGTTCATCACGTTAACGTTCAGGTGCTGacca                                     >  1:1106761/1‑36 (MQ=255)                                                                 tCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatggatg                        >  1:3355235/1‑49 (MQ=255)                                                                 tCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCgatggatg                        >  1:3546612/1‑49 (MQ=255)                                                                 tCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCgtggt               >  1:1086844/1‑58 (MQ=255)                                                                 tCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatc  >  1:1453375/1‑71 (MQ=255)                                                                  cACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAAccatc  >  1:924633/1‑70 (MQ=255)                                                                   |                                                                     GCCAGATACACGGATGGTCAGCTGCGGATATTTTTCCGGGTTTTCCATCGCGTCGAGCAGCATTTCACGGTTCATCACGTTAACGTTCAGGTGCTGACCACCTTCGATGGATGCTTCGTGGTGGAAGTAACCATC  >  minE/608049‑608183

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