How To Create A Read Density Profile Within A Interval?

HI!

I need some help: I have to create density profile with a window specific of 1kb (how many time a sequence is detected after NGS method). I have to use SAM and BEDtools, I think I can use genomeCov in BEDtools but I don't have genome reference.

So, if anybody is abble to help me...

Thanks

bedtools

• 3,436 views • link

updated 13.5 years ago by Istvan Albert

• written 13.5 years ago by Ophélie •

What programming langauge(s) are you familiar with?

• 0 views • link

written 13.5 years ago by Steve Lianoglou

I have to use 'BEDtools' to do that it's really specific.

• 0 views • link

written 13.5 years ago by Ophélie •

1 answer

There are many similar questions on BioStar, although I admit these are not that straightforward to find as the terminology is somewhat fuzzy. Perhaps this search helps:

http://www.biostars.org/search/?q=density+plots&t=top&submit=

In general the answer is to use a tool like HTSeq like so:

http://www-huber.embl.de/users/anders/HTSeq/doc/tss.html

• 0 views • link

written 13.5 years ago by Istvan Albert

It appears the top link is generating a 404, bummer.

• 0 views • link

written 11.7 years ago by Zev.Kronenberg

More posts like this

Sam To Bam
written 12.4 years ago by freddy •

I sequenced a mito genome using paired end method on the mi-seq I mapped all the reads back to reference gnome using <a href='http://bio-bwa.sourceforge.net/'>BWA</a> and …
How To Sort Sam File
written 12.5 years ago by Chen Sun

I want to process the <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file created by <a href='http://bio-bwa.sourceforge.net/'>BWA</a>, if it is sorted by the sequence coordinate on reference, it will be …
Sam To Bam - Loss Of Data Or Just Great Compression?
written 12.5 years ago by Clare

I use picard tools to take a <a href='http://bio-bwa.sourceforge.net/'>BWA</a> alignment <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file and convert it into a sorted <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> file. Normally this works …
Conversion From Sam To Bed With Max Information
written 12.6 years ago by roll

I am trying to convert my <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file into bed file in such a way that i will keep some information from the <a …
Calculate Short Read Mapping Sensitivity And Ppv
written 13.2 years ago by juniper-

Say I use ART to simulate some short reads with the purpose of evaluating the performance of a mapper like bowtie or <a href='http://bio-bwa.sourceforge.net/'>bwa</a>. There …
Merge The Sam Files For Dexseq
written 13.2 years ago by camelbbs

I have RNAseq data for one sample, with two technical replicates. As I know, I need to merge them to get the gene or exon …
How To Rearrange Paired End Bam File?
written 13.2 years ago by deepthithomaskannan

Hello all, I have a paired end <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> file and I want to use <a href='https://code.google.com/p/bedtools/'>bedtools</a> for them. After merging, the paired end read …
Filtering A Sam File For Quality Scores
written 14.3 years ago by Varun Gupta

Hi I have a <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file. The 5th column of the <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file tells us the quality score. In my <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> alignment …
Getting Unmapped Reads: Comparing Fastq To Bam
written 14.7 years ago by User 9996

given a FASTQ file and a <a href='http://samtools.sourceforge.net/SAM1.pdf'>BAM</a> file of aligned reads, is there an efficient way to get all FASTQ reads that are in …
Soap Result Converted To Bam Not Working In Samtool
written 14.7 years ago by Sarwar •

I converted the output of soap2 which is map file to <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> file using map2sam.pl and then converted this <a href='http://samtools.sourceforge.net/SAM1.pdf'>sam</a> to <a href='http://samtools.sourceforge.net/SAM1.pdf'>bam</a> …