当前位置: 首页 > news >正文

python脚本统计fastq数据的GC含量

 

001、

(base) [b20223040323@admin1 test]$ ls
SRR1770413_1.fastq  test.py
(base) [b20223040323@admin1 test]$ cat test.py              ## 脚本
#!/usr/bin/env python
# -*- coding:utf-8 -*-in_file = open("SRR1770413_1.fastq", "r")
out_file = open("result.txt", "w")gc_count = 0
total_base = 0lines = 0
for i in in_file:lines += 1if lines % 4 == 2:i = i.strip().upper()gc_count += i.count("G") + i.count("C")total_base += len(i)gc_percent = gc_count/total_base
out_file.write("SRR1770413_1.fastq: " + str(gc_percent) + "\n")in_file.close()
out_file.close()
(base) [b20223040323@admin1 test]$ python test.py
(base) [b20223040323@admin1 test]$ ls
result.txt  SRR1770413_1.fastq  test.py
(base) [b20223040323@admin1 test]$ cat result.txt     ## 结果文件
SRR1770413_1.fastq: 0.5009478013778936

 

image

 。

 

http://www.hskmm.com/?act=detail&tid=22826

相关文章:

  • 第一次软件工程作业
  • 关于数颜色
  • 2025仿石漆厂家 TOP 企业品牌推荐排行榜,真石漆,水包砂,冠晶石,外墙,多彩,批刮,别墅,批发,出口仿石漆推荐这十家公司!
  • CF1615F LEGOndary Grandmaster
  • 2025 年仿石漆品牌最新推荐排行榜:聚焦真石漆仿石漆,水包砂仿石漆,冠晶石仿石漆,外墙仿石漆,多彩仿石漆供采购参考
  • 诚信液压阀块定制厂家推荐榜:实力制造与精准工艺口碑之选
  • 24 两两交换链表中的节点
  • 算法导论图论部分总结
  • 19 删除链表的倒数第 N 个结点
  • 浅谈 Bakas Trick / 不带删尺取 / 对顶栈
  • 聚变堆:中国BEST装置全面开建
  • web
  • 如何用pivotby函数实现数据透视(2)
  • 2025 年彩钢板厂家 TOP 企业品牌推荐排行榜,复合彩钢板,保温彩钢板,耐腐蚀彩钢板,净化彩钢板推荐这十家公司!
  • 251002
  • AT_agc020_d [AGC020D] Min Max Repetition
  • 题解:P7810 [JRKSJ R2] Upper
  • 记录自己被AWS坑了6刀
  • 如何用pivotby函数实现数据透视(1)
  • 10/1
  • 详细介绍:【数据结构】图论核心应用:关键路径算法详解——从AOE网到项目管理实战​
  • tnkstat3e-merge-0
  • JavaScript零基础入门速通(完整) - 指南
  • @RequestParam 什么时候可以省略?
  • 段页式管理方式
  • 推进电子设计革新:为什么模拟仿真正是核心助力?
  • 深入解析:前端开发,iframe 相关经验总结
  • list 容器 listr容器与vector容器 list 示例代码
  • advance 函数
  • *和 指针与地址 ++i和 i++