使用C语言快速统计fastq文件q20、q30、GC含量

二代测序的分析过程中，经常需要统计原始下机数据的数据量，看数据量是否符合要求；另外还需要统计reads数量、base数据量、q20,q30,GC含量等反应测序质量的指标；

在kseq.h 的基础上稍加改造，就可以实现从fastq 文件中统计这些指标的功能，而且速度非常的快
kseq.h下载地址：

fastq_stat.c文件内容

#include <zlib.h>
#include <stdio.h>
#include <string.h>  #include "kse

使用C语言快速统计fastq文件q20、q30、GC含量相关推荐

Linux下使用wc命令快速统计所有文件的行数
语法:wc [选项] 文件- 说明:该命令统计给定文件中的字节数.字数.行数.如果没有给出文件名,则从标准输入读取.wc同时也给出所有指定文件的总统计数.字是由空格字符区分开的最大字符串. 该命令各选 ...
fastq序列质量Q30统计软件
FaQCs 使用perl编写,R画图.不仅可以统计质量,Q30,还可以trim序列.安装简单准备好两个需要的perl 模块,下载github上的文件,然后直接cd lib;./INSTALL.sh. ...
R语言导入其他统计软件文件
使用R中的foreign包可以较方便地读取其他统计软件的数据文件,比如SPSS.SAS等数据文件等函数描述 read.arff 从ARFF文件中读取文件,著名的数据挖掘开源软件weka的数据就是这 ...
Linux查询压缩文件行数,linux命令行快速统计文件（压缩文件）的行数
统计(文件|压缩文件)的行数 zcat file.gz | sed -n '$=' #迅速.直接打印出多少行.-n 取消 ...
易语言取c盘文件夹中的文件被占用,win7系统怎么快速找出文件被什么进程占用|win7解除文件占用的方法...
‍‍ 很多时候我们想删除文件或文件夹的时候,会出现文件正在使用的提示导致文件无法删除,如果知道是什么程序在使用这个文件好,但多数时候我们并不知道是被什么程序占用,那么在win7系统怎么快速找出文件被什 ...
linux快速统计目录大小,linux下统计文件夹、文件的大小--du
Linux下统计文件夹大小 du -sh ./ 统计文件夹占用的空间 find ./ -type f xargs ls -l awk 'BEGIN { size=0;}{size+=$5};END{p ...
linux fastq.gz文件如何解压,科学网—批量gzip，ungzip快速压缩和解压fastq文件 - 陈明杰的博文...
二代测序fastq文件太大,解压,压缩很浪费时间. 常规:使用gzip gunzip单线程压缩进阶1:python调用threading,多线程压缩,相当于开了N个窗口进行解压和压缩,略. 进阶2: ...
fastp: 极速全能的FASTQ文件自动质控+过滤+校正+预处理软件
软件作者介绍陈实富博士,海普洛斯联合创始人 / CTO 海普洛斯是全球领先的精准医疗和基因大数据国家高新技术企业,拥有 Illumina NovaSeq. HiSeq X10.NextSeq等全系列 ...
C 语言快速入门，21 个小项目足矣！「不走弯路就是捷径」
C 语言作为大学理工科专业的必修,是很多同学走进编程世界的第一课.那么怎样才能更好的入门 C 语言呢? 下面整理了 21 个 C 语言练手项目,从基础语法开始,逐步深入,通过一个个练手项目,让你轻松驰 ...
python向服务器上传fq文件,用python-fas读取大型fastq文件
我有几个fastq文件,平均有500.000.000行(125000.000个序列).有没有一种快速读取这些fastq文件的方法. 我想做的是,读取每个序列并使用前16个序列作为条形码.然后统计每个文 ...

使用C语言快速统计fastq文件q20、q30、GC含量

使用C语言快速统计fastq文件q20、q30、GC含量相关推荐

最新文章

热门文章