BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,788	T→A	21.8%	K489I (AAA→ATA) ‡	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,788	0	T	A	21.8%	223.5 / 54.4	105	K489I (AAA→ATA) ‡	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base T (42/40); new base A (11/12); total (53/52)
Fisher's exact test for biased strand distribution p-value = 8.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 111 total aligned reads displayed.

CATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTT  >  minE/608718‑608855

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

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

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

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

  tCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGtt                                                                     >  1:233172/1‑69 (MQ=255)

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

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

       gagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagccag                                                              >  1:3099829/1‑71 (MQ=255)

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

          gCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagccagc                                                             >  1:1012830/1‑69 (MQ=255)

          gCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagccagc                                                             >  1:3910658/1‑69 (MQ=255)

            ttCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGt                                                         >  1:1300864/1‑71 (MQ=255)

             tCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTa                                                                   >  1:2267359/1‑60 (MQ=255)

              cGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagcca                                                               >  1:2269684/1‑63 (MQ=255)

              cGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGt                                                         <  1:972323/69‑1 (MQ=255)

               gTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCATCGCGCTGATGTACTGCACCGCCAGCCAGTCCa                                                      >  1:3113024/1‑71 (MQ=255)

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

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

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

                   cTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATg                                                    >  1:1187658/1‑69 (MQ=255)

                     gTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGt                                                >  1:691787/1‑71 (MQ=255)

                      tACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGt                                                <  1:3668348/70‑1 (MQ=255)

                      tACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTg                                               <  1:2744294/71‑1 (MQ=255)

                      tACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGt                                                >  1:3348768/1‑70 (MQ=255)

                      tACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGt                                                >  1:723764/1‑70 (MQ=255)

                        cTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGAt                                             >  1:1699360/1‑71 (MQ=255)

                          tGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGAt                                             >  1:2902988/1‑69 (MQ=255)

                           gTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagccag                                                              >  1:2458232/1‑51 (MQ=255)

                           gTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATc                                            <  1:2918680/69‑1 (MQ=255)

                            tCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCa                                          <  1:2677928/70‑1 (MQ=255)

                            tCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAg                                                 >  1:1007084/1‑63 (MQ=255)

                            tCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAg                                                 >  1:3728046/1‑63 (MQ=255)

                            tCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATcgag                                         >  1:2502151/1‑68 (MQ=255)

                             cGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTATCCAGCCAGTCCATGAAGTGATCCATg                                        >  1:3262567/1‑71 (MQ=255)

                                 cATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagccagc                                                             <  1:2457636/46‑1 (MQ=255)

                                 cATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTa                                                                   >  1:2013503/1‑40 (MQ=255)

                                 cATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCt                                    <  1:2889267/71‑1 (MQ=255)

                                 cATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCt                                    >  1:1615469/1‑71 (MQ=255)

                                  aTGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTc                                   <  1:1118954/71‑1 (MQ=255)

                                  aTGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTc                                   >  1:787339/1‑71 (MQ=255)

                                  aTGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTc                                   >  1:79017/1‑71 (MQ=255)

                                  aTGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAg                                                 >  1:1027596/1‑57 (MQ=255)

                                    gTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCa                                 >  1:2264940/1‑71 (MQ=255)

                                    gTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGAGATCCATGCGCTcc                                  >  1:2305669/1‑70 (MQ=255)

                                     tAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCAt                                <  1:2956175/71‑1 (MQ=255)

                                     tAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCAt                                <  1:2640792/71‑1 (MQ=255)

                                       gTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCAt                                >  1:3067653/1‑69 (MQ=255)

                                        tGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTcc                                                       <  1:1230465/45‑1 (MQ=255)

                                        tGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCAc                             <  1:2996861/71‑1 (MQ=255)

                                        tGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCAc                             >  1:2600382/1‑71 (MQ=255)

                                         ggATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACt                            >  1:1342334/1‑71 (MQ=255)

                                          aatgataTTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCAc                             >  1:1476818/2‑69 (MQ=255)

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

                                           atgataTTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATc                               <  1:2189933/66‑1 (MQ=255)

                                            tgatgTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTTATCCATGCGCTCCATCACTtct                         >  1:2892586/1‑70 (MQ=255)

                                            tgatgTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTtca                         >  1:2406421/1‑71 (MQ=255)

                                            tgataTTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACtttg                         <  1:3172816/71‑3 (MQ=255)

                                            tgataTTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACtttg                         <  1:179449/71‑3 (MQ=255)

                                            tgataTTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACtttg                         >  1:3276593/1‑69 (MQ=255)

                                            tgataTTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACtttg                         <  1:2592587/71‑3 (MQ=255)

                                            tgataTTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACtttg                         <  1:1538895/71‑3 (MQ=255)

                                               tgTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTtcatca                      >  1:757589/1‑71 (MQ=255)

                                                gTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTtcatca                      <  1:3148674/70‑1 (MQ=255)

                                                 ttCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCAc                             <  1:508204/62‑1 (MQ=255)

                                                 ttCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTcatcat                     <  1:280813/70‑1 (MQ=255)

                                                 ttCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTcatcat                     <  1:674810/70‑1 (MQ=255)

                                                 ttCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATACATGCGCTCCATCACTTcatcat                     <  1:572039/70‑1 (MQ=255)

                                                  tCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATg                                        >  1:3772060/1‑50 (MQ=255)

                                                   cagtgcagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTtc                          >  1:2811083/1‑63 (MQ=255)

                                                   cagtgcagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAg                   <  1:521722/70‑1 (MQ=255)

                                                   cagtgcagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGAGCTCCATCACTTCATCATAGt                  >  1:1753131/1‑71 (MQ=255)

                                                    agtgcagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGt                  <  1:2595087/70‑1 (MQ=255)

                                                    agtgcagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGt                  <  1:3743181/70‑1 (MQ=255)

                                                    agtgcagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGt                  <  1:3873760/70‑1 (MQ=255)

                                                      tgcagtgATGTACTGTTTAGCCAGCCAGTCCATGaa                                                  <  1:752392/36‑1 (MQ=255)

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

                                                        cagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAgg             <  1:964238/71‑1 (MQ=255)

                                                        cagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAg              >  1:1938460/1‑70 (MQ=255)

                                                        cagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAg              >  1:2603018/1‑70 (MQ=255)

                                                        cagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAg              <  1:1024179/70‑1 (MQ=255)

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

                                                              tGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTtcatca                      <  1:2715161/56‑1 (MQ=255)

                                                              tGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATc        >  1:3524696/1‑70 (MQ=255)

                                                               gTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCa                                          >  1:3436955/1‑35 (MQ=255)

                                                               gTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACtt                           <  1:2552218/50‑1 (MQ=255)

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

                                                                tACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGc      <  1:773468/70‑1 (MQ=255)

                                                                 aCTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCt    >  1:536133/1‑71 (MQ=255)

                                                                 aCTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCt    >  1:2696081/1‑71 (MQ=255)

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

                                                                  cTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCt    <  1:2466096/70‑1 (MQ=255)

                                                                  cTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCt    <  1:169843/70‑1 (MQ=255)

                                                                  cTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCt    <  1:447392/70‑1 (MQ=255)

                                                                  cTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCt    <  1:2239444/70‑1 (MQ=255)

                                                                  cTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCt    <  1:2239621/70‑1 (MQ=255)

                                                                  cTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCt    <  1:3328611/70‑1 (MQ=255)

                                                                  cTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCt    <  1:2316637/70‑1 (MQ=255)

                                                                  cTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCt    <  1:1595246/70‑1 (MQ=255)

                                                                   tGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCAt                                <  1:3108930/41‑1 (MQ=255)

                                                                   tGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCttt  >  1:2347368/1‑71 (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: `‑`