Huge difference in assembly stats between metaSPAdes and metaviralSPAdes when assembling gut viromes #1100
-
Hi, I performed assembly of my gut viromes uisng both metaSPAdes and metaviralSPAdes. It seems that the two pipelines give very different results. Below is the summary stats of the contigs.fasta files I got.
There are huge differences, except for the maximum length of contigs. Could there be huge differences like these between metaSPAdes and metaviralSPAdes? My virome sequences were obtained from fecal samples through procedures including viral particle enrichment. Taxonomic profiling using Phanta (https://www.biorxiv.org/content/10.1101/2022.08.05.502982v1, https://github.com/bhattlab/phanta) showed that ~90% of reads are from viruses. Thanks. |
Beta Was this translation helpful? Give feedback.
Replies: 2 comments 1 reply
-
Tagging @Dmitry-Antipov |
Beta Was this translation helpful? Give feedback.
-
Hi, |
Beta Was this translation helpful? Give feedback.
Hi,
metaviralSPAdes was designed to extract viral part of regular metagenomes, so for viral enriched datasets I'd suggest to use metaspades.
Difference between metaviral and meta is expected, since metaviral does not output anything but some paths, suspcicious to be complete viruses (using coverage and graph topology), and metaviral outputs all the assembled sequence, including also fragmented viral genomes, bacterial chromosomes and whatever you sequenced.