Skip to content
Open
Changes from all commits
Commits
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
8 changes: 4 additions & 4 deletions README.md
Original file line number Diff line number Diff line change
@@ -1,6 +1,6 @@
![](http://i.giphy.com/LwhhZsEHFQgSs.gif)

#Identifying recurrent mutations in cancer[<img target="_blank" align="right" border="0" alt="" src="http://www.cbioportal.org/images/cbioportal_logo.png" width="125" height="30">](http://www.cbioportal.org/) <img align="right" border="0" alt="" src="http://www.cbioportal.org/images/oncokb-flame.svg" width="30" height="30">
# Identifying recurrent mutations in cancer[<img target="_blank" align="right" border="0" alt="" src="http://www.cbioportal.org/images/cbioportal_logo.png" width="125" height="30">](http://www.cbioportal.org/) <img align="right" border="0" alt="" src="http://www.cbioportal.org/images/oncokb-flame.svg" width="30" height="30">



Expand All @@ -21,7 +21,7 @@ source("http://bioconductor.org/biocLite.R")
biocLite("IRanges","BSgenome.Hsapiens.UCSC.hg19")
```

####Usage:
#### Usage:
```
./hotspot_algo.R
--input-maf=[REQUIRED: mutation file]
Expand All @@ -42,7 +42,7 @@ Command to run hotspot algorithm on genes listed in file genes_of_interest.txt:
--output-file=sig_hotspots.txt
```

####Contents:
#### Contents:
**[ Required ]** `hotspot_algo.R` - R script to execute hotspot detection algorithm

**[ Required ]** `hotspot_algo.Rdata` - Rdata object with necessary files for algorithm (mutability, expression filters, etc)
Expand All @@ -53,7 +53,7 @@ Command to run hotspot algorithm on genes listed in file genes_of_interest.txt:

`minimalist_test_maf.txt` - minimalist MAF needed from maf2maf. [mskcc/maf2maf](https://github.com/mskcc/vcf2maf)

####Notes:
#### Notes:
`--align100mer` and `--align24mer` are optional filters based on how uniquely k-mer sequences align to a region of the hg19 genome. Note, both filters were used as part of this analysis. See more information at [ENCODE Mapability](http://genome.ucsc.edu/cgi-bin/hgFileUi?db=hg19&g=wgEncodeMapability).

The use of these filters will require downloading the 100-mer and 24-mer alignability tracks from UCSC that are not included here:
Expand Down