diff --git a/biopet-framework/pom.xml b/biopet-framework/pom.xml
index d735e4bbd703a6dad97a6768a515bd7a6044e270..8385a13dd70e2982641e41a7c61bc26fc370943d 100644
--- a/biopet-framework/pom.xml
+++ b/biopet-framework/pom.xml
@@ -49,7 +49,7 @@
<dependency>
<groupId>org.broadinstitute.gatk</groupId>
<artifactId>gatk-queue-package-distribution</artifactId>
- <version>3.2</version>
+ <version>3.3</version>
</dependency>
<dependency>
<groupId>io.argonaut</groupId>
diff --git a/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/extensions/gatk/HaplotypeCaller.scala b/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/extensions/gatk/HaplotypeCaller.scala
index 6400354553b70f25984db39a5ec5070239017991..a83ee965e86207e1a0e73d48cdeb8ba86426ce8b 100644
--- a/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/extensions/gatk/HaplotypeCaller.scala
+++ b/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/extensions/gatk/HaplotypeCaller.scala
@@ -27,12 +27,10 @@ class HaplotypeCaller(val root: Configurable) extends org.broadinstitute.gatk.qu
if (config("inputtype", default = "dna").getString == "rna") {
dontUseSoftClippedBases = config("dontusesoftclippedbases", default = true)
- recoverDanglingHeads = config("recoverdanglingheads", default = true)
stand_call_conf = config("stand_call_conf", default = 5)
stand_emit_conf = config("stand_emit_conf", default = 0)
} else {
dontUseSoftClippedBases = config("dontusesoftclippedbases", default = false)
- recoverDanglingHeads = config("recoverdanglingheads", default = false)
stand_call_conf = config("stand_call_conf", default = 5)
stand_emit_conf = config("stand_emit_conf", default = 0)
}
diff --git a/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/tools/CheckAllelesVcfInBam.scala b/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/tools/CheckAllelesVcfInBam.scala
index 2d822e2c662ad2075b09a0e70e47d6fc9dca1aa1..315a8db8c8b2743c81cdb618685a6eabf381891f 100644
--- a/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/tools/CheckAllelesVcfInBam.scala
+++ b/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/tools/CheckAllelesVcfInBam.scala
@@ -1,6 +1,7 @@
package nl.lumc.sasc.biopet.tools
import htsjdk.samtools.SAMFileReader
+import htsjdk.samtools.QueryInterval
import htsjdk.samtools.SAMRecord
import htsjdk.variant.variantcontext.VariantContext
import htsjdk.variant.variantcontext.VariantContextBuilder
@@ -84,7 +85,7 @@ object CheckAllelesVcfInBam extends ToolCommand {
val countReports: Map[String,CountReport] = bamReaders.map(x => (x._1, new CountReport))
val refAllele = vcfRecord.getReference.getBaseString
for ((sample, bamReader) <- bamReaders) {
- val queryInterval = new SAMFileReader.QueryInterval(bamHeaders(sample).getSequenceIndex(vcfRecord.getChr),
+ val queryInterval = new QueryInterval(bamHeaders(sample).getSequenceIndex(vcfRecord.getChr),
vcfRecord.getStart, vcfRecord.getStart + refAllele.size - 1)
val bamIter = bamReader.query(Array(queryInterval), false)
diff --git a/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/tools/FindRepeatsPacBio.scala b/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/tools/FindRepeatsPacBio.scala
index 06ae6b64f4e556c3dd9e01d2ee14c5231161b226..add0849bc4e60e4df8662c5432dad0c197b26631 100644
--- a/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/tools/FindRepeatsPacBio.scala
+++ b/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/tools/FindRepeatsPacBio.scala
@@ -16,7 +16,9 @@
package nl.lumc.sasc.biopet.tools
+import htsjdk.samtools.QueryInterval
import htsjdk.samtools.SAMFileReader
+import htsjdk.samtools.ValidationStringency
import htsjdk.samtools.SAMRecord
import java.io.File
import nl.lumc.sasc.biopet.core.ToolCommand
@@ -41,7 +43,7 @@ object FindRepeatsPacBio extends ToolCommand {
val argsParser = new OptParser
val commandArgs: Args = argsParser.parse(args, Args()) getOrElse sys.exit(1)
val bamReader = new SAMFileReader(commandArgs.inputBam)
- bamReader.setValidationStringency(SAMFileReader.ValidationStringency.SILENT)
+ bamReader.setValidationStringency(ValidationStringency.SILENT)
val bamHeader = bamReader.getFileHeader
val header = List("chr", "startPos", "stopPos","Repeat_seq", "repeatLength",
@@ -51,7 +53,7 @@ object FindRepeatsPacBio extends ToolCommand {
for (bedLine <- Source.fromFile(commandArgs.inputBed).getLines;
val values = bedLine.split("\t"); if values.size >= 3) {
- val interval = new SAMFileReader.QueryInterval(bamHeader.getSequenceIndex(values(0)), values(1).toInt, values(2).toInt)
+ val interval = new QueryInterval(bamHeader.getSequenceIndex(values(0)), values(1).toInt, values(2).toInt)
val bamIter = bamReader.query(Array(interval), false)
val results = for (samRecord <-bamIter) yield procesSamrecord(samRecord, interval)
val chr = values(0)
@@ -100,7 +102,7 @@ object FindRepeatsPacBio extends ToolCommand {
}
}
- def procesSamrecord(samRecord:SAMRecord, interval:SAMFileReader.QueryInterval): Option[Result] = {
+ def procesSamrecord(samRecord:SAMRecord, interval:QueryInterval): Option[Result] = {
val readStartPos = List.range(0, samRecord.getReadBases.length)
.find(x => samRecord.getReferencePositionAtReadPosition(x) >= interval.start)
var readPos = if (readStartPos.isEmpty) return None else readStartPos.get
diff --git a/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/tools/WipeReads.scala b/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/tools/WipeReads.scala
index 4e3d118cea246117ec9f7724cda6a9b177334fdd..a990bab67fb7dc2ed67ad91304e74814fb9ba5cc 100644
--- a/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/tools/WipeReads.scala
+++ b/biopet-framework/src/main/scala/nl/lumc/sasc/biopet/tools/WipeReads.scala
@@ -10,7 +10,8 @@ import scala.collection.JavaConverters._
import htsjdk.samtools.AlignmentBlock
import htsjdk.samtools.SAMFileReader
-import htsjdk.samtools.SAMFileReader.QueryInterval
+import htsjdk.samtools.QueryInterval
+import htsjdk.samtools.ValidationStringency
import htsjdk.samtools.SAMFileWriterFactory
import htsjdk.samtools.SAMRecord
import htsjdk.tribble.AbstractFeatureReader.getFeatureReader
@@ -161,7 +162,7 @@ object WipeReads extends ToolCommand {
new SAMFileReader(inBAM, inBAMIndex)
else {
val sfr = new SAMFileReader(inBAM)
- sfr.setValidationStringency(SAMFileReader.ValidationStringency.LENIENT)
+ sfr.setValidationStringency(ValidationStringency.LENIENT)
if (!sfr.hasIndex)
throw new IllegalStateException("Input BAM file must be indexed")
else
@@ -314,7 +315,7 @@ object WipeReads extends ToolCommand {
.setUseAsyncIo(async)
val templateBAM = new SAMFileReader(inBAM)
- templateBAM.setValidationStringency(SAMFileReader.ValidationStringency.LENIENT)
+ templateBAM.setValidationStringency(ValidationStringency.LENIENT)
val targetBAM = factory.makeBAMWriter(templateBAM.getFileHeader, true, outBAM)
diff --git a/mvn_install_queue.sh b/mvn_install_queue.sh
index 3e5d12fdd6d40cd63c8d03cc29f23a5033a63cdc..3bd9be7ca612b7e5f124afb8e0f47dc7a411cce5 100755
--- a/mvn_install_queue.sh
+++ b/mvn_install_queue.sh
@@ -1,6 +1,6 @@
-JAR=/data/DIV5/SASC/common/programs/Queue-3.2-2/Queue.jar
+JAR=/data/DIV5/SASC/common/programs/Queue-3.3-0/Queue.jar
ID=gatk-queue-package-distribution
-VERSION=3.2
+VERSION=3.3
GROUP=org.broadinstitute.gatk
mvn install:install-file -Dfile=$JAR -DgroupId=$GROUP -DartifactId=$ID -Dversion=$VERSION -Dpackaging=jar