This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Nearest genes to a sequence?

Hi all,

Quick query, one can obviously blat/blast short sequences (125bp) to see if they align to a protein coding region - but for those that don't align, is there a way to identify the nearest gene to the inputted sequence?

Cheers all

Jamie

gene blast snp next-gen genome
• 1,945 views • link
written by jt358 • jt358

If you have an alignment and a bed file of the genes then you could use bedtools closest. Link for closest doc.

• 0 views • link
written by GenoMax GenoMax

Hey genomax2 - helpful as always

I can create a BED file in PLINK I assume, but perhaps a silly question, what does bedtools run short sequences against? I see you mention 'if I have an alignment' - does that mean I need an annotated indexed genome of my species?

Cheers

Jamie

• 0 views • link
written by jt358 • jt358

Please use ADD REPLY to answer to earlier comments, as such this thread remains logically structured and easy to follow.

Also note that a bed file from plink isn't the same as a bed file in the rest of bioinformatics. What genomax2 means is an interval file containing chromosome, begin and end of every gene. You can probably convert a gtf/gff file for that.

• 0 views • link
written by WouterDeCoster WouterDeCoster

0 answers

No answers yet.

Log in to answer this question.

More posts like this