User Tools

Site Tools


mkatari-bioinformatics-august-2013-gatknotes

Differences

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

Link to this comparison view

Next revisionBoth sides next revision
mkatari-bioinformatics-august-2013-gatknotes [2014/06/11 12:49] – created mkatarimkatari-bioinformatics-august-2013-gatknotes [2014/06/11 12:52] mkatari
Line 15: Line 15:
 bowtie2-build PTC_Human.fasta PTC_Human bowtie2-build PTC_Human.fasta PTC_Human
 samtools faidx PTC_Human.fasta samtools faidx PTC_Human.fasta
-java -jar /export/apps/picard-tools/1.112/CreateSequenceDictionary.jar R=PTC_Human.fasta O=PTC_Human.dict+java -jar /export/apps/picard-tools/1.112/CreateSequenceDictionary.jar 
 +   R=PTC_Human.fasta 
 +   O=PTC_Human.dict
  
 </code> </code>
Line 40: Line 42:
 The picard method to sort is preferred by GATK The picard method to sort is preferred by GATK
 <code> <code>
-java -jar /export/apps/picard-tools/1.112/SortSam.jar INPUT=Cohen.bam OUTPUT=Cohen.sorted.bam SORT_ORDER=coordinate +java -jar /export/apps/picard-tools/1.112/SortSam.jar 
-java -jar /export/apps/picard-tools/1.112/SortSam.jar INPUT=Sherman.bam OUTPUT=Sherman.sorted.bam SORT_ORDER=coordinate+   INPUT=Cohen.bam 
 +   OUTPUT=Cohen.sorted.bam 
 +   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 47: Line 56:
  
 <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 
 +   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 INPUT=Cohen.sorted.bam OUTPUT=CohenRG.bam RGLB=Cohen RGPL=IonTorrent RGPU=None RGSM=Cohen+java -jar /export/apps/picard-tools/1.112/AddOrReplaceReadGroups.jar 
 +   INPUT=Cohen.sorted.bam 
 +   OUTPUT=CohenRG.bam 
 +   RGLB=Cohen 
 +   RGPL=IonTorrent 
 +   RGPU=None 
 +   RGSM=Cohen
 </code> </code>
  
 This will remove any reads that map to the same exact place. It is helpful to get rid of artifacts.  This will remove any reads that map to the same exact place. It is helpful to get rid of artifacts. 
 <code> <code>
-java -jar /export/apps/picard-tools/1.112/MarkDuplicates.jar INPUT=CohenRG.bam OUTPUT=Cohen.dedup.bam METRICS_FILE=Cohen.dedup.metrics REMOVE_DUPLICATES=TRUE ASSUME_SORTED=TRUE+java -jar /export/apps/picard-tools/1.112/MarkDuplicates.jar 
 +   INPUT=CohenRG.bam 
 +   OUTPUT=Cohen.dedup.bam 
 +   METRICS_FILE=Cohen.dedup.metrics 
 +   REMOVE_DUPLICATES=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+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 98: Line 129:
  
 <code> <code>
-java -jar /export/apps/picard-tools/1.112/MergeSamFiles.jar INPUT=Sherman.dedup.realign.bam INPUT=Cohen.dedup.realign.bam OUTPUT=ShermanCohenMerged.bam  +java -jar /export/apps/picard-tools/1.112/MergeSamFiles.jar 
 +   INPUT=Sherman.dedup.realign.bam 
 +   INPUT=Cohen.dedup.realign.bam 
 +   OUTPUT=ShermanCohenMerged.bam   
 samtools sort ShermanCohenMerged.bam ShermanCohenMerged.sorted samtools sort ShermanCohenMerged.bam ShermanCohenMerged.sorted
 +
 samtools index ShermanCohenMerged.sorted.bam  samtools index ShermanCohenMerged.sorted.bam 
 </code> </code>
mkatari-bioinformatics-august-2013-gatknotes.txt · Last modified: 2016/08/17 08:37 by mkatari