BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,789	T→C	20.9%	K489E (AAA→GAA) ‡	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,789	0	T	C	20.9%	315.2 / 74.9	158	K489E (AAA→GAA) ‡	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base T (61/64); new base C (14/19); total (75/83)
Fisher's exact test for biased strand distribution p-value = 5.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 164 total aligned reads displayed.

CAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATC  >  Exported/608721‑608859

 agagagGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGttt                                                                        >  1:2964557/1‑68 (MQ=255)

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

 agagagGCTTCGTAGCTGTACTTGTCGAGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagc                                                                     >  1:3157256/1‑71 (MQ=255)

          tCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTc                                                            >  1:4314953/1‑71 (MQ=255)

          tCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTc                                                            >  1:4167962/1‑71 (MQ=255)

             tAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTa                                                                       <  1:5263774/57‑1 (MQ=255)

              agatGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGGCAGCCAGTCCAt                                                         >  1:2373097/4‑70 (MQ=255)

               gCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagccagc                                                                 <  1:824590/61‑1 (MQ=255)

               gCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagccagc                                                                 <  1:5332163/61‑1 (MQ=255)

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

               gCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGa                                                       <  1:4621071/71‑1 (MQ=255)

               gCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATg                                                        <  1:1020106/70‑1 (MQ=255)

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

                  gTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGt                                                    >  1:1201001/1‑71 (MQ=255)

                   tACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagc                                                                     <  1:565187/53‑1 (MQ=255)

                   tACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagc                                                                     <  1:2415906/53‑1 (MQ=255)

                   tACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTg                                                   <  1:4662897/71‑1 (MQ=255)

                   tACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTg                                                   <  1:5146292/71‑1 (MQ=255)

                    aCTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGa                                                  >  1:4961221/1‑71 (MQ=255)

                    aCTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGa                                                  >  1:5543885/1‑71 (MQ=255)

                    aCTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGa                                                  >  1:252003/1‑71 (MQ=255)

                     cTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagc                                                                     <  1:4455116/51‑1 (MQ=255)

                     cTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagc                                                                     <  1:4376442/51‑1 (MQ=255)

                     cTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGAt                                                 <  1:628207/71‑1 (MQ=255)

                     cTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGAt                                                 <  1:80889/71‑1 (MQ=255)

                      ttGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagcca                                                                   <  1:284552/52‑1 (MQ=255)

                      ttGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagcca                                                                   <  1:3823966/52‑1 (MQ=255)

                      ttGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATc                                                >  1:5420974/1‑71 (MQ=255)

                       tGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagcca                                                                   >  1:317845/1‑51 (MQ=255)

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

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

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

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

                       tGTCGTGCATGTAGTGGATGATGTTCACTGCAGTGATGTACTGTTTagcca                                                                   >  1:395611/1‑51 (MQ=255)

                        gTCTTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATc                                                <  1:1123840/69‑1 (MQ=255)

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

                        gTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCa                                              >  1:3310956/1‑71 (MQ=255)

                         tCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATcga                                              <  1:2636195/70‑3 (MQ=255)

                          cGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATg                                            >  1:5374400/1‑71 (MQ=255)

                          cGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATg                                            >  1:5019602/1‑71 (MQ=255)

                            tGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGa                                                       <  1:3970828/58‑1 (MQ=255)

                             gCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagccagcc                                                                <  1:1966901/48‑1 (MQ=255)

                              cATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATgcgc                                         <  1:35490/70‑1 (MQ=255)

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

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

                                tGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagccagc                                                                 <  1:919093/44‑1 (MQ=255)

                                tGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagcaagc                                                                 <  1:558436/44‑1 (MQ=255)

                                tGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTc                                       >  1:4003165/1‑70 (MQ=255)

                                tGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTc                                       >  1:3015956/1‑70 (MQ=255)

                                  tAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGaa                                                      >  1:5577554/1‑53 (MQ=12)

                                  tAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGaa                                                      >  1:770928/1‑53 (MQ=12)

                                  tAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGaa                                                      >  1:1852218/1‑53 (MQ=12)

                                  tAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATc                                                <  1:1008395/59‑1 (MQ=255)

                                  tAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATc                                                <  1:1843464/59‑1 (MQ=255)

                                    gTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCa                                  <  1:3008225/71‑1 (MQ=255)

                                    gTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAg                                                     <  1:195627/52‑1 (MQ=12)

                                    gTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCa                                  >  1:2323210/1‑71 (MQ=255)

                                     tgaatgataTTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCAc                                 <  1:159215/68‑1 (MQ=255)

                                     tGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGa                                                       <  1:879306/49‑1 (MQ=255)

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

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

                                       gatgatGTTCAGTGCAGTGATGTACTGTTTagccag                                                                  <  1:921637/36‑1 (MQ=255)

                                       gatgatGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGAt                                                 >  1:4689700/1‑53 (MQ=255)

                                       gatgatGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACt                                <  1:2768294/70‑1 (MQ=255)

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

                                        atgataTTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACttt                              >  1:4764583/1‑70 (MQ=255)

                                        atgataTTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACttt                              <  1:4779975/71‑2 (MQ=255)

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

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

                                            tgTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATg                                                        <  1:5018171/41‑1 (MQ=255)

                                            tgTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTtcatca                          <  1:472241/71‑1 (MQ=255)

                                              ttCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTtcatca                          >  1:4852286/1‑69 (MQ=255)

                                              ttCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTtcatca                          >  1:5059907/1‑69 (MQ=255)

                                              ttCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCa                                  <  1:5063728/61‑1 (MQ=255)

                                               tCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTtcat                            >  1:1206162/1‑66 (MQ=255)

                                               tCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAg                       >  1:1383646/1‑71 (MQ=255)

                                               tCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAg                       >  1:4597266/1‑71 (MQ=255)

                                                cagtgcagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAg                       >  1:2147734/1‑70 (MQ=255)

                                                cagtgcagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAg                       >  1:528916/1‑70 (MQ=255)

                                                cagtgcagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAg                       >  1:5662568/1‑70 (MQ=255)

                                                  gtgcagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTc                    <  1:1142408/71‑1 (MQ=255)

                                                  gtgcagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTc                    <  1:5513822/71‑1 (MQ=255)

                                                   tgcagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATcc                                               <  1:660559/43‑1 (MQ=255)

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

                                                       gtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGa                                                  >  1:1914363/1‑36 (MQ=255)

                                                       gtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGAt                                                 >  1:1336284/1‑37 (MQ=255)

                                                       gtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATcc                                               <  1:2905721/39‑1 (MQ=255)

                                                        tgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCa                   <  1:2086920/66‑1 (MQ=255)

                                                        tgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGAc               <  1:1278528/70‑1 (MQ=255)

                                                        tgATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCa                                  <  1:5043685/51‑1 (MQ=255)

                                                          aTGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAg                  >  1:3079151/1‑65 (MQ=255)

                                                             tACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCa                                     <  1:4835358/43‑1 (MQ=255)

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

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

                                                                 gTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCa                                     >  1:524488/1‑39 (MQ=255)

                                                                 gTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCa                                     >  1:4808103/1‑39 (MQ=255)

                                                                 gTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCa                                     >  1:1775835/1‑39 (MQ=255)

                                                                    tAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCa                                  >  1:2599519/1‑39 (MQ=255)

                                                                    tAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATc  >  1:2914802/1‑71 (MQ=255)

                                                                    tAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATc  >  1:5209155/1‑71 (MQ=255)

CAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTTGATC  >  Exported/608721‑608859

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