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,"=");