Hello,
I have Illumina paired end data (phage/bacteria+phage) which I have assembled using CLC genomics workbench. I did a blastn of the resulting scaffolds and found hits for both the bacteria and the phage. Now I am doubtful if Blast is a true test that my data contains sequences from both the bacteria in question and it's phage.
What other method can I use to check the same and how can remove reads belonging to the phage from that belonging to bacteria?
2 answers
Take a look at this thread: BBSplit syntax for generating builds for the reference genome and how to call different builds. You should be able to use BBSplit to separate your reads into genome specific bins (with the exception of multi-mappers which you will need to choose how to handle). I recommend that you start with the raw data. Let BBSplit do the binning and then do your assembly in CLC.
You know for sure that the phage is not integrated in the bacterial genome (or are there 2 phage)?
Using the latest version of CLC Genomics Workbench there are plenty of options for this
- use the Microbial Genomics Module's binning tools to split your contigs from your assembly into bins either by taxonomic similarity or by sequence/kmer similarity (or both) - "pure bins" are probably good to go, but those that are could be investigated further for contamination.
- you could download UniVec from NCBI as a FASTA, import it into CLC and use read mapper with stringent settings to filter your FASTQ
for common contaminatiing vector sequences.
- use MGM to do a full taxonomic profiling of your raw fastq data using a custom database of UniVec, Phage, Plasmid and bacterial genomes. Pull out reads in the clades you want to keep, and the repeat your assembly with just those and with/without any "unclassified" reads.
Log in to answer this question.