BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	608,788	T→A	20.8%	K489I (AAA→ATA) ‡	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	608,788	0	T	A	20.8%	355.8 / 82.4	158	K489I (AAA→ATA) ‡	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base T (61/64); new base A (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 162 total aligned reads displayed.

CATCAGAGAGGCTTCGTAGCTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAGGACATCGCCTTT  >  Exported/608718‑608855

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                   cTGTACTTGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATg                                                    >  1:3284711/1‑69 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

                         ttGTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTagcca                                                               <  1:531762/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:5567874/1‑51 (MQ=255)

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

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

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

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

                           gTCGTGCATGTAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATcc                                           <  1:2236535/70‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                     tAGTGGATGATGTTCAGTGCAGTGATGTACTGTTTAGCCAGCCAGTCCa                                                      <  1:1974467/49‑1 (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)

                                     tAGTGAATGATATTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATc                                            <  1:1008395/59‑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)

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

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

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

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

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

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

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

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

                                            tgataTTCAGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACtttg                         >  1:3939805/1‑69 (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)

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

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

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

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

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

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

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

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

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

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

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

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

                                                         agtgATGTACTGTTTAGCCAGCCAGTCCATGAAGTGATCCATGCGCTCCATCACTTCATCATAGTTCAgg             <  1:3987384/70‑1 (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)

                                                           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)

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