导读

很多程序输出的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文件的换行符相关推荐

  1. python 写文件 换行_python中写入txt文件需要换行,以及\r 和\n

    原博文 2019-09-25 16:16 − 在Python中,用open()函数打开一个txt文件,写入一行数据之后需要一个换行 如果直接用 f.write('\n')只会在后面打印一个字符串'\n ...

  2. python 写txt 换行_python中写入txt文件需要换行,以及\r 和\n

    在Python中,用open()函数打开一个txt文件,写入一行数据之后需要一个换行 如果直接用 f.write('\n') 只会在后面打印一个字符串'\n',而不是换行' 需要用 f.write(' ...

  3. python处理fasta文件_Python脚本:fasta文件单序列信息提取

    使用Python对fasta格式的序列进行基本信息统计 预期设计输出文件中包括fasta文件名,序列长度,GC含量以及ATCG各自的含量. Python脚本编辑 使用的文件 test.fasta st ...

  4. python文本提取序列信息_从fasta文件中通过头中的ID号提取序列

    accessionids.txt是否只包含四位数代码?在 如果是,请将accessorID更改为:accessorID = accessorIDWithArrow[1:5] 一些方法可以让这更像Pyt ...

  5. python创建空文本文件_Python干货:「文件处理整合」 创建、打开、写入和删除...

    什么是Python文件处理? Python文件处理就是如何在python中创建.打开.读取.写入和删除文件的示例.从文件(EXCELSheet.doc文件.文本文件.cv等)或Add中读取数据.是一个 ...

  6. python服务端编程_Python实现的简单文件传输服务器和客户端

    还是那个题目(题目和流程见java版本),感觉光用java写一点新意也没有,恰巧刚学习了python,何不拿来一用,呵呵: 服务器端: import SocketServer, time class ...

  7. python控制台清屏_Python Shell 怎样清屏?

    启动Python有两种方式,分别为"Windows命令行窗口"和"IDLE" "命令行窗口"下可以通过如下两种方法: 1. import s ...

  8. python mysql日期输出格式_Python/Shell/MySQL时间获取与格式转换

    一.说明 时间的获取及时间各格式间的转换是比较常用的操作,但一是多种语言经常容易弄混,二是同一种语言同一个功能可能有不同的实现函数,导致每次处理时间经常要百度所以来记录一下. 另外个人真不喜欢同样功能 ...

  9. python去重复元素_Python实现去除列表中重复元素的方法总结【7种方法】

    这里首先给出来我很早之前写的一篇博客,Python实现去除列表中重复元素的方法小结[4种方法],感兴趣的话可以去看看,今天是在实践过程中又积累了一些方法,这里一并总结放在这里. 由于内容很简单,就不再 ...

最新文章

  1. 阿里巴巴的五大平台野心,让“连接”论成为过去式
  2. Android AutoFocusCallback is not being called or not returning
  3. Java开发WebService(使用Java-WS)
  4. HashMap和Hashtable的区别 2
  5. 「PowerBI」使用TabularEditor进行PowerBIDeskTop模型开发最佳实践
  6. Loj #6089. 小 Y 的背包计数问题
  7. 牛客网 【每日一题】4月15日 Treepath
  8. Shell数组操作(添加、删除元素、清空、遍历)
  9. H3C S1526交换机端口镜像配置
  10. Java Web学习总结(12)——使用Session防止表单重复提交
  11. 聪明的ITPRO之二IT人做事要“圆”
  12. python is beautiful_python自动化报告BeautifulReport用法
  13. FeelYourSound Chillout Engine Pro for Mac - MIDI文件生成插件
  14. 28.卷1(套接字联网API)---原始套接字
  15. 安全工具系列 -- 信息收集(二)
  16. lisp visual 开根号_分享Visual Lisp程序第二波|文件的读与写-lsp文件
  17. oracle定时任务按照小时,ORACLE定时任务时间间隔设置
  18. ubuntu安装WPS字体缺失的解决办法
  19. 分享电脑中截图的五种方法(包括截长图)
  20. Linux 返回根目录

热门文章

  1. Chrome浏览器官方离线安装包下载
  2. centos llvm安装_在CentOS上编译安装llvm-3.8.1详细教程
  3. ac9560不支持承载_IPFS 最新进展:0.5 版本改进较大,已出现系统级支持 | 火星号精选...
  4. gns3中两个路由器分别连接主机然后分析ip数据转发报文arp协议_TCP/IP协议知识总结...
  5. tf.Session().as_default的作用
  6. 逆向建模软件介绍_逆向技术在2代机修理中的应用
  7. python简单的爬虫教程中召唤小海龟_Python练习,网络小爬虫(初级)
  8. 【AWS】DynamoDB扫描操作获取表全部数据
  9. android p新功能_android 11今天发布这是新功能
  10. 打包无法识别lombok