Ron Levine
9d4b876ccd
Process X and = CIGAR operators
...
Add simple BaseRecalibrator integration test for CIGAR = and X operators
2015-01-29 17:00:00 -05:00
Valentin Ruano Rubio
e26e55efe1
Merge pull request #802 from broadinstitute/pd_selectvariants_subset
...
Added -trimAlternates argument to SelectVariants
2015-01-22 05:05:42 -05:00
Phillip Dexheimer
72f76add71
Added -trimAlternates argument to SelectVariants
...
* PT 84021222
* -trimAlternates removes all unused alternate alleles from variants. Note that this is pretty aggressive for monomorphic sites
2015-01-21 21:33:35 -05:00
jmthibault79
d209699485
Merge pull request #804 from broadinstitute/jt_read_filters
...
Allows users to disable the DuplicateReadFilter
2015-01-21 15:45:48 -05:00
Joel Thibault
5ce34d81b8
Allows users to disable specific read filters from the command line
...
- enable this for DuplicateReadFilter only
- enable the @DisabledReadFilters annotation to do this at the Walker level
2015-01-21 13:17:29 -05:00
Geraldine Van der Auwera
7102477168
Merge pull request #798 from broadinstitute/rhl_split_n_cigar_reads_exception
...
Rhl split n cigar reads exception
2015-01-15 00:20:49 -05:00
Ron Levine
804b2a36b7
Fix SplitNCigar reads exception by making the list of RNAReadTransformer non-abstract, add test for -fixNDN
...
Includes documentation changes for -fixNDN argument and the read transformer documentation.
Documentation changes to CombineVariants
2015-01-14 22:22:05 -05:00
rpoplin
0292d49842
Merge pull request #801 from broadinstitute/pd_gatkvcfconstants
...
Collected VCF IDs and header lines into one place
2015-01-14 09:43:48 -05:00
rpoplin
8599132713
Merge pull request #800 from broadinstitute/ldg_CGPpedError
...
Address user-reported bug featuring "trio" family with two children, one...
2015-01-14 09:35:33 -05:00
Phillip Dexheimer
6190d660e0
Edits to work with the latest htsjdk release:
...
* TextCigarCodec.decode() is now static, and the getSingleton() method is gone
* MergingSamRecordIterator now wants a Collection<SamReader> rather than Collection<SAMFileReader> in the constructor
* SeekableBufferedStream now correctly reads the requested number of bytes, removed workaround in GATKBAMIndex
2015-01-13 21:32:10 -05:00
Phillip Dexheimer
b73e9d506a
Added GATKVCFConstants and GATKVCFHeaderLines to consolidate the GATK-specific VCF annotations
...
* Removed unused annotations (CCC and HWP)
* Renamed one of the two GC annotations to "IGC" (for Interval GC)
* Revved picard & htsjdk (GATK constants are now removed from htsjdk)
* PT 82046038
2015-01-13 21:32:09 -05:00
Laura Gauthier
6b2bd5ed09
Address user-reported bug featuring "trio" family with two children, one parent
...
Add test to cover case with family of one parent, two children
2015-01-13 18:35:44 -05:00
rpoplin
cf4ac4a388
Merge pull request #799 from broadinstitute/rp_fix_active_region_downsampling
...
Raising per-sample limits on the number of reads in ART and HC.
2015-01-09 12:45:58 -05:00
Ryan Poplin
2e5f9db758
Raising per-sample limits on the number of reads in ART and HC.
...
-- Active Region Traversal was using per sample limits on the number of reads that were too low, especially now that we are running one sample at a time. This caused issues with high confidence variants being dropped in high coverage data.
-- HaplotypeCallerGVCFIntegrationTest PL/annotation changes due to using more reads in those tests
-- Removed a CountReadsInActiveRegionsIntegrationTest test for excessive coverage because the read coverage no longer goes over the limits in ART
2015-01-09 11:21:42 -05:00
rpoplin
03203e249e
Merge pull request #792 from broadinstitute/rhl_pairhmm_log_stderr
...
Rhl pairhmm log stderr
2015-01-07 12:41:10 -05:00
rpoplin
99ff5fb75b
Merge pull request #797 from broadinstitute/vrr_ploidy_limitation_removal
...
MERGE AFTER TEST PASS: Fixes explicit limitation of the maximum ploidy of the reference-confide...
2015-01-07 12:38:19 -05:00
Valentin Ruano-Rubio
aae04b6122
Fixes explicit limitation of the maximum ploidy of the reference-confidence model
...
Story:
=====
- https://www.pivotaltracker.com/story/show/83803796
Changes:
=======
- From a fix maximum ploidy indel RCM likelihood cache to a
dynamically resizable one.
- Used the occassion to removed an unused and deprecated method from ReferenceConfidenceModel
Testing:
=======
- Added integration test to check on ploidies larger than the previous limit of 20.
2015-01-07 10:43:22 -05:00
Ron Levine
7d58544f17
Do not use logger, write to stderr, could not get the correct logger dependency in pom.xml
2015-01-06 10:32:11 -05:00
rpoplin
ef32c44688
Merge pull request #795 from broadinstitute/ldg_VQSRcmdlineOutput
...
VQSR VCF header command line now contains annotations and tranche levels
2015-01-05 14:59:04 -05:00
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
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
Laura Gauthier
e28ee24fc6
Update private single sample pipeline script
2015-01-05 12:11:21 -05:00
Laura Gauthier
88b6f3aa50
Change []-type arrays to lists so argument parsing works in VCF header commandline output
2015-01-05 10:21:06 -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
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
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