gpt4 book ai didi

将 SNP 映射到 CDS 的 Python 脚本

转载 作者:行者123 更新时间:2023-11-30 22:39:00 25 4
gpt4 key购买 nike

我正在尝试制作一个程序来识别 CDS包含SNP 。它从两个文件填充两个字典,一个包含 SNP,另一个包含 GFF3文件。从 GFF3 文件填充的字典之一包含 CDS 名称及其作为元组的位置。

一个例子:

cds_pos = {'PRELSG_0019500_6': ('2091320', '2092988'), 'PRELSG_1338600_3': ('1542760','1542853'), 'PRELSG_0013000_1': ('1275531', '1275568')}

另一个字典包含染色体和 SNP 的位置,例如。

chrom_pos = {'PRELSG_13_v1': {'272093', '964287', '844454', '65770', '336211', '36660'}, 'PRELSG_12_v1': {'1270177', '1368630'}}

我的想法是迭代这两个字典并进行成对比较,看看是否在 CDS 的区间内找到 SNP 位置。我尝试了下面的代码,但它似乎不起作用。

for chrom, snp_pos in chrom_pos.items():
for cds, pos in cds_pos.items():
if pos[0] <= str(snp_pos) <= pos[1]:
print(cds)
print(snp_pos)

对于我发现不起作用的事情。首先,没有任何东西满足间隔 if 语句。其次,由于 SNP 的位置可以在多个染色体上找到,因此需要考虑到这一点,我用 chrom ==gene 语句尝试了这一点。但这似乎不起作用。

如果有任何想法和评论,我们将很高兴继续进行。谢谢

编辑:

到目前为止,我的脚本如下所示:

cds_snp = defaultdict(set)
for chrom, snp_pos in chrom_pos.items():
for cds, pos in cds_pos.items():
ed_chrom = chrom[:9]
ed_cds = cds[:9]
if ed_cds == ed_chrom:
for i in snp_pos: # Iterate through the set of snp positions
if int(pos[0]) <= int(i) <= int(pos[1]):
cds_snp[cds].add(i)

for i,j in sorted(cds_snp.items()):
print(i)
print('\n'.join(j))

我必须找到一种方法来评估输出是否正确,但它似乎是合理的。

最佳答案

为了让此代码正常工作,需要解决一些问题:

since the position of the SNP can be found on several chromosomes, this needs to be taken into account, which I tried with the chrom == gene statement.

您正在尝试将 SNP 与 CDS 中的位置关联起来。但是,您似乎没有 CDS 位置的染色体(只有名称)。从您发布的示例中,您的 CDS 名称为 'PRELSG_0019500_6', 'PRELSG_1338600_3', 'PRELSG_0013000_1'你的染色体名称是 'PRELSG_10_v1', 'PRELSG_12_v1' 。这些匹配的情况都不存在,说实话,它们看起来有不同的格式,并且您永远不会遇到 chrom == gene 的情况。

CDS 名称中是否有一些识别信息可以告诉您它位于哪条染色体上?如果有,您可以从染色体名称中提取染色体编号(即“PRELSG_12_v1”中的 12),并将其与从 CDS 名称中提取的染色体编号进行比较。

nothing satisfy the interval if statement

我认为您指的是 if pos[0] <= str(snp_pos) <= pos[1]:

对于为什么这不起作用有一个简单的解释。当您第一次从字典中提取特征时:

for chrom, snp_pos in chrom_pos.items():
for gene, pos in gene_pos.items():

您提取snp_pos 。这不是一个单独的立场,而是一组立场。要迭代各个位置,您可以添加另一个循环:

for snpPos in snp_pos:

最后,为了正确性起见,最好将间隔语句保留为整数。即写成

if int(pos[0]) <= int(snpPos) <= int(pos[1]): 

关于将 SNP 映射到 CDS 的 Python 脚本,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/43325602/

25 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com