BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,760	G→A	14.0%	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	14.0%	311.2 / 34.3	114	I498I (ATC→ATT)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base G (51/47); new base A (4/12); total (55/59)
Fisher's exact test for biased strand distribution p-value = 5.93e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.07e-01

Only 100 of 132 total aligned reads displayed.

GGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTT  >  minE/608693‑608830

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

  aTAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgatgat                                                                   <  1:3531184/71‑1 (MQ=255)

    aaCGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgatgat                                                                   >  1:348213/1‑69 (MQ=255)

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

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

     acatCACGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGaatg                                                                     <  1:2497046/63‑1 (MQ=255)

      cGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGtt                                                                >  1:413721/1‑70 (MQ=255)

      cGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTc                                                               <  1:837458/71‑1 (MQ=255)

      cGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTc                                                               <  1:1147769/71‑1 (MQ=255)

      cGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTc                                                               <  1:1321775/71‑1 (MQ=255)

       gTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTgg                                                                        >  1:1133633/1‑61 (MQ=255)

       gTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgatga                                                                    >  1:3446339/1‑65 (MQ=255)

        tCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcag                                                             >  1:2991471/1‑71 (MQ=255)

         cACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcag                                                             <  1:811859/70‑1 (MQ=255)

          aCGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgatg                                                                     <  1:308759/61‑1 (MQ=255)

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

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

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

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

          aCGGTAGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtg                                                           <  1:3826744/71‑1 (MQ=255)

          aCGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGaatgata                                                                  >  1:1691388/1‑63 (MQ=255)

          aCGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGaatgata                                                                  >  1:766778/1‑63 (MQ=255)

          aCGATCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcagcg                                                           >  1:1016592/1‑69 (MQ=255)

           cGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcag                                                             <  1:2037560/68‑1 (MQ=255)

            ggTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGgatgat                                                                   >  1:2564569/1‑61 (MQ=255)

              tCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgcagt                                                       <  1:3079021/71‑1 (MQ=255)

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

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

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

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

                   cAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATg                                                   <  1:218100/70‑1 (MQ=255)

                      cgcCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACt                                               <  1:2627648/71‑1 (MQ=255)

                       gcCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgca                                                         >  1:320981/1‑60 (MQ=255)

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

                          aTCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGtt                                            >  1:1522671/1‑70 (MQ=255)

                           tCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGt                                                  <  1:3051200/63‑1 (MQ=255)

                           tCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGtt                                                                <  1:4145358/49‑1 (MQ=255)

                           tCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTa                                          <  1:2708923/71‑1 (MQ=255)

                             agagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgcagtg                                                      <  1:3079110/57‑1 (MQ=255)

                             agagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgcagtg                                                      <  1:1099994/57‑1 (MQ=255)

                             agagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGa                                                     >  1:673111/1‑58 (MQ=255)

                             agagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGAt                                                    >  1:1827846/1‑59 (MQ=255)

                             agagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATg                                                   <  1:1085544/60‑1 (MQ=255)

                             agagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGttt                                           >  1:3886515/1‑68 (MQ=255)

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

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

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

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

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

                             agagagGCTTCGTAGCTGTACTTGTCGTGAATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagc                                        >  1:2619530/1‑71 (MQ=255)

                              gagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcag                                                             <  1:3372781/49‑1 (MQ=255)

                               agagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGtt                                                                <  1:1384714/45‑1 (MQ=255)

                               agagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGtt                                                                <  1:2126694/45‑1 (MQ=255)

                               agagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGtt                                                                <  1:861542/45‑1 (MQ=255)

                               agagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgcag                                                        <  1:184330/53‑1 (MQ=255)

                               agagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTa                                                 >  1:2418589/1‑60 (MQ=255)

                                gagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagcca                                      >  1:4136021/1‑70 (MQ=255)

                                 agGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTca                                                              >  1:4096061/1‑45 (MQ=255)

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

                                  agCTTCGTAGCTGTACTTGTCGTGCATGTAGTGaatga                                                                    >  1:2778588/2‑38 (MQ=39)

                                   gCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCgccagcca                                  <  1:4113818/71‑1 (MQ=255)

                                   gCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCgccagcca                                  <  1:1034688/71‑1 (MQ=255)

                                   gCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCgccagcca                                  <  1:3823915/71‑1 (MQ=255)

                                    cTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGtt                                                                <  1:1259773/40‑1 (MQ=255)

                                    cTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGccagccag                                 >  1:1368862/1‑71 (MQ=255)

                                        gTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCa                             >  1:2380215/1‑71 (MQ=255)

                                        gTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTc                               <  1:1198291/69‑1 (MQ=255)

                                         tAGCTGTACTTGTCGTGCATGTAGTGGATGATGtt                                                                <  1:2990876/35‑1 (MQ=255)

                                          aGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACt                                               >  1:4136330/1‑51 (MQ=255)

                                           gatGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGa                          <  1:2609227/69‑1 (MQ=255)

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

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

                                           gCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATg                           <  1:2546473/70‑1 (MQ=255)

                                           gCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATg                           <  1:3559966/70‑1 (MQ=255)

                                             tGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAg                        >  1:4140331/1‑71 (MQ=255)

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

                                              gTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGt                       >  1:1471044/1‑71 (MQ=255)

                                               tACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTg                      <  1:1098142/71‑1 (MQ=255)

                                                 cTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACt                                               <  1:1753592/44‑1 (MQ=255)

                                                 cTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGa                     <  1:2422392/70‑1 (MQ=255)

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

                                                 cTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTg                      >  1:134661/1‑69 (MQ=255)

                                                 cTTGTCCTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGAt                    >  1:5771/1‑71 (MQ=255)

                                                    gTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATc                   <  1:1579295/69‑1 (MQ=255)

                                                    gTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATc                   <  1:2385370/69‑1 (MQ=255)

                                                    gTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATc                   <  1:717764/69‑1 (MQ=255)

                                                       gTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATg               <  1:402324/70‑1 (MQ=255)

                                                           aTGTAGTGGATGATGTTCAGTGCAGTGATGTAGTGTTTAGCGAGCCAGTCCATGAAGTGATCCATGCGCTc          <  1:316212/71‑1 (MQ=255)

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

                                                            tGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTc          >  1:2153744/1‑70 (MQ=255)

                                                              tAGTGGATGATGTTCAGTGCATTGATTTACTGTTTAGCCAGCCAGTCCATGAAg                        >  1:920581/1‑54 (MQ=38)

                                                              tAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATgc              >  1:218226/1‑64 (MQ=255)

                                                              tAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCAt       >  1:1540462/1‑71 (MQ=255)

                                                               aGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCAATGAAGTGATc                   >  1:3660021/1‑58 (MQ=255)

                                                                gTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATg                           <  1:561172/49‑1 (MQ=255)

                                                                 tGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATgcgc            <  1:1005433/63‑1 (MQ=255)

                                                                  ggATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGa                          <  1:2923111/48‑1 (MQ=255)

                                                                   gatgatGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACtt  <  1:153771/71‑1 (MQ=255)

                                                                   gatgatGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACtt  >  1:510941/1‑71 (MQ=255)

                                                                   gatgatGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCAATCACtt  <  1:3595939/71‑1 (MQ=255)

GGATAACGTCACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTT  >  minE/608693‑608830

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