Announcement

Collapse
No announcement yet.

Making a bedGraph from a Bowtie file

Collapse
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • lbthrice
    replied
    Hi,

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

    Code:
    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

    Leave a comment:


  • quinlana
    replied
    Originally posted by gez View Post
    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:

    Code:
    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

    Leave a comment:


  • gez
    started a topic Making a bedGraph from a Bowtie file

    Making a bedGraph from a Bowtie file

    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!
Working...
X