Changes between Version 20 and Version 21 of SOPs/miningSAMBAM
- Timestamp:
- 11/17/14 10:43:37 (10 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
SOPs/miningSAMBAM
v20 v21 97 97 === QC to get a (visual) summary of mapping statistics. For eg. coverage/distribution of mapped reads across the genome or transcriptome. === 98 98 99 ==== RSeQC: RNA-Seq quality control package for getting mapping statistics (eg. unique/multi-mapped reads) ====99 ==== [http://rseqc.sourceforge.net/ | RSeQC]: RNA-Seq quality control package for getting mapping statistics (eg. unique/multi-mapped reads) ==== 100 100 {{{ 101 101 bam_stat.py -i myFile.bam 102 102 }}} 103 103 104 ==== Picard: CollectRnaSeqMetrics.jar to find coverage across gene body for 5' or 3' bias ====104 ==== [http://broadinstitute.github.io/picard/ | Picard]: CollectRnaSeqMetrics.jar to find coverage across gene body for 5' or 3' bias ==== 105 105 {{{ 106 106 java -jar /usr/local/share/picard-tools/CollectRnaSeqMetrics.jar INPUT=accepted_hits.bam REF_FLAT=refFlat.txt STRAND_SPECIFICITY=NONE OUTPUT=Out_RnaSeqMetrics.txt REFERENCE_SEQUENCE=hg19.fa CHART_OUTPUT=Out_RnaSeqMetrics.pdf 107 107 }}} 108 108 109 ==== QualiMap: can be used on DNA or RNA-Seq to get summary of mapping and coverage/distribution ====109 ==== [http://qualimap.bioinfo.cipf.es/ | QualiMap]: can be used on DNA or RNA-Seq to get summary of mapping and coverage/distribution ==== 110 110 {{{ 111 111 # Graphical interface: enter 'qualimap' on the command line