该分析流程主要根据2016年发表在Nature Protocols上的一篇名为Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown的文章撰写的,主要用到以下三个软件: 1. 软件介绍 1. HISAT (ccb.jhu.edu/software/hi)利用大量FM索引,以覆盖整个基因组,能够将RNA-Seq的读取与基因组进...
perl CountToFPKM.pl 04_Result/Stringtie/featureCounts/01.all.count.txt > 04_Result/Stringtie/featureCounts/02.all.FPKM.txt 一、写在前面 今天分享一个转录组上游分析的流程(Hisat2-Stringtie-Count),此流程的操作依旧是非常简单的。我们的流程主要使用软件的安装、数据下载、过滤、比对、Count、Count To FPKM...
stringtie -e -B -p 6 -G black_stringtie_merged.gtf -o sec_gtf/sample_1C/str_1C.gtf 1C.sorted.bam 第二种情况: #不经过融合多个样本的结果这一步骤后,可以很好的避免差异分析中大量的STRING.xxx stringtie -e -B -p 6 -G sheep.gtf -o third_gtf/th_1C.gtf 1C.sorted.bam #后续用dese...
StringTie是一个用于从RNA-seq数据中组装转录本的工具,其主要功能包括: 转录本组装:StringTie 可以从比对后的reads中重新构建转录本,包括已知的和新发现的转录本。这意味着它可以识别出新的剪接变体或未知转录本。 定量表达:除了组装之外,StringTie 还可以估计每个转录本的表达量,通常使用FPKM(fragments per kilobase o...
在短读组装器中，StringTie的转录本水平精度平均比Cufflinks高11％，转录本水平灵敏度比Cufflinks高25％。
stringtie -e -B -p 4 wt_Rep1.sorted.bam -G stringtie_merged.gtf -o athaliana_wt_Rep1/athaliana_wt_Rep1 stringtie -e -B -p 4 wt_Rep2.sorted.bam -G stringtie_merged.gtf -o athaliana_wt_Rep2/athaliana_wt_Rep2.count -A athaliana_wt_Rep2/wt_Rep2_gene_abun.out stringtie -e...
在linux服务器上sra数据的有参转录组定量分析,使用sratoolkit进行sra数据转换,使用trimmomatic进行read修剪,再使用hisat2进行比对,samtools进行sorting,最后使用stringtie进行定量。 欢迎大家批评指正,转载请标明出处:https://www.ivistang.com/articles/312 (opens new window) # 准备工作 # 安装miniconda wget -c https...
被称为表观转录组（epitranscriptome）的RNA修饰正成为基因调控的广泛调控机制。由于绘制转录组范围RNA修饰...
转录组分析---Hisat2+StringTie+Ballgown使用 1.Hisat2建立基因组索引: First, using the python scripts included in the HISAT2 package, extract splice-site and exon information from the gene annotation file: $ extract_splice_sites.py gemome.gtf >genome.ss#得到剪接位点信息...
conda install -c bioconda sra-tools hisat2 stringtie -y ###下载trimmomatic并解压 wget http://www.usadellab.org/cms/uploads/supplementary/Trimmomatic/Trimmomatic-0.39.zip unzip Trimmomatic-0.39.zip 注:如果调用samtools出现如下报错:samtools: error while loading shared libraries: libcrypto.so.1.0.0:...