BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,775	A→G	16.1%	T493T (ACT→ACC) ‡	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,775	0	A	G	16.1%	218.5 / 32.7	90	T493T (ACT→ACC) ‡	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base A (42/33); new base G (3/11); total (46/44)
Fisher's exact test for biased strand distribution p-value = 2.13e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.02e-01

Only 100 of 109 total aligned reads displayed.

TCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGA  >  minE/608707‑608845

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

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

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

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

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

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

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

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

         gcCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgca                                                                        >  1:1868695/1‑60 (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)

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

              cAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTag                                                        >  1:3120778/1‑71 (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)

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

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

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

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

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

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

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

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

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

                     gCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGa                                                                    >  1:438177/1‑52 (MQ=255)

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

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

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

                      cTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCCGTGATGTACTGTTTAGccagccag                                                >  1:2053782/1‑71 (MQ=255)

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

                       ttCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGccagccag                                                >  1:2181979/1‑70 (MQ=255)

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

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

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

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

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

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

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

                             gCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATg                                          >  1:1107939/1‑70 (MQ=255)

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

                             gCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATg                                          <  1:2546473/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)

                                 tACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAg                                       <  1:2537923/69‑1 (MQ=255)

                                 tACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTg                                     <  1:1339310/71‑1 (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)

                                    ttGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATg                                                                  <  1:3508920/39‑1 (MQ=255)

                                     tGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTa                                                                >  1:3057423/1‑40 (MQ=255)

                                      gTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGAt                                   >  1:3876769/1‑68 (MQ=255)

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

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

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

                                      gTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCa                                >  1:79584/1‑71 (MQ=255)

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

                                          tGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATgcg                            >  1:3656334/1‑71 (MQ=255)

                                          tGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATg                                          <  1:659908/57‑1 (MQ=255)

                                           gCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAgccagcca                                                 >  1:239678/1‑49 (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)

                                                 aGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATc                                  >  1:26146/1‑58 (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)

                                                     aatgataTTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACtt                 >  1:3186688/2‑71 (MQ=255)

                                                      atgatgTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACtt                 <  1:2269558/70‑1 (MQ=255)

                                                        gatgTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCAt                               <  1:1056947/54‑1 (MQ=255)

                                                        gatgTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCAt                               <  1:3695324/54‑1 (MQ=255)

                                                        gatgTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCAt                               <  1:1682426/54‑1 (MQ=255)

                                                        gatgTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACt                  <  1:1927204/67‑1 (MQ=255)

                                                        gatgTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTtcat              >  1:1997222/1‑71 (MQ=255)

                                                          tgTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTtcatca            <  1:2765494/71‑1 (MQ=255)

                                                          tgTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTtcatca            <  1:4129407/71‑1 (MQ=255)

                                                          tgTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTtcatca            <  1:3868023/71‑1 (MQ=255)

                                                              cagtgcagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTtcat              >  1:3256504/1‑65 (MQ=255)

                                                              cagtgcagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGt        >  1:4129483/1‑71 (MQ=255)

                                                              cagtgcagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGt        >  1:1321380/1‑71 (MQ=255)

                                                              cagtgcagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGt        >  1:948826/1‑71 (MQ=255)

                                                               agtgcagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTcc                        <  1:113244/54‑1 (MQ=255)

                                                               agtgcagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTcc                        <  1:3116493/54‑1 (MQ=255)

                                                                 tgcagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTc                         <  1:114212/51‑1 (MQ=255)

                                                                 tgcagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTtc                >  1:107161/1‑60 (MQ=255)

                                                                    agtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGa  <  1:363974/71‑1 (MQ=255)

TCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGA  >  minE/608707‑608845

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