BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	minE	608,819	G→T	11.2%	R479S (CGC→AGC)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	minE	608,819	0	G	T	11.2%	333.4 / 15.4	106	R479S (CGC→AGC)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base G (36/58); new base T (5/7); total (41/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

Only 100 of 108 total aligned reads displayed.

CATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTT  >  minE/608751‑608888

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                ttCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATACATGCGCTCCATCACTTcatcat                                                      <  1:572039/70‑1 (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: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)

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

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

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

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

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

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

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

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

                               tACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGc                                       <  1:2925037/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:2316637/70‑1 (MQ=255)

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

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

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

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

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

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

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

                                  tGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCAt                                                                 <  1:3108930/41‑1 (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:2571067/71‑1 (MQ=255)

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

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

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

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

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

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

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

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

                                                   cATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCg                  <  1:576513/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)

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

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

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

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

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

                                                          tGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCttt                                   >  1:74818/1‑47 (MQ=255)

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

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

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

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

                                                                   cgcTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGAt                                >  1:1506823/1‑41 (MQ=255)

                                                                   cgcTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCAtttt  <  1:1198758/71‑1 (MQ=255)

CATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATCGGTTCAGACTTCGGACCAACCTGCATTTT  >  minE/608751‑608888

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