用户工具

站点工具


stat:rnaseq

这是本文档旧的修订版!


RNA-Seq

教程

流程

https://bioconductor.org/packages/release/BiocViews.html#___GeneExpressionWorkflow
Raw Data: Fastq (.gz) →

  1. QC: FastQC
  2. Filter
  3. Alignment: STAR
  4. Count
  5. Cluster
  6. Heatmap
  7. Differential

1 QC

FastQC

apt install fastqc

fastqc --noextract RawData/I409/I409_1.fq.gz -o results/1_initial_qc/ 

2 Alignment

基因组注释数据

Alignment Indexing

对齐工具

HISAT2

STAR

2.7.10b

apt install rna-star

# 创建索引,索引文件创建一次即可. 需要从Ensembl下载对应物种的Fasta文件和GTF文件。
STAR --runMode genomeGenerate --genomeDir star_index --genomeFastaFiles fasta/* --sjdbGTFfile gtf/* --runThreadN 14

# 运行分析
STAR --genomeDir star_index --readFilesIn filtered/sample_filtered.fq --outSAMtype BAM SortedByCoordinate --quantMode GeneCounts --runThreadN 14

STAR --genomeDir star_index --readFilesIn rna4/RawData/I409/I409_1.fq --outSAMtype BAM SortedByCoordinate --quantMode GeneCounts --runThreadN 14

kallisto

# bat
kallisto bus [arguments] FASTQ-files

kallisto quant -i rattus_index_ki/transcriptome.idx -o reads.kallisto_quant -t 64 --fusion --pseudobam --genomebam --gtf gtf\Rattus_norvegicus.mRatBN7.2.109.gtf rna4\RawData\I409\I409_1.fq.gz rna4\RawData\I409\I409_2.fq.gz

3 Count

stat/rnaseq.1687791642.txt.gz · 最后更改: 2023/06/26 15:00 由 inkit