SuperAPI AI 工具集

bio-read-alignment-hisat2-alignment

通用与其它

查看项目

更新于 10-15 21:55

技能介绍

用途:HISAT2 是内存高效的剪接感知比对器,用于将 RNA-seq 读段比对到参考基因组,适合 STAR 内存需求过高或常规 RNA-seq 比对场景。

主要能力:支持单端与双端比对;可构建带剪接位点与 exon 注释的索引;支持链特异性文库(FR/RF/F/R);可输出新剪接位点并用于后续比对;支持两遍比对;可设置 read group;结果可管道输出为排序 BAM 供下游定量。

使用方式:构建索引时可用 hisat2-build -p 8 reference.fa hisat2_index,推荐先用 hisat2_extract_splice_sites.py 与 hisat2_extract_exons.py 从 GTF 提取注释,再以 --ss 和 --exon 建索引。双端比对示例:hisat2 -p 8 -x hisat2_index -1 reads_1.fq.gz -2 reads_2.fq.gz -S aligned.sam;单端用 -U。可直接管道到 samtools sort 生成排序 BAM 并 samtools index。链特异性用 --rna-strandness,双端 RF 常见于 dUTP/TruSeq。新剪接位点用 --novel-splicesite-outfile 与 --novel-splicesite-infile。两遍比对先对全部样本输出位点,合并去重后再比对。read group 用 --rg-id 与 --rg SM/PL/LB。下游定量可输出 name-sorted BAM 供 htseq-count,或 coordinate-sorted 供 featureCounts;StringTie/Cufflinks 场景加 --dta 或 --dta-cufflinks。

关键参数:-p 线程数默认 1;-x 索引前缀;--rna-strandness 默认 unstranded;--min-intronlen 默认 20;--max-intronlen 默认 500000;-k 默认 5。比对摘要输出到 stderr,可重定向保存。

适用场景与限制:适合内存有限的人类基因组比对(资料给出 STAR 约 30GB、HISAT2 约 8GB 的对比)。具体版本、安装方式、平台支持及收费信息暂未提供。