d8a4d78854Bugfix for exceptions with unknown source whose error was not being shown in tableau
Mark DePristo
2012-01-31 10:49:06 -0500
faba3dd530Merge branch 'master' of ssh://chartl@ni.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable
Christopher Hartl
2012-01-31 10:25:29 -0500
17dbe9a95dA few cleanups in the LocusIteratorByState
Mauricio Carneiro
2012-01-30 21:37:02 -0500
e7ace8efc4Fix NullPointerException caused in cases with too many ALT alleles
Menachem Fromer
2012-01-30 21:00:16 -0500
f9162ea705Merge branch 'master' of ssh://nickel.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable
Ryan Poplin
2012-01-30 19:45:19 -0500
abb91cf26bIncreasing the size of the active regions that are produced by the active probability integrator, more context is needed to call more complex events
Ryan Poplin
2012-01-30 15:36:12 -0500
f1e07f169eOnly apply filters if there are filters to apply
Menachem Fromer
2012-01-30 02:22:13 -0500
d1aa5204d7Merge branch 'master' of ssh://copper.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable
Menachem Fromer
2012-01-29 23:39:34 -0500
3186f0f1b0Try more memory and fewer ALT alleles so that we don't run out of memory
Menachem Fromer
2012-01-29 23:38:32 -0500
3164c8dee5S3 upload now directly creates the XML report in memory and puts that in S3
Mark DePristo
2012-01-29 15:14:58 -0500
0e17cbbce9Merged bug fix from Stable into Unstable
Menachem Fromer
2012-01-27 16:03:16 -0500
a9671b73caFix to permit proper handling of mapping qualities between 128 to 255 (which get converted to byte values of -128 to -1)
Menachem Fromer
2012-01-27 16:01:30 -0500
f7ac1f4a69Merge branch 'master' of ssh://nickel.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable
Ryan Poplin
2012-01-27 15:12:55 -0500
fc08235ff3Bug fix in active region traversal, locusView.getNext() skips over pileups with zero coverage but still need to count them in the active probability integrator
Ryan Poplin
2012-01-27 15:12:37 -0500
052a4bdb9cTurning off PHONE HOME option in the MDCP
Mauricio Carneiro
2012-01-27 11:13:30 -0500
0de6aca29bMerge branch 'master' of ssh://copper.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable
Menachem Fromer
2012-01-26 18:26:41 -0500
d54e237671Take advantage of Eric's fix for multiAllelic AC calculation, and also add fix to have the original allele's INFO field be passed through for batch merging
Menachem Fromer
2012-01-26 18:25:30 -0500
246e085ec9Unit tests for GATKSAMRecord class
Mauricio Carneiro
2012-01-26 16:59:08 -0500
0d4027104fReduced reads are now aware of their original alignments
Mauricio Carneiro
2012-01-26 16:07:29 -0500
07f72516aeUnsupported platform should be a user error
Eric Banks
2012-01-26 16:14:25 -0500
cdff23269dHaplotypeCaller now uses insertions and softclipped bases as possible triggers. LocusIteratorByState tags pileup elements with the required info to make this calculation efficient. The days of the extended event pileup are coming to a close.
Ryan Poplin
2012-01-26 15:56:33 -0500
44c3a41f67Merge branch 'master' of ssh://nickel.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable
Ryan Poplin
2012-01-26 13:25:33 -0500
dbe9eb70feUpdating HC integration tests after merge
Ryan Poplin
2012-01-26 13:25:22 -0500
673ceadd11While this fix worked for the evaluator module, it could potentially have bad effects in the phasing walkers. Special-case nocalls in the PhasingEvaluator and return AllelePair to previous state.
Christopher Hartl
2012-01-26 13:06:36 -0500
9c6fda7e15Yup. I was right.
Christopher Hartl
2012-01-26 12:54:11 -0500
bd27f3bc9bMerge branch 'master' of ssh://nickel.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable
Guillermo del Angel
2012-01-26 12:52:16 -0500
67c89cadadFixes for pool caller to match UG outputs at certain sites: implement min base qual/min mapping qual read filter so those reads are filtered from pileups, and implemented filter for sites that have a too large a fraction of deletions
Guillermo del Angel
2012-01-26 12:52:00 -0500
7d059540a4Allow segments of genome to be excluded in generating a reference panel. Occasionally targets would contain no variation (typically, in the middle of the centromere), which beagle doesn't particularly like, and errors out rather than producing empty output files. The best way to deal with these is to just exclude the regions on a second-pass, and the remaining bits will be gathered with no additional work.
Christopher Hartl
2012-01-26 12:43:52 -0500
9d4b84f6bdMerge branch 'master' of ssh://chartl@ni.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable
Christopher Hartl
2012-01-26 12:38:24 -0500
25532bdc37Merge branch 'master' of ssh://nickel.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable
Ryan Poplin
2012-01-26 11:43:32 -0500
390d493049Updating ActiveRegionWalker interface to output a probability of active status instead of a boolean. Integrator runs a band-pass filter over this probability to produce actual active regions. First version of HaplotypeCaller which decides for itself where to trigger and assembles those regions.
Ryan Poplin
2012-01-26 11:37:08 -0500
9a63a9ae3cMerge branch 'master' of ssh://nickel.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable
Eric Banks
2012-01-26 00:38:25 -0500
859dd882c9Don't make it standard for now
Eric Banks
2012-01-26 00:38:16 -0500
c5e81be978Adding pairwise AF table. Not polished at all, but usable none-the-less.
Eric Banks
2012-01-26 00:37:06 -0500
774e540042Fixing broken test
Eric Banks
2012-01-26 00:31:41 -0500
5b9c8ab01bAnother quick update missed in the merge
Eric Banks
2012-01-25 21:53:20 -0500
702a2d768fInitial version of multi-allelic summary module in VariantEval
Eric Banks
2012-01-25 19:42:55 -0500
9a60887567Lost an import in the merge
Eric Banks
2012-01-25 19:41:41 -0500
cba5f1a8b1Merge branch 'master' of ssh://nickel.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable
Eric Banks
2012-01-25 19:19:03 -0500
ddaf51a50fUpdated one integration test for indels
Eric Banks
2012-01-25 19:18:51 -0500
add6918f32Cleaner, more efficient way of determining the last dependent set in the queue.
Eric Banks
2012-01-25 16:21:10 -0500
db645a94caAdded options to make the batch-merger more all-inclusive: keep all indels, SNPs (even filtered ones) but maintain their annotations. Also, VariantContextUtils.simpleMerge can now merge variants of all types using the Hidden non-default enum MultipleAllelesMergeType=MIX_TYPES
Menachem Fromer
2012-01-25 16:10:59 -0500
d405ec2a0dMerge branch 'master' of ssh://nickel.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable
Guillermo del Angel
2012-01-25 15:53:32 -0500
4337dcd7e4More pool caller bug fixes: the QUAL field was actually multiplied by 10 (accounting for a lot of singletons that shouldn't have been there), and correct AD output
Guillermo del Angel
2012-01-25 15:53:03 -0500
66772d0ebfNext iteration in the pool caller: more bug fixes, start of big refactoring to clean up interfaces, moved up a lot of attributes that really belonged to a site up from the Pool class, added by default option to filter out a call if there's no reference depth (instead of just skipping the call which makes it hard to figure out what happened afterwards).
Guillermo del Angel
2012-01-25 15:41:08 -0500
ef335a5812Better implementation of the fix; PL index is now traversed in order.
Eric Banks
2012-01-25 15:15:42 -0500
8e2d372ab0Use remove instead of setting the value to null
Eric Banks
2012-01-25 14:41:34 -0500
05816955aaIt was possible that we'd clean up a matrix column too early when a dependent column aborted early (with not enough probability mass) because we weren't being smart about the order in which we created dependencies. Fixed.
Eric Banks
2012-01-25 14:28:21 -0500
2799a1b686Catch exception for bad type and throw as a TribbleException
Eric Banks
2012-01-25 12:15:51 -0500
96b62daff3Minor tweak to the warning message.
Eric Banks
2012-01-25 11:55:33 -0500
fb863dc6a7Warn user when trying to run with EMIT_ALL_SITES with indels; better docs for that option.
Eric Banks
2012-01-25 11:50:12 -0500
e349b4b14bAllow appending with the dbSNP ID even if a (different) ID is already present for the variant rod.
Eric Banks
2012-01-25 11:35:54 -0500
ea3d4d60f2This annotation requires rods and should be annotated as such
Eric Banks
2012-01-25 11:35:13 -0500
7a26fcb86fSetting the max alternate alleles for the exact model in the HaplotypeCaller's copy of the UG engine.
Ryan Poplin
2012-01-25 09:51:13 -0500
bbefe4a272Added option to be able to write out the active regions to an interval list file
Ryan Poplin
2012-01-25 09:47:06 -0500
9818c69df6Can now specify active regions to process at the command line, mainly for debugging purposes
Ryan Poplin
2012-01-25 09:32:52 -0500
8f7d9bff0aFix MultiplyLikelihood logic, minor alterations to scala scripts.
Christopher Hartl
2012-01-25 00:12:19 -0500
22f0caccacMerge branch 'master' of ssh://nickel.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable
Guillermo del Angel
2012-01-24 16:22:39 -0500
97499529c7another small bug with the file extension.
Mauricio Carneiro
2012-01-24 16:13:53 -0500
ffd61f4c1cRefactor the Pileup Element with regards to indels
Mauricio Carneiro
2012-01-17 18:56:50 -0500
4aacaf8916Merge branch 'master' of ssh://gsa1/humgen/gsa-scr1/gsa-engineering/git/unstable
Matt Hanna
2012-01-24 15:32:26 -0500
c312bd5960Weirdly, PicardException inherits from SAMException, which means that our specialty code for reporting malformed BAMs was actually misreporting any error that happened in the Picard layer as a BAM ERROR.
Matt Hanna
2012-01-24 15:30:04 -0500
47f4440aeaMerged bug fix from Stable into Unstable
David Roazen
2012-01-24 15:19:29 -0500
b07fdb1089Rename alltests* targets in build.xml
David Roazen
2012-01-24 14:53:31 -0500
8aaca83e91Next intermediate iteration on pool caller: start using UnifiedArgumentCollection in order to have similar syntax to UG, several numerical bug fixes, add more logging to vcf (not done yet)
Guillermo del Angel
2012-01-24 11:45:19 -0500
7c7ca0d799fixing bug with fastq extension
Mauricio Carneiro
2012-01-24 10:59:25 -0500
0a3172a9f1Fix for ref 0 bases for Chris
Mark DePristo
2012-01-24 10:53:37 -0500
945cf03889IntelliJ ate my import!
Mauricio Carneiro
2012-01-23 21:46:45 -0500
2bb9525e7fDon't set base qualities if fastQ is provided
Mauricio Carneiro
2012-01-23 17:56:27 -0500
b6c816fe12Turn off unnecessary printing in analyzeRunReports
Mark DePristo
2012-01-23 17:04:12 -0500
1172517abbBugfix for version parsing.
Mark DePristo
2012-01-23 16:02:53 -0500
ceca7e0b37Bugfix to now separate completed, sting and user exceptions. Added dry run mode
Mark DePristo
2012-01-23 12:12:38 -0500
1f620c79e6Add busers and bugroup information to queueStatus
Mark DePristo
2012-01-23 12:12:00 -0500
10bc26079dbugfix to actually run correct python script
Mark DePristo
2012-01-23 12:11:38 -0500
4b17fc3cc1Parallel implementation of random forest training. Very cool (and easy) example of parallel processing in R
Mark DePristo
2012-01-23 12:06:46 -0500
bb203ccf0acombined analyses of snps and indels.
Mark DePristo
2012-01-23 10:51:07 -0500
0ec6f86c21Tests for event length, combined snps and indels. Partial infrastructure to train and eval trees.
Mark DePristo
2012-01-21 14:08:25 -0500
c18beadbdbDevice files like /dev/null are now tracked as special by Queue and are not used to generate .out file paths, scattered into a temporary directory, gathered, deleted, etc. Attempted workaround for xdr_resourceInfoReq unsatisfied link during loading of libbat.so.
Khalid Shakir
2012-01-23 16:17:04 -0500
cc4ba7372fWhy is reference_bases even an option anymore?
Christopher Hartl
2012-01-23 15:18:59 -0500
3392d67c1aMaybe a switch to reference bases will fix this
Christopher Hartl
2012-01-23 15:10:03 -0500
15c0c294c1Adding in this walker to try to debug the 0-byte ref bases
Christopher Hartl
2012-01-23 14:51:24 -0500
02450e4b12Merged bug fix from Stable into Unstable
Mark DePristo
2012-01-23 12:08:39 -0500
798596257bEnable the Genotype Phasing Evaluator. Because it didn't have the same argument structure as the base class, update2 of VariantEvaluator was being called, rather than update2 of the actual module.
Christopher Hartl
2012-01-23 10:50:16 -0500
80a4ce0edfBugfix for incorrect error messages for missing BAMs and VCFs
Mark DePristo
2012-01-23 09:52:07 -0500
31d2f04368Merge branch 'master' of ssh://nickel.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable
Guillermo del Angel
2012-01-23 09:23:03 -0500
966387ca0bNext intermediate commit in the pool caller. Lots of bug fixes and now we can emit true vcf's with calls in discovery mode (still of unknown quality) - old validation mode is temporarily broken,will be fixed in next refactoring.
Guillermo del Angel
2012-01-23 09:22:31 -0500
4a08e8ca6eMinor tweaks to T2D-related qscripts. Replacing old md5s from the BeagleIntegrationTest. All differences boiled down either to the accounting of genotypes changed (./. --> 0/0 is no longer a "changed" genotype, and original genotypes that were ./. are represented as OG=. rather than OG=./. .)
Christopher Hartl
2012-01-23 08:25:34 -0500
4d6312d4eaHaplotypeCaller is now an ActiveRegionWalker.
Ryan Poplin
2012-01-22 14:31:01 -0500
3b1aad4f17After a minor and abject freakout, alter the T2D script to seek out truth sensitivities between 80 and 100, rather than between 0.8 and 1. Also, don't consider a genotype "changed by beagle" if the initial genotype is a no-call.
Christopher Hartl
2012-01-20 23:43:51 -0500
9b4f6afa21Alterations to scripts for better performance. Grid search now expands the sens/spec tradeoff (90 was far too aggressive against hapmap chr20), and 20 max gaussians was too many, and caused errors. For consensus genotypes: remember to gunzip the beagle outputs before converting to VCF. Also, beagle can in fact create 'null' alleles in certain circumstances. I'm not sure what exactly those circumstances are, but those sites should be ignored. When it does, all alleles apear to be set to null, so this should not affect the actual phasing in the output VCF.
Christopher Hartl
2012-01-20 23:07:59 -0500
f3564bbf43Ugh. Darn intelliJ not telling me I was missing an import statement.
Christopher Hartl
2012-01-20 13:25:11 -0500
b902d778ca.
Christopher Hartl
2012-01-20 13:22:46 -0500
7c6a9471e8After ensuring MultiplyLikelihoods does what I want it to do, add a quick and simple integration test to ensure I don't break it.
Christopher Hartl
2012-01-20 13:20:13 -0500