Commit Graph

  • 18728ec5bd Updates to the bundle script: 1. Add the symbolic 'current' link for the new bundle dir 2. Don't gzip and copy .out files 3. Don't call chr20 SNPs on the example BAM because it's now just a few reads on chr1 Eric Banks 2012-12-18 11:16:42 -0500
  • 16eb1c5436 Optimization to TraverseReadsNano Mark DePristo 2012-12-17 16:40:04 -0500
  • b33f804cdc Inline increment function in RecalDatum to avoid minor duplication of work and multiple synchronized method calls Mark DePristo 2012-12-17 16:38:10 -0500
  • 66d32f646b Minor cleanup of BAQ calculation (final variables, etc) Mark DePristo 2012-12-17 16:37:20 -0500
  • 67fe81391c ProgressMeter optimization: don't do genome loc formatting, but instead create an object that only formats when printing is actually needed Mark DePristo 2012-12-17 15:55:21 -0500
  • 1de2f527b9 Optimization of recalibrateRead Mark DePristo 2012-12-13 17:21:36 -0500
  • bd6cda7542 Trivial optimization of TraverseReadsNano -- don't format the shard toString if logger isn't debug enabled Mark DePristo 2012-12-12 16:21:43 -0500
  • a481d006f0 Optimizations for applying BQSR table with PrintReads Mark DePristo 2012-12-12 16:14:09 -0500
  • 5ec25797b3 Optimizations for BaseRecalibrator Mark DePristo 2012-12-12 13:54:14 -0500
  • 6cb1a0a698 Avoid db naming collisions when running multiple instances of na12878kb tests concurrently David Roazen 2012-12-17 16:13:03 -0500
  • e6f468b647 Refactored the quasi-useful IndelType annotation into the more useful VariantType. The indels are still annotated as before, but now all other variant types are annotated too. I'm doing this because of requests on the forum but am not making it standard. If we find it to be useful we can turn it on by default later. Eric Banks 2012-12-17 11:54:47 -0500
  • 6171419e6c Merge branch 'master' of github.com:broadinstitute/gsa-unstable Ryan Poplin 2012-12-17 11:27:52 -0500
  • 98f18b5f9e Changing the HC over to using the non-contamination-downsampled read maps for the purposes of annotations. This behavior now matches the UG. There is a new command line option to go back to the older behavior to explore the differences. Ryan Poplin 2012-12-17 11:27:44 -0500
  • 40b3941678 Merge pull request #7 from jsilter/master depristo 2012-12-17 08:05:55 -0800
  • 762f184262 Bug fix for strict validation: rsID checking wasn't working if there were multiple IDs Eric Banks 2012-12-17 10:32:41 -0500
  • 1d18ee26cc Merge remote-tracking branch 'unstable/master' Eric Banks 2012-12-17 09:03:03 -0500
  • ea704d688f chose smaller buffer size for the bufferedStream Yossi Farjoun 2012-12-15 13:01:38 -0500
  • 6da2338ea7 removed comments and uneeded imports Yossi Farjoun 2012-12-15 12:31:37 -0500
  • 19dd2d628a some changes. Yossi Farjoun 2012-12-14 16:58:14 -0500
  • 8198414352 Add removeCall, which removes all calls matching the provided MongoVariantContext. _id and Date fields are excluded. Generally intended te remove only a single call Jacob Silterra 2012-12-14 13:07:44 -0500
  • 5f1afb4136 Fixing an off-by-one clipping error in ReduceReads for reads off the contig Mauricio Carneiro 2012-12-13 21:56:44 -0500
  • 74344a3871 Bringing in the changes from the CMI repo Mauricio Carneiro 2012-12-13 21:59:37 -0500
  • 696bf95fba Fix for PBT bug reported on the forum: the AD is actually output correctly now (rather than with 'null' or some gibberish memory pointer). Eric Banks 2012-12-13 23:28:30 +0000
  • aeab932c63 Actual working version of unflushing VCFWriter Mark DePristo 2012-12-13 12:24:08 -0500
  • 97e880654b Bug fix in setting output file names Guillermo del Angel 2012-12-13 15:08:22 -0500
  • cf44ec1ace Experiment: don't compress picard BAM outputs Guillermo del Angel 2012-12-12 17:45:55 -0500
  • 021929788c Don't compress intermediate BAMs, will save time Guillermo del Angel 2012-12-12 17:40:25 -0500
  • 5e66109268 Replaced a useless getInt with a skipInt to remove 1/4 of the initial seek time in the BAM Index. Yossi Farjoun 2012-12-12 15:26:38 -0500
  • 5780926254 Merge branch 'master' of gsa2:/humgen/gsa-scr1/chartl/dev/unstable Chris Hartl 2012-12-12 15:42:40 -0500
  • 62eaffdf0a Fix docs for ReadBackedPhasing Eric Banks 2012-12-12 20:28:04 +0000
  • bba63a3b0e Fix for GSA-615: UnifiedGenotyperEngine.getGLModelsToUse takes 5% of the runtime of UG, should be optimized away. Eric Banks 2012-12-12 20:25:45 +0000
  • 211a6e78ea Further related bug fixes to GGA mode in the HC: some variants (especially MNPs) were causing problems because they don't have to start at the current location to match the allele being genotyped. Fixed. Ryan Poplin 2012-12-12 14:53:02 -0500
  • 33290bfe0c Added integration test to catch the read off contig in ReduceReads. Mauricio Carneiro 2012-12-12 13:41:01 -0500
  • a52e3c7e15 Revert "Bug fix for RR: don't let the softclip start position be less than 1" this introduced a bug in reduce reads by de-activating it's hard clipping of the out of bounds soft-clips (specially in the MT). DEV-322 #resolve #time 4m Mauricio Carneiro 2012-12-12 13:09:39 -0500
  • 5632c13bf2 Resolves GSA-681 / Compressed VCF.gz output is too big because of unnecessary call to flush(). Mark DePristo 2012-12-12 10:27:07 -0500
  • 216f92276c Disable scatter-gather with PrintReads since we're already setting high nct so it's unnecessary Guillermo del Angel 2012-12-12 09:06:37 -0500
  • 0e5b1093fb initial implementation of contamination estimation, tested on single gene (which doesn't have enough data) waiting to test on exome/chr20 Kristian Cibulskis 2012-12-11 15:59:59 -0500
  • dd52a70d45 Fix AFCalcResult unit test Mark DePristo 2012-12-11 10:39:58 -0500
  • 6bf31065e3 Merge branch 'master' of github.com:broadinstitute/gsa-unstable Ami Levy-Moonshine 2012-12-11 10:34:50 -0500
  • 2f99569dda change the md5 in one of the CV intergration tests, since it wasn't use the priority list when printing the origin of the annotation (the setValue field) Ami Levy-Moonshine 2012-12-10 22:48:15 -0500
  • 2e3284f306 Continue to fix the case where PRIORITIZE is used but no priority list is given. While fixing that case I also removed unnecessary sorting, when the prioeity list is not provied. When the priority list is not provided, it will continue to be null. Thus, the number of original Variant Contexts should be given as a new parameter to simpleMerge (since priority might be null). This new parameter is used for checking if there are filtered VC, when annotationOrigin is true. Ami Levy-Moonshine 2012-12-10 22:23:58 -0500
  • 19372225af Merge Broad's GATK and CMI gatk Mauricio Carneiro 2012-12-10 15:33:38 -0500
  • 8a115edbaf ReduceReads is now scattered by contig Mauricio Carneiro 2012-12-06 17:06:32 -0500
  • bdda63d973 Related bug fixes to GGA mode in the HC: some variants (especially MNPs) were causing problems because they don't have to start at the current location to match the allele being genotyped. Fixed. Eric Banks 2012-12-10 14:47:04 -0500
  • ceb5431dcb Merge branch 'master' of github.com:broadinstitute/gsa-unstable Ryan Poplin 2012-12-10 12:24:47 -0500
  • c84ff9d75e Adding explicit true negative assessment category to the AssessNA12878 walker. Ryan Poplin 2012-12-10 12:24:43 -0500
  • 573ace4403 restore the right version of VariantContextUtils.java in my unstable dir Ami Levy-Moonshine 2012-12-10 10:28:56 -0500
  • 46edab6d6a Use the new downsampling implementation by default David Roazen 2012-11-26 12:44:48 -0500
  • 5460c96137 Merge branch 'master' of github.com:broadinstitute/gsa-unstable Ami Levy-Moonshine 2012-12-09 23:43:57 -0500
  • 3a420d163e (1) changes in catVariants (work still under development) (2) changes to CV to throw an error when GenotypeMergeType is PRIORITIZE but no priority (rod_priority_list) is not given. Reported by TechnicalVault on the forum on Nov 14 2012 Ami Levy-Moonshine 2012-12-09 23:40:03 -0500
  • 2637f512f8 Merge branch 'master' of github.com:broadinstitute/gsa-unstable Eric Banks 2012-12-09 02:09:40 -0500
  • 574d5b467f Bug fix for indel HMM: protect against situation where long reads (e.g. Sanger) in a pileup can lead to a read starting after the haplotype end for a given haplotype. Eric Banks 2012-12-09 02:09:34 -0500
  • 9b6ee0576f Fix bugs in the consensus genotype creation algorithm for the NA12878 KB Mark DePristo 2012-12-08 13:18:07 -0500
  • bf8421eeb7 Fixes GSA-671 / AFCalcResult.log10pNonRefByAllele should really be log10pRefByAllele Mark DePristo 2012-12-07 16:02:58 -0500
  • 3355216366 Merge branch 'master' of github.com:broadinstitute/gsa-unstable Ryan Poplin 2012-12-07 15:35:17 -0500
  • 9573648f85 Changes to count the sites which might be present in some of the input rods but not present at all in other rods. Now loop over the input rod names instead of looping over the tracker results. Ryan Poplin 2012-12-07 15:35:08 -0500
  • 3b0e3767bf Add a test for a read that extends off the end of chr1 Joel Thibault 2012-12-07 13:46:41 -0500
  • 58e39a8468 Enabling 4-way parallel by default in FastQ2BAM DEV-317 Mauricio Carneiro 2012-12-06 17:24:09 -0500
  • cc4e3ec589 Update TODO list Joel Thibault 2012-12-05 12:45:02 -0500
  • abd94b2976 Bugfix for handling invalid records in NA12878 KB Mark DePristo 2012-12-06 10:24:24 -0500
  • 406adb8d44 The allele biased downsampling should not abort if there's a reduced read. Rather it should always keep the RR and downsample only original reads in the pileup. Eric Banks 2012-12-05 23:15:36 -0500
  • 6d22f4f737 Bringing latest performance updates from the GATK to CMI Mauricio Carneiro 2012-12-05 21:40:03 -0500
  • dbf721968d PrintReads large-scale test to protect against another major low-level performance issue Mark DePristo 2012-12-05 21:35:27 -0500
  • 00c23bf704 Merge branch 'master' of github.com:broadinstitute/gsa-unstable Ryan Poplin 2012-12-05 15:53:05 -0500
  • 234ff64556 Changes to AssessNA12878 to allow for 100s of input callsets to assess against the database. Ryan Poplin 2012-12-05 15:52:57 -0500
  • 5d78a61f7a Merge branch 'master' of github.com:broadinstitute/gsa-unstable Ami Levy-Moonshine 2012-12-05 15:07:12 -0500
  • d0cab795b7 Got caught in the middle of a bad integration test, that was fixed in independent push. Moved test bam into testdata. Mark DePristo 2012-12-05 14:49:01 -0500
  • 465694078e Major performance improvement to the GATK engine Mark DePristo 2012-12-04 22:08:01 -0500
  • 2b601571e7 Better error handling in NanoScheduler Mark DePristo 2012-12-04 20:59:07 -0500
  • 51dbb562c9 Reduce amount of debugging information from NA12878KnowledgeBaseServer Mark DePristo 2012-12-02 16:46:07 -0500
  • efe256ec09 binary search implementation to find the minimum coverage Mauricio Carneiro 2012-12-05 12:19:07 -0500
  • 430d6a07f2 Merge branch 'master' of gsa2:/humgen/gsa-scr1/chartl/dev/unstable Chris Hartl 2012-12-05 11:20:28 -0500
  • 0c925856cb Merge branch 'master' of github.com:broadinstitute/gsa-unstable Eric Banks 2012-12-05 02:00:39 -0500
  • ef87b18e09 In retrospect, it wasn't a good idea to have FisherStrand handle reduced reads since they are always on the forward strand. For now, FS ignores reduced reads but I've added a note (and JIRA) to make this work once the RR het compression is enabled (since we will have directionality in reads then). Eric Banks 2012-12-05 02:00:35 -0500
  • 13896356ad Added bootstrapping and fixed the GLM model of the FMCC Mauricio Carneiro 2012-12-05 01:06:19 -0500
  • 30f013aeb0 Added a copy() method for ReadBackedPileups Mauricio Carneiro 2012-12-04 23:56:30 -0500
  • 6feda540a4 Better error message for SimpleGATKReports Mauricio Carneiro 2012-12-04 23:55:35 -0500
  • 726332db79 Disabling the testNoCmdLineHeaderStdout test in UG because it keeps crashing when I run it locally Eric Banks 2012-12-05 00:54:00 -0500
  • 61bde6210b Restored RemoteFile push and pull in base QScript. kshakir 2012-12-04 12:32:00 -0500
  • 8d2d0253a2 introduce a level of indirection for the forum URLs - this new function will allow me a place to morph the URL into something that is supported by Confluence Randal Moore 2012-12-03 12:54:48 -0600
  • 1af41754e3 Merge branch 'master' of github.com:broadinstitute/gsa-unstable Eric Banks 2012-12-03 22:01:11 -0500
  • bca860723a Updating tests to handle bad validation data files (that used the wrong qual score encoding); overrides push from stable. Eric Banks 2012-12-03 22:01:07 -0500
  • 387c0defed don't change md5 here because I am handling it separately from unstable with a better command-line in the test Eric Banks 2012-12-03 21:49:45 -0500
  • ef95757311 Fix MD5 because of a need to fix a busted bam file in our validation directory (it used the wrong quality score encoding...) Eric Banks 2012-12-03 21:46:46 -0500
  • 4ced2e4ffc Merge branch 'develop' of github.com:broadinstitute/cmi-gatk into develop Guillermo del Angel 2012-12-03 20:14:43 -0500
  • c2c6b858e3 Better checks/more flexibility in fastq2bam parsing. Immediate benefit: we can now process normal-only samples, and metadata should be able to specify tumor/normal pairs in any order. Hard-coded hacks removed. DEV-134 #resolve #time 3m Guillermo del Angel 2012-12-03 20:14:37 -0500
  • 472381245a Allow for more refined control of memory and queues to run with Menachem Fromer 2012-12-03 17:07:03 -0500
  • 67932b357d Bug fix for RR: don't let the softclip start position be less than 1 Eric Banks 2012-12-03 15:59:14 -0500
  • d5ed184691 Updating the HC integration test md5s. According to the NA12878 knowledge base this commit cuts down the FP rate by more than 50 percent with no loss in sensitivity. Ryan Poplin 2012-12-03 15:38:59 -0500
  • a47da9bb2f Merge branch 'master' of github.com:broadinstitute/gsa-unstable Ryan Poplin 2012-12-03 14:30:14 -0500
  • 156d6a5e0b misc minor bug fixes to GenotypingEngine. Ryan Poplin 2012-12-03 12:47:35 -0500
  • 5fed9df295 Quick fix: base qual array in the GATKSAMRecord stores the actual phred values (-33) and not the original bytes (duh). Eric Banks 2012-12-03 12:18:20 -0500
  • b6839b3049 Added checking in the GATK for mis-encoded quality scores. Eric Banks 2012-12-03 11:18:41 -0500
  • 18b002c99c Merge branch 'master' of github.com:broadinstitute/gsa-unstable Ryan Poplin 2012-12-03 10:08:56 -0500
  • 6f523a1ea0 Merge branch 'master' of github.com:broadinstitute/gsa-unstable Eric Banks 2012-12-03 08:41:21 -0500
  • 59fc7456cf Updated expectations for novel TiTv in HSP after Mark's fixes to the exact model Eric Banks 2012-12-03 08:41:13 -0500
  • f0a4710247 Callset summary now includes a table for the consensus itself Mark DePristo 2012-12-02 16:40:12 -0500
  • ce9a323c04 NA12878 knowledge base automatically filters duplicate records out in the SiteIterator Mark DePristo 2012-12-02 14:21:29 -0500
  • 1bdf17ef53 Reworking of how the likelihood calculation is organized in the HaplotypeCaller to facilitate the inclusion of per allele downsampling. We now use the downsampling for both the GL calculations and the annotation calculations. Ryan Poplin 2012-12-02 11:58:32 -0500
  • 1828d33a5a Bugfix to AssessNA12878 Mark DePristo 2012-12-02 11:09:15 -0500