Quantcast
Viewing all articles
Browse latest Browse all 12345

MuTect2 with Queue

Hi,

I'm trying to run Mutect2 with matched normal-tumor whole genome samples and I want to use Queue. I know it will take a lot of time with WGS samples and I want to use Queue as mentioned here http://gatkforums.broadinstitute.org/gatk/discussion/6559/mutect2-runs-much-slower-than-mutect-1-17.

Mutect2 runs fine directly with my samples, and Queue works with the example files that came with the .bz2 file I downloaded. However, I cannot get the two to run together.

This is the error I get:

INFO 19:41:50,429 QScriptManager - Compiling 1 QScript
ERROR 19:41:51,062 QScriptManager - Mutect2.scala:20: type mismatch;
found : java.io.File
required: Seq[java.io.File]
ERROR 19:41:51,064 QScriptManager - mutect2.cosmic = new File("mutect/b37_cosmic_v54_120711.vcf")
ERROR 19:41:51,065 QScriptManager - ^
ERROR 19:41:51,142 QScriptManager - two errors found

ERROR --
ERROR stack trace

org.broadinstitute.gatk.queue.QException: Compile of ../Analysis/mutect/Mutect2.scala failed with 2 errors
at org.broadinstitute.gatk.queue.QScriptManager.loadScripts(QScriptManager.scala:79)
at org.broadinstitute.gatk.queue.QCommandLine.org$broadinstitute$gatk$queue$QCommandLine$$qScriptPluginManager$lzycompute(QCommandLine.scala:94)
at org.broadinstitute.gatk.queue.QCommandLine.org$broadinstitute$gatk$queue$QCommandLine$$qScriptPluginManager(QCommandLine.scala:92)
at org.broadinstitute.gatk.queue.QCommandLine.getArgumentSources(QCommandLine.scala:229)
at org.broadinstitute.gatk.utils.commandline.CommandLineProgram.start(CommandLineProgram.java:213)
at org.broadinstitute.gatk.utils.commandline.CommandLineProgram.start(CommandLineProgram.java:157)
at org.broadinstitute.gatk.queue.QCommandLine$.main(QCommandLine.scala:61)
at org.broadinstitute.gatk.queue.QCommandLine.main(QCommandLine.scala)

ERROR ------------------------------------------------------------------------------------------
ERROR A GATK RUNTIME ERROR has occurred (version 3.6-0-g89b7209):
ERROR
ERROR This might be a bug. Please check the documentation guide to see if this is a known problem.
ERROR If not, please post the error message, with stack trace, to the GATK forum.
ERROR Visit our website and forum for extensive documentation and answers to
ERROR commonly asked questions https://www.broadinstitute.org/gatk
ERROR
ERROR MESSAGE: Compile of ../Analysis/mutect/Mutect2.scala failed with 2 errors
ERROR ------------------------------------------------------------------------------------------

INFO 19:41:51,235 QCommandLine - Shutting down jobs. Please wait...

There seems to be some problem with Java (I'm running Java 1.8) This is the scala code I wrote after seeing different discussions here. I haven't really worked with Scala.

   import org.broadinstitute.gatk.queue.QScript
   import org.broadinstitute.gatk.queue.extensions.gatk._
   class Mutect2 extends QScript {

  // Script Arguments passed from command line
    @Input(doc="Normal or unaffected sample", shortName="normal", required=true)
    var normalIn: File = _
    @Input(doc="Tumor or affected sample", shortName="tumor", required=true)
    var tumorIn: File = _
    @Argument(shortName = "L",  required=false, doc = "Intervals file")
    var intervalsFile: List[File] = Nil
    @Argument(shortName = "o",  required=true, doc = "Output file")
    var outputFile: File = _

    // Add functions hard-coded in the script
    def script() {
            val mutect2 = new MuTect2
            //mutect2.jarFile = new File("/software/UHTS/Analysis/GenomeAnalysisTK/3.6/bin/GenomeAnalysisTK")
            mutect2.R = new File("/scratch/cluster/monthly/fsantoni/index/hg19/human_g1k_v37.fasta")
            mutect2.cosmic = new File("/data3/unige/dmg/cmorey/Analysis/mutect/b37_cosmic_v54_120711.vcf")
            mutect2.dbsnp = new File("/data3/unige/dmg/cmorey/Analysis/mutect/dbsnp_132_b37.leftAligned.vcf.gz")
            mutect2.intervalsString = intervalsFile
            mutect2.input_file = List(new TaggedFile(tumorIn, "tumor"), new TaggedFile(normalIn, "normal"))
            mutect2.out = outputFile
    add(mutect2)
    }
    }

And I run it with
java -jar Queue.jar -S ../Analysis/mutect/Mutect2.scala -normal normal.bam -tumor tumor.bam -o Output_mutect2.vcf


Viewing all articles
Browse latest Browse all 12345

Trending Articles