[Bioc-devel] VariantAnnotation: verbose output of readVcf
Julian Gehring
julian.gehring at embl.de
Fri Jul 17 11:02:20 CEST 2015
Hi Valerie,
Thank you for those changes!
Best wishes
Julianl
On Thu, Jul 16, 2015 at 17:12, Valerie Obenchain wrote:
Hi Julian,
Yes, the behavior is intentional though I hadn't thought of the annoying
chatter in the case of chunking. Sorry about that.
The current readVcf() reads/parses fields according to header
multiplicity and type; fields without headers are skipped. It is on the
TODO to be more liberal in reading (especially FORMAT) fields.
Quite a number of people are (a) unaware of header lines and (b) have
vcf files with incomplete headers. It's confusing to those with
incomplete headers why all fields aren't read in. So, printing the
"found" fields was an attempt to communicate which would be read /
parsed by readVcf().
I've made the following changes in 1.15.21:
- added a 'Header lines' section to the man page to explain this further
- added a 'verbose' arg to readVcf(); when TRUE fields found in header
are printed
Valerie
On 07/16/2015 02:35 AM, Julian Gehring wrote:
Hi,
In recent versions of 'VariantAnnotation', the 'readVcf' function prints
information about the header lines:
fl
[[alternative HTML version deleted]]
More information about the Bioc-devel
mailing list