This package hooks into the jrflab/modules
sequencing pipeline for tailored / dynamic analysis.
From root project directory:
sh pascal/init.sh
subsets.txt
file as follows (space delimited):subsetName1 sample1 sample2
subsetName2 sample1 sample3 sample4
Builds a prevelance-ordered, subsetted heatmap + table of recurrent mutations from the master recurrent mutation table.
make recurrent_mutations
Rscript pascal/recurrence/subsetRecurrence.R
Generates heatmap .pdf and table files in recurrent_mutations directory within project root.
Performs base-truth recurrent mutation lookups in the .bam files of all samples within subsets.
make recurrent_mutations
Rscript pascal/recurrence/sufamRecurrent.R
Generates log file in log directory, vcf and sufam output .tsv in recurrent_mutations directory
Use Sufam to arbitrarily search for base-truth genomic events in a sample subset by specifying a vcf file.
#CHROM POS ID REF ALT
Rscript pascal/recurrence/sufamEventSearch.R
Rscript pascal/clonality/pyclone.R
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.