/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. zoj 1610

    Count the Colors Time Limit: 2 Seconds      Memory Limit: 65536 KB Painting some colored segments on ...

  2. Oracle round函数是什么意思?怎么运用?

    如何使用 Oracle Round 函数 (四舍五入) 描述 : 传回一个数值,该数值是按照指定的小数位元数进行四舍五入运算的结果. SELECT ROUND( number, [ decimal_p ...

  3. python的高级特性3:神奇的__call__与返回函数

    __call__是一个很神奇的特性,只要某个类型中有__call__方法,,我们可以把这个类型的对象当作函数来使用. 也许说的比较抽象,举个例子就会明白. In [107]: f = abs In [ ...

  4. wget: unable to resolve host address的解决方法

    摘要:wget:无法解析主机地址.这就能看出是DNS解析的问题. wget:无法解析主机地址.这就能看出是DNS解析的问题. 解决办法: 登入root(VPS).进入/etc/resolv.conf. ...

  5. 【转】TestFlight APP测试(IOS如何让上架前给其他人测试)

    原文网址:http://blog.csdn.net/dexin5195/article/details/43966571 大家都知道, 以前iOS项目要测试只需要上传到testflightapp.co ...

  6. some basic graph theoretical measures

    · mean characteristic path length calculated as the average length of the shortest path between two ...

  7. SQL80001: Incorrect syntax near ':'

    原文连接:http://geekswithblogs.net/tonyt/archive/2010/03/05/138363.aspx   SQL80001: Incorrect syntax nea ...

  8. C# Math类简介

    Math.abs() 计算绝对值. Math.acos() 计算反余弦值. Math.asin() 计算反正弦值. Math.atan() 计算反正切值. Math.atan2() 计算从x 坐标轴到 ...

  9. 关于在线预览word,excel,ppt,pdf的需求处理方法。

    参考文档:http://www.cnblogs.com/wolf-sun/p/3574278.html 我选用的方案:先用office com组件生成pdf,然后使用pdf.js在线预览pdf文档.在 ...

  10. 微软职位内部推荐-B&I Site Lead

    微软近期Open的职位: Job Title: B&I Site Lead Division: Phones Quality, SWIQ Location: Beijing, China Mi ...