Question: de novo aseembly counting
0
gravatar for s.good
3.1 years ago by
s.good • 0
Canada
s.good • 0 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 COMMENT • link • modified 3.1 years ago by Jennifer Hillman Jackson ♦ 25k • written 3.1 years ago by s.good • 0
0
gravatar for Jennifer Hillman Jackson
3.1 years ago by
United States
Jennifer Hillman Jackson ♦ 25k 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 COMMENT • link written 3.1 years ago by Jennifer Hillman Jackson ♦ 25k

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

ADD REPLY • link written 2.2 years ago by Jennifer Hillman Jackson ♦ 25k
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