@Override
protected void doWork(String inputSource, VcfIterator r, VariantContextWriter w)
throws IOException {
VCFHeader header = r.getHeader();
VCFHeader h2 = new VCFHeader(header.getMetaDataInInputOrder(), header.getSampleNamesInOrder());
h2.addMetaDataLine(
new VCFInfoHeaderLine(
TAG,
VCFHeaderLineCount.UNBOUNDED,
VCFHeaderLineType.String,
"metadata added from " + TABIX + " . Format was " + FORMAT));
h2.addMetaDataLine(
new VCFHeaderLine(
getClass().getSimpleName() + "CmdLine", String.valueOf(getProgramCommandLine())));
h2.addMetaDataLine(
new VCFHeaderLine(getClass().getSimpleName() + "Version", String.valueOf(getVersion())));
h2.addMetaDataLine(
new VCFHeaderLine(
getClass().getSimpleName() + "HtsJdkVersion", HtsjdkVersion.getVersion()));
h2.addMetaDataLine(
new VCFHeaderLine(getClass().getSimpleName() + "HtsJdkHome", HtsjdkVersion.getHome()));
SAMSequenceDictionaryProgress progress = new SAMSequenceDictionaryProgress(header);
w.writeHeader(h2);
while (r.hasNext()) {
VariantContext ctx = progress.watch(r.next());
Set<String> annotations = new HashSet<String>();
CloseableIterator<BedLine> iter =
this.bedReader.iterator(ctx.getContig(), ctx.getStart() - 1, ctx.getEnd() + 1);
while (iter.hasNext()) {
BedLine bedLine = iter.next();
if (!ctx.getContig().equals(bedLine.getContig())) continue;
if (ctx.getStart() - 1 >= bedLine.getEnd()) continue;
if (ctx.getEnd() - 1 < bedLine.getStart()) continue;
String newannot = this.parsedFormat.toString(bedLine);
if (!newannot.isEmpty()) annotations.add(VCFUtils.escapeInfoField(newannot));
}
CloserUtil.close(iter);
if (annotations.isEmpty()) {
w.add(ctx);
continue;
}
VariantContextBuilder vcb = new VariantContextBuilder(ctx);
vcb.attribute(TAG, annotations.toArray());
w.add(vcb.make());
incrVariantCount();
if (checkOutputError()) break;
}
progress.finish();
}
public String toString(BedLine tokens) {
String s = tokens.get(index);
if (s == null) s = "";
return s + (next == null ? "" : next.toString(tokens));
}