GenomicAlignments and mapping CIGAR coordinates
1
0
Entering edit mode
@stephenrudd-21782
Last seen 2.3 years ago

I am using the GenomicAlignments package for a bioinformatics tutorial that I am preparing. I would like to map coordinates between query reads and reference genome on the basis of the reported CIGAR string.

GenomicAlignments offers the documentation for e.g. queryLoc2refLoc - but the method (as of version 1.21.7) appears not to have been implemented with a rather brutal stop("NOT IMPLEMENTED YET, SORRY!")

Can anyone recommend any other R packages that might help here - the closest method that I have been able to find if the PERL Bio::Cigar (https://github.com/MullinsLab/Bio-Cigar/blob/master/lib/Bio/Cigar.pm)

Thanks

Stephen

cigar not-implemented-yet • 291 views
ADD COMMENT
0
Entering edit mode
@herve-pages-1542
Last seen 1 minute ago
Seattle, WA, United States

Hi,

Take a look at ?mapToAlignments in the GenomicAlignments package.

The reason queryLoc2refLoc() never got implemented is because (1) it's a low-level utility not really intended to be used directly by the end user, and (2) the purpose of the utility was to support higher level (i.e. user facing) mapping functions like mapToAlignments() and family but these functions were finally implemented without the need for queryLoc2refLoc().

Hope this helps,

H.

ADD COMMENT
0
Entering edit mode

I've tried to make the error message more useful in GenomicAlignments 1.21.8:

> queryLoc2refLoc(11:13, "55M")
Error in queryLoc2refLoc(11:13, "55M") : 
  Not implemented yet, sorry! Please use mapToAlignments() and family if
  you need to map coordinates between query reads and reference genome on
  the basis of the reported CIGAR. See '?mapToAlignments'

Best,

H.

ADD REPLY

Login before adding your answer.

Traffic: 360 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6