awk -v RS='>' 'NR>1{i++}i>=10&&i<=20{print ">"$0}' input.fasta | sed '/^$/d' > output.fasta
# 输出文件名是1.fasta, 2.fasta, 3.fasta...
awk '/^>/{f=++d".fasta"} {print > f}' input.fasta
# 输出文件名是序列ID
awk '/^>/ {fout=sprintf("%s.fasta",substr($0,2));}{print >> fout;}' input.fasta
awk -F'>' 'NR==FNR{ids[$0]; next} NF>1{f=($2 in ids)} f' id.txt seq.fasta
从一个文件(input.fasta)提取第10至20个序列存到另一个文件(output.fasta) 将某一文件(...
将gbk格式文件转换成fasta格式文件 将多个fasta序列中提取含有'Homo sapiens'的序列 请关注...
fastq 转fasta文件： awk '{if(NR%4 == 1){print ">" substr($0, ...
被我用来处理不换行的fasta 文件，很早以前就写了；拿出来分享，perl oneline 和awk 似乎通用，毕...
1.切割染色体fasta： awk '/^>/ { gsub(">","",$1) FILE=$1 "....
在浏览核酸蛋白质数据库的时候会经常遇见不同的文件格式，常见的有Fasta格式文件、NBRF/PIR格式文件、 EM...
在3d-dna流程中，凡是涉及到edit-fasta-according-to-new-cprops.awk的流程...
之间处理fasta或者fastq时总是自己进行写脚本，比如：fasta-fastaq的转化，根据名称调取fasta...
fasta和fastq格式文件的shell小练习 1 统计reads_1.fq 文件中共有多少条序列信息 2 输出...
awk处理过程: 依次对每一行进行处理，然后输出awk命令形式:awk [-F|-f|-v] ‘BEGIN{} /...
本文标题：awk处理FASTA格式文件
本文链接：https://www.haomeiwen.com/subject/izjructx.html
网友评论