Thread: HTseq
View Single Post
Old 09-06-2012, 07:46 AM   #5
narges
Member
 
Location: Finland

Join Date: Aug 2012
Posts: 29
Default

Quote:
Originally Posted by dpryan View Post
@narges: the problem was that you sorted by coordinate rather than read name, which is required for pair-end reads.

Try instead:
Code:
samtools sort -n accepted_hits.bam accepted_hist.sorted
samtools view accepted_hits.sorted.bam | htseq-count - /hg19/Annotation/Genes/genes.gtf > count.txt
Thank you so much, now it works
narges is offline   Reply With Quote