This is a test version of Biostars. For the public version, visit https://www.biostars.org.
sorting a multi-sample (genotype) vcf file

I have a freebayes genotyped multisample vcf file. I want to sort the names of the samples in alphabetical order to make my life a bit easy with post variant calling analysis. I have tried Picared SortVcf which work fine which works fine on a small file but failed on a large file (~4 Gb). If you suggest ways to sort a large multi-sample file without distorting the variants information, it would a great help.

snp
• 8,855 views • link
updated by eli_ss501_midori • eli_ss501_midori • written by nagarsaggi • nagarsaggi

1 answer

Hey, try this:

$ bcftools query -l input.vcf | sort > samples.txt
$ bcftools view -S samples.txt input.vcf > output.vcf

If not already done, I would also suggest to use bcf instead of vcf or vcf.gz. This really improves speed when working with bcftools on large datasets.

fin swimmer

• 0 views • link
written by finswimmer finswimmer

I spent a little bit too much time trying to figure out how to do this just to come here and find this simple solution. Thanks!

• 0 views • link
written by curious curious

I have the same problem. Can you put an example of "samples.txt" please

• 0 views • link
written by eli_ss501_midori • eli_ss501_midori

it should just be a text file with no index or header and a single sample in each line

sample1
sample2
sample3

etc

• 0 views • link
written by curious curious

$ bcftools query -l populations.snps_whitelistRetainedLociSamplesSingletonsHWETags.vcf | sort > sample.txt [main] Unrecognized command.

I can´t not run

Somebody can help me?

• 0 views • link
written by eli_ss501_midori • eli_ss501_midori

$ bcftools query -l populations.snps_whitelistRetainedLociSamplesSingletonsHWETags.vcf | sort > sample.txt [main] Unrecognized command.

I can´t not run I send my database Somebody can help me?

• 0 views • link
written by eli_ss501_midori • eli_ss501_midori

Log in to answer this question.

More posts like this