[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