0398ae9695I hate these disabled unit tests, #2
Eric Banks
2012-06-14 15:19:27 -0400
676a57de7bI hate these disabled unit tests
Eric Banks
2012-06-14 14:03:58 -0400
065fe9a7b9Updating md5 for bug fixes
Eric Banks
2012-06-14 13:03:51 -0400
de5508fceaBug fixes for cycle and context covariates
Eric Banks
2012-06-14 13:01:14 -0400
5c3c6cbc40Long -> long conversions in BQSR
Eric Banks
2012-06-14 09:07:02 -0400
29a74908bbThe next round of BQSR optimizations: no more Long[] array creation
Eric Banks
2012-06-14 00:05:42 -0400
cd2074b1dcMerge branch 'master' of ssh://gsa4.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable
Guillermo del Angel
2012-06-13 20:59:30 -0400
92669a0468Second intermediate commit for indel pool caller - now works (more or less) in reference sample-free mode. Still needs a lot of cleanups/add more tests and not done w/refactoring quite yet
Guillermo del Angel
2012-06-13 20:59:17 -0400
0550b27799Make downsampler classes themselves generic (instead of just the Downsampler interface)
David Roazen
2012-06-13 14:42:10 -0400
fc36177004Accidentally removed the scattering from this test
Joel Thibault
2012-06-13 14:17:31 -0400
67c0569f9cMerge branch 'master' of ssh://gsa4.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable
Guillermo del Angel
2012-06-13 11:50:00 -0400
9a4c39357aFixed bugs in previous commit
Guillermo del Angel
2012-06-13 11:49:20 -0400
81993b08e2Don't put null entries into the key array
Eric Banks
2012-06-13 11:43:44 -0400
bdf5945dccFixed bugs in DiagnoseTargets
Roger Zurawicki
2012-06-12 16:49:52 -0400
538cdf9210Created the FindCoveredIntervals
Roger Zurawicki
2012-06-11 16:19:57 -0400
aee66ab157Big UG refactoring and intermediate commit to support indels in pool caller (not done yet). Lots of code pulled out of long spaghetti-like functions and modularized to be easily shareable. Add functionality in ErrorModel to count indel matches/mismatches (but left part disabled as not to change integration tests in this commit), add computation of pool genotype likelihoods for indels (not fully working yet in more realistic cases, only working in artificial nice pools). Lot's of TBD's still but existing UG and pool SNP functionality should be intact
Guillermo del Angel
2012-06-13 11:14:44 -0400
bb77aa88c3Drat, forgot the unit tests again
Eric Banks
2012-06-12 19:00:47 -0400
bacd25d1a4Merge branch 'master' of ssh://gsa2.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable
Eric Banks
2012-06-12 16:15:36 -0400
37f56ce8fdA couple of minor updates to BQSR
Eric Banks
2012-06-12 16:12:13 -0400
277493dd83Yet more instances of Lists changed over to native arrays
Eric Banks
2012-06-12 15:56:09 -0400
613badc835Very minor optimizations for the context covariate
Eric Banks
2012-06-12 15:47:32 -0400
c8cda209b5Make intersection easy to change
Joel Thibault
2012-06-11 13:39:14 -0400
0f79adb2aaChanging more Java Lists to native arrays in BQSR for performance optimization.
Eric Banks
2012-06-12 15:41:01 -0400
1da3e43679Wow, apparently it's way, way less efficient to iterate over Java Lists than native arrays. With this change and the bit fiddling, Ryan's 10-day test case now runs in 1 day. More to come.
Eric Banks
2012-06-12 13:32:56 -0400
a96c5da884Oops, forgot to push the unit tests
Eric Banks
2012-06-12 11:38:30 -0400
a057cf31b3Merge branch 'master' of ssh://gsa2.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable
Eric Banks
2012-06-12 11:02:50 -0400
fec0bd5e11Fixing UG argument docs
Eric Banks
2012-06-12 09:46:16 -0400
a4defdfb29Adding a GT header line to SomaticIndelDetector output
Eric Banks
2012-06-12 09:39:17 -0400
891ce51908Refactoring of BQSRv2 to use longs (and standard bit fiddling techniques) instead of Java BitSets for performance improvements.
Eric Banks
2012-06-12 09:19:36 -0400
0c98c34f5fFixing bugs caught by Eric
Mauricio Carneiro
2012-06-11 17:32:22 -0400
e0fbffac3aFixing bugs caught by Eric
Mauricio Carneiro
2012-06-11 16:05:15 -0400
ff5749599dMerge branch 'master' of ssh://gsa2.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable
Eric Banks
2012-06-11 15:46:17 -0400
fea625632fDon't use asList because it maintains an iterator to the original list and then the result can't be used to create a new one
Eric Banks
2012-06-11 15:45:58 -0400
41ef22fe7cRe-enabling BQSR integration tests
Eric Banks
2012-06-11 15:44:55 -0400
e4d371dc80Merge branch 'master' of ssh://gsa2.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable
Ryan Poplin
2012-06-11 10:38:50 -0400
683d4b508eBug fix in fragment utils: the read name wasn't being set in the merged read. Misc minor updates to the HaplotypeCaller.
Ryan Poplin
2012-06-11 10:38:35 -0400
4aad7e23efNew ReduceReads v2 with unclipped variant regions and soft-clipped bases
Mauricio Carneiro
2012-05-14 13:39:42 -0400
edf08ca928Updates to ReduceReads
Mauricio Carneiro
2012-05-03 12:23:49 -0400
afa9b2718aMerge branch 'master' of ssh://gsa2.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable
Eric Banks
2012-06-08 13:54:48 -0400
92280b4068BQSR optimization: cache the BitSetUtils.bitSetFrom() calls since they are called over and over again with the same values. Another 10% reduction in runtime.
Eric Banks
2012-06-08 13:54:37 -0400
898a0e6161Minor optimizations
Eric Banks
2012-06-08 12:07:58 -0400
0a37e19998Bug fix in VQSR so that the VCF index will be created for the recalFile.
Ryan Poplin
2012-06-08 11:51:28 -0400
d463ab2cbfBQSR optimization: String manipulation is extremely expensive in Java (accounts for 8% of BQSR runtime). Instead use byte[] and StringBuilder when possible.
Eric Banks
2012-06-08 10:42:42 -0400
2bd48a7351Bad comments made it into the previous commit
Eric Banks
2012-06-07 23:12:56 -0400
31c3a6be48BQSR optimization: getRequiredCovariates() and getOptionalCovariates() were creating a new List every time they were being called, and unfortunately getRequiredCovariates().size() is used as the stop condition in for-loops throughout the code. Just maintaining the original list of covariates results in a 15% reduction in runtime for BQSR.
Eric Banks
2012-06-07 20:04:10 -0400
0fb9179f76BQSR optimization: don't clone the original quals for each read, we can just overwrite the original array
Eric Banks
2012-06-07 19:41:03 -0400
a13d125ba1Split out contig names from Reference .fai file on white space (to support the GATK resource bundle's file human_g1k_v37.fasta.fai.gz, which does not use tab delimiters)
Menachem Fromer
2012-06-07 16:56:32 -0400
32e4fe5c87Enable zero-sample testing
Joel Thibault
2012-06-06 17:41:29 -0400
fe9dc0cc3fExplicitly set intervals in new test
Joel Thibault
2012-06-06 17:36:43 -0400
0ac4ba9ad3unintended change crept in
Ryan Poplin
2012-06-07 10:59:04 -0400
d449f169d3Merge branch 'master' of ssh://gsa2.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable
Ryan Poplin
2012-06-07 10:56:55 -0400
0b4281fdd0misc minor update to HC debug output for when there are a lot of samples
Ryan Poplin
2012-06-07 10:56:41 -0400
bad50a1b05Fix docs
Eric Banks
2012-06-06 22:45:38 -0400
b093ba9dccStabilized NGSPlatform code: don't assume all reads have read groups (e.g. artificial SAM records)
Eric Banks
2012-06-06 15:17:30 -0400
54f682a99cUnify to NGSPlatform framework. TechnologyComposition annotation now generalizes to Illumina and not just SLX.
Eric Banks
2012-06-06 11:44:37 -0400
dd46d843fbIR should skip Ion reads just like it does with 454 reads; Tim has confirmed that official platform name for Ion.
Eric Banks
2012-06-06 11:04:55 -0400
b6a7c3f780Merged bug fix from Stable into Unstable
David Roazen
2012-06-05 17:44:44 -0400
3b2fab9a37Update stable for shptmp -> hptmp migration
David Roazen
2012-06-05 17:38:06 -0400
2cbd6e5f90Merged bug fix from Stable into Unstable
Guillermo del Angel
2012-06-05 15:58:23 -0400
ce4dc2128dAdding minor clarification to -mbq argument documentation
Guillermo del Angel
2012-06-05 15:17:56 -0400
e02ec8c8b6Don't update the record ID unless we are actually going to emit the record
Eric Banks
2012-06-04 14:58:50 -0400
8405156ae1Refactored VariantsToTable so that 1) genotype-level fields can be specified (stabilized and supported code) and 2) the --moltenize argument could be supported to produce molten output of the data. Added tests that cover these capabilities.
Eric Banks
2012-06-04 14:28:32 -0400
f11e7ebc3aFixing the previous fix related to clipping. Adding extra reference padding in the HaplotypeCaller to get those larger alleles during GGA.
Ryan Poplin
2012-06-04 12:49:36 -0400
320956ee4bBug fix in clipping function in ReadUtils for when the read ends at exactly the clipping boundary. Bug fixes in HaplotypeCaller GGA mode for when Smith-Waterman produces a different allele than what was given in the input alleles VCF. GGA mode now works with multiallelic records. Adding min pruning factor argument which is combined with the pruning factor that is determined dynamically by the coverage.
Ryan Poplin
2012-06-04 10:55:36 -0400
7a54baf08cMerged bug fix from Stable into Unstable
Guillermo del Angel
2012-06-03 08:42:08 -0400
47df7bbc14Merge branch 'master' of ssh://gsa4.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/stable
Guillermo del Angel
2012-06-03 08:38:54 -0400
2ddbdee3bcFixed broken VariantEval stratifications VariantType and IndelSize - integration tests to follow
Guillermo del Angel
2012-06-03 08:38:38 -0400
12a8c54f9aFixing VCF header for filter elements (thanks Eric)
Mauricio Carneiro
2012-06-01 15:45:15 -0400
3b20b4615dant fasttest: for when you just can't wait
David Roazen
2012-06-01 15:03:43 -0400
2e3af1a739Enable contracts for the release jar tests (whoops...)
David Roazen
2012-06-01 13:08:46 -0400
fdf607c3a9One-off walker but techdev may find it useful and expand on it: given a BAM, output series of intervals to file (as interval_list [default] or BED [if -bed used]) where there is coverage. Additionally, print out to stdout at the end total covered bases and total covered bases including soft-clips
Guillermo del Angel
2012-05-31 20:47:37 -0400
3a15ba2102Malformed VCF headers should be User Errors
Eric Banks
2012-05-31 16:05:53 -0400
d1822c926cTemporary fake release directories for Bamboo testing purposes
David Roazen
2012-05-31 16:02:57 -0400
1aa9742e55Decrease block size from 1Mb to 100Kb
Joel Thibault
2012-05-30 15:32:32 -0400
c4f7df4dceWhen an underlying exception occurs because of the user error, if the exception instance does not include a message instead of telling the user "because null", tell them "because <exception class name>".
Khalid Shakir
2012-05-30 16:39:06 -0400
421d0d1435Merge branch 'master' of ssh://gsa2.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable
Ryan Poplin
2012-05-30 15:21:35 -0400