Open ndrubins opened 2 years ago
Hi, I'm also interested in getting the CIGAR output when using the -split/-splitD argument in bamtobed.
I'm using long reads from ONT with minimap2. Large deletions are detected and I would like to split the alignment in multiple parts at deletion position. Without the CIGAR, i can not determine the original position of fragments in the original read.
Moreover, it would be great if -splitD can accept a threshold in order to split the alignment upon user defined deletion size. Best,
Hi,
I'm running
bamToBed
on an RNA-seqbam
file with the-splitD
argument in order to view each aligned segment in a different line in thebed
file, but I would like to be able to distinguish between ambiguously mapped reads in thebed
file. I thought that adding the-cigar
to the option would do the trick (as it would serve as a read mapping ID), but it seems that it gets ignored when the-splitD
option is used (I do get the cigar in thebed
output of I exclude the-splitD
option).Any idea how to include the
cigar
field in the splitbed
(or any other bam record identified, such as row number in thebam
file)?Seems like this has been asked before but I didn't find a follow up.
Thanks