導讀

很多程序輸出的fasta序列每60個堿基會換一次行，下面是一個去除堿基序列后的換行符的方法。

一、輸入

cat test

# 一個隨便寫的文件

>sequence_1

AAAA

GGGG

CCCC

TTTT

>sequence_2

aaaa

gggg

cccc

tttt

二、去換行符：python

腳本：python3

#!/usr/bin/env python3

import os

import sys

import re

ms, infile, outfile = sys.argv

with open(infile) as f:

Dict = {}

for line in f:

if line[0] == ">":

key = line.strip()

Dict[key] = []

else:

Dict[key].append(line.strip())

with open(outfile, 'w') as o:

for key, value in Dict.items():

o.write("{}\n{}\n".format(key, ''.join(value)))

運行：python3

python3 trim_enter.py test test2

結果：python3

>sequence_1

AAAAGGGGCCCCTTTT

>sequence_2

aaaaggggcccctttt

三、去換行符：shell

運行：shell單行命令

awk '!/^>/{printf "%s", $0; n="\n"}/^>/{print n $0; n=""}END{printf "%s", n}' test > test2

在不以>開頭的行中，打印不帶換行符的行，并存儲一個換行字符(在變量n中)以供以后使用。

在以>開頭的行中，打印存儲的換行字符(如果有的話)和該行。

重置n，以防這是最后一行。如果需要，以換行結束。

結果：shell

cat test2

# 隨便看看結果

>sequence_1

AAAAGGGGCCCCTTTT

>sequence_2

aaaaggggcccctttt

總結

以上是生活随笔為你收集整理的python shell怎么换行_python shell：去除fasta文件的换行符的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。