code: def mitocompile(args):
      
      from jcvi.formats.vcf import VcfLine
      from six.moves.urllib.parse import parse_qsl
  
      p = OptionParser(mitocompile.__doc__)
      opts, args = p.parse_args(args)
  
      if len(args) < 1:
          sys.exit(not p.print_help())
  
      vcfs = args
      print("\t".join("vcf samplekey depth seqid pos alt svlen pe sr".split()))
      for i, vcf in enumerate(vcfs):
          if (i + 1) % 100 == 0:
              logging.debug("Process `{}` [{}]".
                            format(vcf, percentage(i + 1, len(vcfs))))
          depthfile = vcf.replace(".sv.vcf.gz", ".depth")
          fp = must_open(depthfile)
          chrm, depth = fp.next().split()
          depth = int(float(depth))
          samplekey = op.basename(vcf).split("_")[0]
  
          fp = must_open(vcf)
          for row in fp:
              if row[0] == '#':
                  continue
              v = VcfLine(row)
              info = dict(parse_qsl(v.info))
              print("\t".join(str(x) for x in (vcf,
                                               samplekey, depth,
                                               v.seqid, v.pos, v.alt,
                                               info.get("SVLEN"), info["PE"], info["SR"])))
docs: %prog mitcompile *.vcf.gz
  
      Extract information about deletions in vcf file.