This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How can get very low counts in RNA-seq ?

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

rna-seq igv
• 1,990 views • link
updated by swbarnes2 swbarnes2 • written by sabaghianamir70 • sabaghianamir70

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.

• 0 views • link
written by Pierre Lindenbaum Pierre Lindenbaum

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

• 0 views • link
written by sabaghianamir70 • sabaghianamir70

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

• 0 views • link
written by ATpoint ATpoint

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

• 0 views • link
written by sabaghianamir70 • sabaghianamir70

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.

• 0 views • link
written by swbarnes2 swbarnes2

1 answer

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.

• 0 views • link
written by swbarnes2 swbarnes2

Log in to answer this question.

More posts like this