Changes between Version 21 and Version 22 of SOP/CallingVariantsRNAseq


Ignore:
Timestamp:
08/16/17 12:45:20 (7 years ago)
Author:
krichard
Comment:

--

Legend:

Unmodified
Added
Removed
Modified
  • SOP/CallingVariantsRNAseq

    v21 v22  
    3333}}}
    3434
     352 - '''Index the reference genome for Second Pass.'''
     36        Create folder, "SecondPass" before running these commands.
     37
     38To generate the 2nd pass genome index files for STAR:
     39{{{
     40
     41bsub STAR --runMode genomeGenerate --genomeDir /path/to/GenomeDir2nd --genomeFastaFiles /path/to/genome/fasta1 --sjdbFileChrStartEnd /path/to/first/pass/directory/SJ.out.tab --sjdbGTFfile /path/to/GTF/FileName.gtf --sjdbOverhang 100 --runThreadN 8
     42}}}
     43
     44
    3545The parameters included in the above sample commands are:
    3646  * '''--runMode <alignReads, genomeGenerate>'''   "alignReads" does the actual mapping. "genomeGenerate" generates the genomeDir required for mapping (default = alignReads).
     
    3949  * '''--readFilesIn <read1.fastq read2.fastq> ''' Specifies the fastq files containing the reads, can be single-end or paired-end.
    4050  * '''--runThreadN <n> ''' Specifies the number of threads to use.
     51
     52
     53