python shell怎么换行_python shell:去除fasta文件的换行符
导读
很多程序输出的fasta序列每60个碱基会换一次行,下面是一个去除碱基序列后的换行符的方法。
一、输入
cat test
# 一个随便写的文件
>sequence_1
AAAA
GGGG
CCCC
TTTT
>sequence_2
aaaa
gggg
cccc
tttt
二、去换行符:python
脚本:python3
#!/usr/bin/env python3
import os
import sys
import re
ms, infile, outfile = sys.argv
with open(infile) as f:
Dict = {}
for line in f:
if line[0] == ">":
key = line.strip()
Dict[key] = []
else:
Dict[key].append(line.strip())
with open(outfile, 'w') as o:
for key, value in Dict.items():
o.write("{}\n{}\n".format(key, ''.join(value)))
运行:python3
python3 trim_enter.py test test2
结果:python3
>sequence_1
AAAAGGGGCCCCTTTT
>sequence_2
aaaaggggcccctttt
三、去换行符:shell
运行:shell单行命令
awk '!/^>/{printf "%s", $0; n="\n"}/^>/{print n $0; n=""}END{printf "%s", n}' test > test2
在不以>开头的行中,打印不带换行符的行,并存储一个换行字符(在变量n中)以供以后使用。
在以>开头的行中,打印存储的换行字符(如果有的话)和该行。
重置n,以防这是最后一行。如果需要,以换行结束。
结果:shell
cat test2
# 随便看看结果
>sequence_1
AAAAGGGGCCCCTTTT
>sequence_2
aaaaggggcccctttt
python shell怎么换行_python shell:去除fasta文件的换行符相关推荐
- python 写文件 换行_python中写入txt文件需要换行,以及\r 和\n
原博文 2019-09-25 16:16 − 在Python中,用open()函数打开一个txt文件,写入一行数据之后需要一个换行 如果直接用 f.write('\n')只会在后面打印一个字符串'\n ...
- python 写txt 换行_python中写入txt文件需要换行,以及\r 和\n
在Python中,用open()函数打开一个txt文件,写入一行数据之后需要一个换行 如果直接用 f.write('\n') 只会在后面打印一个字符串'\n',而不是换行' 需要用 f.write(' ...
- python处理fasta文件_Python脚本:fasta文件单序列信息提取
使用Python对fasta格式的序列进行基本信息统计 预期设计输出文件中包括fasta文件名,序列长度,GC含量以及ATCG各自的含量. Python脚本编辑 使用的文件 test.fasta st ...
- python文本提取序列信息_从fasta文件中通过头中的ID号提取序列
accessionids.txt是否只包含四位数代码?在 如果是,请将accessorID更改为:accessorID = accessorIDWithArrow[1:5] 一些方法可以让这更像Pyt ...
- python创建空文本文件_Python干货:「文件处理整合」 创建、打开、写入和删除...
什么是Python文件处理? Python文件处理就是如何在python中创建.打开.读取.写入和删除文件的示例.从文件(EXCELSheet.doc文件.文本文件.cv等)或Add中读取数据.是一个 ...
- python服务端编程_Python实现的简单文件传输服务器和客户端
还是那个题目(题目和流程见java版本),感觉光用java写一点新意也没有,恰巧刚学习了python,何不拿来一用,呵呵: 服务器端: import SocketServer, time class ...
- python控制台清屏_Python Shell 怎样清屏?
启动Python有两种方式,分别为"Windows命令行窗口"和"IDLE" "命令行窗口"下可以通过如下两种方法: 1. import s ...
- python mysql日期输出格式_Python/Shell/MySQL时间获取与格式转换
一.说明 时间的获取及时间各格式间的转换是比较常用的操作,但一是多种语言经常容易弄混,二是同一种语言同一个功能可能有不同的实现函数,导致每次处理时间经常要百度所以来记录一下. 另外个人真不喜欢同样功能 ...
- python去重复元素_Python实现去除列表中重复元素的方法总结【7种方法】
这里首先给出来我很早之前写的一篇博客,Python实现去除列表中重复元素的方法小结[4种方法],感兴趣的话可以去看看,今天是在实践过程中又积累了一些方法,这里一并总结放在这里. 由于内容很简单,就不再 ...
最新文章
- 阿里巴巴的五大平台野心,让“连接”论成为过去式
- Android AutoFocusCallback is not being called or not returning
- Java开发WebService(使用Java-WS)
- HashMap和Hashtable的区别 2
- 「PowerBI」使用TabularEditor进行PowerBIDeskTop模型开发最佳实践
- Loj #6089. 小 Y 的背包计数问题
- 牛客网 【每日一题】4月15日 Treepath
- Shell数组操作(添加、删除元素、清空、遍历)
- H3C S1526交换机端口镜像配置
- Java Web学习总结(12)——使用Session防止表单重复提交
- 聪明的ITPRO之二IT人做事要“圆”
- python is beautiful_python自动化报告BeautifulReport用法
- FeelYourSound Chillout Engine Pro for Mac - MIDI文件生成插件
- 28.卷1(套接字联网API)---原始套接字
- 安全工具系列 -- 信息收集(二)
- lisp visual 开根号_分享Visual Lisp程序第二波|文件的读与写-lsp文件
- oracle定时任务按照小时,ORACLE定时任务时间间隔设置
- ubuntu安装WPS字体缺失的解决办法
- 分享电脑中截图的五种方法(包括截长图)
- Linux 返回根目录
热门文章
- Chrome浏览器官方离线安装包下载
- centos llvm安装_在CentOS上编译安装llvm-3.8.1详细教程
- ac9560不支持承载_IPFS 最新进展:0.5 版本改进较大,已出现系统级支持 | 火星号精选...
- gns3中两个路由器分别连接主机然后分析ip数据转发报文arp协议_TCP/IP协议知识总结...
- tf.Session().as_default的作用
- 逆向建模软件介绍_逆向技术在2代机修理中的应用
- python简单的爬虫教程中召唤小海龟_Python练习,网络小爬虫(初级)
- 【AWS】DynamoDB扫描操作获取表全部数据
- android p新功能_android 11今天发布这是新功能
- 打包无法识别lombok