Kevin2015的个人博客分享 http://blog.sciencenet.cn/u/Kevin2015

博文

全外显子组生信分析流程-8-depth_and_coverage_readfq

已有 2384 次阅读 2019-3-23 11:24 |个人分类:全外显子项目|系统分类:科研笔记

#! /bin/bash

#download https://github.com/billzt/readfq
#purpose: calculate reads from fastq files
#set env
readfq="/home/zhanghl/supporting_softwares/readfq/readfq-master/kseq_fastq_base"
fq_dir=""
cd $fq_dir  #fq所在位置
ls *.fq.gz > fq.list #collect fastq files
for i in $(cat fq.list)
do 
  printf $i "\t" >> fq.reads_num
  $readfq $i     >> fq.reads_num
done




https://blog.sciencenet.cn/blog-2609994-1169178.html

上一篇:全外显子组生信分析流程-7-depth_and_coverage_qualimap
下一篇:全外显子组生信分析流程-9-depth_and_coverage_samtools_flagstat
收藏 IP: 14.204.63.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...
扫一扫,分享此博文

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-7-27 21:24

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部