【无标题】

awk -F '[' '/^>/ { sub(/\..*/, "",$5); print $5 } !/^>/ { print }' GCF_000001405.39_GRCh38.p13_rna_from_genomic.fna

linux中提取包含指定字符的列 

awk '{for(i=1;i<=NF;i++) if(index($i,"transcript_id")) print $i}' GCF_000001405.39_GRCh38.p13_rna_from_genomic.fna

 linux中去掉每行=之前的字符串

sed 's/^[^=]*=//' process.fna > process_1.fna

去掉末尾的]

sed -i 's/[]]$//' process_1.fna

每一行加>

sed -i 's/^/>/' process_1.fna

以transcript_id为分隔符打印,再把=换成> (去掉点后了)

awk -F 'transcript_id' '/^>/ { sub(/\..*/, "",$2); print $2 } !/^>/ { print }' GCF_000001405.39_GRCh38.p13_rna_from_genomic.fna > process_4.fna
sed -i 's/=/>/g' process_4.fna

 保留点

awk -F 'transcript_id' '/^>/ { sub(/\.*/, "",$2); print $2 } !/^>/ { print }' GCF_000001405.39_GRCh38.p13_rna_from_genomic.fna > process_5.fna

#去掉每行第一个]后所有内容
sed -i 's/].*//' process_5.fna

  • 1
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值