python 计算 fastq数据中的reads数目、碱基数目

 

001、

复制代码
(base) [b20223040323@admin1 test]$ ls
SRR1770413_1.fastq  SRR1770413_2.fastq  test.py
(base) [b20223040323@admin1 test]$ cat test.py
#!/usr/bin/env python
# -*- coding:utf-8 -*-

from Bio import SeqIO
import sys

fq1 = list(SeqIO.parse(sys.argv[1], "fastq"))
fq2 = list(SeqIO.parse(sys.argv[2], "fastq"))

total_reads = 0
total_bases = 0

for i in fq1:
        i = str(i.seq)
        total_reads += 1
        total_bases += len(i)
for i in fq2:
        i = str(i.seq)
        total_reads += 1
        total_bases += len(i)

print("total_reads: " + str(total_reads))
print("total_bases: " + str(total_bases))
(base) [b20223040323@admin1 test]$ python test.py SRR1770413_1.fastq SRR1770413_2.fastq
total_reads: 6
total_bases: 1806
复制代码

 

 。

 

posted @   小鲨鱼2018  阅读(9)  评论(0编辑  收藏  举报
相关博文:
阅读排行:
· 震惊!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(编译安装)
点击右上角即可分享
微信分享提示