Commit Graph

13778 Commits (1d5996dd1de1cf195baebfe0a3b7a2f69913bb5d)

Author SHA1 Message Date
rpoplin 1d5996dd1d Merge pull request #796 from broadinstitute/rp_update_bundle_Axiom_poly
Adding Axiom_Exome_Plus.sites_only.all_populations.poly.vcf to the resou...
2015-01-05 14:53:45 -05:00
Ryan Poplin 10b23bfb04 Adding Axiom_Exome_Plus.sites_only.all_populations.poly.vcf to the resource bundle because it is used in the v3.3 best practices 2015-01-05 14:52:31 -05:00
rpoplin 3240b3538a Merge pull request #794 from broadinstitute/rhl_read_backed_phasing
Rhl read backed phasing
2015-01-05 09:47:25 -05:00
rpoplin 6f566a0e66 Merge pull request #793 from broadinstitute/mf_xhmm_add_genotyping
Add option to genotype additional user-defined interval lists
2015-01-05 09:46:06 -05:00
rpoplin ac10b9c3cd Merge pull request #791 from broadinstitute/ldg_fixCGPbug
Add error handling for genotypes that are called but have no PLs
2015-01-05 09:40:38 -05:00
Ron Levine c6840124fe clean up, add final 2015-01-04 23:01:24 -05:00
Ron Levine 85dc703461 Add TestMergeIntoMNP() and TestReallyMergeIntoMNP() 2015-01-01 09:51:20 -05:00
Ron Levine bb94833750 Add more tests 2014-12-30 22:45:44 -05:00
Ron Levine 714d575e3b correct reference file name 2014-12-25 14:00:39 -05:00
Ron Levine a7fba5c209 restructure and add more tests 2014-12-25 13:57:54 -05:00
Menachem Fromer 11cd0080c3 Add option to genotype additional user-defined interval lists
Add Qscript 'ONLY_GENOTYPE_xhmmCNVpipeline' to genotype additional user-defined interval lists

Add Qscript 'ONLY_GENOTYPE_xhmmCNVpipeline' to genotype additional user-defined interval lists (and similar option to Qscript 'xhmmCNVpipeline')
2014-12-21 13:02:17 -05:00
Ron Levine 069398ad46 Added more tests and documentation 2014-12-19 12:57:43 -05:00
Laura Gauthier a9694951d2 Add error handling for genotypes that are called but have no PLs 2014-12-18 15:03:20 -05:00
ldgauthier a4a9e73ec8 Merge pull request #789 from broadinstitute/gg_update_vqsr_example
Updated VQSR tool docs
2014-12-18 13:32:15 -05:00
Geraldine Van der Auwera b0e615251b Updated VQSR tool docs 2014-12-18 12:59:37 -05:00
rpoplin 4a2ac38308 Merge pull request #790 from broadinstitute/rp_nsubtil_fix-snp-detection
BQSR bug fix from @nsubtil
2014-12-18 09:19:53 -05:00
rpoplin 071389b507 Merge pull request #785 from broadinstitute/rhl_variant_multiallelic_annotations
Rhl variant multiallelic annotations
2014-12-18 09:11:42 -05:00
Ron Levine 08790e1dab Fix mmultiallelic info field annotation for VariantAnnotator
Add multi-allele test for info field annotations

Fix to process all types of INFO annotations

roll back to previous version, removes INFO and FORMAT

Correct @return for VariantAnnotatorEngine.getNonReferenceAlleles()

Enhance comments and clean up multi-allelic logic, handle header info number = R

only parse counts of A & R

Add INFO for AC

update MD5

Performance enhancement, only parse multiallelic with a count A or R

Make argument final in getNonReferenceAlleles()

Code cleanup, add exceptions for bad expression/allele size mismatch and missing header info for an expression

Change exception to warning for expression value/number of alleles check

remove adevertised exceptions
2014-12-17 22:21:00 -05:00
Ron Levine ba949389c5 matchHaplotypeAlleles() no longer calls alleleSegregationIsKnown(), added a TODO to investigate 2014-12-17 14:02:24 -05:00
Ryan Poplin d84970ff75 BQSR bug fix from @nsubtil
-- Ignore SNP matches that lie outside the clipped read window
-- This fixes an issue where GATK would skip the entire read if a SNP is entirely
contained within a sequencing adapter.
2014-12-17 10:04:37 -05:00
Ron Levine 56f8e4f9cf Add comments, alleleSegregationIsKnown() check is added to matchHaplotypeAlleles() 2014-12-17 03:25:26 -05:00
Valentin Ruano Rubio 3deaa4832c Merge pull request #788 from broadinstitute/ldg_fixCombineGVCFsBug
Fixed huge bug from 9895005a (CombineGVCFs used to stop after the first ...
2014-12-16 15:17:33 -05:00
Laura Gauthier 011843c569 Fixed huge bug from 9895005a (CombineGVCFs used to stop after the first contig) 2014-12-16 12:43:32 -05:00
rpoplin bcc6b73e9b Merge pull request #786 from broadinstitute/pd_variantstotable_sma
Fix VariantsToTable output of FORMAT record lists when -SMA is specified
2014-12-16 10:37:22 -05:00
rpoplin 3d32babc29 Merge pull request #787 from broadinstitute/vrr_wrong_reference_base_bugfix
Fixing CombineGVCFs that writes out the wrong REF allele
2014-12-16 10:35:47 -05:00
Valentin Ruano-Rubio 736a857e82 Fixing CombineGVCFs that writes out the wrong REF allele
Story:
=====

  - https://www.pivotaltracker.com/story/show/83259038

Changes:
=======

  - Done minimal changes to make the fix after an arduous attempt to understand
    CombineGVCFs code.

Test:
====

  - Added a integration test to explicitly test for the bug.

  - Updated a md5 changes as the bug was actually affecting one of the existing
    integration tests.
2014-12-13 22:38:24 -05:00
Phillip Dexheimer 71bdfbe465 Fix VariantsToTable output of FORMAT record lists when -SMA is specified
* PT 84242218
 * Note that FORMAT fields behave the same as INFO fields - if the annotation has a count of A (one entry per Alt Allele), it is split across the multiple output lines.  Otherwise, the entire list is output with each field
2014-12-10 21:41:15 -05:00
Geraldine Van der Auwera b561febbfd Merge pull request #784 from broadinstitute/gg_update_gsalib_version
Updated gsalib version to 2.1 for resubmitting with updated license to C...
2014-12-10 13:23:47 -05:00
Geraldine Van der Auwera 45eddb4ecb Updated gsalib version to 2.1 for resubmitting with updated license to CRAN 2014-12-09 17:07:48 -05:00
rpoplin bf2911d62c Merge pull request #783 from broadinstitute/pd_splitsamfile
Fix NPE in SplitSamFile
2014-12-08 09:39:03 -05:00
rpoplin 5f2a2642d6 Merge pull request #782 from broadinstitute/vrr_gvcf_gq_inconsistency_fix
Makes GQ of Hom-Ref Blocks in GVCF output to be consistent with PLs
2014-12-08 09:36:23 -05:00
Valentin Ruano-Rubio 385186e11b Makes GQ of Hom-Ref Blocks in GVCF output to be consistent with PLs
Story:
-----

  - https://www.pivotaltracker.com/story/show/83800586

Changes:
-------

  - In GVCFWriter GQ is now recalculated out of the fianl PL array for the block.

Testing:
-------

  - Updated affected integration test md5s
2014-12-07 16:45:32 -05:00
Phillip Dexheimer a5dee8a42e Fix NPE in SplitSamFile
* PT 82892316
  * Added integration test
  * Fixed similar error in debug output of HC
2014-12-07 10:37:30 -05:00
Ron Levine c9175eeee8 Renamed PhasingUtilitiesUnitTest to PhasingUtilsUnitTest 2014-12-02 18:20:12 -05:00
rpoplin bbe78bd181 Merge pull request #779 from broadinstitute/rhl_output_pairhmm
Rhl output pairhmm
2014-12-02 14:30:55 -05:00
rpoplin 4e8297360d Merge pull request #780 from broadinstitute/aw_gatkdoclet_fixes
Aw gatkdoclet fixes
2014-12-02 10:57:57 -05:00
Ron Levine b8f0f3fdd2 Add argument for loading the vector HMM library once 2014-12-02 10:13:56 -05:00
Alec Wysoker 4fe6ccec98 Add -output-file-extension option to GATKDoclet to produce html instead of php. 2014-12-01 18:06:36 -05:00
Alec Wysoker 62e5d42380 Fix code to filter current directory from paths pass to Reflection library. 2014-12-01 17:45:46 -05:00
droazen 2fa8fe2415 Merge pull request #777 from broadinstitute/rhl_specify_pairhmm
Rhl specify pairhmm
2014-11-26 11:34:51 -05:00
Ron Levine 386aeda022 Add HaplotypeCaller argument so integration tests can specify the hardware dependent PairHMM sub-implementation 2014-11-25 21:53:53 -05:00
Ron Levine 34241a62f6 Use a publicly accessible sequence file 2014-11-24 11:18:21 -05:00
Ron Levine 6ff698c556 Added HP and non-HP tests for matchHaplotypeAlleles(), added a nominal test for mergeIntoMNPvalidationCheck() 2014-11-24 11:08:04 -05:00
Ron Levine 61e1a3ecd1 Added the framework for testing the PhasingUtilies methods matchHaplotypeAlleles() and reallyMergeIntoMNP() 2014-11-22 22:01:39 -05:00
rpoplin ba5a32107f Merge pull request #778 from broadinstitute/mf_RBP_MNP_merging
Fix MNP merging bugs
2014-11-21 16:00:42 -05:00
Menachem Fromer 9b73c8a841 Fix MNP merging bugs 2014-11-21 06:42:51 -05:00
rpoplin 00027e1555 Merge pull request #774 from broadinstitute/ldg_makeSelectVariantsTrimAlleles
Add -trim argument to SelectVariants to trim alleles to minimal represen...
2014-11-13 13:58:13 -05:00
ldgauthier bcac930aad Merge pull request #773 from broadinstitute/rhl_annotator_warning
Make annotators emit a warning if they can't be applied
2014-11-12 08:01:19 -05:00
Ron Levine 67656bab23 Resolved conflict during rebasing
Add more logging to annotators, change loggers from info to warn

Add comments to testStrandBiasBySample()

Clarify comments in testStrandBiasBySample

remove logic for not prcossing an indel if strand bias (SB) was not computed

remove per variant warnings in annotate()

Log warnings if using the wrong annotator or missing a pedgree file

Log test failures once in annotate(), because HaplotypeCaller does not call initialize(). Avoid using exceptions

Fix so only log once in annotate(), Hardey-Weinberg does not require pedigree files, fix test MD5s so pass

Check if founderIds == null

Update MD5s from HaplotypeCaller integrations tests and clean up code

Change logic so SnpEff does not throw excpetions, change engine to utils in imports

Update test MD5s, return immediately if cannot annotate in SnpEff.initialization()

Post peer review, add more logging warnings

Update MD5 for testHaplotypeCallerMultiSampleComplex1, return null if PossibleDeNovo.annotate() is not called by VariantAnnotator
2014-11-12 02:45:49 -05:00
Laura Gauthier 783a4fd651 Change default behavior of SelectVariants to trim remaining alleles when samples are subset. -noTrim argument preserves original alleles. Add test for trimming. 2014-11-11 16:32:25 -05:00