FreeBSD Manual Pages
samtools-cat(1) Bioinformatics tools samtools-cat(1) NAME samtools cat - concatenate files together SYNOPSIS samtools cat [-b list] [-h header.sam] [-o out.bam] in1.bam in2.bam [ ... ] DESCRIPTION Concatenate BAMs or CRAMs. Although this works on either BAM or CRAM, all input files must be the same format as each other. The sequence dictionary of each input file must be identical, although this command does not check this. This command uses a similar trick to reheader which enables fast BAM concatenation. OPTIONS -b FOFN Read the list of input BAM or CRAM files from FOFN. These are concatenated prior to any files specified on the command line. Multiple -b FOFN options may be specified to concatenate multi- ple lists of BAM/CRAM files. -h FILE Uses the SAM header from FILE. By default the header is taken from the first file to be concatenated. -o FILE Write the concatenated output to FILE. By default this is sent to stdout. -q [CRAM only] Query the number of containers in the CRAM file. The output is the filename, the number of containers, and the first and last container number as an inclusive range, with one file per line. Note this works in conjunction with the -r RANGE option, in which case the 3rd and 4th columns become useful for identify- ing which containers span the requested range. -r RANGE [CRAM only] Filter the CRAM file to a specific RANGE. This can be the usual chromosome:start-end syntax, or "*" for unmapped records at the end of alignments. If the range is of the form "#:start-end" then the start and end coordinates are interpreted as inclusive CRAM container numbers, starting at 0 and ending 1 less than the number of containers reported by -q. For example -r "#:0-9" is the first 10 CRAM containers of data. All range types filter data in as fast a manner as possible, using operating system read/write loops where appropriate. -p A/B [CRAM only] Filter the CRAM file using a specific fraction. The file is split into B approximately equal parts and returns element A where A is between 1 and B inclusive. If there are more parts specified than CRAM containers then some of the out- put will be empty CRAMs. This can also be combined with the range option above to oper- ate of parts of that range. For example -r chr2 -p 1/10 re- turns the first 1/10th of data aligned against chromosome 2. -f [CRAM only] Enable fast mode. When filtering by chromosome range with -r we normally do careful recoding of any containers that overlap the start and end of the range so the record count precisely matches that returned by a samtools view equivalent. Fast mode does no filtering, so may return additional align- ments in the same container but outside of the requested re- gion. --no-PG Do not add a @PG line to the header of the output file. -@, --threads INT Number of input/output compression threads to use in addition to main thread [0]. EXAMPLES o Extract a specific chromosome from a CRAM file, outputting to a new CRAM. samtools cat -o chr10.cram -r chr10 in.cram o Split a CRAM file up into separate files, each containing at most 123 containers. set -- $(samtools cat -q in.cram); nc=$2; s=0 while [ $s -lt $nc ] do e=`expr $s + 123` if [ $e -ge $nc ] then e=$nc fi r="$s-`expr $e - 1`"; echo $r fn=/tmp/_part-`printf "%08d" $s`.cram samtools cat -o $fn in.cram -r "#:$r" s=$e done o Split any unaligned data from a (potentially aligned) CRAM file into 10 approximately equal sized pieces. for i in `seq 1 10` do samtools cat in.cram -r "*" -p $i/10 -o part-`printf "%02d" $i`.cram done AUTHOR Written by Heng Li from the Sanger Institute. Updated for CRAM by James Bonfield (also Sanger Institute). SEE ALSO samtools(1) Samtools website: <http://www.htslib.org/> samtools-1.21 12 September 2024 samtools-cat(1)
NAME | SYNOPSIS | DESCRIPTION | OPTIONS | EXAMPLES | AUTHOR | SEE ALSO
Want to link to this manual page? Use this URL:
<https://man.freebsd.org/cgi/man.cgi?query=samtools-cat&sektion=1&manpath=FreeBSD+Ports+14.3.quarterly>