VCF header line counting
Hello happy bioinformaticians :)
It can be a very simple question but I want to ask that how can I count line (row) of header of VCF ?
I can be done manually but I want to get accurate result.
Simply count the number of rows starting with a
grep -c '^#' <vcf>
If you want to count the number of sequence headers:
grep -c '##contig' <vcf>
And even more.. if you want to count the number of non-headers:
grep -c '^[^#]' <vcf>
On the command line of course ;)
I think this would be the fastest
bcftools view --header-only <vcf> | wc -l
Traffic: 2109 users visited in the last hour