# HG changeset patch # User kosrou # Date 1425302205 18000 # Node ID b848d0eb3841bea705208c76ad949d7e30f4f441 Uploaded diff -r 000000000000 -r b848d0eb3841 ngs.plot.xml --- /dev/null Thu Jan 01 00:00:00 1970 +0000 +++ b/ngs.plot.xml Mon Mar 02 08:16:45 2015 -0500 @@ -0,0 +1,314 @@ + + +#if $gene_list_source_type.gene_list == "plot_whole_genome": +#if $genomic_region_source_type.genomic_region == "tss": +ngs.plot.tss_tes.sh $genome_name $genomic_region_source_type.genomic_region $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $genomic_region_source_type.flanking_region_size1 $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area +#else if $genomic_region_source_type.genomic_region == "tes": +ngs.plot.tss_tes.sh $genome_name $genomic_region_source_type.genomic_region $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $genomic_region_source_type.flanking_region_size2 $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area +#else if $genomic_region_source_type.genomic_region == "genebody": +#if $genomic_region_source_type.flanking_region_option_source_type1.flanking_region_option3 == "flanking_region_size": +ngs.plot.genebody_exon_cgi.sh $genome_name $genomic_region_source_type.genomic_region $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area $genomic_region_source_type.flanking_region_option_source_type1.flanking_region_option3 $genomic_region_source_type.flanking_region_option_source_type1.flanking_region_size3 $genomic_region_source_type.further_information1 $genomic_region_source_type.interval_size1 +#else if $genomic_region_source_type.flanking_region_option_source_type1.flanking_region_option3 == "flanking_floating_size": +ngs.plot.genebody_exon_cgi.sh $genome_name $genomic_region_source_type.genomic_region $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area $genomic_region_source_type.flanking_region_option_source_type1.flanking_region_option3 $genomic_region_source_type.flanking_region_option_source_type1.specify_flanking_region3 $genomic_region_source_type.further_information1 $genomic_region_source_type.interval_size1 +#end if +#else if $genomic_region_source_type.genomic_region == "exon": +#if $genomic_region_source_type.flanking_region_option_source_type2.flanking_region_option4 == "flanking_region_size": +ngs.plot.genebody_exon_cgi.sh $genome_name $genomic_region_source_type.genomic_region $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area $genomic_region_source_type.flanking_region_option_source_type2.flanking_region_option4 $genomic_region_source_type.flanking_region_option_source_type2.flanking_region_size4 $genomic_region_source_type.further_information2 $genomic_region_source_type.interval_size2 +#else if $genomic_region_source_type.flanking_region_option_source_type2.flanking_region_option4 == "flanking_floating_size": +ngs.plot.genebody_exon_cgi.sh $genome_name $genomic_region_source_type.genomic_region $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area $genomic_region_source_type.flanking_region_option_source_type2.flanking_region_option4 $genomic_region_source_type.flanking_region_option_source_type2.specify_flanking_region4 $genomic_region_source_type.further_information2 $genomic_region_source_type.interval_size2 +#end if +#else if $genomic_region_source_type.genomic_region == "cgi": +#if $genomic_region_source_type.flanking_region_option_source_type3.flanking_region_option5 == "flanking_region_size": +ngs.plot.genebody_exon_cgi.sh $genome_name $genomic_region_source_type.genomic_region $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area $genomic_region_source_type.flanking_region_option_source_type3.flanking_region_option5 $genomic_region_source_type.flanking_region_option_source_type3.flanking_region_size5 $genomic_region_source_type.further_information3 $genomic_region_source_type.interval_size3 +#else if $genomic_region_source_type.flanking_region_option_source_type3.flanking_region_option5 == "flanking_floating_size": +ngs.plot.genebody_exon_cgi.sh $genome_name $genomic_region_source_type.genomic_region $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area $genomic_region_source_type.flanking_region_option_source_type3.flanking_region_option5 $genomic_region_source_type.flanking_region_option_source_type3.specify_flanking_region5 $genomic_region_source_type.further_information3 $genomic_region_source_type.interval_size3 +#end if +#else if $genomic_region_source_type.genomic_region == "bed": +#if $genomic_region_source_type.flanking_region_option_source_type4.flanking_region_option6 == "flanking_region_size": +ngs.plot.bed.sh $genome_name $bed_file $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area $genomic_region_source_type.flanking_region_option_source_type4.flanking_region_option6 $genomic_region_source_type.flanking_region_option_source_type4.flanking_region_size6 $genomic_region_source_type.interval_size4 +#else if $genomic_region_source_type.flanking_region_option_source_type4.flanking_region_option6 == "flanking_floating_size": +ngs.plot.bed.sh $genome_name $bed_file $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area $genomic_region_source_type.flanking_region_option_source_type4.flanking_region_option6 $genomic_region_source_type.flanking_region_option_source_type4.specify_flanking_region6 $genomic_region_source_type.interval_size4 +#end if +#end if +#else if $gene_list_source_type.gene_list == "offer_gene_list": +#if $genomic_region_source_type.genomic_region == "tss": +ngs.plot.tss_tes_genelist.sh $genome_name $genomic_region_source_type.genomic_region $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $genomic_region_source_type.flanking_region_size1 $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area $gene_list_source_type.gene_list_detail +#else if $genomic_region_source_type.genomic_region == "tes": +ngs.plot.tss_tes_genelist.sh $genome_name $genomic_region_source_type.genomic_region $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $genomic_region_source_type.flanking_region_size2 $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area $gene_list_source_type.gene_list_detail +#else if $genomic_region_source_type.genomic_region == "genebody": +#if $genomic_region_source_type.flanking_region_option_source_type1.flanking_region_option3 == "flanking_region_size": +ngs.plot.genebody_exon_cgi_genelist.sh $genome_name $genomic_region_source_type.genomic_region $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area $genomic_region_source_type.flanking_region_option_source_type1.flanking_region_option3 $genomic_region_source_type.flanking_region_option_source_type1.flanking_region_size3 $genomic_region_source_type.further_information1 $genomic_region_source_type.interval_size1 $gene_list_source_type.gene_list_detail +#else if $genomic_region_source_type.flanking_region_option_source_type1.flanking_region_option3 == "flanking_floating_size": +ngs.plot.genebody_exon_cgi_genelist.sh $genome_name $genomic_region_source_type.genomic_region $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area $genomic_region_source_type.flanking_region_option_source_type1.flanking_region_option3 $genomic_region_source_type.flanking_region_option_source_type1.specify_flanking_region3 $genomic_region_source_type.further_information1 $genomic_region_source_type.interval_size1 $gene_list_source_type.gene_list_detail +#end if +#else if $genomic_region_source_type.genomic_region == "exon": +#if $genomic_region_source_type.flanking_region_option_source_type2.flanking_region_option4 == "flanking_region_size": +ngs.plot.genebody_exon_cgi_genelist.sh $genome_name $genomic_region_source_type.genomic_region $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area $genomic_region_source_type.flanking_region_option_source_type2.flanking_region_option4 $genomic_region_source_type.flanking_region_option_source_type2.flanking_region_size4 $genomic_region_source_type.further_information2 $genomic_region_source_type.interval_size2 $gene_list_source_type.gene_list_detail +#else if $genomic_region_source_type.flanking_region_option_source_type2.flanking_region_option4 == "flanking_floating_size": +ngs.plot.genebody_exon_cgi_genelist.sh $genome_name $genomic_region_source_type.genomic_region $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area $genomic_region_source_type.flanking_region_option_source_type2.flanking_region_option4 $genomic_region_source_type.flanking_region_option_source_type2.specify_flanking_region4 $genomic_region_source_type.further_information2 $genomic_region_source_type.interval_size2 $gene_list_source_type.gene_list_detail +#end if +#else if $genomic_region_source_type.genomic_region == "cgi": +#if $genomic_region_source_type.flanking_region_option_source_type3.flanking_region_option5 == "flanking_region_size": +ngs.plot.genebody_exon_cgi_genelist.sh $genome_name $genomic_region_source_type.genomic_region $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area $genomic_region_source_type.flanking_region_option_source_type3.flanking_region_option5 $genomic_region_source_type.flanking_region_option_source_type3.flanking_region_size5 $genomic_region_source_type.further_information3 $genomic_region_source_type.interval_size3 $gene_list_source_type.gene_list_detail +#else if $genomic_region_source_type.flanking_region_option_source_type3.flanking_region_option5 == "flanking_floating_size": +ngs.plot.genebody_exon_cgi_genelist.sh $genome_name $genomic_region_source_type.genomic_region $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area $genomic_region_source_type.flanking_region_option_source_type3.flanking_region_option5 $genomic_region_source_type.flanking_region_option_source_type3.specify_flanking_region5 $genomic_region_source_type.further_information3 $genomic_region_source_type.interval_size3 $gene_list_source_type.gene_list_detail +#end if +#else if $genomic_region_source_type.genomic_region == "bed": +#if $genomic_region_source_type.flanking_region_option_source_type4.flanking_region_option6 == "flanking_region_size": +ngs.plot.bed_genelist.sh $genome_name $bed_file $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area $genomic_region_source_type.flanking_region_option_source_type4.flanking_region_option6 $genomic_region_source_type.flanking_region_option_source_type4.flanking_region_size6 $genomic_region_source_type.interval_size4 $gene_list_source_type.gene_list_detail +#else if $genomic_region_source_type.flanking_region_option_source_type4.flanking_region_option6 == "flanking_floating_size": +ngs.plot.bed_genelist.sh $genome_name $bed_file $bam_file $out_zip $out_avg_png $out_hm_png "$image_title" $gene_database $randomly_sample $gene_order $chunk_size $fragement_length $quality_requirement $standard_error $radius_size $flooding_fraction $smooth_method $shaded_area $genomic_region_source_type.flanking_region_option_source_type4.flanking_region_option6 $genomic_region_source_type.flanking_region_option_source_type4.specify_flanking_region6 $genomic_region_source_type.interval_size4 $gene_list_source_type.gene_list_detail +#end if +#end if +#end if + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + + +**NGS PLOT Overview** + +ngs.plot_ needs an indexed bam file or a configuration file as an input to plot the peak trends and patterns cross the genomic region or users' interested region. ngs.plot.r will generate multiple files including avgprof file, heatmap file, and a zip file for replotting. + +.. _ngs.plot: http://code.google.com/p/ngsplot/ + +------ + +**Input formats** + +ngs.plot needs an indexed bam file or a configuration file as an input to plot the peak trends and patterns cross the genomic region or users' interested region. + +------ + +**Outputs** + +ngs.plot generates multiple files including avgprof file, heatmap file, and a zip file for replotting. + +------ + +**ngs.plot parameter list** + +This is a list of implemented ngs.plot options:: + +-G Genome name, currently supported: mm9, hg19, rn4, sacCer3(ensembl only) +-R Genomic regions to plot: tss, tes, genebody, exon, cgi or *.bed +-C Indexed bam file or a configuration file for multiplot +-E Gene list to subset regions +-T Image title(default=Noname), will be used in figure legend. +-F If you select genebody, exon or cgi, further information can be provided: for genebody: ChIP-seq(default), RNA-seq. RNA-seq is handled differently from ChIP-seq to cope with intronic regions. for exon: canonical(default), variant, promoter, polyA, altAcceptor, altDonor, altBoth. the alt's mean exons with alternative boundaries. for cgi, choose one based on location: ProximalPromoter(default), Promoter1k, Promoter3k, Genebody, Genedesert, OtherIntergenic, Pericentromere +-D Gene database: refseq(default), ensembl +-I Shall interval be larger than flanking in plot?(0 or 1) +-L Flanking region size (Size in bps. By default, when -R=tss, tes, genebody, -L=2000; when -R=exon, cgi, -L=500; when -R=*.bed, -L=1000) +-N Flanking region factor(will override flanking size) -S Randomly sample the regions for plot, must be:(0, 1] +-S Randomly sample the regions for plot, must be:(0, 1] (default=1, i.e. all regions) +-GO Gene order algorithm used in heatmaps: total(default), hc, max, prod, diff, pca and none(according to gene list supplied) +-CS Chunk size for loading genes in batch(default=100) +-FL Fragment length used to calculate physical coverage(default=150) +-MQ Mapping quality cutoff to filter reads(default=20) +-SE Shall standard errors be plotted?(0 or 1) +-RB The fraction of extreme values to be trimmed on both ends (default=0, 0.05 means 5% of extreme values will be trimmed) +-FC Flooding fraction:[0, 1), default=0.02 +-MW Moving window width to smooth avg. profiles, must be integer:1=no(default); 3=slightly; 5=somewhat; 9=quite; 13=super. +-H Opacity of shaded area, suggested value:[0, 0.5] (default=0, i.e. no shading, just curves) + + +