Changes

From Genome Analysis Wiki
Jump to navigationJump to search
302 bytes added ,  11:56, 6 March 2012
no edit summary
Line 8: Line 8:     
It has options to allow for the conversion of the sequence to/from '=' from/to the actual bases by using the reference sequence.  
 
It has options to allow for the conversion of the sequence to/from '=' from/to the actual bases by using the reference sequence.  
 +
 +
It also has an option to left shift indels in the CIGARs before writing the output file.
    
If you want to convert a BAM file to a SAM file, just call:  
 
If you want to convert a BAM file to a SAM file, just call:  
Line 18: Line 20:  
= Usage  =
 
= Usage  =
   −
  ./bam convert --in <inputFile> --out <outputFile.sam/bam/ubam (ubam is uncompressed bam)> [--refFile <reference filename>] [--seqBases|--seqEquals|--seqOrig] [--recover] [--noeof] [--params]
+
  ./bam convert --in <inputFile> --out <outputFile.sam/bam/ubam (ubam is uncompressed bam)> [--refFile <reference filename>] [--useBases|--useEquals|--useOrigSeq] [--lshift] [--noeof] [--params]
      Line 27: Line 29:  
     Optional Parameters:
 
     Optional Parameters:
 
--refFile  &nbsp;: reference file name
 
--refFile  &nbsp;: reference file name
         --noeof    &nbsp;: do not expect an EOF block on a bam file.
+
--lshift  &nbsp;: left shift indels when writing records
 +
         --noeof    &nbsp;: do not expect an EOF block on a bam file
 
         --params  &nbsp;: print the parameter settings
 
         --params  &nbsp;: print the parameter settings
         --recover  &nbsp;: attempt to recover the input bam file.
+
         --recover  &nbsp;: attempt to recover the input bam file
 
     Optional Sequence Parameters (only specify one):
 
     Optional Sequence Parameters (only specify one):
--seqOrig  &nbsp;: Leave the sequence as is (default &amp; used if reference is not specified).
+
--seqOrig  &nbsp;: Leave the sequence as is (default &amp; used if reference is not specified)
--seqBases &nbsp;: Convert any '=' in the sequence to the appropriate base using the reference (requires --ref).
+
--seqBases &nbsp;: Convert any '=' in the sequence to the appropriate base using the reference (requires --ref)
--seqEquals&nbsp;: Convert any bases that match the reference to '=' (requires --ref).
+
--seqEquals&nbsp;: Convert any bases that match the reference to '=' (requires --ref)
 
</pre>  
 
</pre>  
    
{{InBAMInputFile}}
 
{{InBAMInputFile}}
 +
{{OutBAMOutputFile}}
 +
{{refFile}}
 +
 +
== Left Shift Indels in the CIGAR (<code>--lshift</code>) ==
 +
 +
Left shift indels as far as they can go in the read.
 +
 +
 +
{{noeofBGZFParameter}}
 +
{{paramsParameter}}
   −
{{OutBAMOutputFile}}
      
== Recover a corrupted BAM file (<code>--recover</code>) ==
 
== Recover a corrupted BAM file (<code>--recover</code>) ==

Navigation menu