Changes between Version 5 and Version 6 of SOPs/Using
- Timestamp:
- 09/21/20 16:43:04 (4 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
SOPs/Using
v5 v6 1 1 == Using ngsplot to make stacked heatmaps and profiles of genes or genomic regions (like ChIP-seq peaks) == 2 [[https://github.com/shenlab-sinai/ngsplot/wiki | ngsplot]] is a great tool to make heatmaps and profiles of mapped reads (i.e. bam files of RNA-seq, ChIP-seq etc) across genes, transcription starts sites or other regions of interest. 3 The regions used can be genes included on the ngsplot annotation or custom regions included on a bed file. 2 3 [[https://github.com/shenlab-sinai/ngsplot/wiki | ngsplot]] is a great tool to make heatmaps and profiles of mapped reads (i.e. bam files of RNA-seq, ChIP-seq etc) across genes, transcription starts sites or other regions of interest. 4 The regions used can be genes included on the ngsplot annotation or custom regions included on a bed file. 4 5 5 6 * Here is an example of a ngsplot command using custom regions 6 7 {{{ 7 8 bsub ngs.plot.r -G hg19 -R bed -L 3000 -C config.txt -O Peak_Summits 8 9 9 * Example of the tab delimited config.txt file 10 using this example of the tab delimited config.txt file: 10 11 11 TF1.sorted.bampeaks.SUMMIT.bed "TF1"12 TF2.sorted.bampeaks.SUMMIT.bed "TF2"13 Pol2.sorted.bampeaks.SUMMIT.bed "Pol2"12 TF1.sorted.bam peaks.SUMMIT.bed "TF1" 13 TF2.sorted.bam peaks.SUMMIT.bed "TF2" 14 Pol2.sorted.bam peaks.SUMMIT.bed "Pol2" 14 15 }}} 15 16 * Here is an example of a ngsplot command using gene annotations within the ngsplot package 16 17 17 18 {{{ 18 # using transcription start sites (tss)19 # Using transcription start sites (tss) as the features of interest 19 20 bsub ngs.plot.r -G hg19 -R tss -L 3000 -C config.txt -O ngsplot_tss 20 21 21 # using gene body22 # Using gene bodies as the features of interest 22 23 bsub ngs.plot.r -G hg19 -R genebody -L 3000 -C config.txt -O ngsplot_tss 23 24 24 25 }}} 25 26 26 ngsplot uses a two-step normalization: i) coverage vectors' length i.e. gene length and ii) library size.27 [[https://github.com/shenlab-sinai/ngsplot/wiki/HowAreYaxisValsCalculated | ngsplot: How Are Y-axis Vals Calculated]] 27 ngsplot uses a two-step normalization: i) coverage vectors' length i.e., gene length and ii) library size. 28 See [[https://github.com/shenlab-sinai/ngsplot/wiki/HowAreYaxisValsCalculated | ngsplot: how the Y-axis values are calculated]]. 28 29 29 30