# HG changeset patch # User iuc # Date 1390480206 18000 # Node ID 16dc0e3d659bb321392a1f80272dc333708ade53 # Parent da0a1fcf7fe0deb3250e97e922b92525812cab7e Uploaded diff -r da0a1fcf7fe0 -r 16dc0e3d659b macs2_bdgbroadcall.xml --- a/macs2_bdgbroadcall.xml Tue Jan 21 13:14:33 2014 -0500 +++ b/macs2_bdgbroadcall.xml Thu Jan 23 07:30:06 2014 -0500 @@ -1,5 +1,5 @@ - - from bedGraph + + Call broad peaks from bedGraph output diff -r da0a1fcf7fe0 -r 16dc0e3d659b macs2_bdgcmp.xml --- a/macs2_bdgcmp.xml Tue Jan 21 13:14:33 2014 -0500 +++ b/macs2_bdgcmp.xml Thu Jan 23 07:30:06 2014 -0500 @@ -1,5 +1,5 @@ - - by comparing two signal tracks in bedGraph (bdgcmp) + + Deduct noise by comparing two signal tracks in bedGraph diff -r da0a1fcf7fe0 -r 16dc0e3d659b macs2_bdgdiff.xml --- a/macs2_bdgdiff.xml Tue Jan 21 13:14:33 2014 -0500 +++ b/macs2_bdgdiff.xml Thu Jan 23 07:30:06 2014 -0500 @@ -1,5 +1,5 @@ - - based on paired four bedgraph files + + Differential peak detection based on paired four bedgraph files diff -r da0a1fcf7fe0 -r 16dc0e3d659b macs2_bdgpeakcall.xml --- a/macs2_bdgpeakcall.xml Tue Jan 21 13:14:33 2014 -0500 +++ b/macs2_bdgpeakcall.xml Thu Jan 23 07:30:06 2014 -0500 @@ -1,5 +1,5 @@ - - from bedGraph output + + Call peaks from bedGraph output diff -r da0a1fcf7fe0 -r 16dc0e3d659b macs2_callpeak.xml --- a/macs2_callpeak.xml Tue Jan 21 13:14:33 2014 -0500 +++ b/macs2_callpeak.xml Thu Jan 23 07:30:06 2014 -0500 @@ -1,5 +1,5 @@ - - from alignment results + + Call peaks from alignment results R_3_0_1 @@ -62,8 +62,8 @@ ## TODO ## run R to create pdf from model script - ##if os.path.exists( os.path.join( tmp_dir, "MACS2_PREFIX_model.r" ) ): - ## cmdline = 'Rscript "MACS2_PREFIX_model.r" > "MACS2_PREFIX_model.r.log"' ) + ##if os.path.exists( os.path.join( tmp_dir, "MACS2_model.r" ) ): + ## cmdline = 'Rscript "MACS2_model.r" > "MACS2_model.r.log"' ) ## proc = subprocess.Popen( args=cmdline, shell=True, cwd=tmp_dir ) ## proc.wait() @@ -89,7 +89,7 @@ ##parse xls files to interval files as needed ##TODO is in working dir #if 'peaks_interval' in str($outputs).split(','): - #set $file = os.path.join( $tmp_dir, 'MACS2_PREFIX_peaks.xls' ) + #set $file = os.path.join( $tmp_dir, 'MACS2_peaks.xls' ) #if os.path.exists( $file ): echo '#peaks file' > $output_xls_to_interval_peaks_file; awk '$2-=1' $file >> $output_xls_to_interval_peaks_file; @@ -145,14 +145,13 @@ - + - @@ -181,19 +180,19 @@ - + 'peaks_bed' in outputs - + 'narrow' in outputs - + 'broad' in outputs - + 'gapped' in outputs - + 'summits' in outputs diff -r da0a1fcf7fe0 -r 16dc0e3d659b macs2_filterdup.xml --- a/macs2_filterdup.xml Tue Jan 21 13:14:33 2014 -0500 +++ b/macs2_filterdup.xml Thu Jan 23 07:30:06 2014 -0500 @@ -1,5 +1,5 @@ - - at the same position + + Remove duplicate reads at the same position diff -r da0a1fcf7fe0 -r 16dc0e3d659b macs2_predict.xml --- a/macs2_predict.xml Tue Jan 21 13:14:33 2014 -0500 +++ b/macs2_predict.xml Thu Jan 23 07:30:06 2014 -0500 @@ -1,5 +1,5 @@ - - d or fragment size from alignment results + + Predict 'd' or fragment size from alignment results R_3_0_1 @@ -8,7 +8,7 @@ macs2_macros.xml - macs2 predict + macs2 predictd -i #echo ','.join($infiles)# --tsize $tsize @effective_genome_size@ diff -r da0a1fcf7fe0 -r 16dc0e3d659b macs2_randsample.xml --- a/macs2_randsample.xml Tue Jan 21 13:14:33 2014 -0500 +++ b/macs2_randsample.xml Thu Jan 23 07:30:06 2014 -0500 @@ -1,5 +1,5 @@ - - of total reads + + Randomly sample number or percentage of total reads diff -r da0a1fcf7fe0 -r 16dc0e3d659b macs2_refinepeak.xml --- a/macs2_refinepeak.xml Tue Jan 21 13:14:33 2014 -0500 +++ b/macs2_refinepeak.xml Thu Jan 23 07:30:06 2014 -0500 @@ -1,5 +1,5 @@ - - + + Refine peak summits and give scores measuring balance of forward- backward tags (Experimental)