htgt / CRISPR-Analyser

C++ package for analysing CRISPR off targets
MIT License
21 stars 8 forks source link

Retreive CRISPR sequence and location using ID #6

Open fh-fingolfin opened 9 years ago

fh-fingolfin commented 9 years ago

The search option gives me an ID for a given gRNA sequence. Is there also a way to retreive the CRISPR sequence and location computed in the get_all_crisprs step?

Also, is there a way to retreive the same thing (sequence + location) in the off-target search step? Similar to the Off-Target Sites Table on the WGE CRISPR summary page?

ah19 commented 9 years ago

Unfortunately the location of each crispr isn't stored in the index so there's no way of returning it; the simplest way around this currently is to put the data from the get_all_crisprs step into a database and pull up the locations from there (this is how we do it in WGE). Adding the table also means you can store off targets

The code could also be modified to include the location of each crispr in the index, but this would increase the size of each index by 50% or so. I may add an option to do this in the future, which would mean the server could be run standalone

coronin commented 9 years ago

I modifed the original code to output gRNA sequence, with which you can retrieve the position in genome.

https://github.com/coronin/CRISPR-Analyser/commit/594e41d31a8f022f10ce2082a6abc7a7bb6e5f23

ah19 commented 9 years ago

Coronin: The API method to give a gRNA sequence for a given id that you wrote has now been added to the master branch, thanks for the contribution!

It is now also no longer necessary to use the get_all_crisprs script, it has been integrated into crispr_analyser