SEQanswers

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

Reply
 
Thread Tools
Old 09-21-2021, 01:22 AM   #1
Nanopore_newbie
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
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 7,131
Default

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 online now   Reply With Quote
Reply

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 07:15 AM.


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