arq5x/bedtools2

shuffle -incl: Fall out of -incl include.bed file

chenludi opened this issue · 1 comments

Dear developer,

I want to replace genomic blocks (A.bed) randomly on blocks of specified genomic regions (include.bed). However, the new genomic coordinates fall out of the "include.bed". human.genome

$cat A.bed
chr1 0 100
chr2 0 1000

$cat include.bed
chr1 200 203
chr1 204 400
chr1 300 325
chr2 2000 4000

$ head human.genome
chr1 240000000
chr2 240000000

$bedtools shuffle -i A.bed -g human.genome -chrom -incl include.bed
chr2 2811 2911
chr2 2666 3666

After shuffling, it exceeds the include.bed and the chromosome number was not kept. Did I misunderstand what the command does? Also, what if there is not enough room in include.bed to place the blocks in A.bed?

Thank you!
Best,
Chenlu

In addition, I used betools v2.30.0.