compare samples

[jgSoton]

Anyone know a quick way to compare two samples that have undergone exome-sequencing?

I seem to have two samples that have very similar variant call and I wonder if it is actually the same individual. -I have sam and bam files. I have an edited version of the variant call but not the raw vcf files and the moment.

thanks,
Jane

[sdvie]

Hi Jane,

In case you can get hold of the raw vcf files, you can go for this:
GATK Combine Variants
using set=Intersection to obtain the intersection of variant positions in two files

Otherwise, you could call the variants again using your bam files as input for the GATK Unified Genotyper.

cheers,
Sophia

[jgSoton]

Hi,

Thanks Sophia. I guessed I would need the vcf files really so I am remaking them. I will try the GATK option you suggest, and I think I will try vcftools to compare files too.

Cheers,
Jane