3c88e6859f
fix 1430 for genotype filters; refactored filter() method; added unit and integration test; more comprehensive fix must be done first on htsjdk side in JEXLMap ( #1456 )
2016-09-06 17:30:18 -04:00
601c26a592
More small refactorings of Mutect2 code
2016-08-29 12:45:13 -04:00
a883905101
Remove -stand_emit_conf argument
2016-08-26 22:09:11 -04:00
bc3b3ac0ec
Cleaned up SomaticGenotypingEngine::callMutations and added some TODOs.
2016-08-26 14:12:05 -04:00
4aede99697
Merge pull request #1453 from broadinstitute/db_mutect2
...
Lots of small improvements to Mutect2 code
2016-08-23 18:35:38 -04:00
cd5c04f806
Make getElementForRead() in RankSumTest robust
2016-08-23 14:56:32 -04:00
01142dfb1c
Lots of small improvements to Mutect2 code
2016-08-23 09:50:48 -04:00
abc4d5b7b3
Bypass spanning deletions in Rank Sum tests
2016-08-17 14:02:22 -04:00
3510906c7f
addresses issue #1280 now interval padding works for exclude intervals
2016-08-10 13:45:45 -04:00
49507faaa3
Changed maximum allowed GQB value to 100.
2016-08-05 13:06:31 -04:00
a465c87ff8
Added support for directly reading SRA runs
2016-08-02 15:21:14 -04:00
832a383acd
Fixed MD5 broken by PR #1440 .
2016-07-27 13:51:20 -04:00
9a6ce7a347
Merge pull request #1440 from broadinstitute/sl_issue_1345
...
Added exception for GQB values greater than MAX_GENOTYPE_QUAL and tests.
2016-07-26 14:55:59 -04:00
fef63ce6a8
Make sure that multi-alleleic uninformative PLs (0,0,...,0) stay uninformative after biallelization.
...
Addresses issue #1439 (thus #1437 ).
Fixes a bug where non informative PLs were not handled appropriatelly when calculating multi-allelic site QUAL values.
This was resulting in long execution times for very large datasets (~200,000 samples in the case of ExAC2).
2016-07-25 17:19:03 -04:00
3daed9e5a1
Added exception for GQB values greater than MAX_GENOTYPE_QUAL and tests.
2016-07-20 16:48:59 -04:00
7392c4d1b0
Removed spanning deletions if the deletion was removed
2016-07-19 12:23:49 -04:00
641382eb8b
Fix BetaTestingAnnotation group
...
Add test
2016-07-13 16:05:21 -04:00
d6d0678b50
Build on Laura's code and finish porting MuTect1 clustered read position filter.
2016-07-11 17:33:08 -04:00
9b32cf5291
Fixed merging of GVCF blocks by fixing rounding of GQ values in ReferenceConfidenceModel.
2016-07-06 10:08:08 -04:00
2c94f74a95
Merge pull request #1404 from broadinstitute/ldg_M2_addM1filters
...
MuTect 2: port strand artifact filter from MuTect 1
2016-07-05 13:27:35 -04:00
45607d1b30
RCM Variant sites merger won't output PL when there are too many alleles in order to avoid memory issues with large cohort runs.
...
Small additional "cosmetic" changes to the code
Addresses issue #1419 .
2016-07-01 11:33:07 -04:00
1ff234e7dd
remove alt alleles, when genotype count is explosively large, based on alleles' highest supporting haplotype score; max tolerable genotype count is controlled by a default value overridable by user
...
remove alt alleles, when genotype count is explosively large, based on alleles' highest supporting haplotype score; max tolerable genotype count is controlled by a default value overridable by user
2016-06-30 22:36:49 -04:00
63e0865491
Built on Laura's code to port the strand bias filter from M1 and refactored code aroud SomaticGenotypingEngine. Added a new integration test.
2016-06-29 22:46:40 -04:00
4066bcd75c
Add new annotator for M1 clustered read position filter and M1 strand bias filter.
2016-06-29 22:46:37 -04:00
1b921666a7
Change to max value of ExcessHet
2016-06-29 16:33:50 -04:00
556cc69185
Fix for int overflow in RankSum calculation
2016-06-29 12:02:13 -04:00
07052ba8ea
Changes to use the median rather than the second best likelihood for the NON_REF allele
...
Addresses issue #1378 following the first proposal using the 'median' rather than the 'mean'.
2016-06-28 13:10:22 -04:00
76bb8fd9e5
Allows GatherBqsrReports to accept a .list file as input.
2016-06-22 12:39:18 -04:00
427645162b
SelectVariants works with non-diploids
2016-06-21 12:26:13 -04:00
857459e420
Silly mistake '<' for a '<='. It was causing the exception when the exact number of alleles to drop was matching MAX_DROPPED_ALTERNATIVE_ALLELES_TO_LOG exactly (fixed to 20).
...
I change the code to impose a maximum allele list message length instead and in the process I fixed the bug.
2016-06-17 15:22:59 -04:00
e119feee61
Added regression test for genotyping of spanning deletions in GenotypeGCVFs.
2016-06-15 09:48:26 -04:00
ba2e7be05b
Add integration test using -maxNumPLValues for GenotypeGVCFs
2016-06-07 14:38:12 -04:00
85dce75f3f
Update pom versions to mark the start of GATK 3.7 development
2016-06-01 17:21:48 -04:00
f185a75e1c
Update pom versions for the 3.6 release
2016-06-01 17:08:17 -04:00
b95b76b0e2
Merge pull request #1394 from broadinstitute/gvda_add_colt_dependency
...
Add colt > cern.jet.normal dependency
2016-06-01 14:16:57 -04:00
bd2626bea2
Add colt > cern.jet.normal dependency
2016-06-01 13:24:50 -04:00
30665c7dbc
Move Move htsjdk and picard to version 2.4.1
2016-05-31 22:36:38 -04:00
a76cb052e2
Ability to retry building VQSR model (contributed by mdp)
2016-05-31 18:57:55 -04:00
d87345cd1d
GATKDocs overhaul
...
- Fixed displaying of default values
- Removed code cruft
- Reorganized tooldoc categories and improved names
- Reorganized tools within categories where applicable
- Touched up various tool docs
- Switched default gatkdocs output to html
- Added parameter in agrregator pom to control output type
- Set gatkdocs publishing script to output php
- Deprecated GenotypeAndValidate walker
- Added back PhoneHome arguments with @Deprecated annotations
2016-05-29 16:35:08 -04:00
efbbbb1bd9
Add M2 to the HC annotations check
2016-05-27 13:49:31 -04:00
c4a06ad20a
Move indel realignment to public
2016-05-27 12:39:58 -04:00
9d32dec9cd
Fix for the sum(AD) > DP bug.
...
Closes issue #1340
2016-05-26 15:04:52 -04:00
25fa25b618
Added option to validate gvcf (for ValidateVariants) ( #1379 )
...
* with option --gvcf CLP will now put extra checks that a gvcf must adhere to (existance of <NON_REF> allele at every variant, and that the variants in total cover the entire requested intervals, or the whole genome if no intervals have been specified)
* works on gvcf produced by HC when using either GVCF or BP_RESOLUTION mode
* added positive and negative tests
2016-05-26 06:42:45 -04:00
e1fadae139
Fix error in InfiniteRandomMatingPopulationModel.getLikelihoodsCalculator
...
Same issue noticed in GATK4 [here](https://github.com/broadinstitute/gatk/issues/1856 )
2016-05-25 17:23:26 -04:00
2c8356519c
Merge pull request #1375 from broadinstitute/gvda_fix_offbyone_maxAltAlls_M2_#1297
...
Fixed M2 max alt alleles threshold evaluation error
2016-05-19 13:59:02 -04:00
7fdd3c2a0c
Merge pull request #1358 from broadinstitute/sl_issue_1327
...
Changed calls with GQ=0 to no-call for HaplotypeCaller in normal mode.
2016-05-19 12:26:35 -04:00
644076b1e1
Add fix and test for finalizing MQ annotation at BP resolution for variant and ref samples
...
Addresses issue #1356
2016-05-19 08:15:30 -04:00
f5456a3761
Fixed M2 max alt alleles threshold evaluation error
...
Also clarified some argument docs
2016-05-18 21:54:30 -04:00
bf4b1a5421
Changed calls for GQ=0 from 0/0 to ./. for HaplotypeCaller in normal mode.
2016-05-18 13:17:27 -04:00
35a06879f1
Move htsjdk and picard to version 2.3.0
2016-05-16 14:50:00 -04:00
Steve Huang
David Benjamin
Ron Levine
Takuto Sato
Peter Fan
Samuel Lee
Andrii Nikitiuk
samuelklee
Valentin Ruano Rubio
Laura Gauthier
meganshand
Geraldine Van der Auwera
Yossi Farjoun