SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
DE analysis with multiple groups NicoBxl Bioinformatics 11 08-06-2012 05:41 AM
picard add read groups HGENETIC Bioinformatics 0 01-25-2012 03:19 AM
Cuffdiff differential expression two groups jamminbeh Bioinformatics 0 12-09-2010 02:13 PM
PubMed: A simple method using PyrosequencingTM to identify de novo SNPs in pooled DNA Newsbot! Literature Watch 0 12-07-2010 02:00 AM
Comparison of experimental groups in Cufflinks ega2d Bioinformatics 1 11-08-2010 06:36 PM

Reply
 
Thread Tools
Old 01-03-2012, 05:54 PM   #1
slowsmile
Member
 
Location: long island

Join Date: May 2011
Posts: 22
Default how to identify different snps between two groups?

Dear seq members

I am fairly new to SNP analysis. In my study, I have 2 groups, treatment and control, each containing 3 biological replicates. My goal is to find different snps between the two groups and I'm using samtools mpileup | bcftools in my analysis.

Shown below are my codes:

samtools mpileup -C 50 -uDS -q 10 -f MYPATH/human_ref_genome.fa output/tr1.bam output/tr2.bam output/tr3.bam output/ctr1.bam output/ctr2.bam output/ctr3.bam | bcftools view -bvcg - > ./snp_output/snp_cmp.bcf

bcftools view -c -1 -3 ./snp_output/snp_cmp.bcf | vcfutils.pl varFilter -D 100 > ./snp_output/snp_cmp.vcf


Did I use the correct command? Why did the output snp_cmp.vcf file contain millions of snps while I'm expecting hundreds of snps of interest? How did I only limit (filter) my scope to the different snps between the two groups and which program should I use to achieve this goal?

Thanks in advance.
slowsmile is offline   Reply With Quote
Old 01-03-2012, 09:58 PM   #2
Xi Wang
Senior Member
 
Location: MDC, Berlin, Germany

Join Date: Oct 2009
Posts: 317
Default

I am also interested in this question. Hope some experts could help.
__________________
Xi Wang
Xi Wang is offline   Reply With Quote
Old 01-03-2012, 10:33 PM   #3
gringer
David Eccles (gringer)
 
Location: Wellington, New Zealand

Join Date: May 2011
Posts: 838
Default

Quote:
Why did the output snp_cmp.vcf file contain millions of snps while I'm expecting hundreds of snps of interest?
See here:

http://seqanswers.com/forums/showthread.php?t=16247

Hundreds of thousands of SNPs are expected and perfectly normal. However, I think millions would be pushing it a bit too much for comparing only a few individuals.

Quote:
In my study, I have 2 groups, treatment and control, each containing 3 biological replicates
Three biological replicates seems quite a low number for SNP frequency studies. I'm more used to hundreds or thousands of individuals compared (e.g. WTCCC). I'd expect you'll have a lot of trouble distinguishing between false and true positive results with only 6 individuals.

Edit: I've just noticed that the thread I linked to was also from you discussing the same issue. Oh well, I guess it must have needed some clarification.

Last edited by gringer; 01-03-2012 at 10:35 PM.
gringer is offline   Reply With Quote
Old 01-04-2012, 02:51 AM   #4
david.tamborero
Member
 
Location: spain

Join Date: Feb 2011
Posts: 60
Default

I'm actually newbie in bioinformatics and maybe I am missing something, but as far as I understand you should be interested in those SNPs which are different among the two groups. If so, you could use some of the tools designed for finding somatic mutations.

FYI, now I am getting a try to Varscan. It uses the group1_pileup and group2_pileup files as input and then outputs some nice statistics about the 'somatic status' of any variant.

hope it helps,

cheers
david
david.tamborero is offline   Reply With Quote
Reply

Tags
bcftools, samtools, snp

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 11:52 PM.


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