Changes between Version 31 and Version 32 of SOPs/SAMBAMqc
- Timestamp:
- 11/19/21 09:02:38 (3 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
SOPs/SAMBAMqc
v31 v32 127 127 Fraction of reads explained by "1+-,1-+,2++,2--": 0.0193 128 128 Fraction of reads explained by other combinations: 0.0000 129 # For gene counting: (feature counts, 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. 130 130 131 131 # sample output on strand-specific PE reads (since the second fraction is much larger than the first fraction): … … 134 134 Fraction of reads explained by "1+-,1-+,2++,2--": 0.9807 135 135 Fraction of reads explained by other combinations: 0.0000 136 # For gene counting: (feature counts, 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. 137 137 138 138 # sample output on non-stranded PE reads (since both fractions are about the same): … … 141 141 Fraction of reads explained by "1+-,1-+,2++,2--": 0.4897 142 142 Fraction of reads explained by other combinations: 0.0000 143 # For gene counting: (feature counts, 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. 144 144 145 145 #sample output on stranded SE reads: … … 148 148 Fraction of reads explained by "++,--": 0.9865 149 149 Fraction of reads explained by "+-,-+": 0.0068 150 # For gene counting: (feature counts, 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. 151 151 152 152 #sample output on stranded SE reads: … … 155 155 Fraction of reads explained by "++,--": 0.0068 156 156 Fraction of reads explained by "+-,-+": 0.9865 157 # For gene counting: (feature counts, 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. 158 158 159 159