BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,776	G→C	27.7%	T493S (ACT→AGT) ‡	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,776	0	G	C	27.7%	282.0 / 121.3	173	T493S (ACT→AGT) ‡	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base G (66/59); new base C (16/32); total (82/91)
Fisher's exact test for biased strand distribution p-value = 2.70e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 203 total aligned reads displayed.

GTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGAC  >  Exported/608706‑608846

gttgTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgcag                                                                        <  1:4200103/68‑1 (MQ=255)

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

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

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

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

     gCAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATg                                                                   >  1:172179/1‑71 (MQ=255)

     gCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATg                                                                   >  1:5244114/1‑71 (MQ=255)

      cAGCGCCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATg                                                                   <  1:4888514/70‑1 (MQ=255)

      cAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATg                                                                   <  1:1620676/70‑1 (MQ=255)

         cgcCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACt                                                               <  1:1975400/71‑1 (MQ=255)

         cgcCATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACt                                                               <  1:1537417/71‑1 (MQ=255)

          gcCATCAGAGAGTCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTg                                                              <  1:3744199/71‑1 (MQ=255)

          gcCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTg                                                              <  1:3488505/71‑1 (MQ=255)

          gcCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTg                                                              <  1:2673486/71‑1 (MQ=255)

            cATCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTgca                                                            <  1:1111821/71‑3 (MQ=255)

            cATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTg                                                              >  1:4474950/1‑69 (MQ=255)

            cATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTg                                                              >  1:2101143/1‑69 (MQ=255)

            cATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTg                                                              >  1:1419436/1‑69 (MQ=255)

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

             aTCAGCGAAGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTgcac                                                           <  1:4290162/71‑4 (MQ=255)

             aTCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTcagtgcagt                                                                       <  1:1812992/59‑1 (MQ=255)

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

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

                          cGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTc                                               <  1:5058023/70‑1 (MQ=255)

                          cGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTc                                               <  1:2279433/70‑1 (MQ=255)

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

                             aGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagccag                                                     >  1:4998409/1‑61 (MQ=255)

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

                              gCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGa                                                                     >  1:1453534/1‑44 (MQ=255)

                              gCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGa                                                                     >  1:2600751/1‑44 (MQ=255)

                              gCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGa                                                                     >  1:4110529/1‑44 (MQ=255)

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

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

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

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

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

                                  tACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAg                                        <  1:648144/69‑1 (MQ=255)

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

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

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

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

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

                                   aCTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTgca                                                            <  1:912092/48‑3 (MQ=255)

                                    cTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAg                                        >  1:5651246/1‑67 (MQ=255)

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

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

                                     ttGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATg                                                                   <  1:986249/39‑1 (MQ=16)

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

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

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

                                      tGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATc                                   >  1:1219189/1‑70 (MQ=255)

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

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

                                        tCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATc                                   <  1:3352880/68‑1 (MQ=255)

                                        tCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATc                                   <  1:2626188/68‑1 (MQ=255)

                                        tCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCAt                                >  1:4930906/1‑71 (MQ=255)

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

                                         cGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATg                               >  1:5374400/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)

                                               tGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATg                                           <  1:3589835/53‑1 (MQ=255)

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

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

                                                 tAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCa                        <  1:4580983/70‑1 (MQ=255)

                                                 tAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCAt                       <  1:4527083/71‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                          atgTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATcc                                  >  1:4782532/1‑51 (MQ=255)

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

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

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

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

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

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

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

                                                                    cagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAg     <  1:5444385/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)

                                                                      gtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACtt                  <  1:4432361/55‑1 (MQ=255)

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

GTCGTGCAGCGCCATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGAC  >  Exported/608706‑608846

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