获取指定位置序列

最新推荐文章于 2022-12-03 20:29:54 发布
最新推荐文章于 2022-12-03 20:29:54 发布
阅读量3.1k 收藏 5
点赞数
分类专栏： 生物信息 python 文章标签： 基因组-生物信息学 python
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/msw521sg/article/details/52562259
版权
本文介绍如何在基因组序列中通过Python获取指定位置的反向互补序列，重点在于理解负号表示反向序列以及start和end位置在反向互补后的变化。
摘要由CSDN通过智能技术生成

负号表示反向序列，并且start and end 的位置是序列反向互补之后的位置

Name	strand	start	end	total_length
UN000226	-2	1088	1411	1431
UN073473	3	861	1082	1158
UN082299	-1	838	1014	1105
UN064320	2	791	1006	1150
UN070736	-3	780	1271	1293
UN070422	2	776	955	1361
UN022140	2	758	1135	1235
UN215439	1	718	891	1044

代码

# encoding: utf-8
from Bio import SeqIO

record_dict = SeqIO.index
最低0.47元/天 解锁文章
msw521sg
关注
专栏目录
java获取机器序列号_Java获取主板序列号、MAC地址、CPU序列号工具类
weixin_33846508的博客
12-23 206
importjava.io.BufferedReader;importjava.io.IOException;importjava.io.InputStreamReader;importjava.net.InetAddress;importjava.net.NetworkInterface;importjava.util.ArrayList;importjava.util.Enumeration;...
linux提取指定行至指定位置
weixin_30916125的博客
11-08 916
grep查找ERROR，定位位置 awk打印到指定行数 sed打印到文本末尾 awk打印到文本末尾 方法一 #!/bin/csh -f if(-f errorlog.rpt) then rm -rf errorlog.rpt endif ls *.log > loglst.lst #将log文件写到指定文件 #loop execute set n=`wc -l ...
根据位置信息提取 fasta 文件中的序列 -- extract fasta sequence by their position
weixin_30552635的博客
08-08 1633
#!/usr/bin/env python # usages: python extract_seq_by_pos.py input.fasta id_start_end > result.fasta import sys import re FASTA= sys.argv[1] BED= sys.argv[2] fasta= open(FASTA, 'U') fasta_dict= {...
根据 基因名、bed 文件的基因位置，提取 DNA 序列 bedtools
dujidan的博客
03-19 3896
根据 基因名、bed 文件的基因位置，提取 DNA 序列 bedtools1、根据 Gene Symbol 查找在序列上的位置2、根据 基因位置 提取参考上的序列 1、根据 Gene Symbol 查找在序列上的位置 从UCSC下载匹配的 文件 链接: http://genome.ucsc.edu/cgi-bin/hgTables. 1）、主要注意版本的信息，按默认就行 2）、点击 get output，至下一页，就可以选择自己勾选自己想要的信息了 3）、下载后的文件，三列分别是 Transcripti
已知基因染色体号，坐标如何快速获取基因序列
llhgj的博客
12-03 2836
bedtools提取序列
django自带serializers序列化返回指定字段的方法
09-18
- `format`：指定序列化格式，例如'json'、'xml'等。 - `queryset`：要序列化的QuerySet对象。 - `fields`（可选）：一个包含要保留的字段名的元组。如果提供此参数，只有这些字段会被序列化。 - `use_natural_...
oracle 重置序列从指定数字开始的方法详解
12-15
以下是一些关于如何重置Oracle序列从指定数字开始的方法。 1. **不删除重建方式** 在不删除重建序列的情况下，可以使用PL/SQL动态SQL语句来实现重置。首先，获取序列的当前`nextval`，然后根据目标起始值计算出...
Django REST Framework序列化外键获取外键的值方法
09-18
这种方法涉及到在序列化器中创建新的字段，并通过 `source` 参数指定该字段应获取的数据源。例如，我们可以在 `QuestionSerializer` 中添加一个名为 `questionnaire_title` 的字段，用于获取 `Questionnaire` 模型的...
NCBI获取指定区域基因序列及其引物设计
aagk6406的博客
01-08 2280
获取指定基因区间序列链接：https://www.ncbi.nlm.nih.gov/nuccore/NC_000017.11 步骤一：进入NCBI主页，左边框中选择Nucleotide，并且在右边框中输入相应的搜索信息，点击Search，如下图所示： 　　　　https://www.ncbi.nlm.nih.gov/ 　　　　 步骤二：出现如下图，点击第一个搜索结果 ...
从fasta中根据序列坐标提取序列的四种方法（seqkit、seqtk、bedtools、gffread）
wanwan的博客
11-09 1万+
链接: 可以从fasta中提取基因序列的4款软件. 看到很好的有关根据坐标信息提取基因序列的总结，mark一下
edger和deseq2_edgeR-DESeq2分析RNA-seq差异表达
weixin_42393272的博客
12-31 722
转录组分析工具edgeR-DESeq2差异表达分析edgeR 包的安装∙edgeR 包是基于Bioconductor平台发布的，所以安装不能直接用install.packages()命令从 CRAN 上来下载∙安装:# try http:// if https:// URLs are not supported>source("https://http://doc.xuehai.net/bi...
提取参考基因组某个位置的碱基
whenfree的博客
12-27 1万+
引言 参考基因组是现代遗传学的基石，更是生物信息学研究的基础。不论你研究的是DNA、RNA，还是蛋白，你都无法避开参考基因组而进行你的研究。在生物信息分析过程中，参考基因组最常用的一个场景就是序列比对---将测序的read比对回参考基因组上。那么，如果我们想知道参考基因组某个位置或者某个区域的碱基序列，我们应该如何快速的解决这个问题呢？本文将介绍三种解决方案。 参考基因组格式 参考基因组通常...
利用python处理dna序列_科学网-简单的Python脚本提取对应位置基因序列（fasta文件）-王彬忠的博文...
weixin_39845430的博客
11-30 1232
最近，用Python脚本提取，在基因号已知，位置已知条件下，相对应位置的基因序列时发现，这样很简单但是很实用的脚本，在网上却比较难找。而且，能被找到的脚本，相对于具有初级编程能力的人而言，有点难。本人写了相对于初学者同样很简单脚本分享给大家。首先，我将fa文件处理为单行（嫌麻烦，没有写成scaffold_x一行，序列一行的样子，如图三），将下面的序列处理（图一）：（补充）经过：import ref...
python爬取genek视频_Python 脚本提取基因组指定位置序列
weixin_39785081的博客
12-09 822
说明在基因组分析中，我们经常会有这么一个需求，就是在一个fasta文件中提取一些序列出来。有时这些序列是一段完整的序列，而有时仅仅为原fasta文件中某段序列的一部分。特别是当数据量很多时，使用肉眼去挑选序列会很吃力，那么这时我们就可以通过简单的编程去实现了。例如此处在网盘附件中给定了某物种的全基因组序列(0-refer/ Bacillus_subtilis.str168.fasta)，及其基因组...
根据SNP的位置从基因组提取上下游序列
weixin_30586257的博客
09-04 2449
代码如下： #!/usr/bin/perl -w use strict; die "perl $0 <vcf> <genome>" if(@ARGV == 0); #Author:yueyao@genomics.cn my $vcf=shift; my $genome=shift; my%hash; my $id;...
如何根据染色体坐标快速得到基因组的 DNA 序列
klcola的专栏
01-18 1万+
http://pythonhosted.org/twobitreader/ 提供了一个方便的小工具 python -m twobitreader hg19.2bit < example.bed 染色体的位置信息在 bed 文件中给出，.2bit 文件格式是 UCSC Genome Browser 的基因组序列文件索引格式，可以在 http://hgdownload.soe.ucsc.edu...
python 获取list特定元素下标
热门推荐
qq_24737639的博客
12-19 15万+
在平时开发过程中，经常遇到需要在数据中获取特定的元素的信息，如到达目的地最近的车站，橱窗里面最贵的物品等等。怎么办？看下面 方法一： 利用数组自身的特性 a.index(target), 其中a是你的目标list，target是你需要的下标对应的值 a=[72, 56, 76, 84, 80, 88] print(a.index(76)) output: 2但是，如果a中有多个76呢？ 我们
mysql 8.0+获取序列函数
最新发布
06-06
MySQL 8.0及以上版本并没有内置的序列函数，但可以通过自定义函数实现类似的功能。以下是一个简单的自定义函数实现序列的示例： 1. 创建自定义函数 ``` DELIMITER $$ CREATE FUNCTION seq_nextval(seq_name VARCHAR(50)) RETURNS BIGINT BEGIN DECLARE seq_val BIGINT DEFAULT 0; SET seq_val = (SELECT seq_val FROM sequence WHERE seq_name = seq_name FOR UPDATE); IF seq_val IS NULL THEN SET seq_val = 1; INSERT INTO sequence (seq_name, seq_val) VALUES (seq_name, seq_val); ELSE SET seq_val = seq_val + 1; UPDATE sequence SET seq_val = seq_val WHERE seq_name = seq_name; END IF; RETURN seq_val; END$$ DELIMITER ; ``` 上述代码中，我们创建了一个名为 seq_nextval 的自定义函数，它接受一个参数 seq_name，表示序列的名称。函数的作用是获取指定名称的序列的下一个值。 函数首先从名为 sequence 的表中获取指定名称的序列的当前值，如果该序列不存在，则将其初始化为1，并返回1。如果序列存在，则将其值加1，并更新 sequence 表中的值，最后返回更新后的值。 2. 创建序列表 在自定义函数中，我们用到了一个名为 sequence 的表，它用于存储所有序列的当前值。我们需要先创建该表： ``` CREATE TABLE sequence ( seq_name VARCHAR(50) PRIMARY KEY, seq_val BIGINT ); ``` 3. 使用自定义函数获取序列的下一个值 现在我们可以使用 seq_nextval 函数获取指定名称的序列的下一个值了： ``` SELECT seq_nextval('my_seq'); ``` 上述代码中，我们调用 seq_nextval 函数，并传入名称为 my_seq 的序列。函数会返回该序列的下一个值。 需要注意的是，由于自定义函数中用到了事务和行级锁，因此在高并发场景下可能会出现性能问题。如果需要生成大量的序列，建议使用自增主键或其他更高效的方法。
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值