I have a .sam file (for single-end reads) and I wish to find the end position of each read. I noticed that I can do that using Rsamtools if the file is a .bam file. However, I wish to do it directly on .sam file. How should I do? (using cigar value?)
Thanks!
One more question: what the use of N.regions.removed argument in the function GenomicAlignments::cigarWidthAlongReferenceSpace() ?