代码之家 › 专栏 › 技术社区 › Jose Ricardo Bustos M.

如何将fasta文件中的序列(几行)重新格式化为单行?[副本]

fasta sed bash

Jose Ricardo Bustos M. · 技术社区 · 9 年前

输入“file.fasta”( 笔记 ,这是一个示例……在fasta文件中,序列可能有三行以上)

>chr1:117223140-117223856 TAG:GTGGG
GTGGgggggcgCATAGCAGCAGCAGCAGCAGCAGCAGCAGCAGCAGAGtt
aGTAGTATCGAATCGCACGACTGACAGCTCAGCATCAGCGACGACTAGTG
GTGGGCGACGACAgCGATATA
>chr2:117223140-117223856 TAG:GGGCT
ACGAGCAGCAGCAGCAGCagCCGATCGACGACTCAAGTACGATACGCGaa
cCCCCCGACGACGACTCACGA

预期产量

>chr1:117223140-117223856 TAG:GTGGG
GTGGgggggcgCATAGCAGCAGCAGCAGCAGCAGCAGCAGCAGCAGAGttaGTAGTATCGAATCGCACGACTGACAGCTCAGCATCAGCGACGACTAGTGGTGGGCGACGACAgCGATATA
>chr2:117223140-117223856 TAG:GGGCT
ACGAGCAGCAGCAGCAGCagCCGATCGACGACTCAAGTACGATACGCGaacCCCCCGACGACGACTCACGA

我的努力: sed 命令

sed ':a;N;$!ba;s/\([actgACGT]\)\n\([actgACGT]\)/\1\2/g' file.fasta

我的错误输出:

>chr1:117223140-117223856 TAG:GTGGGGTGGgggggcgCATAGCAGCAGCAGCAGCAGCAGCAGCAGCAGCAGAGttaGTAGTATCGAATCGCACGACTGACAGCTCAGCATCAGCGACGACTAGTGGTGGGCGACGACAgCGATATA
>chr2:117223140-117223856 TAG:GGGCTACGAGCAGCAGCAGCAGCagCCGATCGACGACTCAAGTACGATACGCGaacCCCCCGACGACGACTCACGA

标头的正则表达式(第一个字母为“>”的行)是 "^>.*$" ,但我不知道如何包含在 已安装的 命令

提前谢谢

2 回复 | 直到 9 年前

potong 9 年前

这可能适用于您(GNU sed):

sed ':a;N;/>/!s/\n//;ta;P;D' file

查看两行,如果其中一行不包含 > 删除它们之间的换行符并重复。如果其中一行包含 > 然后打印并删除其中的第一个,然后重复。

John1024 David C. Rankin 9 年前

$ awk '/^>/ {print (NR>1?"\n":"")$0;; next} {printf "%s",$0;} END{print "";}' file.fasta 
>chr1:117223140-117223856 TAG:GTGGG
GTGGgggggcgCATAGCAGCAGCAGCAGCAGCAGCAGCAGCAGCAGAGttaGTAGTATCGAATCGCACGACTGACAGCTCAGCATCAGCGACGACTAGTGGTGGGCGACGACAgCGATATA
>chr2:117223140-117223856 TAG:GGGCT
ACGAGCAGCAGCAGCAGCagCCGATCGACGACTCAAGTACGATACGCGaacCCCCCGACGACGACTCACGA

它的工作原理

/^>/ {print (NR>1?"\n":"")$0;; next}

如果该行以 > ,即如果正则表达式 /^>/ 如果为真,则打印该行。如果这不是第一行,即 NR>1 ,然后在该行前面打印一个换行符。然后,跳过其余命令并跳转到 next 线
printf "%s",$0;

对于所有其他行,打印时不要使用换行符。
END{print "";}

到达文件末尾后,打印最后一个换行符。

推荐文章

user2954003 · AWK使用正则表达式匹配字符串并与前一个字符串组合

2 年前

Giancarlo D · 在IPv4地址匹配后使用SED删除行尾的冒号

2 年前

John Smith · 在特定行的末尾添加文本

2 年前

Code With Banchi · 在sed命令中使用变量-sed-e异常:“s”的选项未知

2 年前

Waddah Shamroukh · 使用sed将字符串替换为变量

2 年前

DAG · bash使用sed在文件的第一个位置插入单词

2 年前

Aravinth Kalai · 如何使用Linux shell命令[duplicate]对两列求和并将值保存到第三列

2 年前

NicolÃ¡s Abuchar Moya · 如何在Linux中从给定字符串中删除奇数索引字符[closed]

2 年前

nickcrv06 · 在两个常量字符串之间提取单词

2 年前

DevOps_Engg · 使用sed提取路径值子字符串

2 年前