VCF header line counting
2
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.
Thanks,
BG
vcf
header
• 7.8k views
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 ;)
•
link
updated 2.5 years ago by
Ram
44k
•
written 9.7 years ago by
Coryza
▴
430
I think this would be the fastest
bcftools view --header-only <vcf> | wc -l
Login before adding your answer.
Traffic: 1720 users visited in the last hour
The grep command will read the entire file, so if your VCF file is very large, something like this will run faster:
Thanks it was so helpful !!!
Thanks !! :)
Please accept Coryza's answer if it was helpful, otherwise give them feedback on why it did not address your problem.