This is not exactly new (it was fixed in GATK 3.0) but it's come to our attention that many people are unaware of this bug, so we want to spread the word since it might have some important impacts on people's results.

Affected versions: 2.x versions up to 2.8 (not sure when it started)

Affected tool: SelectVariants

Trigger conditions: Extracting a subset of samples with SelectVariants while using multi-threading (-nt)

Effects: Genotype-level fields (such as AD) swapped among samples

This bug no longer affects any tools in versions 3.0 and above, but callsets generated with earlier versions may need to be checked for consistency of genotype-level annotations. Our sincere apologies if you have been affected by this bug, and our thanks to the users who reported experiencing this issue.


ArthurGoldberg


Also SelectVariants reordered PL values in genotypes of variants with multiple alternate alleles.

Fri 11 Apr 2014




At a glance



Follow us on Twitter

GATK Dev Team

@gatk_dev

@skanwal23 @geoffjentry @WDL_dev Google setting up public WDL service; we'll provide WDLs & test data buckets. Max reproducibility, soon!
30 Aug 16
@skanwal23 @wdl_dev Note: the WDL @geoffjentry linked to is the 1st of 2 parts of our WGS pipeline. Second part not yet shared (but soon)
30 Aug 16
RT @WDL_dev: @BroadGenomics Storage pls. How does volume of bio material compare to volume of servers needed to store the data? https://t.…
29 Aug 16
RT @BroadGenomics: After the Samples Lab, we move into Extractions or Storage. What would you like to see next? #BroadGenomicsTour https://…
29 Aug 16
@BrandonDWomack Sure, just ask in the GATK support forum https://t.co/6z2bZ2AeN0
29 Aug 16

Our favorite tweets from others

it's the nightly build owl for GATK :D https://t.co/OwTRrk6KHA https://t.co/rfmAbdIIQp
11 Aug 16
We're going to make an hg38 version of ExAC. And we'll make @dgmacarthur pay for it. #BioinformaticsCampaignPromises
2 Aug 16
You’re right @gatk_dev honesty is key! About variants manual filtering: “In any case you're probably in for a world of pain.” Ha now I know!
11 Jul 16
.@gatk_dev I like the new documentation index page, the subheading has made my day! :D #doge #GeekHumourFTW https://t.co/9RXnDTMoBm
8 Jul 16
There is no NGS, NG is today so should only be called high-throughput sequencing #CSC #GATKworkshop https://t.co/paHcNimD7o
16 Jun 16
See more of our favorite tweets...
Search blog by tag

ad appistry ashg benchmarks best-practices bug bug-fixed cancer catvariants challenge cloud cluster commandline commandlinegatk community competition compute conferences cram cromwell denovo depthofcoverage diagnosetargets error fix forum gatk3 genotype genotype-refinement genotypegvcfs google grch38 gvcf haploid haplotypecaller hg38 holiday hts htsjdk ibm java8 job job-offer jobs joint-discovery license meetings mendelianviolations multisample multithreading mutect mutect2 ngs nt outreach pairhmm parallelism patch performance phone-home picard pipeline plans ploidy polyploid poster presentations printreads profile promote reference-model release release-notes rnaseq runtime saas script search selectvariants sequencing service slides snow speed status sting support syntax talks team terminology third-party-tools topstory trivia troll tutorial unifiedgenotyper variantannotator variantrecalibrator vcf-gz version-highlights versions vqsr wdl webinar workflow workshop