How can I create gene specific bed file from whole genome sequencing data using command line or using software? In case of software what will be the pipeline. We are able to generate bam file and vcf. After that which pipeline should we follow. We have enough data but our aim/objectives are not fixed yet.