Question: raw read mapping for differential gene expression
0
gravatar for Bioinfonext
2.7 years ago by
Bioinfonext160
Korea
Bioinfonext160 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 • 694 views
ADD COMMENTlink written 2.7 years ago by Bioinfonext160

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.

BBMap.sh 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 2.7 years ago • written 2.7 years ago by genomax73k
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 2.3.0
Traffic: 752 users visited in the last hour