This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Tool to check what alternate allele is dominant across samples per line of the VCF file.

Hello All,

I am very new to WGS analysis. I have a multisample VCF file that I have annotated using snpEFF. I wanted to see if I can find what alternate alleles are conserved between samples for each genomic location. For eg: Chr1 pos: 1001, has a reference A and the alternate allele seen are T, AAT, TTT, AA and there are 10 samples. I want to count what samples have TTT, AA and so on. This way I can understand what allele is dominant across samples for each position.

All help is appreciated

Thanks.

next-gen wgs
• 2,025 views • link
updated by chrchang523 chrchang523 • written by halo22 halo22

2 answers

I made very good experiences with the bioconductor vcfR package.

In case you're new to R, too and this is a one-off project, there's nothing wrong to just use Excel and multiple text-to-column operations to split the data (provided your machine is powerful enough to handle it). It's a bit tedious, but the learning curve is less steep

• 0 views • link
written by Carambakaracho Carambakaracho

:-D damn it, I got the excel shame AND didn't realise the thread was 20 days old.

@halo22 you better not use my excel advise and try to hire with Pierre I guess.

@pierre or anyone as this is most likely irrelevant for the OP anyway, does biostars feature strikeout markdown?

Cheers

• 0 views • link
written by Carambakaracho Carambakaracho

Thank you, guys! I wrote my own to get this done.

• 0 views • link
written by halo22 halo22

With plink 2.0,

plink2 --vcf <VCF path> --freq counts

gets this information for you. (Remove 'counts' if you want proportions instead.)

• 0 views • link
written by chrchang523 chrchang523

Log in to answer this question.

More posts like this