From 7dc9a0569ad87a07fbfeea842cb3e16ec80148b2 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Jens=20Preu=C3=9Fner?= Date: Fri, 30 Oct 2015 09:17:55 +0100 Subject: [PATCH] Fixed a wrong reference to a bed file --- src/admire | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/src/admire b/src/admire index 3c76a57..83f2d82 100755 --- a/src/admire +++ b/src/admire @@ -436,7 +436,7 @@ for a; do bedtools sort -i $RESULTS/$N/$g.pvals.bed > $RESULTS/$N/$g.pvals.sorted.bed fi bedtools map -c 5 -o min,median,max,count -a $RESULTS/$N/$r/$g.pvals.passed_fdr_$Q.bed -b $RESULTS/$N/$g.pvals.sorted.bed > $RESULTS/$N/$r/$g.pvals.passed_fdr_$Q.abs.bed - bedtools map -c 5 -o min,median,max,count -a $RESULTS/$N/$r/$p.bed -b $RESULTS/$N/$g.pvals.sorted.bed > $RESULTS/$N/$r/$g.pvals.all.abs.bed + bedtools map -c 5 -o min,median,max,count -a $RESULTS/$N/$r/$g.pvals.bed -b $RESULTS/$N/$g.pvals.sorted.bed > $RESULTS/$N/$r/$g.pvals.all.abs.bed awk -v g="$g" 'BEGIN{FS="\t";OFS=",";split(g,c,"-");print "Feature","Gene","GeneID","Type","Chr","Start","Stop","Strand","P-val "c[1],"Q-val "c[1],"P-val "c[2],"Q-val "c[2],"Higher in","with Q-val","Coverage","Min.abs.diff","Median.abs.diff","Max.abs.diff";} {n = split($4,a,";");for(i=1; i<=n; i++){ split(a[i],b,"=");