3,045
edits
Changes
From Genome Analysis Wiki
Jump to navigationJump to searchUpdated to actually split all chromsomes
The <code>splitChromosome</code> option on the [[bamUtil]] executable splits an indexed BAM file into multiple files based on the Chromosome (Reference Name).
The <code>splitChromosome</code> option on the [[bamUtil]] executable splits an indexed BAM file into multiple files based on the Chromosome (Reference Name).
The files all have the same base name, but with an _# where # corresponds with the associated reference id from the BAM file.
The files all have the same base name, but with the chromosome name ".bam" or ".sam" appended.
= Parameters =
= Parameters =
--in : the BAM file to be split
--in : the BAM file to be split
--out : the base filename for the SAM/BAM files to write into. Does not include the extension.
--out : the base filename for the SAM/BAM files to write into. Does not include the extension.
CHROM.bam or CHROM.sam will be appended to the basename where CHROM is the chromosome name.
Optional Parameters:
Optional Parameters:
--noeof : do not expect an EOF block on a bam file.
--noeof : do not expect an EOF block on a bam file.
--bamout : write the output files in BAM format (default).
--bamout : write the output files in BAM format (default).
--samout : write the output files in SAM format.
--samout : write the output files in SAM format.
= Usage =
= Usage =
./bam splitChromosome --in <inputFilename> --out <outputFileBaseName> [--bamIndex <bamIndexFile>] [--noeof] [--bamout|--samout] [--params]
./bam splitChromosome --in <inputFilename> --out <outputFileBaseName> [--noeof] [--bamout|--samout] [--params]
= Example Output =
= Example Output =
<pre>
<pre>
Reference ID -1 has 2 records
Reference Name: 1 has 5 records
Reference Name: 2 has 2 records
Reference ID 1 has 2 records
Reference Name: 3 has 1 records
Reference ID 2 has 1 records
Reference Name: * has 2 records
Reference ID 3 has 0 records
Number of records = 10
Number of records = 10
Returning: 0 (SUCCESS)
Returning: 0 (SUCCESS)
</pre>
</pre>
