Mapped Read Counts For Rna-Seq Analysisusing Samtools In R
Hi,
Is it possible to use samtools, just as HTSeq count, count the mapped reads for RNA-Seq analysis? This is not clear in the documentation.
Thanks
• 3,490 views
•
link
written
by
robjohn7000
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Dealing with multiple-mapped reads in sRNA-seq
written by joaquinj12 •
-
counting reads over repeats
written by R •
-
rRNA Removal In Rna-Seq Data
written by jockbanan
-
Why The Properly Paired Results From Bwa Mem Is An Odd Number
written by liupfskygre
-
Extract Raw Counts From Sam File
written by mad.cichlids
-
Htseq Count Finding No Feature
written by joaslucas •
-
What Type Of Filtering Is Typically Applied To A Bam File?
written by Ashutosh Pandey
-
Rpkm Calculation
written by Varun Gupta
-
How To Process Soapalign Results Before Call Snps And Indels By Samtools?
written by Haiping •
-
How To Calculated Coverage Of The Genome/Contig Slice?
written by Leszek
No, it's not possible. You could probably hack together a solution (using a BED file as the annotation), but that'd be a waste of time given that there are already tools to do that.
Edit: It's unclear what any of this has to do with R. There's an Rsamtools package, which I have to admit I've never found very useful, but you could use to do counting in combination with other packages (I think there's an example in the
GenomicFeaturesorGenomicRangespackages).I don't think so. at best, you can compute the coverage with the mpileup command. HTSeq-count is usually used for this task, is there a reason why you don't want to use it?
Thanks for your comments everyone. I have had problems in the past with gff which is required by HTSeq-count, so I was just looking for alternative tools.