code: def map_to_precursors(seqs, names, loci, out_file, args):
      
      with make_temp_directory() as temp:
          pre_fasta = os.path.join(temp, "pre.fa")
          seqs_fasta = os.path.join(temp, "seqs.fa")
          out_sam = os.path.join(temp, "out.sam")
          pre_fasta = get_loci_fasta(loci, pre_fasta, args.ref)
          out_precursor_file = out_file.replace("tsv", "fa")
          seqs_fasta = get_seqs_fasta(seqs, names, seqs_fasta)
  
          # print(open(pre_fasta).read().split("\n")[1])
          if find_cmd("razers3"):
              cmd = "razers3 -dr 2 -i 80 -rr 90 -f -o {out_sam} {temp}/pre.fa  {seqs_fasta}"
              run(cmd.format(**locals()))
              out_file = read_alignment(out_sam, loci, seqs, out_file)
              shutil.copy(pre_fasta, out_precursor_file)
      return out_file
docs: map sequences to precursors with razers3