counting reads for genes
0
0
Entering edit mode
10 months ago

Hello, I am counting number of reads assigned to each transcript using htseq-count after aligning antennae rna-seq to the annotated tsetse transcripts with STAR but all counts are zero and reads are specified as without features and alignment not unique .can someone help..

here is my command

htseq-count -m union -f bam -r pos -o sam --stranded=yes -t exon -i gene_id /home/billiah/Documents/Projects/responses_of_G.m.morsitans_to_odors/Star_Mapping_Results/Control_Rep_1_output_pairedAligned.sortedByCoord.out.bam  /home/billiah/Documents/data_sets/Transctipts/glossina-morsitans-yalebasefeaturesgmory1.9.gtf  > /home/billiah/Documents/Projects/responses_of_G.m.morsitans_to_odors/htseq_results/Star_htseq_results/Control_Rep_1_output_paired.counts.bam
alignment • 219 views
ADD COMMENT
0
Entering edit mode

a common cause for these kind of issues is that the gff/gtf file does not use the same sequence names as in your star index for mapping (or the fasta file or such). Can you check that first?

(so make sure that the sequence names are consistent over all files you are using)

ADD REPLY

Login before adding your answer.

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