BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,789	T→C	13.3%	K489E (AAA→GAA) ‡	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,789	0	T	C	13.3%	256.3 / 28.3	98	K489E (AAA→GAA) ‡	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base T (48/37); new base C (3/10); total (51/47)
Fisher's exact test for biased strand distribution p-value = 3.58e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.00e-01

Only 100 of 104 total aligned reads displayed.

TCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATC  >  minE/608720‑608859

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

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

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

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

  agagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagc                                                                     >  1:451347/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)

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

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

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

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

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

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

               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)

                      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)

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

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

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

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

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

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

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

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

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

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

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

                                                  agtgcagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTcc                                      <  1:113244/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)

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

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

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

                                                         tgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACa              >  1:1848361/1‑71 (MQ=255)

                                                         tgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACa              >  1:165123/1‑71 (MQ=255)

                                                           aTGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCTCTCCa                                     >  1:3885380/1‑46 (MQ=255)

                                                           aTGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATc            >  1:4005471/1‑71 (MQ=255)

                                                             gTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGc          >  1:2692319/1‑71 (MQ=255)

                                                              tACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCa                                  <  1:290030/46‑1 (MQ=255)

                                                              tACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTtcat                            <  1:1026029/52‑1 (MQ=255)

                                                              tACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCg           <  1:3411269/69‑1 (MQ=255)

                                                              tACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCg           <  1:2423739/69‑1 (MQ=255)

                                                              tACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGcc         >  1:3009623/1‑71 (MQ=255)

                                                              tACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGcc         >  1:1399521/1‑71 (MQ=255)

                                                                cTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCtt       >  1:2246904/1‑71 (MQ=255)

                                                                  gTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCttt      >  1:2730455/1‑70 (MQ=255)

                                                                    ttAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCg           <  1:3425855/63‑1 (MQ=255)

                                                                     tAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGcc         >  1:4054520/1‑64 (MQ=255)

                                                                     tAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGAt   <  1:1946862/70‑1 (MQ=255)

                                                                     tAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATc  >  1:70148/1‑71 (MQ=255)

TCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATC  >  minE/608720‑608859

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