BRESEQ :: Evidence

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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,787	0	T	C	13.3%	259.6 / 28.2	98	K489K (AAA→AAG) ‡	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.94e-01

Only 100 of 101 total aligned reads displayed.

CATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTT  >  minE/608718‑608855

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

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

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

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

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

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

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

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

          gCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCgccagcca                                                           <  1:1034688/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:1282908/70‑1 (MQ=255)

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

                  gCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATg                                                    <  1:140941/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:2385370/69‑1 (MQ=255)

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

                           gTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATc                                            <  1:717764/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:528735/1‑49 (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)

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

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

                                                   cagtgcagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGt                  >  1:2050008/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:2423739/69‑1 (MQ=255)

                                                                tACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCg       <  1:3411269/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)

CATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTT  >  minE/608718‑608855

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