Question: raw read mapping for differential gene expression
gravatar for Bioinfonext
3.6 years ago by
Bioinfonext220 wrote:

I have found that in some plant species, a large number of genes have similarity in the coding regions so when we do mapping on cDNA, it will not count the appropriate raw read counts.

Do you think best way is to map raw reads to extract whole mRNA sequence including 5' and 3' prime? but it should not have intron region.

Because most of the highly similar genes may have differences in 5' prime and 3 prime regions.

What is the best way to count real raw read in this type of plant species?

rna-seq • 853 views
ADD COMMENTlink written 3.6 years ago by Bioinfonext220

Multi-mapping reads are a common issue with short reads. How you treat them is somewhat up to you. I am sure there are threads on Biostars that have wiser answers. allows you to choose from one of these options. There are two extremes you can choose from.

ambiguous=best          (ambig) Set behavior on ambiguously-mapped reads (with 
                        multiple top-scoring mapping locations).
                            best    (use the first best site)
                            toss    (consider unmapped)
                            random  (select one top-scoring site randomly)
                            all     (retain all top-scoring sites)
ADD REPLYlink modified 3.6 years ago • written 3.6 years ago by genomax87k
Please log in to add an answer.


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