generate fasta file from chipseq peaks

[says_anova]

Hi all, hope someone has tried this before:

I have a collection of peaks from chipseq data. For each peak, I have the chr#, start position, and stop position. I'd like to generate a fasta file that contains the sequences from hg19 that correspond to each peak (one peak per line). Seems like there should be a tool out there to do this but my frantic googling is coming up short.

Advice??

[balaji]

Could you explain, What do you want to do with the fasta sequence of each peak and for how many peaks you want to extract the fasta sequences.

an example http://genome.ucsc.edu/cgi-bin/das/h...r5:91000,91050

or see

http://www.biostars.org/p/7481/

http://www.biostars.org/p/6156/

and probably there are many more

[dariober]

Hi all, hope someone has tried this before:

I have a collection of peaks from chipseq data. For each peak, I have the chr#, start position, and stop position. I'd like to generate a fasta file that contains the sequences from hg19 that correspond to each peak (one peak per line). Seems like there should be a tool out there to do this but my frantic googling is coming up short.

Advice??

Hi- Have a look at fastaFromBed in BEDtools. You should be able to do something like:

fastaFromBed -fi hg19.fasta -bed mypeaks.bed -fo mypeaks.fasta

[says_anova]

Thanks dariober! This looks like it'll do just what I need.

And thanks for the links, balaji, they're also helpful.

Update: yup, fastaFromBed did the trick!