Question: Deblock a VCF
0
gravatar for win
3.8 years ago by
win810
India
win810 wrote:

Hi all,

I am trying to deblock a gVCF file and it requires an input BED file. While I have the BED files for the standard genomes such as HG18/HG19 etc my question is how do I generate a BED file from the 1000 genomes reference fasta file or say from a custom reference fasta file?

Thanks in advance.

vcf • 985 views
ADD COMMENTlink modified 3.7 years ago by Nicola Casiraghi440 • written 3.8 years ago by win810
0
gravatar for Nicola Casiraghi
3.7 years ago by
Trento, IT
Nicola Casiraghi440 wrote:

Hi, getting the lenght of each chromosome of your reference genome (or fasta file), you can use bedtools makewindows to divide it in regions and represent them in bed format.

For example, to divide the human genome into windows of 1MB:

$ bedtools makewindows -g hg19.txt -w 1000000 > hg19_1Mb.bed

where hg19.txt, the file containing chromosome lengths, has the following tab delimited format:

chr1    249250621
chr2    243199373
...

 

ADD COMMENTlink written 3.7 years ago by Nicola Casiraghi440
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 829 users visited in the last hour