Biostar Beta. Not for public use.
Question: How can get very low counts in RNA-seq ?
0
Entering edit mode

Hi every one. i did rna seq analysis. in PER3 gene and PER4 pseudogene i didnt get any counts, but in IGV i saw some small peaks which means they have some counts, even very low, but still they have !!. i want to know is there any way i can get any counts from that ? were my options in my aligner wrong or what else ?

i used blood samples. and used with trimmomatic, STAR, featureCounts

ADD COMMENTlink 8 months ago sabaghianamir70 • 0 • updated 8 months ago swbarnes2 5.7k
Entering edit mode
0

i know these two are similar. but im make my point more clear this time. thank you anyway

ADD REPLYlink 8 months ago
sabaghianamir70
• 0
Entering edit mode
0

Did you try salmon? I think it is still the best option to make use of problematic reads such as multimappers.

ADD REPLYlink 8 months ago
ATpoint
17k
Entering edit mode
0

since i working in galaxy plat form . working with galaxy is a little confusing and different from other tools i used. because it using data without needing to aligner. what the exact steps of Salmon ? i mean after trimming. we use salmon and featureCounts? i dont know what workflow should i use

ADD REPLYlink 8 months ago
sabaghianamir70
• 0
Entering edit mode
0

Why don't you just try it first? It's right there on Galaxy. Try it and then you can come back with specific questions.

ADD REPLYlink 8 months ago
swbarnes2
5.7k
Entering edit mode
0

i want to know is there any way i can get any counts from that ?

may be there is a low MAPQ for those reads . Use samtools view in.bam "chr1:234-567" to explore the reads in the region of the gene.

ADD REPLYlink 8 months ago
Pierre Lindenbaum
120k
1
Entering edit mode

Your read counter might be refusing to count reads that do not align uniquely. Pierre's command will let you count reads in the region, even if your read counter won't assign them to a gene.

ADD COMMENTlink 8 months ago swbarnes2 5.7k

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.0