Entering edit mode
4.1 years ago
Rashedul Islam
▴
480
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.
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.