BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,772	T→C	15.5%	A494A (GCA→GCG)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,772	0	T	C	15.5%	224.0 / 30.0	90	A494A (GCA→GCG)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base T (43/33); new base C (3/11); total (46/44)
Fisher's exact test for biased strand distribution p-value = 2.02e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 116 total aligned reads displayed.

ACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGT  >  minE/608703‑608839

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

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

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

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

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

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

    tCGTGCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTcagcgcgct                                                                <  1:4139864/71‑4 (MQ=255)

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

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

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

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

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

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

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

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

                     agagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgcag                                                                 <  1:184330/53‑1 (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: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)

                                     tACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAg                                 <  1:2537923/69‑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)

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

                                                         gatgatGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACtt           <  1:153771/71‑1 (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)

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

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

                                                                  cagtgcagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGt  >  1:4129483/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)

ACGGTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGT  >  minE/608703‑608839

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 36 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`