cbf290ada0Reving CoverageByRG to the new GATKReport API
Mauricio Carneiro
2012-08-15 16:09:51 -0400
9b84fa20bfFix an issue with the classpath for external builds in build.xml
David Roazen
2012-08-15 15:56:54 -0400
db92671b7fMerge branch 'master' of ssh://gsa4.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable
Guillermo del Angel
2012-08-15 13:32:07 -0400
7599f94296Updating large-scale validation script to new UG syntax
Guillermo del Angel
2012-08-15 13:31:45 -0400
87e41c83c5In AlleleCount stratification, check to make sure the AC (or MLEAC) is valid (i.e. not higher than number of chromosomes) and throw a User Error if it isn't. Added a test for bad AC.
Eric Banks
2012-08-14 15:02:30 -0400
8e3774fb0eFixing behavior of the --regenotype argument in SelectVariants to properly run in GenotypeGivenAlleles mode. Added integration tests to cover recent SV changes.
Eric Banks
2012-08-14 14:21:42 -0400
34b62fa092Two changes to SelectVariants: 1) don't add DP INFO annotation if DP wasn't used in the input VCF (it was adding DP=0 previously). 2) If MLEAC or MLEAF is present in the original VCF and the number of samples decreases, remove those annotations from the VC.
Eric Banks
2012-08-14 12:54:31 -0400
cfb994abd2Trivial removal of ununsed variable (mentioned in resolved JIRA entry)
Eric Banks
2012-08-13 22:55:02 -0400
f809f24afbRemoved SelectHeader's --include_reference_name option since the reference is always included. In SelectHeaders instead of including the path to the file, only include the name of the reference since dbGaP does not like paths in headers.
Khalid Shakir
2012-08-13 16:49:27 -0400
22b4466cf5Added setupRetry() to modify jobs when Queue is run with '-retry' and jobs are about to restart after an error. Implemented a mixin called "RetryMemoryLimit" which will by default double the memory. GridEngine memory request parameter can be selected on the command line via '-resMemReqParam mem_free' or '-resMemReqParam virtual_free'. Java optimizations now enabled by default: - Only 4 GC threads instead of each job using java's default O(number of cores) GC threads. Previously on a machine with N cores if you have N jobs running and java allocates N GC threads by default, then the machines are using up to N^2 threads if all jobs are in heavy GC (thanks elauzier). - Exit if GC spends more than 50% of time in GC (thanks ktibbett). - Exit if GC reclaims lest than 10% of max heap (thanks ktibbett). Added a -noGCOpt command line option to disable new java optimizations.
Khalid Shakir
2012-08-13 15:37:35 -0400
6ad75d2f5cReverting changes to BCF2 ranges
Mark DePristo
2012-08-13 15:06:28 -0400
4d3fad38e9Increase allowable range for BCF2 by -1 on low-end
Mark DePristo
2012-08-13 14:20:26 -0400
4cbd11faf5Fixed spelling error in BQSR.R
Mark DePristo
2012-08-13 10:01:33 -0400
aab417c94dFix missing argument in unittest
Mark DePristo
2012-08-12 13:58:14 -0400
f032e0aba4A bit better output for ContextCovariate context size logging
Mark DePristo
2012-08-12 13:45:52 -0400
243af0adb1Expanded the BQSR reporting script
Mark DePristo
2012-08-12 13:44:07 -0400
458bbdee8fAdd useful logger.info telling us the mismatch and indel context sizes
Mark DePristo
2012-08-12 10:27:05 -0400
6fefdaf428"update integration tests in CombineVariantsIntegrationTest"
Ami Levy Moonshine
2012-08-10 17:00:35 -0400
4968daf0a5update integration tests at CombineVariantsIntegrationTest
Ami Levy Moonshine
2012-08-10 16:58:05 -0400
1a87f67258Merge branch 'master' of ssh://gsa2.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable
Eric Banks
2012-08-10 15:11:12 -0400
40f0320a1cWhen adding a unit test to LIBS for X and = CIGAR operators, I uncovered a bug with the implementation of the ReadBackedPileup.depthOfCoverage() method.
Eric Banks
2012-08-10 14:58:29 -0400
eca9613356Adding support of X and = CIGAR operators to the GATK
Eric Banks
2012-08-10 14:54:07 -0400
949ed207cacapMaxAllelesForIndels -> capMaxAltAllelesForIndels
Joel Thibault
2012-08-08 15:52:24 -0400
b17edaad66Change memoryValues to List[Double]
Joel Thibault
2012-07-31 11:15:54 -0400
32a66b5ae4Add -nt parameter
Joel Thibault
2012-07-31 10:31:17 -0400
2a113977a9Resolving merge conflicts with the new MD5s
Ryan Poplin
2012-08-10 11:47:00 -0400
5f82ffd5d8Adding LowQual filter to the output of the HaplotypeCaller.
Ryan Poplin
2012-08-10 11:25:14 -0400
d7d7ccf789Revert unintentional license change
David Roazen
2012-08-09 17:10:47 -0400
d56a4631dcUpdate cofoja version in build.xml
David Roazen
2012-08-09 17:08:43 -0400
68fb04b8f7Merge branch 'master' of ssh://gsa2.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable into testing
Ami Levy Moonshine
2012-08-09 16:48:22 -0400
3362584014Updating cofoja to the latest version
Mark DePristo
2012-08-09 15:34:20 -0400
06258c8a01BCF2 optimizations
Mark DePristo
2012-08-09 15:34:07 -0400
c6bd9b15ffBCF2 optimizations
Mark DePristo
2012-08-09 11:09:09 -0400
9a0dda71d4BCF2 optimizations
Mark DePristo
2012-08-09 09:00:44 -0400
9887bc4410Merge branch 'master' of ssh://gsa2.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable
Ryan Poplin
2012-08-09 16:31:06 -0400
f4c72a26d5A few quick, minor findbugs fixes.
Ryan Poplin
2012-08-09 16:30:58 -0400
c7f22e410fA few quick, minor findbugs fixes.
Ryan Poplin
2012-08-09 16:22:08 -0400
abb168e1baMerged bug fix from Stable into Unstable
Mauricio Carneiro
2012-08-09 16:09:58 -0400
67d4148b32Fixing but reported by Thomas in the forum where reads were soft-clipped beyond the limits of the contig and ReduceReads was failing with a NoSuchElement exception. Now we hard clip anything that goes beyond the boundaries of the contig.
Mauricio Carneiro
2012-08-09 15:58:18 -0400
58420098acMerged bug fix from Stable into Unstable
Mauricio Carneiro
2012-08-09 13:02:23 -0400
c6132ebe26Fixed divide by zero bug when downsampler goes over regions where reads are all filtered out. Added Guillermo's bug report as an integration test
Mauricio Carneiro
2012-08-09 13:02:11 -0400
def077c4e5There's actually a subtle but important difference between foo++ and ++foo
Eric Banks
2012-08-09 12:42:50 -0400
e48727dae3Merge branch 'master' of ssh://gsa2.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable
Ryan Poplin
2012-08-09 10:31:10 -0400
5be7e0621dMerge branch 'master' of ssh://gsa4.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable
Guillermo del Angel
2012-08-09 09:58:34 -0400
71ee8d87b3Rename per-sample ML allelic fractions and counts so that they don't have the same name as the per-site INFO fields, and clarify wording in VCF header
Guillermo del Angel
2012-08-09 09:58:20 -0400
6230b49a86Merge branch 'master' of ssh://gsa2.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable
Eric Banks
2012-08-08 21:44:43 -0400
35cec8530cMake coverage threshold in FindCoveredIntervals a command-line argument
Eric Banks
2012-08-08 21:44:24 -0400
250ffd2ad7Merged bug fix from Stable into Unstable
Mauricio Carneiro
2012-08-08 15:50:07 -0400
78c1556186Fixing ReduceReads downsampling bug -- downsampled reads were not being excluded from the read window, causing them to trail back and get caught by the sliding window exception
Mauricio Carneiro
2012-08-08 15:49:31 -0400
1223d77546Removing argument from HaplotypeCaller that was made unneccesary by recent improvements to triggering around large events
Ryan Poplin
2012-08-08 15:13:20 -0400
0a2a646a52Other random FindBugs fixes
Eric Banks
2012-08-08 14:56:27 -0400
f652d7806eFindBugs found an infinite loop in the code
Eric Banks
2012-08-08 14:44:49 -0400
4c84cc9486Quick pass of FindBugs 'should be static inner class' fixes.
Eric Banks
2012-08-08 14:42:06 -0400
a0196c9f5bQuick pass of FindBugs 'method invokes inefficient Number constructor' fixes.
Eric Banks
2012-08-08 14:34:16 -0400
4b2e3cec0bQuick pass of FindBugs 'inefficient use of keySet iterator instead of entrySet iterator' fixes for core tools.
Eric Banks
2012-08-08 14:29:41 -0400
3e2752667cIntermediate checkin for ReducedReads with HaplotypeCaller - change min read count over k-mer to average count over k-mer when doing assembly of a reduced read (not optimal, currently trying max and then will decide on best approach), fix merge conflicts
Guillermo del Angel
2012-08-08 12:07:33 -0400
a7811d673fUpdate URL for phone home / GATK key documentation output by the GATK upon error
David Roazen
2012-08-08 09:29:54 -0400
cda8d944b7Bugfixes for BCF with VQSR
Mark DePristo
2012-08-07 17:22:27 -0400
80b94a4f9aAdaptiveContexts implement pruning to a given chi2 p value
Mark DePristo
2012-08-07 16:36:57 -0400
982c735c76VisualizeAdaptiveTree now considers only leaf nodes when computing max/min penalty
Mark DePristo
2012-08-07 08:38:16 -0400
15085bf03eThe UnifiedGenotyper now makes use of base insertion and base deletion quality scores if they exist in the reads.
Ryan Poplin
2012-08-07 13:58:22 -0400
2c76f71a03Update -maxAlleles argument in integration tests
Eric Banks
2012-08-06 22:48:04 -0400
c66a896b8eFix UG integration test broken by new -maxAltAlleles nomenclature
Guillermo del Angel
2012-08-06 21:29:21 -0400
97c5ed4febMerge branch 'master' of ssh://gsa4.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable
Guillermo del Angel
2012-08-06 20:22:31 -0400
238d55cb61Fixes for running HaplotypeCaller with reduced reads: a) minor refactoring, pulled out code to compute mean representative count to ReadUtils, b) Don't use min representative count over kmer when constructing de Bruijn graph - this creates many paths with multiplicity=1 and makes us lose a lot of SNP's at edge of capture targets. Use mean instead
Guillermo del Angel
2012-08-06 20:22:12 -0400
d2feb5d742Merge branch 'master' of ssh://gsa2.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable
Eric Banks
2012-08-06 15:06:45 -0400
4b5e1aed1fWe might completely overhaul the AssignSomaticStatus tool at some point, but for now I've just tweaked it a little so that it runs well with the TALEN data: fixed VCF header bugs and changed it so that any multi-sample VCF can be used and the user specifies which of the samples is the normal and which should be used as the tumor (only 1 T/N pair allowed now).
Eric Banks
2012-08-06 15:06:36 -0400
00858f16a6Deleting empty unit test for AdaptiveContexts
Mark DePristo
2012-08-06 12:57:55 -0400
f1c30c3a59Merge branch 'master' of ssh://gsa2.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable
Ryan Poplin
2012-08-06 12:02:26 -0400
d85b38e4daUpdating HaplotypeCaller integration tests
Ryan Poplin
2012-08-06 12:02:19 -0400
44f160f29findelGOP and indelGCP are now advanced, not hidden arguments
Mark DePristo
2012-08-06 11:42:36 -0400
2f004665fbFixing public -> private dep
Mark DePristo
2012-08-06 11:42:12 -0400
7bf5ca51eeMajor bugfix for adaptive contexts
Mark DePristo
2012-08-04 11:28:17 -0400
b4841548f1Bug fixes and misc. improvements to running the adaptive context tools
Mark DePristo
2012-08-03 17:53:46 -0400
b8709d8c67Merge branch 'master' of ssh://gsa2.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable
Ryan Poplin
2012-08-06 11:41:28 -0400
afa70a13a9Merge branch 'master' of ssh://gsa2.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable
Ryan Poplin
2012-08-06 11:40:27 -0400
973d1d47edMerging together the computeDiploidHaplotypeLikelihoods functions in the HaplotypeCaller's LikelihoodEngine so they both benefit from the ReducedRead's RepresentativeCount
Ryan Poplin
2012-08-06 11:40:07 -0400
210db5ec27Update -maxAlleles argument to -maxAltAlleles to make it more accurate. The hidden GSA production -capMaxAllelesForIndels argument also gets updated.
Eric Banks
2012-08-06 11:31:18 -0400
8f95a03bb6Prevent NumberFormatExceptions when parsing the VCF POS field
Eric Banks
2012-08-06 11:19:54 -0400
b7eec2fd0eBug fixes related to the changes in allele padding. If a haplotype started with an insertion it led to array index out of bounds. Haplotype allele insert function is now very simple because all alleles are treated the same way. HaplotypeUnitTest now uses a variant context instead of creating Allele objects directly.
Ryan Poplin
2012-08-05 12:29:10 -0400
e1bba91836Ready for full-scale evaluation adaptive BQSR contexts
Mark DePristo
2012-08-03 16:01:27 -0400
d2e8eb7b23Fixed 2 haplotype caller unit tests: a) new interface for addReadLikelihoods() including read counts, b) disable test that test basic DeBruijn graph assembly, not ready yet
Guillermo del Angel
2012-08-03 14:26:51 -0400
c3b6e2b143Merge branch 'master' of ssh://gsa2.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable
Ryan Poplin
2012-08-03 13:14:43 -0400
ff80f17721Using PathComparatorTotalScore in the assembly graph traversal does a better job of capturing low frequency branches that are inside high frequnecy haplotypes.
Ryan Poplin
2012-08-03 13:14:37 -0400
6f8e7692d4Merge branch 'master' of ssh://gsa4.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable
Guillermo del Angel
2012-08-03 12:24:37 -0400
9e25b209e0First pass of implementation of Reduced Reads with HaplotypeCaller. Main changes: a) Active region: scale PL's by representative count to determine whether region is active. b) Scale per-read, per-haplotype likelihoods by read representative counts. A read representative count is (temporarily) defined as the average representative count over all bases in read, TBD whether this is good enough to avoid biases in GL's. c) DeBruijn assembler inserts kmers N times in graph, where N is min representative count of read over kmer span - TBD again whether this is the best approach. d) Bug fixes in FragmentUtils: logic to merge fragments was wrong in cases where there is discrepancy of overlaps between unclipped/soft clipped bases. Didn't affect things before but RR makes prevalence of hard-clipped bases in CIGARs more prevalent so this was exposed. e) Cache read representative counts along with read likelihoods associated with a Haplotype. Code can/should be cleaned up and unified with PairHMMIndelErrorModelCode, as well as refactored to support arbitrary ploidy in HaplotypeCaller
Guillermo del Angel
2012-08-03 12:24:23 -0400
8817fc70d1Merged bug fix from Stable into Unstable
Ryan Poplin
2012-08-03 10:45:01 -0400
f40d0a0a28Updating VQSR to work with the MNP and symbolic variants that are coming out of the HaplotypeCaller. Integration tests change because of the MNPs in dbSNP.
Ryan Poplin
2012-08-03 10:44:36 -0400
51bd03cc36Add RemoveProgramRecords annotation to ActiveRegionWalker
Joel Thibault
2012-08-02 14:46:00 -0400
addbfd6437Add a RemoveProgramRecords annotation
Joel Thibault
2012-08-02 14:21:30 -0400
524d7ea306Choose whether to keep program records based on Walker
Joel Thibault
2012-07-31 15:09:36 -0400
e04989f76dBugfix for new PASS position in dictionary in BCF2
Mark DePristo
2012-08-03 09:42:21 -0400
d22b8cf86bVisualizeContextTree now loops over M, I, and D states generating trees and analyzes
Mark DePristo
2012-08-02 16:45:13 -0400
fb5dabce18Update BCF2 to include a minor version number so we can rev (and report errors) with BCF2
Mark DePristo
2012-08-02 16:44:39 -0400
b4f4d86c77Merge branch 'master' of ssh://gsa2.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable
Eric Banks
2012-08-02 11:49:03 -0400
3ece4c4993Merged bug fix from Stable into Unstable
Ryan Poplin
2012-08-02 11:41:36 -0400
6f7a236cfcMerge branch 'master' of ssh://gsa2.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/stable
Ryan Poplin
2012-08-02 11:41:20 -0400
cb8bc18aebFix for error in HaplotypeCaller. HC has a UG argument collection for the UG engine but some of those arguments aren't appropriate to set.
Ryan Poplin
2012-08-02 11:41:06 -0400
e3f89fb054Missing/malformed GATK report files are user errors
Eric Banks
2012-08-02 11:33:21 -0400
cc01f844d4Merged bug fix from Stable into Unstable
Eric Banks
2012-08-02 11:25:28 -0400