tbl2fasta

A way to transform tab delimited to FASTA

#!/bin/bash
awk ‘{
seq = $2
l = length(seq)
i = 1
printf ">%s\n", $1
while (i <= l){
  printf "%s\n", substr(seq,i,80)
  i=i+80
}
}’ "$@"

Tags: AWK, bioinformatics