Using HTSeq for paired-end files
1
0
Entering edit mode
7.7 years ago
Whoknows ▴ 890

Hi friends

I have 3 RNA-Seq paired end samples which I want to generate differentially expressed gene list, I used these commands:

samtools sort -n accepted_hits.bam accepted_hits_bam_sorted
samtools view -h accepted_hits_bam_sorted.bam >accepted_hits_sorted.sam

But I'm not sure about the following line that must create count table for DESeq input, is it fine (for paired-end)?

htseq-counts --stranded=no --mode=intersection-nonempty -t exon -i gene_id my.sort.sam annotation.gtf > output.txt

Thanks

HTSeq DeSeq RNA-Seq Tophat • 5.7k views
ADD COMMENT
2
Entering edit mode
7.7 years ago

That's fine.

For future reference, you can also just use htseq-counts --stranded=no --mode=intersection-nonempty -r pos accepted_hits.bam annotation.gtf > output.txt and save the hassle of sorting and/or conversion to SAM.

ADD COMMENT

Login before adding your answer.

Traffic: 1554 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6