Hello,

you could first use the -name option. And than change the id with a little python script using Biopython like this:

from Bio import SeqIO

bed_file = "./regions.bed"
original_file = "./original.fasta"
renamed_file = "./renamed.fasta"

regions = dict()
records = []

with open(bed_file, "r") as bed:
    for line in bed:
        data = line.strip().split("\t")
        regions[data[3]] = data

for record in SeqIO.parse(original_file, 'fasta'):
    record.id = "{3}:{0}:{1}-{2}({5})".format(*regions[record.id])
    record.description = record.id
    records.append(record)

SeqIO.write(records, renamed_file, 'fasta')

fin swimmer

EDIT: Code correction. The original has overwritten the renamed.fasta in every iteration.

Alternative to bedtools, you can use a subseq command of seqkit, which help fetching particular sequence region and also provides all the information in the header.

You could use samtools faidx with bed2faidxsta.pl:

$ bed2faidxsta.pl < in.bed > out.fa

If you want the header to have a custom format, just edit line 67:

my $header = ">".join(":",($chr, $start, $stop, $id, $score, $strand));

You could move $id to the front of the list, for instance.