Python:如何使用二进制值编码DNA序列? [英] Python: How to encode DNA sequence using binary values?

问题描述

我想将包含少量DNA序列的文件转换为二进制值，如下所示:

I would like to convert a file that contained few DNA sequences into binary values which is as follow:

A=1000
C=0100
G=0010
T=0001

FileA.txt

CCGAT
GCTTA

所需的输出

01000100001010000001
00100100000100011000

我尝试使用此代码解决问题，但是bin输出文件似乎无法输出所需的答案.谁能帮我吗?

I have tried using this code to solve my problem but the bin output file seem failed to output my desired answer. Can anyone help me?

代码

import sys

if len(sys.argv) != 2 :
  sys.stderr.write('Usage: {} <nucleotide file>\n'.format(sys.argv[0]))
  sys.exit()

# assumes the file only contains dna and newlines
sequence = ''
for line in open(sys.argv[1]) :
  sequence += line.strip().upper()

sequence = sequence.replace('A', chr(0b1000))
sequence = sequence.replace('C', chr(0b0100))
sequence = sequence.replace('G', chr(0b0010))
sequence = sequence.replace('T', chr(0b0001))

outfile = open(sys.argv[1] + '.bin', 'wb')

outfile.write(bytearray(sequence, encoding = 'utf-8'))

推荐答案

您要使用ascii输出还是二进制?下面将为您提供您在帖子中显示的内容(尽管在一行上.需要修改代码以保留换行符).

Do you want ascii output or binary? The below will give you what you show in your post (though on a single line. Code needs to be modified to keep newlines).

import sys

if len(sys.argv) != 2 :
  sys.stderr.write('Usage: {} <nucleotide file>\n'.format(sys.argv[0]))
  sys.exit()

# assumes the file only contains dna and newlines
sequence = ''
for line in open(sys.argv[1]) :
  sequence += line.strip().upper()

sequence = sequence.replace('A', '1000')
sequence = sequence.replace('C', '0100')
sequence = sequence.replace('G', '0010')
sequence = sequence.replace('T', '0001')

outfile = open(sys.argv[1] + '.bin', 'wb')

outfile.write(sequence)

编辑，这将创建一个二进制文件，其中每个核苷酸均为一个字节，换行符以二进制格式保存.

EDIT This creates a binary file where each nucleotide is a byte and the newlines are preserved in binary format.

import sys

if len(sys.argv) != 2 :
  sys.stderr.write('Usage: {} <nucleotide file>\n'.format(sys.argv[0]))
  sys.exit()

# assumes the file only contains dna and newlines
newbytearray=bytearray(b'',encoding='utf-8')
dict={'A':0b1000,'C':0b0100,'G':0b0010,'T':0b0001,'\n':0b1010}
with open(sys.argv[1]) as file:
    while True:
        char=file.read(1)
        if not char:
            file.close()
            break
        newbytearray.append(dict[char])
outfile = open(sys.argv[1] + '.bin', 'wb')
outfile.write(newbytearray)
outfile.close()

#Converts the binary file to unicode and prints the result sequence.
testBin = open('fileA.txt.bin','rb')
sequence=''
for line in testBin:
    line = line.replace(chr(0b1000),'1000')
    line = line.replace(chr(0b0100),'0100')
    line = line.replace(chr(0b0010),'0010')
    line = line.replace(chr(0b0001),'0001')
    line = line.replace(chr(0b1010),'\n')
    sequence += line
#outputVerify = open('outputVerify.txt','wb')
#outputVerify.write(sequence)
#outputVerify.close()
print sequence
testBin.close()

#Shows the data of the binary file. Note that byte 6 is the newline character 0b1010.
testBin = open('fileA.txt.bin','rb')
list = ''
i=0
while True:
    b = testBin.read(1)
    i += 1
    if not b:
    break #due to eof
    list += b
    print 'byte: ' + str(i) + ' is '+ '{0:04b}'.format(ord(b)) +' and has decimal representation: ' + str(ord(b))
testBin.close()

这篇关于Python:如何使用二进制值编码DNA序列?的文章就介绍到这了，希望我们推荐的答案对大家有所帮助，也希望大家多多支持IT屋！