Question: Comparison of SNPs between groups of samples
1
gravatar for agata88
2.0 years ago by
agata88790
Poland
agata88790 wrote:

Hi all!

Could you please recommend a tool that can statistically compare SNPs between two groups of vcf's?

To be more clear, I have 30 vcf files for healthy patients, and 30 for sick patients. I need to search for some kind of SNP patterns, statistically significant between those two groups.

Any idea how to perform such analysis?

Thanks in advance,

Agata

snp • 1.1k views
ADD COMMENTlink modified 2.0 years ago by Pierre Lindenbaum127k • written 2.0 years ago by agata88790
1
gravatar for Pierre Lindenbaum
2.0 years ago by
France/Nantes/Institut du Thorax - INSERM UMR1087
Pierre Lindenbaum127k wrote:

run a fisher test for those variants. https://en.wikipedia.org/wiki/Fisher%27s_exact_test

ADD COMMENTlink written 2.0 years ago by Pierre Lindenbaum127k

Thanks, it's one way. I'am thinking about something different.

What is on my mind is to go back to variant detection step and use VarScan Somatic caller. I might use samtools mpileup to generate two vcf's each containing variants for 30 samples. Do you think it's good idea?

ADD REPLYlink written 2.0 years ago by agata88790

what would be the

"kind of SNP patterns,"

here ?

ADD REPLYlink written 2.0 years ago by Pierre Lindenbaum127k

I would like to see whether there are variants that occur only in sick patients and not in healthy group. In addition, I would like to see that in multiple samples comparison. For example, let's say I've found 5 mutations that occur in 20 sick samples and not occur in healthy group at all - this 5 mutation will create some kind of pattern.

ADD REPLYlink written 2.0 years ago by agata88790

and you'll need a statistical test ...

ADD REPLYlink written 2.0 years ago by Pierre Lindenbaum127k

to show significant difference at p< 0,05.

ADD REPLYlink written 2.0 years ago by agata88790
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 2133 users visited in the last hour