Commit Graph

13762 Commits (26c46ae05ee24d93c75c309a4ad74b26b494a1ba)

Author SHA1 Message Date
Ron Levine 26c46ae05e Change logger.info to logger.error 2015-01-05 14:14:02 -05:00
Ron Levine b4fda38922 Use logging system instead of stderr 2015-01-05 14:04:10 -05:00
Ron Levine 64375f6341 Messages that were going to stdout now going to stderr
Make PairHMM outputs go to stderr instead of stdout

Change output from stdout to stderr in close()

Updated lib with output going to stderr
2014-12-23 11:03:29 -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
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
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
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
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
rpoplin dc0cd7625a Merge pull request #775 from broadinstitute/vrr_hc_bamout_wrong_indel_start
Correct wrong left-alignment of reads in HC bamout
2014-11-10 15:10:16 -05:00
Valentin Ruano-Rubio c5977e5c8f Correct wrong left-alignment of reads in HC bamout
Story:
-----

  https://www.pivotaltracker.com/story/show/80684230

Changes:
-------

  - Corrected the bug: AlignmentUtils#createReadAlignedToRef was
    not realigning against the reference but the best haplotype for
    the read.

Test:
----

  - Added integration test in HaplotypeCallerIntegrationTest to check
    that the bug has been fixed.
  - Fixed md5s modified by this change; these are cause due to small
    changes in the state of the random-number generator and read vs
    variant site overlapping.
2014-11-10 10:09:58 -05:00
rpoplin 31cb47b9e6 Merge pull request #763 from broadinstitute/ldg_combineGVCFsFix
Fix bug in CombineGVCFs so now sample 2 variants occurring within sample ...
2014-11-06 12:56:15 -05:00
jmthibault79 8ddcff159e Merge pull request #765 from broadinstitute/rhl_improve_bad_file_err_msg
Catch RuntimeException and add the bad file name to the newly thown exce...
2014-11-05 10:59:38 -05:00
Ron Levine 6f78bf5041 Catch RuntimeException and add the bad file name to the newly thown exception 2014-11-05 09:14:33 -05:00
Laura Gauthier c09667a20d Fix bug in CombineGVCFs so now sample 2 variants occuring within sample 1 deletions get merged properly.
CombineGVCFs now outputs ref conf for the duration of deletions so that SNPs occuring in other samples aligned with those deletions will be genotyped correctly
2014-11-05 09:11:47 -05:00
kshakir c14fccb571 Merge pull request #771 from broadinstitute/ks_faster_builds
Faster builds of GATK executables
2014-11-05 01:54:59 +08:00
Khalid Shakir 0092a0b9eb Faster builds, with updates to documentation generation.
Reading the multiple GATKText files as a single stream, especially with new top level target executable jar files pointing to a lib folder.
Don't dirty the build with a new GATKText.properties if input files are unmodified.
Stop warning on undocumented abstract classes.
Fixed ClassNotFoundException/NoClassDefFoundError by fixing ResourceBundleExtractorDoclet artifact.
Excluding Exceptions from documentation.
Removed custom log4j dependency from ResourceBundleExtractorDoclet.
Stop generating the dependency reduced pom during shade.
Stop regenerating gsalib when the files are already up to date.
Disabled mvn site generation from external-example.
2014-11-05 00:32:23 +08:00
Khalid Shakir 1cb4b99548 Added faster built executable, non-packaged jars.
Moved top level target symlinks to package jar files to under target/package.
Executable jar files are placed under target/executable with the new target[/lib] directories.
Under top level target, symlinks to *either* the package *or* the executable jars replace what was a symlink to the package jar path.
Allow disabling of the shade package.
ant-bridge.sh by default only builds executable jars, and doesn't package by default, as did the old ant build.xml.
Added a new package_path.sh utility script for other scripts to use instead of anything in the target folder.
2014-11-05 00:30:46 +08:00
rpoplin f6af776590 Merge pull request #767 from broadinstitute/gg_fix_gsalib_doc
Updated gsalib package to match version we put in CRAN
2014-11-04 10:46:27 -05:00
rpoplin 65d834fece Merge pull request #769 from broadinstitute/pd_sb_plusplus
Added StrandAlleleCountsBySample annotation
2014-11-04 10:45:46 -05:00
Phillip Dexheimer 10f99cbe04 Added StrandAlleleCountsBySample annotation
This annotation outputs the number of reads supporting each allele, stratified by sample and read strand.
Addresses PT 76958712
2014-11-03 21:35:58 -05:00
droazen 20bab93634 Merge pull request #772 from broadinstitute/ks_disable_lsf_tests
Disabling tests for Lsf706 specific functionality.
2014-11-03 16:24:57 -05:00
Khalid Shakir 8b81031bf8 Disabling tests for Lsf706 specific functionality. 2014-11-04 01:31:18 +08:00
rpoplin f93b5b2c39 Merge pull request #770 from broadinstitute/pd_ngs_platform
Moved platform flow information into NGSPlatform
2014-11-03 10:15:17 -05:00
Phillip Dexheimer bcfd9ce19a Moved platform flow information into NGSPlatform
* Explicitly added a type for rarely used platforms
 * PT 81767718
2014-10-31 22:27:34 -04:00