转录组概述
6篇重要综述
最佳实践
分析策略
转录组测序数据分析公益视频学习笔记分享
https://mp.weixin.qq.com/s/b_NnpjSTLowos_ABsBeMTQ
中心法则
转录组概念
各种RNA在细胞中的绝对数量
普通转录组测序实验流程图-RNA-Seq
质控
2100峰图
建库(富集,打断,反转录,修复,加A, 接头,PCR)
文库
质控
测序
准备工作
linux命令复习
目录管理
FASTQ数据介绍
数据集介绍
fastq 格式
fastq 第一行
质量值
练习
## 1
zless SRR1039510_1.fastq.gz | grep '^@SRR' -c
# sed 版本
$ zless SRR1039510_1.fastq.gz | sed -n '/^@SRR/p' |wc -l
## 2
zless SRR1039510_1.fastq.gz | grep '^@SRR' > SRR1039510_1_id
# sed 版本
zless SRR1039510_1.fastq.gz | sed -n '/^@SRR/p' > SRR1039510_1_id
zless -S SRR1039510_1.fastq.gz | sed -n '1~4p' > SRR1039510_1_id
## 3
zless -S SRR1039510_1.fastq.gz | paste - - - - |cut -f 2 > SRR1039510_1_seq
# sed 版本
zless -S SRR1039510_1.fastq.gz | sed -n '2~4p' > SRR1039510_1_seq
## 4
zless SRR1039510_1.fastq.gz | paste - - - - |cut -f 2 |grep -o [A-Z] |wc -l
# sed 版本
zless -S SRR1039510_1.fastq.gz | sed -n '2~4p' |grep -o [A-Z] |wc -l
## 5
awk 'BEGIN {reads=25000*2;bases= 1575000*2} { print "reads: " reads*63/10^6" M\n" "dase:" bases/10^6 "M"}'