BRESEQ :: Evidence

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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,787	0	T	C	20.9%	316.4 / 75.5	158	K489K (AAA→AAG) ‡	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 = 9.77e-01

Only 100 of 165 total aligned reads displayed.

CATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTT  >  Exported/608718‑608855

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

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

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

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

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

             tCGTAGCTGTACTTGTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTc                                                        >  1:926233/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)

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

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

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

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

                  gCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATg                                                    >  1:4338760/1‑70 (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:1791742/53‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

                          tGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTTATGTACTGTTTagcca                                                               >  1:4660148/1‑51 (MQ=255)

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

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

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

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

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

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

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

                           gTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTc                                                        >  1:5177967/1‑57 (MQ=255)

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

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

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

                           gTCGTGCATGTAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTcc                                                       <  1:481885/58‑1 (MQ=255)

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

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

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

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

                              gTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGtt                                                                     >  1:3791247/1‑41 (MQ=255)

                              gTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGt                                                                      >  1:2618236/1‑40 (MQ=255)

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

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

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

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

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

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

                                   tGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTc                                   >  1:3015956/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:1852218/1‑53 (MQ=12)

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

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

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

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

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

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

                                          gatgatGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAg                                                 >  1:3481864/1‑49 (MQ=255)

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

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

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

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

                                            tgatgTTCAGTGCAGTGATGTACTGTTTACCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTtca                         <  1:2401335/71‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

                                                   cagtgcagtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGt                  >  1:87301/1‑71 (MQ=255)

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

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

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

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

                                                         agtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAgg             <  1:3987384/70‑1 (MQ=255)

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

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

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

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

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

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

                                                           tgATGTACTGCACCGCCAGCCAGTCCATGAAGTGATc                                            <  1:4753590/37‑1 (MQ=21)

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

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

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

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

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

                                                                    gTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCttt  >  1:610017/1‑70 (MQ=255)

CATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTT  >  Exported/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: `‑`