Making a bedGraph from a Bowtie file

[gez]

Hi everyone,

I have a few Bowtie alignment files that I would like to view unsmoothed (I usually use Fseq to find peaks and it smoothes things over) to get an idea of read depth. Is there any way to convert a Bowtie output file to a bedGraph?

Thanks!

[quinlana]

Hi everyone,

I have a few Bowtie alignment files that I would like to view unsmoothed (I usually use Fseq to find peaks and it smoothes things over) to get an idea of read depth. Is there any way to convert a Bowtie output file to a bedGraph?

Thanks!

If you use sam output from Bowtie, you could use the genomeCoverageBed tool from my BEDTools package as follows:

samtools view -Sb bowtieOut.sam | genomeCoverageBed -ibam stdin -g <genomefile> -bg > out.bedgraph

Note that the <genomefile> is a special file in BEDTools denoting the lengths of each chromosome. The "chromInfo" files from UCSC can be used.

Best,
Aaron

[lbthrice]

Hi,

genomeCoverageBed requires a sorted bowtieOut.sam
the following code worked for me:

samtools view -Sbo bowtieOut.bam bowtieOut.sam # convert to bam format
samtools sort bowtieOut.bam sortedbowtieOut.bam # sort the bam file
genomeCoverageBed -bg -ibam sortedbowtieOut.bam -g chrominfo.txt > out.bedgraph # compute coverage

see http://groups.google.com/group/bedto...usage-examples

-Lionel