Go Back   SEQanswers > Applications Forums > RNA Sequencing

Similar Threads
Thread Thread Starter Forum Replies Last Post
GBS-SNP-CROP: a reference-optional pipeline for SNP discovery arthurmelo Bioinformatics 1 03-15-2016 09:02 AM
novel species discovery metagenomics chris_bioinfo Bioinformatics 5 09-27-2012 09:19 AM
PubMed: A combined massively parallel sequencing - indicator species approach reveale Newsbot! Literature Watch 0 08-24-2010 02:40 AM

Thread Tools
Old 06-22-2017, 08:14 AM   #1
Location: Canada

Join Date: May 2017
Posts: 15
Question Confused about how to approach SNP discovery in hybridizing sister species


I'm new to genomics, but I have assembled a de novo transcriptome (Trinity), and have alignments for two sister species (two individuals per species = 4 individuals total) - the same species I used to develop the transcriptome.

My goal is to screen these species for SNPs using samtools, to see what genes of potential ecological importance are different between these closely related species. I am familiar with the samtools variant calling pipeline, but was looking for insight on how I should call the SNPs.

Do I pool my alignments, and obtain one set of SNPs for all species, or do I obtain SNPs separately for each species? I've read literature that discuss both, but was wondering why I should choose one or another. This is a basic question, but it has me stumped. I guess I'm confused whether the SNPs are called between my reference transcriptome and the alignments, or if they are called between the alignments, and the reference is being used as some sort of guide.

Thanks in advance.
MGCBrown 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 03:22 AM.

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