ILRI Research Computing

User Tools

Site Tools

You are here: start » mkatari-bioinformatics-august-2013-gatknotes
mkatari-bioinformatics-august-2013-gatknotes

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revisionPrevious revision
Next revisionBoth sides next revision
mkatari-bioinformatics-august-2013-gatknotes [2014/06/11 13:36] – mkatarimkatari-bioinformatics-august-2013-gatknotes [2014/06/11 13:46] – mkatari
Line 38: Line 38:
 bowtie2 -x PTC_Human -U Cohen.fastq -S Cohen.sam bowtie2 -x PTC_Human -U Cohen.fastq -S Cohen.sam
 samtools view -bS Cohen.sam > Cohen.bam samtools view -bS Cohen.sam > Cohen.bam
-bowtie2 -x PTC_Human -U Sherman.fastq -S Sherman.sam +/code>
-samtools view -bS Sherman.sam > Sherman.bam +
-</code>+
  
 The picard method to sort is preferred by GATK The picard method to sort is preferred by GATK
Line 49: Line 47:
    SORT_ORDER=coordinate    SORT_ORDER=coordinate
        
-java -jar /export/apps/picard-tools/1.112/SortSam.jar \ +
-   INPUT=Sherman.bam \ +
-   OUTPUT=Sherman.sorted.bam \ +
-   SORT_ORDER=coordinate+
 </code> </code>
  
Line 58: Line 53:
  
 <code> <code>
-java -jar /export/apps/picard-tools/1.112/AddOrReplaceReadGroups.jar \ 
-   INPUT=Sherman.sorted.bam \ 
-   OUTPUT=ShermanRG.bam \ 
-   RGLB=Sherman \ 
-   RGPL=IonTorrent \ 
-   RGPU=None \ 
-   RGSM=Sherman 
  
 java -jar /export/apps/picard-tools/1.112/AddOrReplaceReadGroups.jar \ java -jar /export/apps/picard-tools/1.112/AddOrReplaceReadGroups.jar \
Line 84: Line 72:
    ASSUME_SORTED=TRUE    ASSUME_SORTED=TRUE
  
-java -jar /export/apps/picard-tools/1.112/MarkDuplicates.jar \ 
-   INPUT=ShermanRG.bam \ 
-   OUTPUT=Sherman.dedup.bam \ 
-   METRICS_FILE=Sherman.dedup.metrics \ 
-   REMOVE_DUPLICATES=TRUE \ 
-   ASSUME_SORTED=TRUE 
 </code> </code>
  
Line 95: Line 77:
 <code> <code>
 samtools index Cohen.dedup.bam  samtools index Cohen.dedup.bam 
-samtools index Sherman.dedup.bam  
  
 #identifying indels #identifying indels
Line 104: Line 85:
    -o CohenforIndelRealigner.intervals    -o CohenforIndelRealigner.intervals
    
- #identifying indels 
-java -Xmx2g -jar /export/apps/GenomeAnalysisTK/GenomeAnalysisTK-2.3-9-ge5ebf34/GenomeAnalysisTK.jar \ 
-   -T RealignerTargetCreator \ 
-   -R PTC_Human.fasta \ 
-   -I Sherman.dedup.bam \ 
-   -o ShermanforIndelRealigner.intervals 
  
    
Line 118: Line 93:
   -targetIntervals CohenforIndelRealigner.intervals \   -targetIntervals CohenforIndelRealigner.intervals \
    -o Cohen.dedup.realign.bam    -o Cohen.dedup.realign.bam
- 
- java -Xmx4g -jar /export/apps/GenomeAnalysisTK/GenomeAnalysisTK-2.3-9-ge5ebf34/GenomeAnalysisTK.jar \ 
-   -T IndelRealigner \ 
-   -R PTC_Human.fasta \ 
-   -I Sherman.dedup.bam \ 
-  -targetIntervals ShermanforIndelRealigner.intervals \ 
-   -o Sherman.dedup.realign.bam 
  
 </code> </code>
mkatari-bioinformatics-august-2013-gatknotes.txt · Last modified: 2016/08/17 08:37 by mkatari