Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
BLAST Database error: No alias or index file found for nucleotide database [nt] chris_s Bioinformatics 14 03-24-2019 08:11 AM
Find The Database for HIV-pol,env and gag regions separately AnushaC Bioinformatics 1 10-20-2015 02:25 AM
BLAST Database error: No alias or index file found for protein database gt259 Bioinformatics 2 05-25-2015 12:31 AM
How to find all gi's for an organism in the nucleotide database? hlyates Bioinformatics 12 04-14-2015 01:56 AM
What is the difference between the GEO database and TSA database in genbank ? anuj2054 Bioinformatics 0 03-03-2015 01:59 PM

Thread Tools
Old 09-21-2021, 01:22 AM   #1
Junior Member
Location: Denmark

Join Date: Sep 2021
Posts: 1
Default Find contamination in database

Hi everyone.

I'm quite new to NGS and is currently working on creating a custom database of reference genomes.
I want to check if theres any human contamination in the database and therefore I'm trying to map the database to the human genome (GRCh38).

As I understand it, i have to create pseudoreads/contigs from each genome in my database before I can map them to the human genome.
Is that the right way to proceed? And if so, is there a smart way of doing this?

Any help is really appreciated
Nanopore_newbie is offline   Reply With Quote
Old 09-21-2021, 06:24 AM   #2
Senior Member
Location: East Coast USA

Join Date: Feb 2008
Posts: 7,138

You could use "minimap2" to align against the human genome. That said you will need to be judicious about results as you determine what is "human" contamination. Depending on the genomes you are working with there can be sequence similarities.

There are other aligners like "lastz" that can align chromosome sized sequences as well.
GenoMax is offline   Reply With Quote

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off

All times are GMT -8. The time now is 10:39 PM.

Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2022, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO