The current GATK version is 3.7-0
Examples: Monday, today, last week, Mar 26, 3/26/04

Howdy, Stranger!

It looks like you're new here. If you want to get involved, click one of these buttons!

Get notifications!

You can opt in to receive email notifications, for example when your questions get answered or when there are new announcements, by following the instructions given here.

Did you remember to?

1. Search using the upper-right search box, e.g. using the error message.
2. Try the latest version of tools.
3. Include tool and Java versions.
4. Tell us whether you are following GATK Best Practices.
5. Include relevant details, e.g. platform, DNA- or RNA-Seq, WES (+capture kit) or WGS (PCR-free or PCR+), paired- or single-end, read length, expected average coverage, somatic data, etc.
6. For tool errors, include the error stacktrace as well as the exact command.
7. For format issues, include the result of running ValidateSamFile for BAMs or ValidateVariants for VCFs.
8. For weird results, include an illustrative example, e.g. attach IGV screenshots according to Article#5484.
9. For a seeming variant that is uncalled, include results of following Article#1235.

Did we ask for a bug report?

Then follow instructions in Article#1894.

Formatting tip!

Wrap blocks of code, error messages and BAM/VCF snippets--especially content with hashes (#)--with lines with three backticks ( ``` ) each to make a code block as demonstrated here.

Jump to another community
Picard 2.9.0 is now available. Download and read release notes here.
GATK 3.7 is here! Be sure to read the Version Highlights and optionally the full Release Notes.

an error I confronted when i run the code below

zhangrui9xzhangrui9x ChinaMember Posts: 8
edited December 2013 in Ask the GATK team
import org.broadinstitute.sting.gatk.walkers.variantrecalibration
import org.broadinstitute.sting.gatk.walkers.variantrecalibration.VariantRecalibratorArgumentCollection
import org.broadinstitute.sting.gatk.walkers.variantrecalibration.VariantRecalibratorArgumentCollection.Mode
import org.broadinstitute.sting.queue.extensions.gatk.CountLoci
import org.broadinstitute.sting.queue.extensions.gatk._
import org.broadinstitute.sting.queue.function.JavaCommandLineFunction

class MyScript extends org.broadinstitute.sting.queue.QScript {
  @Input(doc="The reference file.", shortName="R")
  val reference:File= null
  @Input(doc="The left reads.", shortName="lr")
  val leftreads:File= null
  @Input(doc="The right reads.", shortName="rr")
  val rightreads:File= null

 def script ={
   class leftsaiFormationFunction extends CommandLineFunction{
      val leftsai:File=new File("leftsai.sai")
      def commandLine=required("bwa")+required("aln")+required(reference)+required(leftreads)+required(">")+
    val leftsaiFormation=new leftsaiFormationFunction()

    class rightsaiFormationFunction extends CommandLineFunction{
      val rightsai:File=new File("rightsai.sai")
      def commandLine=required("bwa")+required("aln")+required(reference)+required(rightreads)+required(">")+
    val rightsaiFormation=new rightsaiFormationFunction()

    class combinedSamFormationFunction extends CommandLineFunction{
      val leftsai:File=leftsaiFormation.leftsai
      val rightsai:File=rightsaiFormation.rightsai
      val combinedSam:File=new File("combined.sam")
      def commandLine=required("bwa")+required("sampe")+required(reference)+required(leftsai)+
    val combinedSamFormation=new combinedSamFormationFunction()



the error message says as follow:
lxd@lxd-System-Product-Name:~/Documents/queue$ java -jar Queue.jar -S 111111.scala -R human_g1k_v37.fasta -lr C166_Modified_1.fq -rr C166_Modified_2.fq -run -startFromScratch 
INFO  16:05:35,810 QScriptManager - Compiling 1 QScript 
INFO  16:05:39,696 QScriptManager - Compilation complete 
INFO  16:05:39,759 HelpFormatter - ---------------------------------------------------------------------- 
INFO  16:05:39,759 HelpFormatter - Queue v2.7-2-g6bda569, Compiled 2013/08/28 16:33:34 
INFO  16:05:39,759 HelpFormatter - Copyright (c) 2012 The Broad Institute 
INFO  16:05:39,759 HelpFormatter - For support and documentation go to 
INFO  16:05:39,760 HelpFormatter - Program Args: -S 111111.scala -R human_g1k_v37.fasta -lr C166_Modified_1.fq -rr C166_Modified_2.fq -run -startFromScratch 
INFO  16:05:39,760 HelpFormatter - Date/Time: 2013/12/06 16:05:39 
INFO  16:05:39,760 HelpFormatter - ---------------------------------------------------------------------- 
INFO  16:05:39,760 HelpFormatter - ---------------------------------------------------------------------- 
INFO  16:05:39,766 QCommandLine - Scripting MyScript 
INFO  16:05:39,800 QCommandLine - Added 3 functions 
INFO  16:05:39,801 QGraph - Generating graph. 
INFO  16:05:39,811 QGraph - Running jobs. 
INFO  16:05:39,812 QGraph - Removing outputs from previous runs. 
INFO  16:05:39,840 FunctionEdge - Starting:  'bwa'  'aln'  'human_g1k_v37.fasta'  'C166_Modified_1.fq'  '>'  '/home/lxd/Documents/queue/leftsai.sai'  
INFO  16:05:39,840 FunctionEdge - Output written to /home/lxd/Documents/queue/leftsai.sai.out 
INFO  16:05:39,858 FunctionEdge - Starting:  'bwa'  'aln'  'human_g1k_v37.fasta'  'C166_Modified_2.fq'  '>'  '/home/lxd/Documents/queue/rightsai.sai'  
INFO  16:05:39,858 FunctionEdge - Output written to /home/lxd/Documents/queue/rightsai.sai.out 
INFO  16:05:39,862 QGraph - 1 Pend, 2 Run, 0 Fail, 0 Done 
ERROR 16:06:09,834 FunctionEdge - Error:  'bwa'  'aln'  'human_g1k_v37.fasta'  'C166_Modified_1.fq'  '>'  '/home/lxd/Documents/queue/leftsai.sai'  
ERROR 16:06:09,839 FunctionEdge - Contents of /home/lxd/Documents/queue/leftsai.sai.out:
/home/lxd/Documents/queue/.queue/tmp/.exec6752111377935879956: 2: /home/lxd/Documents/queue/.queue/tmp/.exec6752111377935879956: bwa: not found 
ERROR 16:06:09,842 FunctionEdge - Error:  'bwa'  'aln'  'human_g1k_v37.fasta'  'C166_Modified_2.fq'  '>'  '/home/lxd/Documents/queue/rightsai.sai'  
ERROR 16:06:09,842 FunctionEdge - Contents of /home/lxd/Documents/queue/rightsai.sai.out:
/home/lxd/Documents/queue/.queue/tmp/.exec7253750678399969841: 2: /home/lxd/Documents/queue/.queue/tmp/.exec7253750678399969841: bwa: not found 
INFO  16:06:09,843 QGraph - Writing incremental jobs reports... 
INFO  16:06:09,844 QJobsReporter - Writing JobLogging GATKReport to file /home/lxd/Documents/queue/111111.jobreport.txt 
INFO  16:06:09,869 QGraph - 1 Pend, 0 Run, 2 Fail, 0 Done 
INFO  16:06:09,870 QCommandLine - Writing final jobs report... 
INFO  16:06:09,870 QJobsReporter - Writing JobLogging GATKReport to file /home/lxd/Documents/queue/111111.jobreport.txt 
INFO  16:06:09,876 QJobsReporter - Plotting JobLogging GATKReport to file /home/lxd/Documents/queue/111111.jobreport.pdf 
WARN  16:06:10,262 RScriptExecutor - RScript exited with 1. Run with -l DEBUG for more info. 
INFO  16:06:10,263 QCommandLine - Done with errors 
INFO  16:06:10,264 QGraph - ------- 
INFO  16:06:10,266 QGraph - Failed:   'bwa'  'aln'  'human_g1k_v37.fasta'  'C166_Modified_1.fq'  '>'  '/home/lxd/Documents/queue/leftsai.sai'  
INFO  16:06:10,266 QGraph - Log:     /home/lxd/Documents/queue/leftsai.sai.out 
INFO  16:06:10,266 QGraph - ------- 
INFO  16:06:10,268 QGraph - Failed:   'bwa'  'aln'  'human_g1k_v37.fasta'  'C166_Modified_2.fq'  '>'  '/home/lxd/Documents/queue/rightsai.sai'  
INFO  16:06:10,268 QGraph - Log:     /home/lxd/Documents/queue/rightsai.sai.out 
INFO  16:06:10,269 QCommandLine - Script failed with 3 total jobs 

It will be really appreciated for your attention and answer!!!!!

Post edited by Geraldine_VdAuwera on

Best Answer


Sign In or Register to comment.