code: def bwasw(args, opts):
      
      dbfile, readfile = args
      dbfile = check_index(dbfile)
  
      samfile, _, unmapped = get_samfile(readfile, dbfile,
                                         bam=opts.bam, unmapped=opts.unmapped)
      if not need_update(dbfile, samfile):
          logging.error("`{0}` exists. `bwa bwasw` already run.".format(samfile))
          return "", samfile
  
      cmd = "bwa bwasw " + " ".join(args)
      cmd += " -t {0}".format(opts.cpus)
      cmd += " " + opts.extra
      return cmd, samfile
docs: %prog bwasw database.fasta long_read.fastq
  
      Wrapper for `bwa bwasw`. Output will be long_read.sam.