统计fastq文件中读段的数量

mycount=`cat SRR108114_new_1.fastq | wc -l`
echo ‘Number of clean reads, SRR108114_new_1.fastq: ‘$((mycount/4)) >> SRR108114.stat
mycount=`cat SRR108114_new_2.fastq | wc -l`
echo ‘Number of clean reads, SRR108114_new_2.fastq: ‘$((mycount/4)) >> SRR108114.stat

时间： 2024-04-27 04:30:04

统计fastq文件中读段的数量的相关文章

Java学习（4）：统计一个文件中的英文，中文，数字，其他字符以及字符总数

要求:统计一个文件中的英文,中文,数字,其他字符以及字符总数(此随笔以txt文件为例) import java.io.BufferedReader; import java.io.File; import java.io.FileInputStream; import java.io.FileNotFoundException; import java.io.IOException; import java.io.InputStreamReader; /** * 将一个文件中英文,中文,数字,其

统计一个文件中出现字符'a'的次数

# -*- coding: utf-8 -*- #python 27 #xiaodeng #统计一个文件中出现字符'a'的次数 #http://www.cnblogs.com/hongten/p/hongten_python_count.html import os number=0 def getNumber(filePath,c): 'c---->the word numbers' #统计一个文件中出现字符'a'的次数 if os.path.exists(filePath): global

Python 统计yaml文件中数字出现的次数

背景需求: 公司开发了一个抽奖系统,大概功能是可以设置抽奖号码,然后设置抽奖的等级及数量,再从设置的号码中抽取. 由于是抽奖系统,需要评估一下数字中奖的概率.我对这个系统进行了以下处理 1.编写初始化抽选号码,即根据初始化接口清空已有号码数据 2.根据设置号码池 3.根据抽奖接口返回数据取出中奖号码 4.将数据写入yaml文件 5.设置100次循环抽奖 6.对yaml文件进行处理,读取数据 7.按照条件取出数据 8.取出文件中出现频率最高的数字最后结果为附yaml文件读取代码 import

关于从文件中读数据，将字符串写到文件中的操作

import java.io.*; import java.util.ArrayList; import java.util.List; /** * Created by guan on 11/1/16. */ public class ReadFromFile { /** * 随机读取文件内容 */ public static void readFileByRandomAccess(String fileName) { RandomAccessFile randomFile = null; t

统计代码文件中的实际有效行数，去掉空行、单行注释、多行注释

#coding=gbk import os #rootdir='f:\\pylianxi' def count_line_core(file_name): ##传入单个文件,统计行数,之后返回该文件的实际代码行数:区分utf-8.gbk有待优化 print('core_file_name:',file_name) lines_count=0 flag=True try: with open(file_name,'r',encoding='gbk') as fp: print('gbk file_

[linux] shell脚本编程-统计日志文件中的设备号发通知邮件

1.日志文件列表比如:/data1/logs/2019/08/15/ 10.1.1.1.log.gz 10.1.1.2.log.gz 2.统计日志中的某关键字shell脚本 zcat *.gz|grep 关键字 |grep -oP "deviceid=[^=]+"|uniq|sort -u > /tmp/20190815.log date 格式化出年月等信息,拼接成路径 wc -l /tmp/20190815.log , 获取到行数 php /xxxxx/sendmail.ph

统计python文件中的代码,注释,空白对应的行数

其实代码和空白行很好统计,难点是注释行 python中的注释分为以#开头的单行注释或者以'''开头以'''结尾或以"""开头以"""结尾的文档注释,如: ''' hello world '''和 ''' hello world''' 思路是用is_comment记录是否存在多行注释,如果不存在,则判断当前行是否以'''开头,是则将is_comment设为True,否则进行空行.当前行注释以及代码行的判断,如果is_comment已经为True

mothur summary.seqs 统计fasta文件中每条序列的长度

在介绍summary.seqs的用法之前,我们首先需要搞清楚两个概念: 1)ambiguous bases 中文叫做模糊碱基,对于DNA序列来说,只有ATCG 4种碱基,在IUPAC定义的碱基标准中,出了上述4种碱基之外,还包括其他的碱基,可以代表不同类型的碱基代码英文含义中文含义 G Guanine 鸟嘌啉 A Adenine 腺嘌啉 T (U) Thymine (Uracil) 胸腺嘧啶 (尿嘧啶) C Cytosine 胞嘧啶 R (A or G) PuRine 嘌啉 Y

利用MySQL统计一列中不同值的数量方法示例

前言本文实现的这个需求其实十分普遍,举例来说,我们存在一个用户来源表,用来标记用户从哪个渠道注册进来.表结构如下所示- 其中 origin 是用户来源,其中的值有 iPhone .Android .Web 三种,现在需要分别统计由这三种渠道注册的用户数量. 解决方案1 ? 1 2 3 4 5 6 7 8 9 SELECT count(*) FROM user_operation_log WHERE origin = 'iPhone'; SELECT count(*) FROM user_ope