python 计算 fastq数据中的GC 含量
001、
(base) [b20223040323@admin1 test]$ ls SRR1770413_1.fastq SRR1770413_2.fastq test.py (base) [b20223040323@admin1 test]$ cat test.py #!/bin/env python # -*- coding:utf-8 -*- from Bio import SeqIO fq1 = list(SeqIO.parse("SRR1770413_1.fastq", "fastq")) fq2 = list(SeqIO.parse("SRR1770413_2.fastq", "fastq")) total_gc = 0 total_bases = 0 for i in fq1: i = str(i.seq) total_gc += i.count('G') + i.count('C') + i.count('g') + i.count('c') total_bases += len(i) for i in fq2: i = str(i.seq) total_gc += i.count('G') + i.count('C') + i.count('g') + i.count('c') total_bases += len(i) print("GC_count:" + str(total_gc/total_bases)) (base) [b20223040323@admin1 test]$ python test.py GC_count:0.6710963455149501
。
更正一下,python脚本的解释器写错了,正确的应该为:
#!/usr/bin/env python
# -*- coding: uft-8 -*-
分类:
python
【推荐】国内首个AI IDE,深度理解中文开发场景,立即下载体验Trae
【推荐】编程新体验,更懂你的AI,立即体验豆包MarsCode编程助手
【推荐】抖音旗下AI助手豆包,你的智能百科全书,全免费不限次数
【推荐】轻量又高性能的 SSH 工具 IShell:AI 加持,快人一步
· 震惊!C++程序真的从main开始吗?99%的程序员都答错了
· 【硬核科普】Trae如何「偷看」你的代码?零基础破解AI编程运行原理
· 单元测试从入门到精通
· 上周热点回顾(3.3-3.9)
· winform 绘制太阳,地球,月球 运作规律
2024-02-01 linux 中 awk统计 文本指定列 各个类别的均值
2022-02-01 R语言中判断向量是否排序
2022-02-01 linux中如何关闭\开启SElinux
2021-02-01 linux系统redhat8.3安装R-4.0.3(编译安装)
2021-02-01 linux系统redhat 8.3 安装R
2021-02-01 linux系统centos8.0 安装R
2021-02-01 linux系统redhat7.9 安装R(编译安装)