forked from samtools/bcftools
-
Notifications
You must be signed in to change notification settings - Fork 0
/
NEWS
1537 lines (977 loc) · 54.5 KB
/
NEWS
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
## Release a.b
## Release 1.17 (21st February 2023)
Changes affecting the whole of bcftools, or multiple commands:
* The -i/-e filtering expressions
- Error checks were added to prevent incorrect use of vector arithmetics. For example,
when evaluating the sum of two vectors A and B, the resulting vector could contain
nonsense values when the input vectors were not of the same length. The fix introduces
the following logic:
- evaluate to C_i = A_i + B_i when length(A)==B(A) and set length(C)=length(A)
- evaluate to C_i = A_i + B_0 when length(B)=1 and set length(C)=length(A)
- evaluate to C_i = A_0 + B_i when length(A)=1 and set length(C)=length(B)
- throw an error when length(A)!=length(B) AND length(A)!=1 AND length(B)!=1
- Arrays in Number=R tags can be now subscripted by alleles found in FORMAT/GT. For example,
FORMAT/AD[GT] > 10 .. require support of more than 10 reads for each allele
FORMAT/AD[0:GT] > 10 .. same as above, but in the first sample
sSUM(FORMAT/AD[GT]) > 20 .. require total sample depth bigger than 20
* The commands `consensus -H` and `+split-vep -H`
- Drop unnecessary leading space in the first header column and newly print `#[1]columnName`
instead of the previous `# [1]columnName` (#1856)
Changes affecting specific commands:
* bcftools +allele-length
- Fix overflow for indels longer than 512bp and aggregate alleles equal or larger than
that in the same bin (#1837)
* bcftools annotate
- Support sample reordering of annotation file (#1785)
- Restore lost functionality of the --pair-logic option (#1808)
* bcftools call
- Fix a bug where too many alleles passed to `-C alleles` via `-T` caused memory
corruption (#1790)
- Fix a bug where indels constrained with `-C alleles -T` would sometimes be missed (#1706)
* bcftools consensus
- BREAKING CHANGE: the option `-I, --iupac-codes` newly outputs IUPAC codes based on FORMAT/GT
of all samples. The `-s, --samples` and `-S, --samples-file` options can be used to subset
samples. In order to ignore samples and consider only the REF and ALT columns (the original
behavior prior to 1.17), run with `-s -` (#1828)
* bcftools convert
- Make variantkey conversion work for sites without an ALT allele (#1806)
* bcftool csq
- Fix a bug where a MNV with multiple consequences (e.g. missense + stop_gained)
would report only the less severe one (#1810)
- GFF file parsing was made slightly more flexible, newly ids can be just 'XXX'
rather than, for example, 'gene:XXX'
- New gff2gff perl script to fix GFF formatting differences
* bcftools +fill-tags
- More of the available annotations are now added by the `-t all` option
* bcftools +fixref
- New INFO/FIXREF annotation
- New -m swap mode
* bcftools +mendelian
- The +mendelian plugin has been deprecated and replaced with +mendelian2. The
function of the plugin is the same but the command line options and the output
format has changed, and for this was introduced as a new plugin.
* bcftools mpileup
- Most of the annotations generated by mpileup are now optional via the
`-a, --annotate` option and add several new (mostly experimental) annotations.
- New option `--indels-2.0` for an EXPERIMENTAL indel calling model. This model aims
to address some known deficiencies of the current indel calling algorithm, specifically,
it uses diploid reference consensus sequence. Note that in the current version it
has the potential to increase sensitivity but at the cost of decreased specificity.
- Make the FS annotation (Fisher exact test strand bias) functional and remove it
from the default annotations
* bcftools norm
- New --multi-overlaps option allows to set overlapping alleles either to the
ref allele (the current default) or to a missing allele (#1764 and #1802)
- Fixed a bug in `-m -` which does not split missing FORMAT values correctly and
could lead to empty FORMAT fields such as `::` instead of the correct `:.:` (#1818)
- The `--atomize` option previously would not split complex indels such as C>GGG.
Newly these will be split into two records C>G and C>CGG (#1832)
* bcftools query
- Fix a rare bug where the printing of SAMPLE field with `query` was incorrectly
suppressed when the `-e` option contained a sample expression while the formatting
query did not. See #1783 for details.
* bcftools +setGT
- Add new `--new-gt X` option (#1800)
- Add new `--target-gt r:FLOAT` option to randomly select a proportion of genotypes (#1850)
- Fix a bug where `-t ./x` mode was advertised as selecting both phased and unphased
half-missing genotypes, but was in fact selecting only unphased genotypes (#1844)
* bcftools +split-vep
- New options `-g, --gene-list` and `--gene-list-fields` which allow to prioritize
consequences from a list of genes, or restrict output to the listed genes
- New `-H, --print-header` option to print the header with `-f`
- Work around a bug in the LOFTEE VEP plugin used to annotate gnomAD VCFs. There the
LoF_info subfield contains commas which, in general, makes it impossible to parse the
VEP subfields. The +split-vep plugin can now work with such files, replacing the offending
commas with slash (/) characters. See also https://github.com/Ensembl/ensembl-vep/issues/1351
- Newly the `-c, --columns` option can be omitted when a subfield is used in `-i/-e` filtering
expression. Note that `-c` may still have to be given when it is not possible to infer the
type of the subfield. Note that this is an experimental feature.
* bcftools stats
- The per-sample stats (PSC) would not be computed when `-i/-e` filtering options and
the `-s -` option were given but the expression did not include sample columns (1835)
* bcftools +tag2tag
- Revamp of the plugin to allow wider range of tag conversions, specifically all combinations
from FORMAT/GL,PL,GP to FORMAT/GL,PL,GP,GT
* bcftools +trio-dnm2
- New `-n, --strictly-novel` option to downplay alleles which violate Mendelian
inheritance but are not novel
- Allow to set the `--pn` and `--pns` options separately for SNVs and indels and make
the indel settings more strict by default
- Output missing FORMAT/VAF values in non-trio samples, rather than random nonsense values
* bcftools +variant-distance
- New option `-d, --direction` to choose the directionality: forward, reverse, nearest (the default)
or both (#1829)
## Release 1.16 (18th August 2022)
* New plugin `bcftools +variant-distance` to annotate records with distance to the
nearest variant (#1690)
Changes affecting the whole of bcftools, or multiple commands:
* The -i/-e filtering expressions
- Added support for querying of multiple filters, for example `-i 'FILTER="A;B"'`
can be used to select sites with two filters "A" and "B" set. See the documentation
for more examples.
- Added modulo arithmetic operator
Changes affecting specific commands:
* bcftools annotate
- A bug introduced in 1.14 caused that records with INFO/END annotation would
incorrectly trigger `-c ~INFO/END` mode of comparison even when not explicitly
requested, which would result in not transferring the annotation from a tab-delimited
file (#1733)
* bcftools merge
- New `-m snp-ins-del` switch to merge SNVs, insertions and deletions separately (#1704)
* bcftools mpileup
- New NMBZ annotation for Mann-Whitney U-z test on number of mismatches within
supporting reads
- Suppress the output of MQSBZ and FS annotations in absence of alternate allele
* bcftools +scatter
- Fix erroneous addition of duplicate PG lines
* bcftools +setGT
- Custom genotypes (e.g. `-n c:1/1`) now correctly override ploidy
## Release 1.15.1 (7th April 2022)
* bcftools annotate
- New `-H, --header-line` convenience option to pass a header line on command line,
this complements the existing `-h, --header-lines` option which requires a file
with header lines
* bcftools csq
- A list of consequence types supported by `bcftools csq` has been added to
the manual page. (#1671)
* bcftools +fill-tags
- Extend generalized functions so that FORMAT tags can be filled as well, for example:
bcftools +fill-tags in.bcf -o out.bcf -- -t 'FORMAT/DP:1=int(smpl_sum(FORMAT/AD))'
- Allow multiple custom functions in a single run. Previously the program would silently
go with the last one, assigning the same values to all (#1684)
* bcftools norm
- Fix an assertion failure triggered when a faulty VCF file with a '-'
character in the REF allele was used with `bcftools norm --atomize`. This
option now checks that the REF allele only includes the allowed characters
A, C, G, T and N. (#1668)
- Fix the loss of phasing in half-missing genotypes in variant atomization (#1689)
* bcftools roh
- Fix a bug that could result in an endless loop or incorrect AF estimate when
missing genotypes are present and the `--estimate-AF -` option was used (#1687)
* bcftools +split-vep
- VEP fields with characters disallowed in VCF tag names by the specification (such as '-'
in 'M-CAP') couldn't be queried. This has been fixed, the program now sanitizes the field
names, replacing invalid characters with underscore (#1686)
## Release 1.15 (21st February 2022)
* New `bcftools head` subcommand for conveniently displaying the headers
of a VCF or BCF file. Without any options, this is equivalent to
`bcftools view --header-only --no-version` but more succinct and memorable.
* The `-T, --targets-file` option had the following bug originating in HTSlib code:
when an uncompressed file with multiple columns CHR,POS,REF was provided, the
REF would be interpreted as 0 gigabases (#1598)
Changes affecting specific commands:
* bcftools annotate
- In addition to `--rename-annots`, which requires a file with name mappings,
it is now possible to do the same on the command line `-c NEW_TAG:=OLD_TAG`
- Add new option --min-overlap which allows to specify the minimum required
overlap of intersecting regions
- Allow to transfer ALT from VCF with or without replacement using
bcftools annotate -a annots.vcf.gz -c ALT file.vcf.gz
bcftools annotate -a annots.vcf.gz -c +ALT file.vcf.gz
* bcftools convert
- Revamp of `--gensample`, `--hapsample` and `--haplegendsample` family of options
which includes the following changes:
- New `--3N6` option to output/input the new version of the .gen file format,
see https://www.cog-genomics.org/plink/2.0/formats#gen
- Deprecate the `--chrom` option in favor of `--3N6`. A simple `cut` command
can be used to convert from the new 3*M+6 column format to the format printed
with `--chrom` (`cut -d' ' -f1,3-`).
- The CHROM:POS_REF_ALT IDs which are used to detect strand swaps are required
and must appear either in the "SNP ID" column or the "rsID" column. The column
is autodetected for `--gensample2vcf`, can be the first or the second for
`--hapsample2vcf` (depending on whether the `--vcf-ids` option is given), must be
the first for `--haplegendsample2vcf`.
* bcftools csq
- Allow GFF files with phase column unset
* bcftools filter
- New `--mask`, `--mask-file` and `--mask-overlap` options to soft filter
variants in regions (#1635)
* bcftools +fixref
- The `-m id` option now works also for non-dbSNP ids, i.e. not just `rsINT`
- New `-m flip-all` mode for flipping all sites, including ambiguous A/T and C/G sites
* bcftools isec
- Prevent segfault on sites filtered with -i/-e in all files (#1632)
* bcftools mpileup
- More flexible read filtering using the options
--ls, --skip-all-set .. skip reads with all of the FLAG bits set
--ns, --skip-any-set .. skip reads with any of the FLAG bits set
--lu, --skip-all-unset .. skip reads with all of the FLAG bits unset
--nu, --skip-any-unset .. skip reads with any of the FLAG bits unset
The existing synonymous options will continue to function but their use
is discouraged
--rf, --incl-flags STR|INT Required flags: skip reads with mask bits unset
--ff, --excl-flags STR|INT Filter flags: skip reads with mask bits set
* bcftools query
- Make the `--samples` and `--samples-file` options work also in the `--list-samples`
mode. Add a new `--force-samples` option which allows to proceed even when some of
the requested samples are not present in the VCF (#1631)
* bcftools +setGT
- Fix a bug in `-t q -e EXPR` logic applied on FORMAT fields, sites with all
samples failing the expression EXPR were incorrectly skipped. This problem
affected only the use of `-e` logic, not the `-i` expressions (#1607)
* bcftools sort
- make use of the TMPDIR environment variable when defined
* bcftools +trio-dnm2
- The --use-NAIVE mode now also adds the de novo allele in FORMAT/VA
## Release 1.14 (22nd October 2021)
Changes affecting the whole of bcftools, or multiple commands:
* New `--regions-overlap` and `--targets-overlap` options which address
a long-standing design problem with subsetting VCF files by region.
BCFtools recognize two sets of options, one for streaming (`-t/-T`) and
one for index-gumping (`-r/-R`). They behave differently, the first
includes only records with POS coordinate within the regions, the other
includes overlapping regions. The two new options allow to modify the
default behavior, see the man page for more details.
* The `--output-type` option can be used to override the default compression
level
Changes affecting specific commands:
* bcftools annotate
- when `--set-id` and `--remove` are combined, `--set-id` cannot use
tags deleted by `--remove`. This is now detected and the program
exists with an informative error message instead of segfaulting
(#1540)
- while non-symbolic variation are uniquely identified by POS,REF,ALT,
symbolic alleles starting at the same position were undistinguishable.
This prevented correct matching of records with the same positions and
variant type but different length given by INFO/END (samtools/htslib@60977f2).
When annotating froma VCF/BCF, the matching is done automatically. When
annotating from a tab-delimited text file, this feature can be invoked
by using `-c INFO/END`.
- add a new '.' modifier to control wheter missing values should be carried
over from a tab-delimited file or not. For example:
-c TAG .. adds TAG if the source value is not missing. If TAG
exists in the target file, it will be overwritten
-c .TAG .. adds TAG even if the source value is missing. This
can overwrite non-missing values with a missing value
and can create empty VCF fields (`TAG=.`)
* bcftools +check-ploidy
- by default missing genotypes are not used when determining ploidy.
With the new option `-m, --use-missing` it is possible to use the
information carried in the missing and half-missing genotypes
(e.g. ".", "./." or "./1")
* bcftools concat:
- new `--ligate-force` and `--ligate-warn` options for finer control
of `-l, --ligate` behavior in imperfect overlaps. The new default is
to throw an error when sites present in one chunk but absent in the
other are encountered. To drop such sites and proceed, use the new
`--ligate-warn` option (previously this was the default). To keep such
sites, use the new `--ligate-force` option (#1567).
* bcftools consensus:
- Apply mask even when the VCF has no notion about the chromosome. It
was possible to encounter this problem when `contig` lines were not
present in the VCF header and no variants were called on that chromosome
(#1592)
* bcftools +contrast:
- support for chunking within map/reduce framework allowing to collect
NASSOC counts even for empty case/control sample sets (#1566)
* bcftools csq:
- bug fix, compound indels were not recognised in some cases (#1536)
- compound variants were incorrectly marked as 'inframe' even when
stop codon would occur before the frame was restored (#1551)
- bug fix, FORMAT/BCSQ bitmasks could have been assigned incorrectly
to some samples at multiallelic sites, a superset of the correct
consequences would have been set (#1539)
- bug fix, the upstream stop could be falsely assigned to all samples in
a multi-sample VCF even if the stop was relevant for a single sample
only (#1578)
- further improve the detection of mismatching chromosome naming
(e.g. "chrX" vs "X") in the GFF, VCF and fasta files
* bcftools merge:
- keep (sum) INFO/AN,AC values when merging VCFs with no samples (#1394)
* bcftools mpileup:
- new --indel-size option which allows to increase the maximum considered
indel size considered, large deletions in long read data are otherwise
lost.
* bcftools norm:
- atomization now supports Number=A,R string annotations (#1503)
- assign as many alternate alleles to genotypes at multiallelic sites
in the`-m +` mode, disregarding the phase. Previously the program
assumed to be executed as an inverse operation of `-m -`, but when
that was not the case, reference alleles would have been filled
instead of multiple alternate alleles (#1542)
* bcftools sort:
- increase accuracy of the --max-mem option limit, previously the limit
could be exceeded by more than 20% (#1576)
* bcftools +trio-dnm:
- new `--with-pAD` option to allow processing of VCFs without FORMAT/QS.
The existing `--ppl` option was changed to the analogous `--with-pPL`
* bcftools view:
- the functionality of the option --compression-level lost in 1.12
has been restored
## Release 1.13 (7th July 2021)
This release brings new options and significant changes in BAQ parametrization
in `bcftools mpileup`. The previous behavior can be triggered by providing
the `--config 1.12` option. Please see https://github.com/samtools/bcftools/pull/1474
for details.
Changes affecting the whole of bcftools, or multiple commands:
* Improved build system
Changes affecting specific commands:
* bcftools annotate:
- Fix rare a bug when INFO/END is present, all INFO fields are removed
with `bcftools annotate -x INFO` and BCF output is produced. Then the
removed INFO/END continues to inform the end coordinate and causes
incorrect retrieval of records with the -r option (#1483)
- Support for matching annotation line by ID, in addition to CHROM,POS,REF,
and ALT (#1461)
bcftools annotate -a annots.tab.gz -c CHROM,POS,~ID,REF,ALT,INFO/END input.vcf
* bcftools csq:
- When GFF and VCF/fasta use a different chromosome naming convention
(e.g. chrX vs X), no consequences would be added. Newly the program
attempts to detect these differences and remove/add the "chr" prefix
to chromosome name to match the GFF and VCF/fasta (#1507)
- Parametrize brief-predictions parameter to allow explicit number of
aminoacids to be printed. Note that the `-b, --brief-predictions` option
is being replaced with `-B, --trim-protein-seq INT`
* bcftools +fill-tags:
- Generalization and better support for custom functions that allow
adding new INFO tags based on arbitrary `-i, --include` type of
expressions. For example, to calculate a missing INFO/DP annotation
from FORMAT/AD, it is possible to use:
-t 'DP:1=int(sum(FORMAT/AD))'
Here the optional ":1" part specifies that a single value will be
added (by default Number=. is used) and the optional int(...) adds
an integer value (by default Type=Float is used).
- When FORMAT/GT is not present, the INFO/AF tag will be newly calculated
from INFO/AC and INFO/AN.
* bcftools gtcheck:
- Switch between FORMAT/GT or FORMAT/PL when one is (implicitly) requested
but only the other is available
- Improve diagnostics, printing warnings when a line cannot be matched and
the number of lines skipped for various reasons (#1444)
- Minor bug fix, with PLs being the default, the `--distinctive-sites` option
started to require explicit `--error-probability 0`
* bcftools index:
- The program now accepts both data file name and the index file name. This
adds to user convenience when running index statistics (-n, -s)
* bcftools isec:
- Always generate sites.txt with isec -p (#1462)
* bcftools +mendelian:
- Consider only complete trios, do not crash on sample name typos (#1520)
* bcftools mpileup:
- New `--seed` option for reproducibility of subsampling code in HTSlib
- The SCR annotation which shows the number of soft-clipped reads now
correctly pools reads together regardless of the variant type. Previously
only reads with indels were included at indel sites.
- Major revamp of BAQ. Please see https://github.com/samtools/bcftools/pull/1474
for details. The previous behavior can be triggered by providing the `--config 1.12`
option.
- Thanks to improvements in HTSlib, the removal of overlapping reads (which can
be disabled with the `-x, --ignore-overlaps` options) is not systematically biased
anymore (https://github.com/samtools/htslib/pull/1273)
- Modified scale of Mann-Whitney U tests. Newly INFO/*Z annotations will be printed,
for example MQBZ replaces MQB.
* bcftools norm:
- Fix Type=Flag output in `norm --atomize` (#1472)
- Atomization must not discard ALT=. records
- Atomization of AD and QS tags now correctly updates occurrences of duplicate
alleles within different haplotypes
- Fix a bug in atomization of Number=A,R tags
* bcftools reheader:
- Add `-T, --temp-prefix` option
* bcftools +setGT:
- A wider range of genotypes can be set by the plugin by allowing
specifying custom genotypes. For example, to force a heterozygous
genotype it is now possible to use expressions like:
c:'m|M'
c:0/1
c:0
* bcftools +split-vep:
- New `-u, --allow-undef-tags` option
- Better handling of ambiguous keys such as INFO/AF and CSQ/AD. The
`-p, --annot-prefix` option is now applied before doing anything else
which allows its use with `-f, --format` and `-c, --columns` options.
- Some consequence field names may not constitute a valid tag name, such
as "pos(1-based)". Newly field names are trimmed to exclude brackets.
* bcftools +tag2tag:
- New --QR-QA-to-QS option to convert annotations generated by Freebays
to QS used by BCFtools
* bcftools +trio-dnm:
- Add support for sites with more than four alleles. Note that only the
four most frequent alleles are considered, the model remains unchanged.
Previously such sites were skipped.
- New --use-NAIVE option for a naive DNM calling based solely on FORMAT/GT
and expected Mendelian inheritance. This option is suitable for prefiltering.
- Fix behavior to match the documentation, the `--dnm-tag DNG` option now
correctly outputs log scaled values by default, not phred scaled.
- Fix bug in VAF calculation, homozygous de novo variants were incorrectly
reported as having VAF=50%
- Fix arithmetic underflow which could lead to imprecise scores and improve
sensitivity in high coverage regions
- Allow combining --pn and --pns to set the noise trehsholds independently
## Release 1.12 (17th March 2021)
Changes affecting the whole of bcftools, or multiple commands:
* The output file type is determined from the output file name suffix, where
available, so the -O/--output-type option is often no longer necessary.
* Make F_MISSING in filtering expressions work for sites with multiple
ALT alleles (#1343)
* Fix N_PASS and F_PASS to behave according to expectation when reverse
logic is used (#1397). This fix has the side effect of `query` (or
programs like `+trio-stats`) behaving differently with these expressions,
operating now in site-oriented rather than sample-oriented mode. For
example, the new behavior could be:
bcftools query -f'[%POS %SAMPLE %GT\n]' -i'N_PASS(GT="alt")==1'
11 A 0/0
11 B 0/0
11 C 1/1
while previously the same expression would return:
11 C 1/1
The original mode can be mimicked by splitting the filtering into two steps:
bcftools view -i'N_PASS(GT="alt")==1' | \
bcftools query -f'[%POS %SAMPLE %GT\n]' -i'GT="alt"'
Changes affecting specific commands:
* bcftools annotate:
- New `--rename-annots` option to help fix broken VCFs (#1335)
- New -C option allows to read a long list of options from a file to
prevent very long command lines.
- New `append-missing` logic allows annotations to be added for each ALT
allele in the same order as they appear in the VCF. Note that this is
not bullet proof. In order for this to work:
- the annotation file must have one line per ALT allele
- fields must contain a single value as multiple values are appended
as they are and would break the correspondence between the alleles
and values
* bcftools concat:
- Do not phase genotypes by mistake if they are not already phased
with `-l` (#1346)
* bcftools consensus:
- New `--mask-with`, `--mark-del`, `--mark-ins`, `--mark-snv` options
(#1382, #1381, #1170)
- Symbolic <DEL> should have only one REF base. If there are multiple,
take POS+1 as the first deleted base.
- Make consensus work when the first base of the reference genome is
deleted. In this situation the VCF record has POS=1 and the first
REF base cannot precede the event. (#1330)
* bcftools +contrast:
- The NOVELGT annotation was previously not added when requested.
* bcftools convert:
- Make the --hapsample and --hapsample2vcf options consistent with each
other and with the documentation.
* bcftools call:
- Revamp of `call -G`, previously sample grouping by population was not
truly independent and could still be influenced by the presence of other
sample groups.
- Optional addition of INFO/PV4 annotation with `call -a INFO/PV4`
- Remove generation of useless HOB and ICB annotation;
use `+fill-tags -- -t HWE,ExcHet` instead
- The `call -f` option was renamed to `-a` to (1) make it consistent with
`mpileup` and (2) to indicate that it includes both INFO and FORMAT
annotations, not just FORMAT as previously
- Any sensible Number=R,Type=Integer annotation can be used with -G,
such as AD or QS
- Don't trim QUAL; although usefulness of this change is questionable for
true probabilistic interpretation (such high precision is unrealistic),
using QUAL as a score rather than probability is helpful and permits more
fine-grained filtering
- Fix a suspected bug in `call -F` in the worst case, for certain improve
readability
- `call -C trio` is temporarily disabled
* bcftools csq:
- Fix a bug wich caused incorrect FORMAT/BCSQ formatting at sites with too
many per-sample consequences
- Fix a bug which incorrectly handled the --ncsq parameter and could clash
with reserved BCF values, consequently producing truncated or even incorrect
output of the %TBCSQ formatting expression in `bcftools query`. To account
for the reserved values, the new default value is --ncsq 15 (#1428)
* bcftools +fill-tags:
- MAF definition revised for multiallelic sites, the second most common
allele is considered to be the minor allele (#1313)
- New FORMAT/VAF, VAF1 annotations to set the fraction of alternate reads
provided FORMAT/AD is present
* bcftools gtcheck:
- support matching of a single sample against all other samples in the file
with `-s qry:sample -s gt:-`. This was previously not possible, either
full cross-check mode had to be run or a list of pairs/samples had to
be created explicitly
* bcftools merge:
- Make `merge -R` behavior consistent with other commands and pull in
overlapping records with POS outside of the regions (#1374)
- Bug fix (#1353)
* bcftools mpileup:
- Add new optional tag `mpileup -a FORMAT/QS`
* bcftools norm:
- New `-a, --atomize` functionality to decompose complex variants,
for example MNVs into consecutive SNVs
- New option `--old-rec-tag` to indicate the original variant
* bcftools query:
- Incorrect fields were printed in the per-sample output when subset
of samples was requested via -s/-S and the order of samples in the
header was different from the requested -s/-S order (#1435)
* bcftools +prune:
- New options --random-seed and --nsites-per-win-mode (#1050)
* bcftools +split-vep:
- Transcript selection now works also on the raw CSQ/BCSQ annotation.
- Bug fix, samples were dropped on VCF input and VCF/BCF output (#1349)
* bcftools stats:
- Changes to QUAL and ts/tv plotting stats: avoid capping QUAL to
predefined bins, use an open-range logarithmic binning instead
- plot dual ts/tv stats: per quality bin and cumulative as if threshold
applied on the whole dataset
* bcftools +trio-dnm2:
- Major revamp of +trio-dnm plugin, which is now deprecated and replaced by
+trio-dnm2.
The original trio-dnm calling model used genotype likelihoods (PLs) as the
input for calling. However, that is flawed because PLs make assumptions
which are unsuitable for de novo calling: PL(RR) can become bigger than
PL(RA) even when the ALT allele is present in the parents. Note that
this is true also for other programs such as DeNovoGear which rely on
the same samtools calculation.
The new recommended workflow is
bcftools mpileup -a AD,QS -f ref.fa -Ou proband.bam father.bam mother.bam |
bcftools call -mv -Ou |
bcftools +trio-dnm -p proband,father,mother -Oz -o output.vcf.gz
This new version also implements the DeNovoGear model. The original
behavior of trio-dnm is no longer supported.
For more details see http://samtools.github.io/bcftools/trio-dnm.pdf
## Release 1.11 (22nd September 2020)
Changes affecting the whole of bcftools, or multiple commands:
* Filtering -i/-e expressions
- Breaking change in -i/-e expressions on the FILTER column. Originally
it was possible to query only a subset of filters, but not an exact match.
The new behavior is:
FILTER="A" .. exact match, for example "A;B" does not pass
FILTER!="A" .. exact match, for example "A;B" does pass
FILTER~"A" .. both "A" and "A;B" pass
FILTER!~"A" .. neither "A" nor "A;B" pass
- Fix in commutative comparison operators, in some cases reversing sides
would produce incorrect results (#1224; #1266)
- Better support for filtering on sample subsests
- Add SMPL_*/S* family of functions that evaluate within rather than across
all samples. (#1180)
* Improvements in the build system
Changes affecting specific commands:
* bcftools annotate:
- Previously it was not possible to use `--columns =TAG` with INFO tags
and the `--merge-logic` feature was restricted to tab files with BEG,END
columns, now extended to work also with REF,ALT.
- Make `annotate -TAG/+TAG` work also with FORMAT fields. (#1259)
- ID and FILTER can be transferred to INFO and ID can be populated from
INFO. However, the FILTER column still cannot be populated from an INFO
tag because all possible FILTER values must be known at the time of
writing the header (#947; #1187)
* bcftools consensus:
- Fix in handling symbolic deletions and overlapping variants.
(#1149; #1155; #1295)
- Fix `--iupac-codes` crash on REF-only positions with `ALT="."`. (#1273)
- Fix `--chain` crash. (#1245)
- Preserve the case of the genome reference. (#1150)
- Add new `-a, --absent` option which allows to set positions with no
supporting evidence to "N" (or any other character). (#848; #940)
* bcftools convert:
- The option `--vcf-ids` now works also with `-haplegendsample2vcf`. (#1217)
- New option `--keep-duplicates`
* bcftools csq:
- Add `misc/gff2gff.py` script for conversion between various flavors of
GFF files. The initial commit supports only one type and was contributed
by @flashton2003. (#530)
- Add missing consequence types. (PR #1203; #1292)
- Allow overlapping CDS to support ribosomal slippage. (#1208)
* bcftools +fill-tags:
- Added new annotations: INFO/END, TYPE, F_MISSING.
* bcftools filter:
- Make `--SnpGap` optionally filter also SNPs close to other variant types.
(#1126)
* bcftools gtcheck:
- Complete revamp of the command. The new version is faster and allows
N:M sample comparisons, not just 1:N or NxN comparisons.
Some functionality was lost (plotting and clustering) but may be added
back on popular demand.
* bcftools +mendelian:
- Revamp of user options, output VCFs with mendelian errors annotation,
read PED files (thanks to Giulio Genovese).
* bcftools merge:
- Update headers when appropriate with the '--info-rules *:join' INFO rule.
(#1282)
- Local alleles merging that produce LAA and LPL when requested, a draft
implementation of https://github.com/samtools/hts-specs/pull/434 (#1138)
- New `--no-index` which allows to merge unindexed files. Requires the input
files to have chromosomes in th same order and consistent with the order
of sequences in the header. (PR #1253; samtools/htslib#1089)
- Fixes in gVCF merging. (#1127; #1164)
* bcftools norm:
- Fixes in `--check-ref s` reference setting features with non-ACGT bases.
(#473; #1300)
- New `--keep-sum` switch to keep vector sum constant when splitting
multiallelics. (#360)
* bcftools +prune:
- Extend to allow annotating with various LD metrics: r^2,
Lewontin's D' (PMID:19433632), or Ragsdale's D (PMID:31697386).
* bcftools query:
- New `%N_PASS()` formatting expression to output the number of samples
that pass the filtering expression.
* bcftools reheader:
- Improved error reporting to prevent user mistakes. (#1288)
* bcftools roh:
- Several fixes and improvements
- the `--AF-file` description incorrectly suggested "REF\tALT" instead
of the correct "REF,ALT". (#1142)
- RG lines could have negative length. (#1144)
- new `--include-noalt` option to allow also ALT=. records. (#1137)
* bcftools scatter:
- New plugin intended as a convenient inverse to `concat`
(thanks to Giulio Genovese, PR #1249)
* bcftools +split:
- New `--groups-file` option for more flexibility of defining desired
output. (#1240)
- New `--hts-opts` option to reduce required memory by reusing one
output header and allow overriding the default hFile's block size
with `--hts-opts block_size=XXX`. On some file systems (lustre) the
default size can be 4M which becomes a problem when splitting files
with 10+ samples.
- Add support for multisample output and sample renaming
* bcftools +split-vep:
- Add default types (Integer, Float, String) for VEP subfields and make
`--columns -` extract all subfields into INFO tags in one go.
## Release 1.10.2 (19th December 2019)
This is a release fix that corrects minor inconsistencies discovered in
previous deliverables.
## Release 1.10 (6th December 2019)
* Numerous bug fixes, usability improvements and sanity checks were added
to prevent common user errors.
* The -r, --regions (and -R, --regions-file) option should never create
unsorted VCFs or duplicates records again. This also fixes rare cases where