fasta文件拆分與合並


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


免責聲明!

本站轉載的文章為個人學習借鑒使用,本站對版權不負任何法律責任。如果侵犯了您的隱私權益,請聯系本站郵箱yoyou2525@163.com刪除。



 
粵ICP備18138465號   © 2018-2025 CODEPRJ.COM