2018-10-12
组装或定量 eXpress

https://pachterlab.github.io/eXpress/downloads/express-1.5.1/express-1.5.1-linux_x86_64.tgz

使用方法

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
express [options]* <target_seqs.fasta> <aligned_reads.(sam/bam)>

# <target_seqs.fasta>: 参考序列的fasta文件
# <aligned_reads.(sam/bam)>: 比对结果文件
# -o path: 指定输出路径
# -B n: 指定进行batch EM rounds计算的次数,以时间为代价提高定量准确性
# -O n: 指定进行online EM rounds计算的次数,以时间为代价提高定量准确性
# -m n: 指定fragment的平均长度
# -s n:指定fragment长度的标准差
# -H str: 以逗号分隔的位置坐标,指定haplotypes的位置。利于等位基因的特异表达定量
# --output-align-samp:
# --output-align-prob:
# --fr-stranded: 指定链特异性比对的方向,read1比对到forward链,read2比对到reverse链
# --rf-stranded:指定链特异性比对的方向,read1比对到reverse链,read2比对到forward链
# --f-stranded: single-end, read比对到forward链
# --r-stranded: single-end, read比对到reverse链
#

输入文件

eXpress需要一个multi-fasta格式的输入文件,计算该文件序列的转录丰度。

Read More

2018-10-12
组装或定量 kallisto_sleuth

免比对的定量:kallisto

kallisto是一个align-free的测序结果定量工具。速度快得很.他可以在10min内完成index的构建,然后花3min完成30 million的human reads的定量。它不仅快速而且准确。

Pseudoalignment保存了用于定量的关键信息。An important feature of kallisto is that it outputs bootstraps along with the estimates of transcript abundances

Read More

2018-10-12
组装或定量 stringTie

stringTie 是用于 RNA-seq 的转录本组装和定量软件

Read More