This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to quantify the viral sequencing reads from RNA-seq data?

Hi everyone. I'm having problem with quantifying viral reads from my RNA-seq data, which is part of my Final Year Project. Basically, I have filtered the host reads by aligning to host organism genome. The unmapped reads was aligned to viral database with BLASTn, which illustrated the present of some specific virus. Then, I tried to align the unmapped reads to the genome of these virus that I have found using HISAT2, but the output is 0% matched reads. I need some suggestion to solve my problem. Thank you very much!

rna-seq viral alignments quantification reads
• 1,497 views • link
updated by Thind amarinder Thind amarinder • written by Huy Duong • Huy Duong

1 answer

There are few tool, which are desined for this purpose such as 'Decontaminer'. take a look here https://bmcbioinformatics.biomedcentral.com/articles/10.1186/s12859-019-2684-x

Here is the Github link for decontaminer tool. https://github.com/amarinderthind/decontaminer

• 0 views • link
written by Thind amarinder Thind amarinder

Thank you so much! This article is useful to me.

• 0 views • link
written by Huy Duong • Huy Duong

Log in to answer this question.

More posts like this