Hi,In order to do some parallelization, i wanted to get the list of chromosomes from my bcf file, but I couldn't find a command to do that.
My only guess is :bcftools view -h my_file.bcf | perl -lne 'if( /^##contig=<ID=(.+)>/ ){ print $1; }'
any better solution ? Sébastien L. _______________________________________________ Samtools-help mailing list Samtools-help@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/samtools-help
