This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How To Perform An Exhaustive Search (Report All Hits) Using Blat

I am using Blat to map long contigs to annotated transcripts. But Blat does not report all feasible hits exhaustively. For example, there is a transcript T with length 2280. If I aligned all contigs to it, Blat reports one hit with range [770, 2280]:

Transcript T:   ===============================================================
Hit (contig A):                           -------------------------------------
                                          ^770                                ^2280
               |<-    missing region    <-|

But if I truncate the transcript (keep the first 770 bases only) and blat again, there is an exact hit: 770 matches, 0 mismatch, 0 gap.

Partial T:      ==========================
Hit (contig B): --------------------------
                ^0                       ^770                                                      
                |<-    An exact hit    <-|

Somehow Blat does not report Contig B at first (The missing region). This is not what I want.

I do not want to use Blast, because my program is all using Blat format.

I tried the options -fine and -oneOff=1 to do the same mapping, but the results are the same.

Does anybody know how to enforce exhaustive searching using Blat? Or any script to convert blastn results to PSL format? Thanks.

blat alignment mapping

I don't understand the question. It is not clear from your description what are the differences in the results that you are getting.

@istvan_albert modified. Please help

@istvan_albert modified. Please help

@istvan_albert modified. Please help

Sounds like it might be a multiple hit issue. Please post the raw psl output to your question.

0 answers

No answers yet.

Log in to answer this question.