low alignment rate for RNA-Seq data
Dear scientists,
I have RNA-Seq data of rice samples. I analyzed two rice datasets from two different labs (Fastp, Hisat2, Samtools, FeatureCounts). Alignment rate of one rice dataset is 98% on average however alignment rate of other rice dataset ranges between 26% to 63%. Both datasets were sequenced together at the same time on same run.
Ideally alignment rate should be above 80% but i can not discard the dataset with low alignment rate. How can i utilize it. Do i need to extract mapped reads and perform quantification on mapped reads only. Please guide me.
Thankyou so much!
• 1,088 views
•
link
updated
by
GenoMax
•
written
by
Bushra Bibi
•
written
by
Bushra Bibi
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
structural variation calling through delly is stuck at Split-read assembly step
written by Bushra Bibi
-
is there any way to normalize the counts and compare it for the splicing quantification?
written by PK •
-
Good alignment rate for DNA-seq data
written by Bushra Bibi
-
Possible rRNA contamination in one sample
written by bioinfo355 •
-
featureCounts has low successfully assigned reads
written by tomas4482 •
-
Data transformation/normalization/scaling for machine learning using both RNA-Seq and Microarray da…
written by Felipe Colombelli •
-
extract unique sequences on bam file
written by SeaStar •
-
Salmon Quantification for RNA-seq Read Pairs with Different Lengths
written by ATpoint
-
Tophat vs Tophat-fusion concordant alignment rate
written by IP
-
Finding alignment score for multi-read frombam/sam file
written by EVR
Hi,
Check for rRNA reads, if they are still included in the read set (issues during rRNA depletion, or similar). These reads map usually to too many locations and are removed from the alignment. Other than that check for other sources of contamination. If you found it, you may want to remove the reads mapping to that source (e.g. using bbsplit).
Thankyou michael.ante! I will extract unmapped reads and blast to check contamination. I will let you know then. How to check for rRNA reads?
Align to rDNA repeat for Rice. Or use a tool like SortMeRNA: https://github.com/sortmerna/sortmerna