This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Contigs have different mutations than raw reads

Hi,

I ran two assembly softwares and called mutations from the contigs. However, when I compared the contigs with raw reads, I found some instances where raw reads does not have mutations but the contigs have. I don't know how to explain this. Please see the attached.

igv-snapshot

mutation-calling assembly

Which assemblers did you run and how?

The first one is megahit ($megahit -1 fastq1 -2 fastq2 -o /path/) and second one is metaspades ($metaspades.py -1 fastq1 -2 fastq2 -o /path/). Reads were trimmed with Trimmomatic. Its an RNA virus genome.

0 answers

No answers yet.

Log in to answer this question.