I have a bunch of viruses I've assembled from metagenomic data and I have the host genome as well. I want to see what proportion of each virus aligns to the host genome. I ran minimap2 with the host genome as the reference and can calculate the percent of the genome covered using samtools coverage but how can I do this with the reverse? Would it make sense to use the viral genome as the reference? If not, is there a way to do this with the SAM files I currently have already?
1 answer
I assume you used the chromosomes/scaffold/contigs of your host against the viral genomes, which I think oversimplifies the results as whole chromosomes are likely not integrated in a single event. I think you would get more appropriate results (depending on your question) using a k-mer based synteny tool. I am only familiar with Satsuma2, but it is a little old so unsure if there are better options available now.
You could use the satsuma_summary.chained.out output file would be easy to parse to identify which regions of the genome are similar and how similar. Some thresholds for inclusion would be sensible to remove putatively spurious hits with low sequence identity.
Log in to answer this question.
-
alignement of human cells infected with a virus
written by hajar •
-
RNA virus (HIV, HCV, Influenza, or SARS-CoV-2) datasets with corresponding fitness values.
written by jimmysibandah •
-
How can one identify genome segments and mutations in Influenza viruses?
written by eennadi •
-
How to quantify the viral sequencing reads from RNA-seq data?
written by Huy Duong •
-
How to create a merged human host viral RefSeq?
written by MatStat
-
How to assemble viral genomes when my data contains host DNA as well
written by GBC_Zonatos •
-
How to obtain a novel virus' complete genome from de novo assembly
written by brianna.flynn •
-
Virus RNA Seq reads mapping with salmon
written by lokraj2003
-
RNA-seq analysis of virally infected cells
written by Adrian Pelin
-
BAM file statics and heterozygoty rate
written by kk.mahsa
