Help with a reading count command using HTSeq
1
0
Entering edit mode
7.3 years ago

I am working on bacterial RNAseq data analysis, here is what I have down:

Bowtie:

bowtie -S reference ~/RNAseq/QC/decompressed_mybacteria.fastq > mybacteria.sam

Htseq:

htseq-count -m intersection-nonempty --stranded=no mybacteria.sam template.gtf > count.sam

My original seq was created by using Hiseq Illumina, paired end.

I know it could be htseq variable I used here is not good.

Anyone can help with this command line? Thank you

RNA-Seq • 2.1k views
ADD COMMENT
0
Entering edit mode

You didn't ask a question or say what the problem is.

ADD REPLY
0
Entering edit mode
7.3 years ago
Manvendra Singh ★ 2.2k

Firstly, I am not sure that this is right command for paired end fastq files , they should be "prefix_1.fastq" and "prefix_2.fastq"

next, I would use featureCount instead of htseq-count, but if you want to use htseq-count, then this is right command.

ADD COMMENT

Login before adding your answer.

Traffic: 896 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