统计fasta序列条数

1.统计大于号开始的行数或seqkit 工具

# 通过搜索>的数量
grep -c '^>' myFasta.fasta
1397492
#seqkit统计提取，速度也是很快的
seqkit stats t.fa -T | grep -v file | cut -f 4
1397492
# 统计 1-100bp 范围长的序列数
cat t.fa | seqkit seq -m 1 -M 100 | seqkit stat -T | grep -v file | cut -f 4

Total sequence length 5,759,798,599
Total ungapped length 5,759,798,599
Number of contigs 1,397,492
Contig N50 9,587
Contig L50 174,483
Total number of chromosomes and plasmids 0
Number of component sequences (WGS or clone) 1,397,492

2.fastq序列条数统计

压缩格式解压，统计行数除以4

# 通常以fastq.gz格式压缩
zcat  input.fastq.gz | awk 'NR%4==2{c++} END{print c}'# 推荐下面的方法 pigz 会比gzip快10倍
pigz -dc input.fastq.gz | awk 'NR%4==2{c++} END{print c}'# 如果不是压缩格式
cat input.fastq | awk 'NR%4==2{c++} END{print c}'

转载于:https://www.cnblogs.com/huangyinger/p/10420774.html

统计fasta序列条数相关推荐

Excel统计QQ聊天记录条数
将QQ记录从QQ中导出为Txt文件,再导入的Excel表格中,会发现消息记录格式: 2012-06-23 10:58:34 小明你好! 2012-06-23 10:58:34 小红你好! 消息记录 ...
【总结】最短路径条数问题
例题:NOI2007 社交网络[难度不大,主要就是考察了这个知识点] 本文进过原作者同意转载,原博主为: erosun 原地址:http://www.cnblogs.com/acxblog/p/75 ...
去重之后统计条数_BOPET：12的普通包装膜到底去哪了？
导语近期BOPET市场多数客户反映,12μ的普通包装膜现货一货难求,前期订单交货紧张,新订单交期较长,12μ的普通包装膜到底去哪了呢? 对于BOPET市场来说,12μ普通包装膜交货紧.交货慢的现象从 ...
php中统计记录条数,使用GROUP BY的时候如何统计记录条数 COUNT(*) DISTINCT
例如这样一个表,我想统计email和passwords都不相同的记录的条数 CREATE TABLE IF NOT EXISTS `test_users` ( `email_id` int(11) u ...
mysql统计去重记录数量_MySQL分组后，如何统计记录条数
MySQL分组后,统计记录条数的方法:1.统计记录条数,代码为[SELECT num,count(*) AS counts from test_a GROUP BY num]:2.对num去重后的数量 ...
MyBatis Plus 的 Service 统计数据条数
MyBatis Plus 的 Service 统计数据条数 1. 简单介绍 2. 接口说明 3. 参数说明 4. 实例代码 4.1 无条件,查询整个数据表记录总数 4.2 有条件,查询整个数据表记录总 ...
sql 统计记录条数后打印出所有记录_用SQL完成购买行为分析（下篇II）
(接<用SQL完成购买行为分析(下篇I)>内容) 12)查询首条记录为fav,总记录条数为14的记录.将前面getNum(3)红框处替换为12,运行getNum(14)得到第2条记录的数量 ...
mongodb 高级查询统计记录条数
使用count()方法查询表中的记录条数,例如,下面的命令查询表users的记录数量: db.users.find().count();
百度新闻资讯类信息爬虫--统计一年内关键词新闻的条数
背景通过百度词条搜索,来查找300个关键词,在一年内发布新闻的条数. 最终效果实现如下: 实现思路实现思路依然是: 先根据多页的url,来找到规律,构建起一页的url: def format_ur ...
CAD库中统计PBN运行航路条数和总距离
select 'PBN运行航路' 类型, fb.b 总条数, fa.a 总距离 from ( select sum(s) a fro ...

统计fasta序列条数

1.统计大于号开始的行数或seqkit 工具

2.fastq序列条数统计

统计fasta序列条数相关推荐

最新文章

热门文章