너도 AB1파일가지고 pdf 만들수 있어 (1)

국내 회사에서 일반 시퀀싱 즉 Sanger Sequencing을 맡기시면

시퀀싱되어 바이너리 정보로 저장되어 있는 ab1파일, 서열정보가 들어있는 fasta파일과 함께 Applied Biosystems 로고가 찍혀있는 pdf파일 혹은 시퀀싱을 요청한 회사의 로고가 찍혀있는 pdf파일을 받아 보실 수 있으실겁니다. 

자세히 보다보시면 서열과 plot의 위치가 맞는거 같으면서도 안맞는거 같기도하고..

대놓고 틀린건 아닌데 그렇다고 안대놓고 맞는건 아닌거같고..

그러면 이 파일은 대체 어떻게 만들어지는걸까?

그럼 나는 저런거 못만드나?

그럴리가 있겠습니까?

그래서 이 글의 제목이 "너도 AB1파일로 pdf 만들 수 있어" 아니겠습니까?

예전에는 ab1.py? ab1_parse.py? 어떤 용자께서 제작하신 3rd party파일을 python에 import하여 사용하였었는데

이번에는 biopython에서 제공되는 ab1파일을 핸들링 하는 기능을 사용해서 작업해보려고 합니다.

ab1파일을 parsing하려면 ab1파일이 어떻게 생겼는지 알아야 하겠죠?

그래서 ab1파일에 어떤 항목에 어떤 내용들이 담겨져 있는지 확인을 해보시기 바랍니다.

>ABIF_File_Format.pdf<

다음번에는 ab1파일안에 들어있는 내용들을 확인해보는 시간을 가지도록 하겠습니다.

진짜로~

제발~

출처: @sana_twice.09

AB1 파일 작업하기

출처: interactive-biosoftware

Python이 설치된 폴더의 Lib 폴더안에 ABIFReader.py를
복사하여 넣어주고 사용하면 됨.

from ABIFReader import *


# creates an instance of ABIFReader
reader = ABIFReader(<filename>)
# version of ABIF file
reader.version
# print all entries of ABIF file "<name> (<num>) / <type> (<size>)"
reader.showEntries()
# read data for entry named <name> with number <num>, by default <num> is 1
data = reader.getData(<name>[, <num>])
# close the file, since it is kept open
reader.close()

본 파일을 이용하여 ABI3730xl에서 나온 ab1파일을 분석하는 경우 중간중간
몇몇 data type은 3500xl 버전을 참고해야되는 경우가 있었음.
ABIF Type의 Description이 잘 맞지 않은 것인지 AB1파일을 로드할때 중간에 missing되는 것이 있는지 정확히는 모르겠음.