[Bioc-devel] VariantAnnotation: verbose output of readVcf
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