Question: aligner for intersection of chromosomic regions
gravatar for xcalle91
4.0 years ago by
European Union
xcalle9120 wrote:

I have a list of 10 million elements ( each element is 2 base long) in a bed format. I also have a bed files with all the coding genes (each row been the start and ending points of the genes) I want to check how many elements fall into coding protein genes, an also see how many coding genes have been multiple times hit. For doing so, I wrote a script where using "bedtools: coverage" I get the number of elements in each region.

This works find but the problem is that I have many lists of elements, and bedtools takes forever....

I had the idea of using an alinger for doing so (as they are really fast at mapping). From my experience with STAR I know that is a mapping software but, is it anyway so I could build a reference genome just with coordinates (not with the actual sequence) of the coding genes, and then try to map these bed files.

Longs story short: Can I use an aligner such as STAR to align regions instead of sequences?

sequence alignment • 947 views
ADD COMMENTlink written 4.0 years ago by xcalle9120
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1420 users visited in the last hour