Changes between Version 41 and Version 42 of SOPs/rna-seq-diff-expressions


Ignore:
Timestamp:
05/17/17 10:55:39 (8 years ago)
Author:
thiruvil
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • SOPs/rna-seq-diff-expressions

    v41 v42  
    5959      * Sample commands:
    6060{{{
     61#default: unstranded
    6162#single-end reads
    6263featureCounts -a gene_anotations.gtf -o MySample.featureCounts.txt MySample.bam
    6364#PE reads
    6465featureCounts -p -a gene_anotations.gtf -o MySample.featureCounts.txt MySample.bam
     66
     67#stranded (fwd)
     68featureCounts -p -s 1 -a gene_anotations.gtf -o MySample.featureCounts.txt MySample.bam
     69
    6570}}}
    6671    * For some analyses (or for visualization), you can add a pseudocount (such as 1 or another small number) to all genes in all samples to prevent log2 ratios that require dividing by 0 and reduce background count noise -- BUT be aware that some statistical methods (like DESeq) require raw input values without any pseudocounts or normalization.