求助，将基因组fasta格式的每条序列切割成100kb的片段，相邻片段间以50kb重叠

primary assemblies were transformed into very long
overlapping sequences with a maximum of 100 kb (50 kb overlap)

写回答
好问题 0 提建议
追加酬金
关注问题
分享
邀请回答
编辑收藏删除结题
收藏举报

2条回答默认最新

关注

码龄粉丝数原力等级 --

被采纳

被点赞

采纳率
Lawfe 2023-04-04 09:25
关注
你好，你有这个代码了吗？能分享一下吗？

解决无用
评论打赏
分享
举报

评论

按下Enter换行，Ctrl+Enter发表内容

报告相同问题？

关注问题

怎么用python修改fasta序列的ID名字 python
2022-04-19 08:29

回答 1 已采纳 import re def run(): fasta_str = """WP_018731760.1 NAD(P)-binding domain-containing protein [Sal
比较两个fasta序列的碱基个数并提取ATG及之后的序列 python 有问必答
2021-05-19 15:51

回答 2 已采纳参考一下：https://blog.csdn.net/david10251025/article/details/102093869
怎么用Python将Excel转化为fasta？ python 有问必答
2022-04-11 12:04

回答 3 已采纳在github上有处理代码，可参考一下：https://github.com/kuod/cctf/blob/master/parser.py
如何使用MUMmer比对大片段序列
2021-04-08 11:20

徐洲更hoptop的博客如何使用MUMmer比对大片段序列 ...但是1999年后，越来越多的物种全基因组出现，比如说在1999年出现了_Helicobacter pylori_的第二类菌株的基因组序列，就需要研究同一物种不同品系进化过程的基因...
Python获取基因库序列并排序 python 有问必答
2021-05-29 00:51

回答 2 已采纳 list index out of range 错误出现的原因主要有两个，一个可能是下标超出范围，一个可能是list是空的，没有一个元素比如 a = random.randint(0,
re.sub()对A文件中每个以>开头的字符串，在B文件查找并以该字符串在B所在列名替换A中的该字符串 python 大数据有问必答正则表达式
2021-10-26 10:35

回答 1 已采纳是不是类似这样操作，循环替换。 import re s=['gene>ggcctt>tggc','gene>ctct'] s1=[] for x in s: x1=re.s
如何实现读取多个文件夹下的fasta文件，大脑缺氧，求救命！ python
2019-07-03 21:26

回答 2 已采纳使用os.listdir('目录路径') 列表形式返回该目录下的所有文件名，然后你再依次读取和处理就ok了
宏基因组组装神器-MEGAHIT使用及常见问题
2021-02-01 19:05

Neptuneyut的博客 std::bad_alloc/Exit code -6 从中断点继续运行输出结果预设参数对内存需求样本实际组装时间参考简介宏基因组测序获得海量短片段测序数据，这些数据混合着环境中各种各样的微生物基因组序列，如何恢复出这些...
初学python，批量选取文件中某一部分内容，然后将内容逐条存于一个文件中 python
2019-06-27 20:46

回答 2 已采纳 ``` # -*- coding:utf-8 -*- import os from Bio import SeqIO # root_dir为要读取文件的根目录 root_dir =
sed查找与替换具体代码搞不懂 linux
2022-04-19 12:13

回答 2 已采纳 [0-9]匹配数字的意思。+表示上一个表达式出现一次以上。
shell或python脚本批量处理txt文件 bash linux python 有问必答
2022-03-20 04:01

回答 6 已采纳 # -*- utf8 -*- """ pyCode : test.py Created by: Taro Description: 2022/3/20 0020 """ import os impor
病毒组学数据分析 -03 CheckV病毒序列识别
2022-02-21 22:19

道曼曼的博客 CheckV 是一个完全自动化的命令行管道，用于评估单重叠群病毒基因组的质量，包括识别集成原病毒的宿主污染、估计基因组片段的完整性以及识别封闭基因组。
如何用python写一个循环，或者其他方式，完成批量选取内容并保存。 python
2019-06-27 21:38

回答 3 已采纳 from Bio import SeqIO for xx in range(1, 68): xn = "split__" + str(xx) + ".fasta" yn = "selec
NAR：查询未培养病毒基因组的综合生态和进化框架IMG/VR v3
2022-01-04 07:00

刘永鑫Adam的博客 IMG/VR v3：一个用于查询未培养病毒基因组的综合生态和进化框架IMG/VR v3: an integrated ecological and evolutionary framewo...
这个处理不同基因组区域关系的工具集很不错！
2021-08-03 22:00

生信宝典的博客 Bedtools是处理基因组信息分析的强大工具集合，本文列出自己学习其官方文档的几个点，对后面计算不同样品peak相似性的脚本做了下更新和调整，使用起来更为简单方便。内容摘要区域注释，如p...
宏转录组组装软件rnaSPAdes输出文件解读
2024-01-18 11:07

WDPLA的博客 rnaSPAdes作为宏转录组数据组装软件，数据结果表现优秀——共输出5个文件夹及14个单独的文件，transcripts.fasta是主要结果文件。
CONTING EXPRESS
2015-01-31 21:27

4. **拼接**："CONTING EXPRESS"就在此环节发挥作用，它通过比对结果，利用算法寻找相邻片段之间的重叠部分，将这些片段连接成更长的连续序列，即contig。软件可能采用了如De Bruijn图、Overlap-Layer-Graph等方法。...
High-Quality Genome-Scale Models From Error-Prone, Long-Read Assemblies高质量的基因组尺度模型来自易出错的，长时间读取的程序集
2021-01-24 14:08

wangchuang2017的博客高质量的基因组尺度模型来自易出错的，长时间读取的程序集误差杆，长读装配的高质量基因组规模模型贾里德T. Broddrick1 *†，理查德舒宾2，查尔斯J. Norsigian2，乔纳森M.僧2，哈德O.保尔松2和玛丽Parenteau的N....
seqkit安装与使用 v2.5.1（生物信息学工具-003）
2024-02-23 23:56

让学习成为一种生活方式的博客 SeqKit - 跨平台且极速的FASTA/Q文件操作工具包，生信胶水之一。
生信知识与软件小结：BWA、wtdbg
2019-10-05 14:57

BaldStrong的博客三代组装软件:canu、smartdenovo、wtdbg 三代矫正软件:Racon、Nanopolish Falcon：一个实验性的二倍体组装工具，测试...blast ：的全称是 basic local alignment search tool 是一种极其常见的序列比对工具。其...
没有解决我的问题, 去提问

悬赏问题

¥15 metadata提取的PDF元数据，如何转换为一个Excel
¥15 关于arduino编程toCharArray()函数的使用
¥100 vc++混合CEF采用CLR方式编译报错
¥15 coze 的插件输入飞书多维表格 app_token 后一直显示错误，如何解决？
¥15 vite+vue3+plyr播放本地public文件夹下视频无法加载
¥15 c#逐行读取txt文本，但是每一行里面数据之间空格数量不同
¥50 如何openEuler 22.03上安装配置drbd
¥20 ING91680C BLE5.3 芯片怎么实现串口收发数据
¥15 无线连接树莓派，无法执行update，如何解决？（相关搜索：软件下载）
¥15 Windows11, backspace, enter, space键失灵

求助，将基因组fasta格式的每条序列切割成100kb的片段，相邻片段间以50kb重叠

2条回答 默认 最新

悬赏问题

2条回答默认最新