使用python脚本结合Excel筛选功能对BSA结果进行候选基因筛选

本文介绍如何利用Python将BSA测序文件的SNPindex值与基因注释结合,通过Excel筛选功能找出候选基因。通过pd.merge进行数据合并,并设定筛选条件:SNPindex_mt=1且deltaSNPindex>=0.5,同时排除特定变异类型和同义突变,筛选后保存为m_12.xlsx。

拟解决问题-对BSA返回结果数据进行筛选获得候选基因

BSA返回的测序文件,将两个文件的SNPindex值与基因注释结合起来,这样就能直接使用Excel的筛选功能找到候选基因ID。两个.xlsx文件,比如Delta-SNPindex.xlsx和exonic_variant_function.xlsx

pd.merge(df1, df2, how = ‘inner’, on = [‘col1’, ’col2’ ‘col3’], indicator = True )

使用代码如下

该处使用的url网络请求的数据

import pandas as pd

#导入excel文件
filepath1 = "E:\yB\Delta_SNPindex.xlsx"
filepath2 = "E:\yB\exonic_variant_function.xlsx"
df1 = pd.read_excel(filepath1, sheet_name 
评论
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值