BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,760	G→A	27.9%	I498I (ATC→ATT)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,760	0	G	A	27.9%	174.2 / 77.7	111	I498I (ATC→ATT)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base G (44/36); new base A (16/15); total (60/51)
Fisher's exact test for biased strand distribution p-value = 8.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

Only 100 of 125 total aligned reads displayed.

TGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACT  >  minE/608690‑608829

tGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTgg                                                                       <  1:3700027/71‑1 (MQ=255)

  cGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCTTGCATGTAGTGga                                                                      <  1:180471/70‑1 (MQ=255)

  cGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGga                                                                      <  1:3045120/70‑1 (MQ=255)

  cGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGga                                                                      <  1:115108/70‑1 (MQ=255)

   ggATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgatg                                                                    <  1:1051507/71‑1 (MQ=255)

   ggATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgat                                                                     >  1:1712621/1‑70 (MQ=255)

   ggATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgat                                                                     >  1:2498510/1‑70 (MQ=255)

    gATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgatga                                                                   <  1:1510895/71‑1 (MQ=255)

      tAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGatgatg                                                                 <  1:2111386/71‑1 (MQ=255)

       aaCGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGt                                                                >  1:2848019/1‑71 (MQ=255)

        acatCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATAtt                                                               <  1:2451649/68‑1 (MQ=255)

        aCGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGt                                                                <  1:284503/70‑1 (MQ=255)

          gTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTca                                                             <  1:3765323/71‑1 (MQ=255)

          atCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATAtt                                                               >  1:1375221/2‑69 (MQ=255)

            cACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcagc                                                           >  1:1683013/1‑70 (MQ=255)

             aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg                                                          <  1:583284/71‑1 (MQ=255)

             aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg                                                          <  1:2616968/71‑1 (MQ=255)

             aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg                                                          <  1:3659170/71‑1 (MQ=255)

             aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg                                                          >  1:2448796/1‑71 (MQ=255)

             aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg                                                          >  1:3690028/1‑71 (MQ=255)

             aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagt                                                           >  1:1746528/1‑70 (MQ=255)

               ggTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGtt                                                               <  1:1271516/64‑1 (MQ=255)

               ggTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGtt                                                               <  1:923735/64‑1 (MQ=255)

               ggTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg                                                          >  1:3745091/1‑69 (MQ=255)

               ggTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTc                                                              >  1:2678580/1‑65 (MQ=255)

                gTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTca                                                             <  1:827209/65‑1 (MQ=255)

                atCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcagcgcgc                                                       <  1:1986661/70‑3 (MQ=255)

                atCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcagcgcgc                                                       <  1:2285052/70‑3 (MQ=255)

                 tCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgcagt                                                      >  1:1361555/1‑71 (MQ=255)

                  cGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcagcgcgct                                                      >  1:3892463/1‑67 (MQ=255)

                  cGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcagcgcgct                                                      >  1:3818077/1‑67 (MQ=255)

                   gTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcag                                                            <  1:3051711/63‑1 (MQ=255)

                   gTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgcagtg                                                     <  1:3364346/70‑1 (MQ=255)

                   gTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGa                                                    >  1:3172857/1‑71 (MQ=255)

                   gTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGa                                                    >  1:524081/1‑71 (MQ=255)

                    tGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgcagtg                                                     <  1:1989274/69‑1 (MQ=255)

                     gCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGAt                                                   <  1:1474046/70‑1 (MQ=255)

                      cAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGaatgat                                                                  >  1:3949500/1‑54 (MQ=255)

                      cAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATg                                                  >  1:1270476/1‑70 (MQ=255)

                      cAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGt                                                 <  1:1276183/71‑1 (MQ=255)

                         cgcCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACt                                              >  1:1974860/1‑71 (MQ=255)

                           ccATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGt                                            <  1:2932575/71‑1 (MQ=255)

                            cATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGtt                                           <  1:3283633/71‑1 (MQ=255)

                            cATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGtt                                           <  1:1634670/71‑1 (MQ=255)

                             aTCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcagcg                                                          <  1:2296146/55‑3 (MQ=255)

                             aTCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTgcac                                          >  1:1489148/1‑68 (MQ=255)

                              tCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTAc                                               <  1:1963599/65‑1 (MQ=255)

                              tCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTgg                                                                       <  1:7523/41‑1 (MQ=255)

                               cAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGga                                                                      <  1:3129243/41‑1 (MQ=255)

                               cAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTag                                        >  1:258021/1‑71 (MQ=255)

                                agagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcag                                                            >  1:925830/1‑50 (MQ=255)

                                agagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagc                                       >  1:638018/1‑71 (MQ=255)

                                 gagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgatga                                                                   <  1:2817946/42‑1 (MQ=255)

                                  agagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGt                                                                >  1:213058/1‑44 (MQ=255)

                                  agagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGt                                                                >  1:3791007/1‑44 (MQ=255)

                                   gagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGatgatg                                                                 <  1:326538/42‑1 (MQ=255)

                                    agGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagccagc                                   >  1:3524790/1‑71 (MQ=255)

                                      gCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagccagc                                   >  1:3910658/1‑69 (MQ=255)

                                      gCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagccagc                                   >  1:1012830/1‑69 (MQ=255)

                                        ttCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgca                                                        <  1:1460678/46‑1 (MQ=255)

                                         tCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTa                                         >  1:2267359/1‑60 (MQ=255)

                                          cGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGt                               <  1:972323/69‑1 (MQ=255)

                                           gTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCATCGCGCTGATGTACTGCACCGCCAGCCAGTCCa                            >  1:3113024/1‑71 (MQ=255)

                                            tAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCAt                           <  1:2457256/71‑1 (MQ=255)

                                            tAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCAt                           <  1:3509691/71‑1 (MQ=255)

                                             aGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCAt                           >  1:3492343/1‑70 (MQ=255)

                                              gCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATg                          <  1:321491/70‑1 (MQ=255)

                                               cTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTg                                             <  1:563680/50‑1 (MQ=255)

                                               cTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTg                                             <  1:129318/50‑1 (MQ=255)

                                               cTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATg                          >  1:1187658/1‑69 (MQ=255)

                                                 gTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGt                      >  1:691787/1‑71 (MQ=255)

                                                  tACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGt                                            >  1:3491816/1‑48 (MQ=255)

                                                  tACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGt                      <  1:3668348/70‑1 (MQ=255)

                                                  tACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGt                      >  1:1843002/1‑70 (MQ=255)

                                                  tACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTg                     <  1:2744294/71‑1 (MQ=255)

                                                  tACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGt                      >  1:3348768/1‑70 (MQ=255)

                                                  tACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGt                      >  1:723764/1‑70 (MQ=255)

                                                    cTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACt                                              >  1:1219141/1‑44 (MQ=25)

                                                     ttGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATc                  >  1:2524114/1‑71 (MQ=255)

                                                      tGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGAt                   >  1:2902988/1‑69 (MQ=255)

                                                       gTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATc                  <  1:2918680/69‑1 (MQ=255)

                                                        tCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGt                                                 <  1:1544398/37‑1 (MQ=255)

                                                        tCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGt                                                 <  1:1612665/37‑1 (MQ=255)

                                                        tCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGtt                                           <  1:2311670/43‑1 (MQ=255)

                                                        tCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCa                <  1:2677928/70‑1 (MQ=255)

                                                        tCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAg                       >  1:1007084/1‑63 (MQ=255)

                                                        tCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATcgag               >  1:2502151/1‑68 (MQ=255)

                                                        tCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATcga                <  1:726165/70‑3 (MQ=255)

                                                             cATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagccagc                                   <  1:2457636/46‑1 (MQ=255)

                                                             cATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCt          >  1:1615469/1‑71 (MQ=255)

                                                              aTGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTc         >  1:787339/1‑71 (MQ=255)

                                                              aTGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTc         <  1:1118954/71‑1 (MQ=255)

                                                              aTGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAg                       >  1:1027596/1‑57 (MQ=255)

                                                                gTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCa       >  1:2264940/1‑71 (MQ=255)

                                                                gTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGAGATCCATGCGCTcc        >  1:2305669/1‑70 (MQ=255)

                                                                 tAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCAt      <  1:2640792/71‑1 (MQ=255)

                                                                 tAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCAt      <  1:2956175/71‑1 (MQ=255)

                                                                   gTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCAt      >  1:3067653/1‑69 (MQ=255)

                                                                    tGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCAc   <  1:2996861/71‑1 (MQ=255)

                                                                     ggATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACt  >  1:1342334/1‑71 (MQ=255)

TGCGGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACT  >  minE/608690‑608829

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