Commit Graph

8218 Commits (32adbd614f9dc39b09fb3aa3fd9a2e910e3b4e24)

Author SHA1 Message Date
Ryan Poplin 0534149708 Merged bug fix from Stable into Unstable 2011-11-07 14:07:08 -05:00
Ryan Poplin 2d1e385ca4 Adding note to VQSR docs about Rscript being needed in the environment PATH. 2011-11-07 14:04:13 -05:00
Mark DePristo bb448b27d2 Merge branch 'master' of ssh://gsa1/humgen/gsa-scr1/gsa-engineering/git/unstable 2011-11-07 13:26:37 -05:00
Mark DePristo 63eeb5a56c MultiSampleCalling script update 2011-11-07 13:23:34 -05:00
Eric Banks 759f4fe6b8 Moving unclaimed walker with bad integration test to archive 2011-11-07 13:16:38 -05:00
Eric Banks aa0c8c3600 Revving Tribble jar to v40. Our last jar was busted. 2011-11-07 11:30:08 -05:00
Eric Banks c1986b6335 Add notes to the GATKdocs as to when a particular annotation can/cannot be calculated. 2011-11-07 11:06:19 -05:00
Eric Banks 724e3f3b0d Merged bug fix from Stable into Unstable 2011-11-06 22:23:22 -05:00
Eric Banks cdd40d1222 Removing contracts for the SimpleTimer 2011-11-06 22:22:49 -05:00
Ryan Poplin 2f2c81dccb Adding argument to control scatter count 2011-11-06 12:56:22 -05:00
Ryan Poplin 4a1e8ecbb7 Updating HaplotypeCaller qscript to work with most recent Queue 2011-11-06 12:53:54 -05:00
Ryan Poplin 5c565d28b9 Merge branch 'master' of ssh://nickel.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable 2011-11-06 10:26:19 -05:00
Ryan Poplin ebdced412c misc cleanup 2011-11-06 10:26:07 -05:00
Eric Banks 3517489a22 Better --sample selection integration test for VE. The previous one would return true even if --sample was not working at all. 2011-11-06 01:07:49 -04:00
Eric Banks 1c4e429a1c Merge branch 'master' of ssh://nickel.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable 2011-11-06 00:05:56 -04:00
Eric Banks a12bc63e5c Get rid of support for bams without sample information in the read groups. This hidden option wasn't being used anyways because it wasn't hooked up properly in the AlignmentContext. 2011-11-05 23:54:28 -04:00
Eric Banks ad57bcd693 Adding integration test to cover using expressions with IDs (-E foo.ID) 2011-11-05 23:53:15 -04:00
Eric Banks e89fe7770e Merged bug fix from Stable into Unstable 2011-11-05 22:41:21 -04:00
Eric Banks 90a053ea93 Don't change the mapping quality of MQ=255 reads in IR 2011-11-05 22:40:45 -04:00
Ryan Poplin 611a395783 Now properly extending candidate haplotypes with bases from the reference context instead of filling with padding bases. Functionality in the private Haplotype class is no longer necessary so removing it. No need to have four different Haplotype classes in the GATK. 2011-11-05 12:18:56 -04:00
Mauricio Carneiro fd09d92801 Added contracts to SyntheticRead class 2011-11-04 16:27:11 -04:00
Mauricio Carneiro 839cab7427 Generalizing RunningConsensus to SyntheticRead
the class structure of the RunningConsensus is now going to represent both consensus and missing data objects.
2011-11-04 16:27:11 -04:00
Mauricio Carneiro 56674e3518 Removing SlidingReads
Since we're not doing '=' conversion inside the Reduce Reads walker, there is no real need for this class.
2011-11-04 16:27:11 -04:00
Mauricio Carneiro 18837e3e32 Moving SlidingRead functionality into SlidingWindow 2011-11-04 16:27:11 -04:00
Mauricio Carneiro bf822172c7 Renaming the compressor classes
Naming scheme was confusing for the consensus classes. Now that we'll have multiple types of running consensus, I felt like a new naming scheme was necessary.
2011-11-04 16:27:11 -04:00
Mauricio Carneiro 56a6bb6e98 Fixing BAMRecord illegal access
Can't use getReadName on a read that has been completely hardclipped. The BAMRecord doesn't like it.
2011-11-04 16:27:11 -04:00
Eric Banks de07e06cbc Merge remote-tracking branch 'unstable/master' 2011-11-04 15:57:13 -04:00
Ryan Poplin 888d3b4fdc Initial graph pruning algorithm for the assembler. 2011-11-04 14:02:50 -04:00
Mark DePristo e99871f587 Bug fix for decode loc
-- decodeLoc() wasn't skipping input header lines, so the system blew up when there was an = line being split.
2011-11-04 13:20:54 -04:00
Mark DePristo a340a1aeac Bug fix. decodeLoc() should update lineNo so you get meaningful line no when indexing
due to malformed VCF files.
2011-11-04 11:44:24 -04:00
Mark DePristo 849c0757f2 Bug fix for LocusScatterFunction when no intervals are provided
-- Now correctly grabs reference contigs and cuts them all up, rather than NPE as intervalString == null.
2011-11-04 10:55:09 -04:00
Mark DePristo 9f260c0dc1 Zero byte index bug fix for RandomlySplitVariants + cleanup
-- vcfWriter2 was never being closed in onTraversalDone(), so the on the fly index file was being created but never actually properly written to the file.

-- This bug is ultimately due to the inability of the GATK to allow multiple VCF output writers as @Output arguments, though

-- Removed the unnecessary local variable iFraction, = 1000 * the input fraction argument.  Now the system just uses a double random number and compares to the input fraction at all.  Is there some subtle reason I don't appreciate for this programming construct?
2011-11-04 09:45:20 -04:00
Mark DePristo 5a47c3c8a0 Merge branch 'master' of ssh://gsa1/humgen/gsa-scr1/gsa-engineering/git/unstable 2011-11-04 09:36:42 -04:00
Mauricio Carneiro 020b8b88ef GATKSAMRecord refactor in the tools
No tools should create SAMRecords internally. This commit should move all internals of the current tools to GATKSAMRecord.
2011-11-03 17:33:42 -04:00
Mauricio Carneiro e89ff063fc GATKSAMRecord refactor
The GATK engine will now provide a GATKSAMRecord to all tools which incorporates the functionality used by the GATK to the bam file (ReadGroups, Reduced Reads, ...).

* No tools should create SAMRecord anymore, use GATKSAMRecord instead *
2011-11-03 15:43:26 -04:00
Ryan Poplin f1df6c0c81 Misc cleanup in haplotype caller after incorporating Mark's FragmentCollection to merge overlapping read pairs. 2011-11-03 13:51:19 -04:00
Mark DePristo 748f8f1edc Merge branch 'master' of ssh://gsa1/humgen/gsa-scr1/gsa-engineering/git/unstable 2011-11-03 10:06:48 -04:00
Laurent Francioli 385a6abec1 Fixed a bug that wrongly swapped the mother and father genotypes in case the child genotype missing. 2011-11-03 13:04:53 +01:00
Laurent Francioli 893787de53 Functions getAsMap and getNegLog10GQ now handle missing genotype case. 2011-11-03 13:04:11 +01:00
Mark DePristo c7f51e92a0 Mostly working version of multi-sample analysis qscript 2011-11-02 23:00:17 -04:00
Eric Banks e8bceb1eaa Merge branch 'master' of ssh://nickel.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable 2011-11-02 21:13:54 -04:00
Eric Banks 78a00d2ddc Updating UG integration tests (needed updating only because the -mbq default is different from the old -mmq one). 2011-11-02 21:13:44 -04:00
Eric Banks 52b16bf739 Must check whether there's a normal vs. extended pileup before asking for it. 2011-11-02 20:45:24 -04:00
Eric Banks e1edd6bd12 Removing the min mapping quality argument since it wasn't being used in the normal processing of the pileups in UG - only for indel pileups. Instead, we apply the min base quality to the reads in the pileup for indels and define it to be the min 'confidence' of the base. Docs are updated but I didn't rename the argument as I don't want people to complain. 2011-11-02 20:32:58 -04:00
Mauricio Carneiro c22a14ee3b Merged bug fix from Stable into Unstable 2011-11-02 17:53:56 -04:00
Mauricio Carneiro e4a583a53f Fixing docs: No -I in this walker 2011-11-02 17:53:32 -04:00
Ryan Poplin e94fcf537b Merge branch 'master' of ssh://nickel.broadinstitute.org/humgen/gsa-scr1/gsa-engineering/git/unstable 2011-11-02 16:29:19 -04:00
Ryan Poplin 4d35272916 Bug fixes with Mauricio to functions in ReadUtils used by reduced reads and the haplotype caller. 2011-11-02 16:29:10 -04:00
Mark DePristo 8a2929c1dd Merge branch 'master' of ssh://gsa1/humgen/gsa-scr1/gsa-engineering/git/unstable 2011-11-02 16:21:00 -04:00
Mark DePristo e2f40da27f Scala script to run multi-sample analysis 2011-11-02 16:20:57 -04:00