Exemplo n.º 1
0
def doHeader(fd):
  vcf = Vcf(fd)
  vcf.load_meta_header()
  ot_already_there = str(vcf.check_info('OT'))
  if ot_already_there == True:
    error("This vcf seems to have an OT INFO already. Bailing out.")
  vcf.add_info('OT', '0', 'Flag', 'The site is on target.')
  # print "Checking if INFO id=OT is there: " + str(vcf.check_info('OT'))
  print vcf.get_meta()
  return vcf
Exemplo n.º 2
0
def doHeader(fd):
    vcf = Vcf(fd)
    vcf.load_meta_header()
    ot_already_there = str(vcf.check_info('OT'))
    if ot_already_there == True:
        error("This vcf seems to have an OT INFO already. Bailing out.")
    vcf.add_info('OT', '0', 'Flag', 'The site is on target.')
    # print "Checking if INFO id=OT is there: " + str(vcf.check_info('OT'))
    print vcf.get_meta()
    return vcf
Exemplo n.º 3
0
def doHeader(fd):
  vcf = Vcf(fd)
  vcf.load_meta_header()
  ot_already_there = str(vcf.check_info('RDP'))
  if ot_already_there == True:
    error("This vcf seems to have an RDP INFO already. Bailing out.")
  vcf.add_info('RDP', '1', 'Integer', 'Raw read coverage at locus.')
  # print "Checking if INFO id=OT is there: " + str(vcf.check_info('OT'))
  print vcf.get_meta()
  return vcf