More posts like this
-
Nanopore specific assembly tool of Mitochondria
written by adarsh •
-
merging of two assemblies from different short-reads techniques
written by Andrew Yu •
-
Assembly at lineage level from shotgun data ?
written by v.berriosfarias •
-
Check if nanopores reads is from the same species as the one used in an assembly
written by Olivier •
-
How can I align a specific squence against fastq reads
written by Juliana •
-
extract pathogens from metaphlan3 results
written by serene.s •
-
Assembly of 454 raw data (in fastq)
written by bird77 •
-
Predict genes from megahit assembly sequence for virome tools
written by fanglujing •
-
Getting species counts from Illumina metagenome data
written by lotus28
-
significant testing for species abundancy
written by Prasad
Do you have an idea of what genomes to expect? Would there be one or more Depending on that you will need to select/create a database to search against.
If this is a single genome them simply blasting a few reads at NCBI should also give you an idea of what you have.
If your target species is included in a WGS database, you could use one of the metagenomic tools like Kraken2 and compare all your reads against a database of known sequences. Then you would also have a kind of confidence metric given the proportion of reads that map to the correct species. Another benefit would be identification and quantification of any potential contamination.
If your target species is not in a WGS database but sister taxa are, then you could still use the genus level mapping information.
I tend to use Kraken2 for QC of DNA and RNA reads and it is a really useful step.