Question: de novo aseembly counting
0
gravatar for s.good
3.1 years ago by
s.good0
Canada
s.good0 wrote:

Hi,

 

I have assembled de novo reference gene and mapped it to tophat now i want to calculate the read count for mapped reads.

For this, i sorted my bam files and used htseq count for counting the reads . I used the gtf file for this.

Ques 1. Htseq is giving 0 for every read count? Why ? But i can see a lot of reads in IGB. Please suggest a way to count the number of mapped reads

 

Ques 2. When i performed cufflinks with the same gtf file it gave a lot ensemble ids and FPKM but why is HTSeq not working in this case, it is identifying the ensemble ids but giving 0 for every ensemble ids.

 

 

What kind  of programs can be used for this?

 

Thanks

 

 

counting de novo • 890 views
ADD COMMENTlink modified 3.1 years ago by Jennifer Hillman Jackson25k • written 3.1 years ago by s.good0
0
gravatar for Jennifer Hillman Jackson
3.1 years ago by
United States
Jennifer Hillman Jackson25k wrote:

Hello,

These problems are often associated with a reference genome mis-match problem between the inputs. In this case, the GTF annotation dataset and the fasta genome used for mapping (either built-in or custom). 

Start by double checking the chromosome identifiers in all inputs. Here is how: Support#Reference_genomes.

Best, Jen, Galaxy team

ADD COMMENTlink written 3.1 years ago by Jennifer Hillman Jackson25k

Also see: https://biostar.usegalaxy.org/p/19322/#19350

ADD REPLYlink written 2.2 years ago by Jennifer Hillman Jackson25k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 169 users visited in the last hour