samtools depad - convert padded BAM to unpadded BAM
samtools depad [-SsCu1] [-T ref.fa] [-o
output] in.bam
Converts a BAM aligned against a padded reference to a BAM aligned against the
depadded reference. The padded reference may contain verbatim "*"
bases in it, but "*" bases are also counted in the reference
numbering. This means that a sequence base-call aligned against a reference
"*" is considered to be a cigar match ("M" or
"X") operator (if the base-call is "A", "C",
"G" or "T"). After depadding the reference "*"
bases are deleted and such aligned sequence base-calls become insertions.
Similarly transformations apply for deletions and padding cigar operations.
- -S
- Ignored for compatibility with previous samtools versions. Previously this
option was required if input was in SAM format, but now the correct format
is automatically detected by examining the first few characters of
input.
- -s
- Output in SAM format. The default is BAM.
- -C
- Output in CRAM format. The default is BAM.
- -u
- Do not compress the output. Applies to either BAM or CRAM output
format.
- -1
- Enable fastest compression level. Only works for BAM or CRAM output.
- -T FILE
- Provides the padded reference file. Note that without this the @SQ line
lengths will be incorrect, so for most use cases this option will be
considered as mandatory.
- -o FILE
- Specifies the output filename. By default output is sent to stdout.
- --no-PG
- Do not add a @PG line to the header of the output file.
Written by Heng Li from the Sanger Institute with extensions by Peter Cock from
the James Hutton Institute.
samtools(1)
Samtools website: <http://www.htslib.org/>