Link to CGHub manifest describing the contents of its database.
In [11]:
cghub_manifest = 'https://cghub.ucsc.edu/reports/SUMMARY_STATS/LATEST_MANIFEST.tsv'
CGHUB download url
In [12]:
CGHUB = 'https://cghub.ucsc.edu/cghub/data/analysis/download'
There paths are on the local file-system (cellar).
In [13]:
BCBIO = '/cellar/users/agross/sources/bcbio_nextgen_local/'
DBSNP = BCBIO + 'genomes/Hsapiens/GRCh37/variation/dbsnp_138.vcf.gz'
COSMIC = BCBIO + 'genomes/Hsapiens/GRCh37/variation/cosmic-v68-GRCh37.vcf.gz'
REFERENCE = BCBIO + 'genomes/Hsapiens/GRCh37/seq/GRCh37.fa'
Path to key for downloading data from CGHub. You need to contact CGHub to get a key.
In [14]:
KEY = '/cellar/users/agross/.ssh/cghub.key'
Path to MuTect. I was getting bugs with version 1.1.4 so make sure you use this version.
In [15]:
MUTECT_JAR = BCBIO + 'share/java/mutect/muTect-1.1.5.jar'
Path to SomaticIndelDetector jar. This is in a specific version of GATK, so make sure you are using version 2.2-2.
In [16]:
SID_JAR = '/cellar/users/hcarter/programs/GenomeAnalysisTK-2.2-2/GenomeAnalysisTK.jar'
Path to cache directory. This is only used in conjuntion with GT-Fuse. This is important because the bam files can pile up in the cash and eat up all of the space on a hard drive very quickly.
In [17]:
#CACHE = '/home/centos/cache/fusecache'
Number of processes you want to spawn with the variant calling. I'm taking a quick and dirty pass here and just running a bunch of bash scripts simultaniously. At some point we will probably switch to a scheduler, but we don't have one on our Annai VM currently and this seems to be working ok for now.
In [18]:
NUM_PROCESSES = 16
Directory to store the data on whatever machine you are running the scripts on (we use a VM).
In [9]:
#VM_DIRECTORY = '/home/centos/projects'
Local directory to spit out the bash scripts.
In [10]:
LOCAL_DIRECTORY = '/cellar/users/agross/scripts'