Contigs have different mutations than raw reads
0
0
Entering edit mode
3.5 years ago

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 • 745 views
ADD COMMENT
0
Entering edit mode

Which assemblers did you run and how?

ADD REPLY
0
Entering edit mode

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.

ADD REPLY

Login before adding your answer.

Traffic: 3146 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6