首页 > 其他分享 >批量修改Fasta文件中序列的名称

批量修改Fasta文件中序列的名称

时间:2023-10-30 11:02:09浏览次数:33  
标签:rename pysam name 批量 名称 序列 Fasta

比如一个Fasta文件的内容如下:

seq001|aaa
ATCGGGG
seq002|bbb
AAAATTTT

删除序列名称中“|”后的内容,只保留seq001, seq002这样的名称

点击查看代码
#!/usr/bin/env python
import sys
import pysam

with pysam.FastxFile(sys.argv[1]) as fh:
    for r in fh:
        new_name = r.name.split('|')[0]
        print(">"+new_name)
        print(r.sequence)

python rename_fasta.py test.fa > rename.fa

标签:rename,pysam,name,批量,名称,序列,Fasta
From: https://www.cnblogs.com/-bonjour/p/17797276.html

相关文章

  • 关于使用ChatGPT的API开发批量处理工具的总结
    关于使用ChatGPT的API开发批量处理工具的总结由于批量处理,输入和输出最好使用json格式,并在提示词中要求以相同格式返回;以便于工具的解析和处理;返回的结果中要求chatgpt,说明其详细的处理过程及原因,方便对提示词进行调优,否则,面对问题时,只能瞎猜;提示词中尽量说明什么情况要怎么做......
  • Prufer序列 学习笔记
    2023.10.29晚,在随机做AtCoder的时候见到了[ABC303Ex]ConstrainedTreeDegree。然后开始思考DP,未果后看题解,发现是Prufer序列->尝试学习Prufer序列。什么是Prufer序列Prufer序列是一种将带标号的树用一个唯一的整数序列表示的方法,是解决树计数问题的工具。给一棵有根树......
  • python vtk读取dicom序列+鼠标键盘交互
    目标:vtk+pyqt实现四视图。之前不了解vtk,也不了解鼠标键盘交互。网上搜索了资料,发现博客里大都是C++的例子。困扰几天,今天终于做出来一部分,分享一下。参考官方教程:examples.vtk.org/site/Python/IO/ReadDICOM/examples.vtk.org/site/Python/IO/ReadDICOMSeries/第一步:py......
  • session伪造+pickle反序列化
    来自[CISCN2019华北赛区Day1Web2]ikun小黑子真没有树脂!!!进去就是坤哥的舞蹈视频,意思就是让我们找到LV6的号然后买了,翻了好几页没找到:url上有个page参数,直接写个python脚本帮我找:importrequestsimportre#页面URL和GET参数base_url="http://ec000824-34eb-4762-8......
  • 使用FastAPI部署Ultralytics YOLOv5模型
    YOLO是YouOnlyLookOnce(你只看一次)的缩写,它具有识别图像中的物体的非凡能力,在日常应用中会经常被使用。所以在本文中,我们将介绍如何使用FastAPI的集成YOLOv5,这样我们可以将YOLOv5做为API对外提供服务。Python有几个web框架,其中最突出的两个是Flask和FastAPI。Flask是一个轻......
  • 用友反序列化漏洞综合
    用友NC序列化漏洞综合利用工具分享下载地址https://github.com/wgpsec/YongYouNcTooljava11以上的环境可使用集成了以下漏洞BshServletrcejsInvokerceDeleteServletcc6反序列化DownloadServletcc6反序列化FileReceiveServletcc6反序列化fsDownloadServletc......
  • PowerShell批量改名
    $files=Get-Location|Get-ChildItem$current=Get-LocationWrite-Output$current.Pathfor($i=0;$i-lt$files.Length;$i++){$file=$files[$i]if($file.FullName.EndsWith(".png")){Write-Output"--------------......
  • 最长上升子序列
    importjava.util.*;publicclassSolution{/***代码中的类名、方法名、参数名已经指定,请勿修改,直接返回方法规定的值即可**给定数组的最长严格上升子序列的长度。*@paramarrint整型一维数组给定的数组*@returnint整型*/......
  • 2、关于网络中接受的数据如何序列化和反序列化的思考以及实现
    1、背景介绍因工作接触到半导体行业,主要负责EAP相关的东西,其中需要实现SECS/GEM协议,消息协议使用的是SECS-II,其中有一种数据类型是A类型,表示字符串类型。需要将接收到的SECS指令记录在日志中,以及反解析SECS指令。我们知道,网络中接受到的数据都是byte,需要自己根据规......
  • php-SER-libs-main反序列化靶场部分wp
    什么是序列化/反序列化,为什么要进行序列化/反序列化序列化:有时需要把一个对象在网络上传输,为了方便传输,可以把整个对象转化为二进制串,等到达另一端时,再还原为原来的对象,这个过程称之为串行化(也叫序列化)。反序列化:将序列化的结果还原PHP序列化:把对象转化为二进制的字符......