Changes between Version 26 and Version 27 of SOPs/SAMBAMqc


Ignore:
Timestamp:
08/20/20 10:58:04 (5 years ago)
Author:
gbell
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • SOPs/SAMBAMqc

    v26 v27  
    7272
    7373==== [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)
     77qualimap
     78
     79# Full command on the command line:
    7780unset DISPLAY  #needed for submitting to cluster
    7881bsub "qualimap bamqc -bam myFile.bam -outdir output_qualimap"
    7982# For huge data, you can increase memory with --java-mem-size="4800M" to avoid OutOfMemoryError: Java heap space
    8083
    81 
    82 #rnaseq qc
     84# For rnaseq QC
    8385bsub "qualimap rnaseq -bam myFile.bam -gtf Homo_sapiens.GRCh37.72.canonical.gtf -outdir output_qualimap_rnaseq -p non-strand-specific"
    8486
    85 #counts qc (after using htseq-count or similar 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)
    8688qualimap counts -d countsqc_input.txt -c -s HUMAN -outdir counts_qc
    8789#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.
     
    9395N82     N8      totalCounts.txt 6
    9496N83     N8      totalCounts.txt 7
    95 
    9697}}}
    9798