SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Samtools SNP Caller versus MAQ SNP caller BertieWooster Bioinformatics 5 05-18-2012 02:41 PM
GATK unified SNP caller error dakl Bioinformatics 4 11-23-2010 10:27 PM
Aligner/SNP caller for Illumina data sklages Bioinformatics 2 10-06-2010 05:48 AM
SNP Caller for BAM Files AnamikaDarwin Bioinformatics 0 01-22-2010 06:53 AM
SNP caller for RNA-Seq data shurjo Bioinformatics 0 12-08-2009 03:58 PM

Reply
 
Thread Tools
Old 12-19-2011, 06:20 AM   #1
brentp
Member
 
Location: salt lake city, UT

Join Date: Apr 2010
Posts: 72
Default best-practices for assigning read-groups for snp-caller

Hi, We have samples from 4 individuals from 2 tissues each == 8 total:

tissue is blood/lung
individuals are control/disease

so, let's (ignore that 2 samples is too few and) say
lc1,lc2,bc1,bc2,ld1,ld2,bd1,bd2
where l == lung, b==blood, c == control d == disease, and the last number is the individual of that disease category.

If we want to call lung-specific, disease-related SNPs, is it better to send to the caller with ld2,ld1 vs all others by giving ld1,ld1 a read-group of disease and the rest a read-group of control?
Or, is it better to just call snps where each of the 8 groups has it's own read-group and pull out the SNPs by intersecting post-hoc?

I'm not asking wether it's better to send in subsets, but how to divvy up the subsets.
thanks.
brentp is offline   Reply With Quote
Old 12-19-2011, 08:01 AM   #2
Dameon
Member
 
Location: St. Louis, MO - USA

Join Date: Dec 2011
Posts: 14
Default

I would use GATK available from the Broad's Genome Analysis Toolkit and indel realign, recalibrate, and call and filter genotypes following the Broad's recommended best practices. Then I would take the finished and annotated vcf file of just the SNPs (-glm of SNP in the Unified Genotyper) and use vcftools to convert the vcf file into Plink or Beagle format. Plink or Beagle would be better tools to conduct your association testing rather than try to perform a series of intersects. If you are familiar with SAS and JMPGenomics, I believe the latest version of JMPGenomics can readily import vcf format SNP files and SAS is an excellent platform to conduct any number of test you wish to perform. Good luck.
Dameon is offline   Reply With Quote
Old 12-20-2011, 08:34 AM   #3
brentp
Member
 
Location: salt lake city, UT

Join Date: Apr 2010
Posts: 72
Default

Thanks Dameon, do you know of an example of this in Plink or Beagle?
brentp is offline   Reply With Quote
Reply

Tags
mpileup, readgroups, snp calling

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 04:15 PM.


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