用tophat 将clear reads mapping 到植物基因组上
2016-07-07 11:47
381 查看
nohup /home/lixiangyong/software/tophat/tophat2 -p 8 -G /data/plant_genome/Lj3.0/Lj3.0_gene.gtf -o tophat_output /data/plant_genome/Lj3.0/Lj3.0_genome \*_1/R_R1.fastq_filtered_trimmed \*_2/R_R2.fastq_filtered_trimmed &
详情参见
http://blog.sciencenet.cn/u/nsaa001
结果会在tophat_output 文件夹里面生成 accept_hits.bam 文件将mapping上的reads合成为转录组数据
nohup cufflinkss -g /data/plant_genome/Lj3.0/Lj3.0_gene.gtf -o cufflinks_sample -p 16 accepted_hits.bam &
会在cufflink_sample 中输出4个文件
genes.fpkm_tracking isoforms.fpkm_tracking skipped.gtf transcripts.gtf
相关文章推荐
- GWAS分析软件—— 一、Plink的安装
- 监控linux系统某文件的生成,并进行另一个脚本
- 小RNA的测序技术路线以及分析流程
- 如何统计id很复杂的fasta文件的长度?
- 基因组变异检测概述(SNP、InDel、SV)
- centOS集群repeatmarker软件安装
- 生物信息——连锁不平衡 Linkage Disequilibrium
- 国际人类基因组单体型图计划
- 生物信息常用数据库
- 基因序列分析(生物信息学论坛)
- fasta格式 图解
- blast 详解
- (有参考基因组)植物转录组分析之一数据处理
- 生信 求旁系同源基因
- 生信 批量替换文件夹 内 某一字符
- 生信 Fastq 文件讲解
- 分子进化之一,序列的获得
- 转录组分析之 Trimming对reads进行处理
- 如何去除测序数据中的接头和低质量的reads? 软件fastx
- 关于各种转录因子的数据库