BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,766	G→A	28.2%	N496N (AAC→AAT)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,766	0	G	A	28.2%	170.3 / 77.3	110	N496N (AAC→AAT)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base G (45/34); new base A (16/15); total (61/49)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

Only 100 of 124 total aligned reads displayed.

TAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCA  >  minE/608696‑608835

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

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

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

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

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

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

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

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

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

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

        cGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgc                                                               <  1:1547791/71‑1 (MQ=255)

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

         ggTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGtt                                                                     <  1:1271516/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:2285052/70‑3 (MQ=255)

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

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

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

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

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

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

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

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

                 aGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTa                                                      >  1:1871441/1‑71 (MQ=255)

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

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

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

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

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

                        tCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGtt                                                 >  1:233172/1‑69 (MQ=255)

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

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

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

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

                           gagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagcc                                            <  1:1464048/71‑1 (MQ=255)

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

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

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

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

                             gagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagccag                                          >  1:3099829/1‑71 (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)

                                  ttCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGt                                     >  1:1300864/1‑71 (MQ=255)

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

                                    cGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagcca                                           >  1:2269684/1‑63 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                              cTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGAt                         >  1:1699360/1‑71 (MQ=255)

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

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

                                                 gTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagccag                                          >  1:2458232/1‑51 (MQ=255)

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

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

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

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

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

                                                  tCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCAt                     <  1:2686874/71‑1 (MQ=255)

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

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

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

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

                                                   cGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTATCCAGCCAGTCCATGAAGTGATCCATg                    >  1:3262567/1‑71 (MQ=255)

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

                                                       cATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTa                                               >  1:2013503/1‑40 (MQ=255)

                                                       cATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCt                <  1:2889267/71‑1 (MQ=255)

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

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

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

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

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

                                                              tGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTcc                                   <  1:1230465/45‑1 (MQ=255)

                                                              tGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCAc         >  1:2600382/1‑71 (MQ=255)

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

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

                                                                aatgataTTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCAc         >  1:1476818/2‑69 (MQ=255)

                                                                 atgataTTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATc           <  1:2189933/66‑1 (MQ=255)

                                                                  tgatgTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTTATCCATGCGCTCCATCACTtct     >  1:2892586/1‑70 (MQ=255)

                                                                  tgatgTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTtca     >  1:2406421/1‑71 (MQ=255)

                                                                  tgataTTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACtttg     <  1:1538895/71‑3 (MQ=255)

                                                                  tgataTTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACtttg     <  1:2592587/71‑3 (MQ=255)

                                                                  tgataTTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACtttg     <  1:3172816/71‑3 (MQ=255)

                                                                  tgataTTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACtttg     >  1:3276593/1‑69 (MQ=255)

                                                                     tgTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTtcatca  >  1:757589/1‑71 (MQ=255)

                                                                      gTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTtcatca  <  1:3148674/70‑1 (MQ=255)

TAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCA  >  minE/608696‑608835

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