I would like to use blobtools to filter out contaminating sequences from by genome assemblies. Inputs are assembly fasta, coverage file (bam), and a hits file. How do we make the hits file from a BLAST database?
See the documentation here: Hits file. For example, using blast:
-query $ASSEMBLY \
-db nt \
-outfmt ’6 qseqid staxids bitscore std’ \
-max_target_seqs 10 \
-max_hsps 1 \
Login before adding your answer.
Use of this site constitutes acceptance of our User Agreement and Privacy