跟着Nature Communications学习Hisat-Trinity-PASA等分析流程

小杜的生信筆記

已于 2023-12-05 16:12:34 修改

阅读量415

点赞数 5

分类专栏： R语言精美图形绘制教程文章标签：数据库开发语言学习网络图富集功能富集

于 2023-12-05 15:48:27 首次发布

本文链接：https://blog.csdn.net/kanghua_du/article/details/134809471

版权

博客分享了如何跟着Nature Communications学习Hisat2-Samtools、Trinity de novo和PASA等转录组分析流程，适合初学者，提供详细教程和脚本，帮助理解本地适应和气候变化对森林树木的影响。

摘要由CSDN通过智能技术生成

一边学习，一边总结，一边分享！

详细教程请访问：
组学分析流程

本期分析流程

Hisat2-Samtools
Trinity_GG_denovo
PASA

…

本期教程文章

题目：Genomic insights into local adaptation and future climate-induced vulnerability of a keystone forest tree in East Asia

Hisat2-samtools分析流程

#!/bin/bash

genome=$1
index=${genome%.*}
rna_1_fq=`cat $2|grep 1P|sed ":a;N;s/\n/,/g;ta"` #1.fq path list
rna_2_fq=`cat $2|grep 2P|sed ":a;N;s/\n/,/g;ta"` #2.fq path list

#echo $index
hisat2-build -p 20 $genome $index

hisat2 -x $index \
           -1 $rna_1_fq\
           -2 $rna_2_fq\
           --threads 20 \
           --min-intronlen 20 \
           --max-intronlen 20000 \
           --dta \
           --score-min L,0.0,-0.4 \
           -S ${index}.sam


samtools sort -@ 20 \
                  -o ${index}.sorted.bam \
                      -O BAM \
                ${index}.sam

PSSA_align

#!/bin/bash

export PATH="$PATH:/usr_storage/jcf/.conda/envs/PASA"
source  /pub_storage2/new_PASA/.bashrc

#cat $Trinity_GG $Trinity_denovo >transcripts.fasta #
transcripts_fasta="$1" # transcripts.fasta generated from merging fasta file of Trinity denovo and Trinity genome guided mod