#1 seqret "/P/y21/term3/pr9/to_join/*" full.fasta #2 seqretsplit "/P/y21/term3/pr9/to_split/coding2.fasta" #4 transeq "/P/y21/term3/pr9/coding.fasta" prot.fasta -frame 1 -table 0 #10 shuffleseq "/P/y21/term3/pr9/coding.fasta" output.fasta #12 cusp "/P/y21/term3/pr9/coding.fasta" freq.txt