数据处理 第一天 准备工作

数据分析 开始

1 上传数据(.gz)完成

2 解压数据

for i in *.gz; do gzip -d $i; done

du *.fq:查看文件大小

3:数据质控

/baicai2/weiyx/local/FastQC/fastqc -q -t 18 ./Resequencing/cleandata/cleandata.fq/*.fq -o ./Resequencing/fastQC/ &

-q 安静 -t 线程 -o 输出 

4:数据比对

nohup bowtie2 -p 5 -S /baicai2/weiyx/map_bra/755.sam -x /baicai2/weiyx/genome_bra/genome -1 /baicai2/weiyx/Resequencing/cleandata/cleandata.fq/S755_NDHG07043_1.clean.fq -2 /baicai2/weiyx/Resequencing/cleandata/cleandata.fq/S755_NDHG07043_2.clean.fq 2> bowtie2_bra.log &

-p 调用核心数目 -S 输出 -x index 

挺快

参考 https://zhuanlan.zhihu.com/p/31432993

©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成,浏览时请结合常识与多方信息审慎甄别。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

推荐阅读更多精彩内容

友情链接更多精彩内容