BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,776	G→C	15.6%	T493S (ACT→AGT) ‡	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,776	0	G	C	15.6%	218.2 / 28.7	90	T493S (ACT→AGT) ‡	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base G (42/33); new base C (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 = 9.96e-01

Only 100 of 110 total aligned reads displayed.

TCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGAC  >  minE/608707‑608846

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

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

  gTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgcagtg                                                                      <  1:50034/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)

           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:1099994/57‑1 (MQ=255)

               agagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgcagtg                                                                      <  1:3079110/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:1732663/1‑71 (MQ=255)

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

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

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

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

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

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

                    ggCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagccagcc                                                   <  1:3326718/71‑1 (MQ=255)

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

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

                     gCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCgccagcca                                                  <  1:4113818/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)

                          gTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTc                                               <  1:1198291/69‑1 (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:1107939/1‑70 (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)

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

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

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

                                  aCTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTg                                      <  1:2965741/70‑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:2385370/69‑1 (MQ=255)

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

                                      gTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATc                                   <  1:1579295/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)

                                           gCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAgccagcca                                                  >  1:528735/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)

                                                tAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCAt                       >  1:1540462/1‑71 (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:3695324/54‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                     gtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGAt                                    >  1:450180/1‑37 (MQ=255)

                                                                     gtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGAc  <  1:50470/71‑1 (MQ=255)

                                                                     gtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGAc  <  1:161491/71‑1 (MQ=255)

TCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGAC  >  minE/608707‑608846

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