Genome Contamination Check with Blobtools
0
0
Entering edit mode
9 hours ago
saamhasan55 ▴ 30

I generated a genome assembly and wanted to check for contaminant sequences using blobtools. I carried out a genome alignment between the assembly and fastq reads using bwa-mem, used to diamond blastx to search my sequences against the blast database and then used blobtools. But blobtools finds no taxonomic hits, I only get undefined or no-hits. The diamond output has taxonomic IDs and when I search some sample predicted proteins from my assembly they map. Why is blobtools not finding any taxonomic hits?

diamond blobtools genome blast • 441 views
ADD COMMENT
0
Entering edit mode

I generated a genome assembly and wanted to check for contaminant sequences

Why do you think you have contaminant sequences in your assembly? Generally when libraries are made from an organismm there should be no extraneous sequence present.

ADD REPLY

Login before adding your answer.

Traffic: 3363 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