Merge pull request #75 from loosolab/JannikHamp-patch-4

corrected the wrong stats output of number of footprints after filter

bin/1.2_filter_motifs/compareBed.sh

@@ -276,20 +276,18 @@ first_line=`sed -n 1p $data | sed "s/$/\tcontains_maxpos\tsequence/"`
 if [[ ${first_line:0:1} == "#" ]]
 then
 	echo "$first_line" > $output
-	# add some final values to the log file
+	# add initial number of footprints to the log file
 	fp_initial=`cat $data | wc -l`
 	fp_initial=`expr $fp_initial - 1`
-	fp_final=`cat "$workdir"/filtered.bed | wc -l`
-	fp_final=`expr $fp_final - 1`
 	echo $fp_initial | sed 's/^/initial number of footprints: /g' >> "$workdir"/compareBed.stats
-	echo $fp_final | sed 's/^/number of footprints after subtract: /g' >> "$workdir"/compareBed.stats
 else
 	# output will be overwritten if it exists
 	rm -f $output
-	# add some final values to the log file
+	# add initial number of footprints to the log file
 	cat $data | wc -l | sed 's/^/initial number of footprints: /g' >> "$workdir"/compareBed.stats
-	cat "$workdir"/filtered.bed | wc -l | sed 's/^/number of footprints after subtract: /g' >> "$workdir"/compareBed.stats
 fi
+# add number of footprints after filtering to the log file
+cat "$workdir"/filtered_flagged.bed | wc -l | sed 's/^/number of footprints after subtract: /g' >> "$workdir"/compareBed.stats
 
 # add fasta sequences to bed and create fasta file
 out_fasta=`echo $output | sed "s|.bed$|.fasta|g"`