Go Back   SEQanswers > Bioinformatics > Bioinformatics

Similar Threads
Thread Thread Starter Forum Replies Last Post
16s v4 region library preparation with negative control trong.wised Illumina/Solexa 6 03-13-2017 03:40 AM
Best 16S region for 300PE Miseq ZanderH22 Metagenomics 2 03-13-2015 08:39 AM
Miseq Custom primer for 16s V1V3 region lmma Illumina/Solexa 12 02-26-2015 01:13 PM
lots of short read between 16s rDNA Region V1-V3 ? stare024 454 Pyrosequencing 4 11-25-2013 03:08 AM
best 16s region to amplify with NGS ? dph Metagenomics 6 07-02-2012 04:26 PM

Thread Tools
Old 11-13-2018, 04:30 PM   #1
Location: St. Louis

Join Date: Dec 2009
Posts: 74
Default mapping V1V3 region 16s to a database of V1V3+V3V5 region 16s

There is a table of 'most wanted genomes' (MWGs) from the Human Microbiome project that lists the 2698 organisms that the community generally thinks would be useful to have genome sequence. The table includes the consensus sequence of the V1V3 and/or V3V5 region of the 16s for each of these wanted organisms.

Locally I have a collection of 16s V1V3 region sequences that I would like to map to those MWG 16s regions but I am having trouble with the mapping. My initial attempt was to use ncbi blast to align the 16s regions of the MWGs to a database of my own 16s V1V3 regions. But I am seeing many of my local 16s sequences appearing as the strong, tophit for multiple of these MWGs. And I am getting 'strong' (by evalue) hits to almost every single MWG in the HMP collection, which is not what I had expected. Even when I filter the alignments to require (I think?) species level similarity of 97% across 90% of the query length I am still getting a fair number of my local 16s sequences mapping as the best hit to multiple MWGs. So basically I think my approach is too simplistic.

Could someone suggest a good way to do this? My goal is to get an idea of how many of these MWGs I have represented in my local collection of 16s sequences.

John Martin
jmartin 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 12:58 PM.

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