Linux中fasta文件的拆分與合並
FASTA文件的拆分:
(1)如果從一個文件a提取第11至20個序列存到另一個文件b:
awk -v RS='>' 'NR>1{i++}i>=10&&i<=21{print ">"$0}' a.fasta|sed '/^$/d'>b.fasta
(2)將某一文件a中每一條序列保存到一個文件中:
awk '/^>/{f=++d".fasta"} {print > f}' input.fasta
FASTA文件合並:
cat *.fasta > output.fasta
cat `ls -tr *.fasta` > output.fasta #按時間反序的順序合並文件
- 轉載自:http://www.yelinsky.com/blog/archives/380.html