login1$ module load bioperl login1$ bp_genbank2gff3.pl input.gbk
bp_seqconvert --from genbank --to fasta < input.gbk > output.fasta
samtools index input.fasta