BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	609,051	C→T	17.4%	E401E (GAG→GAA)	pflB ←	pyruvate formate lyase I

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	609,051	0	C	T	17.4%	259.1 / 40.1	115	E401E (GAG→GAA)	pflB	pyruvate formate lyase I
Reads supporting (aligned to +/- strand): ref base C (33/62); new base T (11/9); total (44/71)
Fisher's exact test for biased strand distribution p-value = 1.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 129 total aligned reads displayed.

TCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTT  >  Exported/608981‑609120

tcatcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtc                                                                       <  1:3016969/71‑1 (MQ=255)

tcatcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtc                                                                       <  1:2113837/71‑1 (MQ=255)

 catcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                                                                      >  1:5542797/1‑71 (MQ=255)

 catcGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                                                                      >  1:4750025/1‑71 (MQ=255)

   ttgggCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCa                                                                    >  1:5684195/3‑71 (MQ=255)

   ttgggCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCa                                                                    >  1:3290511/3‑71 (MQ=255)

   ttgggCTGACGCAGCAGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCa                                                                    >  1:4556021/3‑71 (MQ=255)

     gggCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATa                                                                  >  1:230515/1‑71 (MQ=255)

         ttACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTg                                                               <  1:2953341/70‑1 (MQ=255)

            cGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCaga                                                           <  1:2738835/71‑1 (MQ=255)

             gcagcaGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCaagg                                                          <  1:4350307/71‑4 (MQ=255)

             gcagcaGGCAATCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCaagg                                                          <  1:3895082/71‑4 (MQ=255)

              cagcaAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAgaga                                                         <  1:2118396/71‑1 (MQ=255)

               agcaagcaATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGATCTCATACTGCAGagaa                                                        <  1:1570117/71‑1 (MQ=255)

                gcaggcaaTCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAggaag                                                       <  1:2366258/71‑1 (MQ=255)

                gcaagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                                                                    <  1:553177/58‑1 (MQ=255)

                gcaagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                                                                    <  1:1710426/58‑1 (MQ=255)

                gcaagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                                                                    <  1:564540/58‑1 (MQ=255)

                gcaagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                                                                    <  1:3215043/58‑1 (MQ=255)

                 caagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCaga                                                           <  1:154232/66‑1 (MQ=255)

                 caagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCaga                                                           <  1:312338/66‑1 (MQ=255)

                  aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaa                                                        >  1:3710429/1‑68 (MQ=255)

                  aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaa                                                        >  1:515667/1‑68 (MQ=255)

                  aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAagag                                                     >  1:3078264/1‑71 (MQ=255)

                  aagcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAagag                                                     <  1:316184/71‑1 (MQ=255)

                   ggcaaTCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCa                                                                    >  1:905889/2‑55 (MQ=255)

                   agcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAagag                                                     <  1:3330683/70‑1 (MQ=255)

                    gcaaTCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGgt                                                   <  1:4835971/71‑1 (MQ=255)

                    gcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGgt                                                   <  1:1435471/71‑1 (MQ=255)

                    gcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGgt                                                   <  1:419717/71‑1 (MQ=255)

                    gcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGgt                                                   >  1:4337394/1‑71 (MQ=255)

                    gcaaTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGgt                                                   <  1:1461097/71‑1 (MQ=255)

                      aaTCGCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGgtga                                                 <  1:4618897/71‑2 (MQ=255)

                       aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                                <  1:1082422/71‑1 (MQ=255)

                       aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                                >  1:2302400/1‑71 (MQ=255)

                       aTAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                                <  1:1551045/71‑1 (MQ=255)

                        tAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCa                                                                    <  1:960653/50‑1 (MQ=255)

                        tAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCaga                                                           <  1:1820379/59‑1 (MQ=255)

                        tAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCaga                                                           <  1:1803017/59‑1 (MQ=255)

                        tAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCg                                               <  1:5383566/71‑1 (MQ=255)

                          gCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGAt                                             <  1:2751295/71‑1 (MQ=255)

                          gCGTAATCGTCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGa                                              <  1:940828/70‑1 (MQ=255)

                              aGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATAc                                                                 <  1:5126632/47‑1 (MQ=255)

                               gTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCAt                                                                   <  1:2293534/44‑1 (MQ=255)

                                tcgtCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGacac                                       <  1:1115463/68‑1 (MQ=255)

                                tcgtCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGa                                          >  1:1970741/4‑68 (MQ=255)

                                tcgtCGCTGTTGAAGTCAGTACGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGa                                          >  1:2152665/4‑68 (MQ=255)

                                tCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTtct                                                                      <  1:3550239/40‑1 (MQ=255)

                                tCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCagag                                                          <  1:5667610/52‑1 (MQ=255)

                                  aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATa                                                                  >  1:770698/1‑42 (MQ=255)

                                  aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCa                                                             <  1:4137198/47‑1 (MQ=255)

                                  aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                                     >  1:3379225/1‑71 (MQ=255)

                                  aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                                     >  1:705335/1‑71 (MQ=255)

                                  aTCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                                     >  1:445539/1‑71 (MQ=255)

                                   tCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaag                                                       <  1:2456467/52‑1 (MQ=255)

                                   tCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaag                                                       <  1:5640415/52‑1 (MQ=255)

                                   tCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaag                                                       <  1:1334793/52‑1 (MQ=255)

                                   tCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTAATACTGCAGAGAAGAGGTGTCGATGGACACtt                                     >  1:3116851/1‑70 (MQ=255)

                                    cGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaag                                                       <  1:2280545/51‑1 (MQ=255)

                                     gttgttGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAg                                  >  1:5548385/1‑71 (MQ=255)

                                       tgttgAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAgcg                                >  1:456444/1‑71 (MQ=255)

                                         ttgAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaaga                                                      <  1:5177378/47‑1 (MQ=255)

                                          tgAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCagag                                                          <  1:4411910/42‑1 (MQ=255)

                                          tgAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCagag                                                          <  1:2901027/42‑1 (MQ=255)

                                            aaGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGagaa                                                        >  1:2942945/1‑42 (MQ=255)

                                                ccGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAAt                          <  1:1150458/68‑1 (MQ=255)

                                                 cGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTc                       <  1:5337889/70‑1 (MQ=255)

                                                 cGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTc                       <  1:2198931/70‑1 (MQ=255)

                                                   gACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTg                    <  1:3077275/71‑1 (MQ=255)

                                                   gACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGAAACTTTAGCGGCGAATTTCtt                     >  1:249047/1‑70 (MQ=255)

                                                    aCGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGa                   <  1:3909367/71‑1 (MQ=255)

                                                     cGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTg                    >  1:1036772/1‑69 (MQ=255)

                                                     cGCATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGACACCTGCGCGGCATATTTTTTg                    >  1:5175063/1‑69 (MQ=255)

                                                       cATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAgcg                                >  1:1500114/1‑55 (MQ=255)

                                                       cATCAGATCGTCATTTTCATACTGCAAGGAAGAGGTGACGATCGACACCTGCGCGGCATATTTTTTGAAgg                <  1:1224327/71‑2 (MQ=255)

                                                         tCAGTTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc              <  1:3847743/71‑1 (MQ=255)

                                                         tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                                <  1:16782/37‑1 (MQ=255)

                                                         tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                                <  1:611434/37‑1 (MQ=255)

                                                         tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTg                    <  1:2152486/65‑1 (MQ=255)

                                                         tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc              <  1:4503201/71‑1 (MQ=255)

                                                         tCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTc              <  1:2073074/71‑1 (MQ=255)

                                                          cAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGtt               >  1:4818782/1‑69 (MQ=255)

                                                           aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTc                                                <  1:3710638/35‑1 (MQ=255)

                                                           aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGa                                              >  1:5544640/1‑37 (MQ=255)

                                                           aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGaa                  >  1:2111126/1‑65 (MQ=255)

                                                           aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGaa                  >  1:1385118/1‑65 (MQ=255)

                                                           aGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGa                   >  1:2873580/1‑64 (MQ=255)

                                                             gTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTa                                   <  1:5636241/46‑1 (MQ=255)

                                                               cATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACtt                                     >  1:3452629/1‑42 (MQ=255)

                                                               cATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCgg         >  1:5694392/1‑70 (MQ=255)

                                                                 tCGTTCTCATTCTGCAGAGAAGAGGTGTCGATGGACACTTTAGGGGCGAATttgttg                    <  1:3695301/57‑1 (MQ=255)

                                                                 tCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAg      <  1:1416829/71‑1 (MQ=255)

                                                                 tCGTTCTCATACTGCAGAGAAGAGCTGTCGATGGa                                          <  1:4657500/35‑1 (MQ=255)

                                                                   gTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGGGGCAGtt    >  1:145938/1‑71 (MQ=255)

                                                                   gTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGtt    >  1:3218500/1‑71 (MQ=255)

                                                                    ttCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGttt   <  1:701365/71‑1 (MQ=255)

                                                                      ctcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGaa                           >  1:3134255/1‑45 (MQ=255)

                                                                      ctcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGaa                           >  1:2721395/1‑45 (MQ=255)

                                                                      ctcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGtttt  >  1:1640732/1‑70 (MQ=255)

                                                                      ctcATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGtttt  <  1:5499503/70‑1 (MQ=255)

TCATCGGGCTTACGCAGCAAGCAATAGCGTAGTCATCGTTGTTGAAGTCCGGACGCATCAGGTCATCGTTCTCATACTGCAGAGAAGAGGTGTCGATGGACACTTTAGCGGCGAATTTCTTGAAGTTCAGCGGCAGTTTT  >  Exported/608981‑609120

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