openvax / neoantigen-vaccine-pipeline

Bioinformatics pipeline for selecting patient-specific cancer neoantigen vaccines
Apache License 2.0
75 stars 25 forks source link

Add strelka2 to pipeline #136

Open julia326 opened 5 years ago

iskandr commented 5 years ago

I'm currently running Strelka2 in these two steps:

1) Create a Python file to run Strelka2:

configureStrelkaSomaticWorkflow.py \
            --normalBam $NORMAL_BAM \
            --tumorBam $TUMOR_BAM \
            --referenceFasta $REFERENCE_FASTA_PATH \
            --exome \
        --runDir $STRELKA_DIR";

2) Run the generated Python file

python runWorkflow.py \
                -m local \
                -j $NUMBER_PROCESSORS \
                -g $MEMORY_LIMIT_AS_INTEGER
iskandr commented 5 years ago

There's also an equivalent configureStrelkaGermlineWorkflow.py

iskandr commented 5 years ago

The --exome flag gets rid of any filtering assumptions about how the depth of a site should relate to the average depth on a chromosome. If you drop the flag then it assumes you have WGS data with relatively uniform coverage.

iskandr commented 5 years ago

The germline pipeline can take multiple BAMs and then the generated gVCF will have a sample column for each BAM. The gVCF is gigantic btw.