BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,816	T→G	11.3%	M480L (ATG→CTG) ‡	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,816	0	T	G	11.3%	237.1 / 26.8	106	M480L (ATG→CTG) ‡	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base T (37/57); new base G (5/7); total (42/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 110 total aligned reads displayed.

TCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATT  >  minE/608746‑608886

tCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCAt                                                                        <  1:2686874/71‑1 (MQ=255)

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

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

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

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

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

      aTGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTc                                                                  >  1:787339/1‑71 (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)

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

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

             ggATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACt                                                           >  1:1342334/1‑71 (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:179449/71‑3 (MQ=255)

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

                tgataTTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACtttg                                                        <  1:3172816/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)

                     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)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                       tGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCttt                                 >  1:2347368/1‑71 (MQ=255)

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

                                            agccagccAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCg                            >  1:1895768/1‑71 (MQ=255)

                                              ccagccagTCCATGAAGTGATCGAGGCTGTCCATc                                                              >  1:2043614/1‑35 (MQ=19)

                                              ccagccagTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGt                          <  1:1578591/71‑1 (MQ=255)

                                               cagccagTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGtt                         <  1:3027864/71‑1 (MQ=255)

                                               cagccagTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGtt                         <  1:2571067/71‑1 (MQ=255)

                                                agccagTCCATGAAGTGATCCATGCGCTCCATCACTtc                                                         >  1:1722884/1‑38 (MQ=255)

                                                agccagTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGtt                         >  1:1966427/1‑70 (MQ=255)

                                                agccagTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGtt                         >  1:724778/1‑70 (MQ=255)

                                                agccagTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGt                          >  1:775429/1‑69 (MQ=255)

                                                  ccagTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTg                                <  1:2881458/61‑1 (MQ=255)

                                                  ccagTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAg                      <  1:3546067/71‑1 (MQ=255)

                                                   cagTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGa                     <  1:741024/71‑1 (MQ=255)

                                                   cagTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGa                     >  1:2264303/1‑71 (MQ=255)

                                                      tCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACtt                  <  1:1281181/71‑1 (MQ=255)

                                                      tCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACtt                  <  1:1563431/71‑1 (MQ=255)

                                                      tCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACtt                  <  1:1962357/71‑1 (MQ=255)

                                                        cATGAAGTGATCCTTGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCgg                           <  1:3915860/60‑1 (MQ=255)

                                                        cATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTc                 <  1:3688635/70‑1 (MQ=255)

                                                        cATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCg                <  1:576513/71‑1 (MQ=255)

                                                        cATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCg                <  1:1274164/71‑1 (MQ=255)

                                                         aTGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTc                                               <  1:624194/39‑1 (MQ=255)

                                                          tGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTc                                               <  1:918301/38‑1 (MQ=255)

                                                          tGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGa              <  1:20433/71‑1 (MQ=255)

                                                           gAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCgg                           <  1:1953278/57‑1 (MQ=255)

                                                            aaGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGt                          >  1:3755521/1‑57 (MQ=255)

                                                            aaGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGAcc            <  1:2599987/71‑1 (MQ=255)

                                                             aGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCt                                   <  1:57228/47‑1 (MQ=255)

                                                             aGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCt                                   <  1:3326833/47‑1 (MQ=255)

                                                             aGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGtt                         >  1:161966/1‑57 (MQ=255)

                                                             aGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCa           <  1:3434313/71‑1 (MQ=255)

                                                             aGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCa           <  1:377369/71‑1 (MQ=255)

                                                             aGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCa           <  1:1952164/71‑1 (MQ=255)

                                                              gTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACtt                  >  1:1624204/1‑63 (MQ=255)

                                                              gTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCATTTGATCGGTTCAGACtt                  >  1:1082188/1‑63 (MQ=255)

                                                                  tCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACtt                  <  1:54063/59‑1 (MQ=255)

                                                                  tCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTg      >  1:1403227/1‑71 (MQ=255)

                                                                  tCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTg      <  1:2934115/71‑1 (MQ=255)

                                                                  tCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTg      >  1:1093962/1‑71 (MQ=255)

                                                                    cATGCGCTCCATCACTTCATCATAGTTAAGGACATCGCCttt                                 <  1:774878/42‑1 (MQ=255)

                                                                     aTGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCAt   >  1:1594974/1‑71 (MQ=255)

                                                                      tGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCAtt  <  1:2717356/71‑1 (MQ=255)

                                                                      tGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCAtt  >  1:3445172/1‑71 (MQ=255)

TCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATT  >  minE/608746‑608886

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