/home/liuhui/bin/trinityrnaseq_r20140413p1/util/support_scripts/get_Trinity_gene_to_trans_map.pl Trinity.fasta  > Trinity.fasta_trans_map

bowtie-build Trinity.fasta Trinity.fasta.bowtie

rsem-prepare-reference --no-polyA --no-bowtie --transcript-to-gene-map Trinity.fasta_trans_map Trinity.fasta Trinity.fasta.RSEM

mkdir -p lib

cd lib

bowtie -q --all --best --strata -m 300 --chunkmbs 512 -X 800 -S -p 4 Trinity.fasta.bowtie -1 left.fq -2 right.fq | samtools view -S -b -o bowtie.bam -

rsem-calculate-expression --paired-end -p 4 --keep-intermediate-files --forward-prob 0 --bam bowtie.bam Trinity.fasta.RSEM RSEM

samtools sort bowtie.bam bowtie.csorted

samtools index bowtie.sorted.bam

gene_abundance_estimation的更多相关文章

随机推荐

  1. gitlab两种连接方式:ssh和http配置介绍

    gitlab环境部署好后,创建project工程,在本地或远程下载gitlab代码,有两种方式:ssh和http (1)ssh方式:这是一种相对安全的方式 这要求将本地的公钥上传到gitlab中,如下 ...

  2. Java8简单的本地缓存实现

    原文出处:lukaseder         Java8简单的本地缓存实现 这里我将会给大家演示用ConcurrentHashMap类和lambda表达式实现一个本地缓存.因为Map有一个新的方法,在 ...

  3. Linux下Mongodb安装和启动配置

    1.下载安装包 wget http://fastdl.mongodb.org/linux/mongodb-linux-i686-1.8.2.tgz 下载完成后解压缩压缩包 tar zxf mongod ...

  4. scala 学习笔记(05) OOP（中）灵活的trait

    trait -- 不仅仅只是接口! 接上回继续,scala是一个非常有想法的语言,从接口的设计上就可以发现它的与众不同.scala中与java的接口最接近的概念是trait,见下面的代码: packa ...

  5. 利用ActiveX实现web页面设置本地默认打印机、纸张大小

    通常web技术无法设置本地计算机的默认打印机,包括用代码设置纸张大小,如果业务系统中真遇到这种需求,只能通过其它辅助手段(比如ActiveX)实现.下面这段代码,出自网上被广泛使用的"泥人张 ...

  6. （十）装饰器模式详解（与IO不解的情缘）

    作者:zuoxiaolong8810(左潇龙),转载请注明出处,特别说明:本博文来自博主原博客,为保证新博客中博文的完整性,特复制到此留存,如需转载请注明新博客地址即可. LZ到目前已经写了九个设计模 ...

  7. 如何批量删除虚拟机及其关联的存储（Windows Azure）

    可以通过运行附件中PowerShell脚本文件RemoveVMandDisk.ps1批量删除VM和Disk,详细代码如下: param($serviceName) echo "Startin ...

  8. EF增删改查操作

    增 using (StudentEntities ent = new StudentEntities()) { User aNewUser = new User() { Age = , Name = ...

  9. [BZOJ1188][HNOI2007]分裂游戏（博弈论）

    题目:http://www.lydsy.com:808/JudgeOnline/problem.php?id=1188 分析: 设SG[i]表示一个石子在位置i上的SG值 这个很容易暴力求,因为i的后 ...

  10. 【BZOJ1001】狼抓兔子

    1001: [BeiJing2006]狼抓兔子 Time Limit: 15 Sec  Memory Limit: 162 MBSubmit: 7530  Solved: 1724[Submit][S ...