Question: Remove contamination before assembly
0
gravatar for bruseq
21 months ago by
bruseq40
bruseq40 wrote:

Hello everyone,

I am trying to draft a fungal genome, but the raw sequence R1 and R2 contain contamination, so before doing assembly I mapped Raw reads using bowtie2 against refseq Bacterial db and mitochondrial db, then exact unmapped reads using samtools {the unmapped reads are contamination free } then did assembly using MaSuRCA and get final.genome.sacffolds.fasta file. But again I did ncbi online blast against NR db with assembled fasta file, they show approx. 93% similarity with bacteria.

So, please guide why assembled fasta show similarity with bacteria after removing these(bacteria) in the previous mapping step using refseq db.

Is there any other way to remove contamination from Raw reads before assembly. Please guide me.

Thank you Divya

draft genome assembly genome • 1.3k views
ADD COMMENTlink modified 21 months ago by h.mon30k • written 21 months ago by bruseq40
1
gravatar for genomax
21 months ago by
genomax85k
United States
genomax85k wrote:

But again I did ncbi online blast against NR db with assembled fasta file, they show approx. 93% similarity with bacteria.

Did you by chance remove the reads you need and leave the bacterial reads in?

Is your data so badly contaminated that you feel the need to align against RefSeq bacterial DB? Can't you use a specific bacterium (or two)? Otherwise use a related fungal genome try to get the reads you need and leave the rest behind. I suggest using bbsplit.sh from BBMap suite for that purpose.

ADD COMMENTlink written 21 months ago by genomax85k
1
gravatar for h.mon
21 months ago by
h.mon30k
Brazil
h.mon30k wrote:

Did you try assembling and then removing the contaminant contigs? BlobTools is a good tool for this.

ADD COMMENTlink written 21 months ago by h.mon30k
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1531 users visited in the last hour