Best practice to find break point in fusion genes
0
0
Entering edit mode
2.7 years ago
pablo ▴ 310

Hello,

I work on Pacbio Isoseq data, on A. thalina genome. I need to find the fusion transcripts and the break point of the fusion. I used the SQANTI3 pipeline to detect them : it needs a alignment step with minimap2 , without secondary alignements .

Also, I would like to get the break point (where there is the "link") in these fusion transcripts. It is not implemented yet with these tools.

Do you know a tool which could solve that problem, or do you have an idea?

I thought If I create a sliding window (n=1000nt for example) at the beginning and at the end of the fasta sequence of these transcripts with a python script.

So I will have, for example, a sliding window which creates the fasta sequences (1st nt->1000th nt , 2nd nt ->1001th nt ...) at the beginning, and same way for the end (last nt -1000nt -> last nt , last nt -1001nt -> last nt -1nt) etc.

Then, I blast these window sequences (the target should be different for the beginning and the ending sequences) and I get the breaking point where the beginning target becomes the ending target, and vice versa.

Any help?

fusiongenes isoseq breakpoint pacbio rnaseq • 405 views
ADD COMMENT

Login before adding your answer.

Traffic: 2314 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