Changes between Version 26 and Version 27 of SOPs/SAMBAMqc
- Timestamp:
- 08/20/20 10:58:04 (5 years ago)
Legend:
- Unmodified
- Added
- Removed
- Modified
-
SOPs/SAMBAMqc
v26 v27 72 72 73 73 ==== [http://qualimap.bioinfo.cipf.es/ QualiMap] can be used on DNA or RNA-Seq to get summary of mapping and coverage/distribution ==== 74 {{{ 75 # Graphical interface: enter 'qualimap' on the command line 76 # Command line: 74 75 {{{ 76 # For Graphical interface (with GUI access to all tools) 77 qualimap 78 79 # Full command on the command line: 77 80 unset DISPLAY #needed for submitting to cluster 78 81 bsub "qualimap bamqc -bam myFile.bam -outdir output_qualimap" 79 82 # For huge data, you can increase memory with --java-mem-size="4800M" to avoid OutOfMemoryError: Java heap space 80 83 81 82 #rnaseq qc 84 # For rnaseq QC 83 85 bsub "qualimap rnaseq -bam myFile.bam -gtf Homo_sapiens.GRCh37.72.canonical.gtf -outdir output_qualimap_rnaseq -p non-strand-specific" 84 86 85 # counts qc (after using htseq-count orsimilar program to generate a matrix of counts)87 # For counts QC (after using htseq-count or a similar program to generate a matrix of counts) 86 88 qualimap counts -d countsqc_input.txt -c -s HUMAN -outdir counts_qc 87 89 #Format of countsqc_input.txt (below), totalCounts.txt is a matrix of counts; header lines must be commented "#" and species is human or mouse only. … … 93 95 N82 N8 totalCounts.txt 6 94 96 N83 N8 totalCounts.txt 7 95 96 97 }}} 97 98