BedFrag
From genomewiki
Jump to navigationJump to search
#!/usr/bin/gawk -f # extract part of a bed file that lie completely between two limits: # example: "cat test.bed | bedFrag chr1 1000 1500" BEGIN { OFS="\t"; chrom = ARGV[1]; min = ARGV[2]; max = ARGV[3]; ARGV[1]=""; ARGV[2]=""; ARGV[3]=""; ARGC=1; } ($1 == chrom) && ($2 >= min) && ($3 <= max) { print $0 }