Changes between Version 8 and Version 9 of SOPs/rna-seq-diff-expressions
- Timestamp:
- 07/25/13 09:50:47 (12 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
SOPs/rna-seq-diff-expressions
v8 v9 38 38 * Note the "-" before the gtf filename; this indicates that the SAM file is coming from standard input (in this case, piped from samtools). 39 39 * For paired-end reads the sam file has to be sorted by read name. You can sort first the BAM file and then run htseq-count like this: 40 * ''bsub "samtools sort -n -m 5G accepted_hits.bam accepted_hitsSortedByname"'' 41 * To request a certain amount of memory and a specific node use ''bsub -R "rusage[mem=50000]" -m NodeName '' 40 42 * ''bsub "samtools sort -n -m 5G accepted_hits.bam accepted_hitsSortedByname"'' 41 43 * ''bsub "samtools view accepted_hitsSortedByname.bam | htseq-count -m intersection-strict --stranded=no - gene_models.gtf >| gene_model.counts.txt"''