I have sequenced and assembled a eukaryotic genome using 10X genomics technology. I used their supernova software package. Now, I want to inspect if there is any bacterial contig among my assembled sequences.
I know that one way is to blast the contigs against NCBI bacterial refseq (ftp://ftp.ncbi.nih.gov/refseq/release/bacteria) and remove the contigs that have a certain percentage matches. I am wondering if there is any other way? Any software package/pipeline?
Thanks!
1 answer
sourmash is perfect for detecting contamination in an assembly:
https://angus.readthedocs.io/en/2019/sourmash.html#
If you do have contamination, I would use bbmap to remove those reads and reassemble without them.
bbmap.sh in1=R1.fq.gz in2=R2.fq.gz ref=contam.fa outu1=R1.clean.fq.gz outu2=R2.clean.fq.gz
Where R1 is your forward read fastq file, R2 is the reverse, contam.fa is the contaminate's genome, outu1 is the now "cleaned" forward fastq file without the contaminate reads, and R2 is the reverse.
https://jgi.doe.gov/data-and-tools/bbtools/bb-tools-user-guide/bbmap-guide/
Maybe look at GC content in the genome bins and separate by that? Check out the tool anvio, there could be some options:
Log in to answer this question.
-
Basic guidance on Using 10x Genomics Data for De Novo Assembly
written by aleponce •
-
ncbi-genome-download - bacterial and viral genomes ?
written by sunnykevin97
-
How to calculate distance matrix and phylogenetic tree for genome alignment
written by Kumar
-
Determine number and length of plasmids in bacterial genome (using something like Prokka)
written by suzuBell •
-
Using 10X Chromium linked reads for genome gap-filling
written by daren.card •
-
How to circularise the bacterial genome?
written by Kumar
-
Thousands of contigs in E.coli assembly
written by vitorgomesbio •
-
Sequence date of Refseq genomes
written by alyosamah •
-
Alignment of mutli-contig fasta files
written by L Ron Hubbard •
-
accession to taxid for fasta files downloaded from refseq
written by badribio
Can you comment on the size of the assembly and the length and number of sequences? You could do a naive blast search as stated.
You could try
sketchfrom BBMap suite.If you suspected that your data had contamination, it would have been much better to have identified that before assembly using
kraken2/centrifuge.Thanks. My sample is not contaminated. I just want to get the endosymbiont sequences. The genome is 600Mb, 260million 150bp reads. Thanks