Changes between Version 31 and Version 32 of SOPs/SAMBAMqc


Ignore:
Timestamp:
11/19/21 09:02:38 (2 years ago)
Author:
gbell
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • SOPs/SAMBAMqc

    v31 v32  
    127127Fraction of reads explained by "1+-,1-+,2++,2--": 0.0193
    128128Fraction of reads explained by other combinations: 0.0000
    129 # For gene counting: (featurecounts, use -p -s 1; htseq-count, use --stranded=yes); mapping with TopHat should have been performed with --library-type=fr-secondstrand.
     129# For gene counting: (featureCounts, use -p -s 1; htseq-count, use --stranded=yes); mapping with TopHat should have been performed with --library-type=fr-secondstrand.
    130130
    131131# sample output on strand-specific PE reads (since the second fraction is much larger than the first fraction):
     
    134134Fraction of reads explained by "1+-,1-+,2++,2--": 0.9807
    135135Fraction of reads explained by other combinations: 0.0000
    136 # For gene counting: (featurecounts, use -p -s 2; htseq-count, use --stranded=reverse); mapping with TopHat should have been performed with --library-type=fr-firststrand.
     136# For gene counting: (featureCounts, use -p -s 2; htseq-count, use --stranded=reverse); mapping with TopHat should have been performed with --library-type=fr-firststrand.
    137137
    138138# sample output on non-stranded PE reads (since both fractions are about the same):
     
    141141Fraction of reads explained by "1+-,1-+,2++,2--": 0.4897
    142142Fraction of reads explained by other combinations: 0.0000
    143 # For gene counting: (featurecounts, use -p -s 0; htseq-count, use --stranded=no); mapping with TopHat should have been performed with --library-type=fr-unstranded.
     143# For gene counting: (featureCounts, use -p -s 0; htseq-count, use --stranded=no); mapping with TopHat should have been performed with --library-type=fr-unstranded.
    144144
    145145#sample output on stranded SE reads:
     
    148148Fraction of reads explained by "++,--": 0.9865
    149149Fraction of reads explained by "+-,-+": 0.0068
    150 # For gene counting: (featurecounts, use -s 1; htseq-count, use --stranded=yes; mapping with TopHat should have been performed with --library-type=fr-secondstrand.
     150# For gene counting: (featureCounts, use -s 1; htseq-count, use --stranded=yes; mapping with TopHat should have been performed with --library-type=fr-secondstrand.
    151151
    152152#sample output on stranded SE reads:
     
    155155Fraction of reads explained by "++,--": 0.0068
    156156Fraction of reads explained by "+-,-+": 0.9865
    157 # For gene counting: (featurecounts, use -s 2; htseq-count, use --stranded=reverse; mapping with TopHat should have been performed with --library-type=fr-firststrand.
     157# For gene counting: (featureCounts, use -s 2; htseq-count, use --stranded=reverse; mapping with TopHat should have been performed with --library-type=fr-firststrand.
    158158
    159159