从 RNA-seq 比对结果组装转录本并估计表达量的工具。Tool for transcript assembly and expression estimation from RNA-seq alignments.
mamba install -c bioconda stringtie
按参考注释估计表达:
stringtie sample.sorted.bam -p 8 -G annotation.gtf \
-e -B -o sample/stringtie.gtf
无参考组装:
stringtie sample.sorted.bam -p 8 -o sample.assembled.gtf
合并多个样本组装结果:
stringtie --merge -p 8 -G annotation.gtf \
-o merged.gtf mergelist.txt
生成 gene abundance 表:
stringtie sample.sorted.bam -p 8 -G annotation.gtf \
-A sample.gene_abund.tsv -o sample.gtf
-G:参考注释 GTF/GFF。-e:仅估计参考转录本表达,不组装新转录本。-B:生成 Ballgown 兼容输出。-A:输出 gene abundance 表。--merge:合并多个 GTF。-p:线程数。