
note that the 't' should start with a $
while IFS=$'t' read -r gene old_locus new_locus contig start end strand
do
efetch -db=nuccore -format=fasta -id=$contig -seq_start=$start -seq_stop=$end -strand=$strand > "$gene"__"$new_locus".fasta
done < VF_table




近期评论