generate fasta file from chipseq peaks

says_anova · 01-30-2013, 09:01 PM

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 · 01-30-2013, 11:56 PM

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 · 01-31-2013, 12:52 AM

Quote:

Originally Posted by 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??

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

Code:

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

says_anova · 01-31-2013, 08:40 AM

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!

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Indexing ChIPseq libraries using Illumina's TruSeq and ChIPseq kits	Alex Clop	Epigenetics	6	11-08-2012 11:07 AM
ChIPSeq: Unusual narrow peaks	aritakum	Bioinformatics	5	04-11-2012 10:50 AM
samtools sort bam file error: generate non-existent file	mediator	Bioinformatics	0	03-05-2012 08:42 PM
ChIPseq: Strange Peaks after Gel Extraction	cmawhinney	General	1	12-13-2010 04:25 AM
ChIPseq: Strange Peaks after Gel Extraction	cmawhinney	Sample Prep / Library Generation	3	12-10-2010 08:23 AM

01-30-2013, 09:01 PM	#1
says_anova Junior Member Location: California Join Date: May 2011 Posts: 7	generate fasta file from chipseq peaks 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??

01-30-2013, 11:56 PM	#2
balaji Junior Member Location: Europe Join Date: Feb 2011 Posts: 8	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 Last edited by balaji; 01-31-2013 at 12:13 AM.

01-31-2013, 08:40 AM	#4
says_anova Junior Member Location: California Join Date: May 2011 Posts: 7	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! Last edited by says_anova; 01-31-2013 at 09:09 PM.