This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Bedtools BEDMAP usage

Can bedtools bedmap be used in a for loop? Can we overwrite the output file using >> ?

assembly genome next-gen

Standard rules of shell usage apply, nothing special with these tools. If you give a representative example and show what you've tried we can try to debug.

 foreach ($v1 as $value) {
$outputh = system("/home/uday/anaconda3/bin/bedmap --chrom ". $value ." --mean ppCkrtYdyJ_3.bed pc.bed >> ppCkrtYdyJ_f10.txt");
}

v1 contains unique chromosomes. I'm getting an empty file as output.

0 answers

No answers yet.

Log in to answer this question.