Hi everybody,
I'd like to replace the sequence id for a huge alignment file in .bam format.
I'd rather avoid using the sam format as it would be around 250GB of data, instead I'd like to use samtools c api (writesam/writebam).
I couldn't find much documentation on how to use this library, does anyone know if there is any quick way to replace the read id? I know that the read id is stored in a bam1_t->data array but I cannot understand how to set/replace it.
Any help is much appreciated. Thanks in advance
I'd like to replace the sequence id for a huge alignment file in .bam format.
I'd rather avoid using the sam format as it would be around 250GB of data, instead I'd like to use samtools c api (writesam/writebam).
I couldn't find much documentation on how to use this library, does anyone know if there is any quick way to replace the read id? I know that the read id is stored in a bam1_t->data array but I cannot understand how to set/replace it.
Any help is much appreciated. Thanks in advance
Comment