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

Both sides previous revisionPrevious revision
Next revisionBoth sides next revision
mkatari-bioinformatics-august-2013-gatknotes [2014/06/13 08:28] mkatarimkatari-bioinformatics-august-2013-gatknotes [2014/07/02 15:19] mkatari
Line 96: Line 96:
 </code> </code>
  
-Now we merge the bam files and then sort and index them+In some cases there may be a need to clean the sam/bam file(s) (soft-trimming the coordinates). To do this use CleanSam in Picard tools. You may want to just do it to all to avoid the error in a workflow, but it may not be necessary.
  
 <code> <code>
-java -jar /export/apps/picard-tools/1.112/MergeSamFiles.jar \+java -jar /export/apps/picard-tools/1.112/CleanSam.jar \
    INPUT=Sherman.dedup.realign.bam \    INPUT=Sherman.dedup.realign.bam \
 +   OUTPUT=Sherman.clean.dedup.realign.bam
 +</code>
 +
 +Now we merge the bam files and then sort and index them. If you cleaned the bam file, remember to use the cleaned ones.
 +
 +<code>
 +java -jar /export/apps/picard-tools/1.112/MergeSamFiles.jar \
 +   INPUT=Sherman.clean.dedup.realign.bam \
    INPUT=Cohen.dedup.realign.bam \    INPUT=Cohen.dedup.realign.bam \
    OUTPUT=ShermanCohenMerged.bam      OUTPUT=ShermanCohenMerged.bam  
Line 110: Line 118:
  
  
-Finall !! run gatk+Finally !! run gatk
  
 <code> <code>
mkatari-bioinformatics-august-2013-gatknotes.txt · Last modified: 2016/08/17 08:37 by mkatari