Разделение заголовка и содержимого с помощью регулярных выражений - PullRequest
0 голосов
/ 14 февраля 2019

У меня есть следующий текст последовательности в файле

каждый заголовок начинается с ">" количество строк является случайным

>XM_024446048.1 PREDICTED: Homo sapiens mannosidase alpha class 2A member 1 (MAN2A1), transcript variant X2, mRNA
CAGCCCCC
TGAGCGAC
TCCCTAATGTG
ACAGTAAAGAA
>NM_001308028.1 Homo sapiens FER tyrosine kinase (FER), transcript variant 2, mRNA
CAGCCC
CCGTGACGC
GGGGTGGTGACT
GGCTC
GGTGGT
GTGAC
>NM_0013082323028.1 H STZ mRSN1A
CAGCCC
CCGTGACGC
GGG
GTGGTGA
CTGGCTCCGGAGT
CTGAGGGGTTCGG

Я хочу создать вложенный dict в следующемформат:


nested_dict = { 'sequence1': {

'header': 'XM_024446048.1 PREDICTED: Homo sapiens mannosidase alpha class 2A member 1 (MAN2A1), transcript variant X2, mRNA',
'content': 'CAGCCCCCTGAGCGACTCCCTAATGTGACAGTAAAGAA'


},
               'sequence2': {

'header': 'NM_001308028.1 Homo sapiens FER tyrosine kinase (FER), transcript variant 2, mRNA',
'content': 'CAGCCCCCGTGACGCGGGGTGGTGACTGGCTCGGTGGTGTGAC'


},

               'sequence3': {

'header': 'NM_0013082323028.1 H STZ mRSN1A',
'content': 'CAGCCCCCGTGACGCGGGGTGGTGACTGGCTCCGGAGTCTGAGGGGTTCGG'


}

Я застрял с выражениями регулярных выражений, кто-нибудь может мне помочь?

1 Ответ

0 голосов
/ 14 февраля 2019

Вы можете использовать

^>
(?P<header>
[^\n()]*
    (?:\((?P<name>[^\n()]+)\))?
[^()\n]*
)
[\n\r]
(?P<content>[\s\S]+?)(?=^>|\Z)

с флагом multiline и verbose, см. демонстрацию на regex101.com .


In Python:
rx = re.compile(r'''
    ^>
    (?P<header>
    [^\n()]*
        (?:\((?P<name>[^\n()]+)\))?
    [^()\n]*
    )
    [\n\r]
    (?P<content>[\s\S]+?)(?=^>|\Z)''', re.M | re.X)

sequences = {"sequence{}".format(idx):
                 {"header": m.group('header'),
                  "content": m.group("content"),
                  "name": m.group('name')
                  }
             for idx, m in enumerate(rx.finditer(data), 1)}
print(sequences)

Это дает

{'sequence1': {'header': 'XM_024446048.1 PREDICTED: Homo sapiens mannosidase alpha class 2A member 1 (MAN2A1), transcript variant X2, mRNA', 'content': 'CAGCCCCC\nTGAGCGAC\nTCCCTAATGTG\nACAGTAAAGAA\n', 'name': 'MAN2A1'}, 'sequence2': {'header': 'NM_001308028.1 Homo sapiens FER tyrosine kinase (FER), transcript variant 2, mRNA', 'content': 'CAGCCC\nCCGTGACGC\nGGGGTGGTGACT\nGGCTC\nGGTGGT\nGTGAC\n', 'name': 'FER'}, 'sequence3': {'header': 'NM_0013082323028.1 H STZ mRSN1A', 'content': 'CAGCCC\nCCGTGACGC\nGGG\nGTGGTGA\nCTGGCTCCGGAGT\nCTGAGGGGTTCGG', 'name': None}}

См. демо на ideone.com .

Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...