Entering edit mode
10.4 years ago
stolarek.ir
▴
700
Hi,
I did alignment with cushaw2, however my sam output is missing the cigar string. Is there a way to add it afterwards, or if someone is better with cushaw2, a way to tell it to add directly?
do you have an MD tag?
is your read aligned ? show us the SAM line please.
the read is aligned, holds the NM and AS as the additional fields. I'm trying the fillmd with samtools, maybe this will help. Convertion to bam and viewing in igv is also working normally. I just wanted CIGAR since it's easy for me to parse it for info
problem solved, stupid mistake
can you add the reason as an answer? I am intrigued as to how can one one get a SAM file with no CIGAR strings