Entering edit mode
9.5 years ago
clear.choi
▴
30
Hello I have blast search results as below,
BLASTN 2.2.29+
Reference: Zheng Zhang, Scott Schwartz, Lukas Wagner, and Webb
Miller (2000), "A greedy algorithm for aligning DNA sequences", J
Comput Biol 2000; 7(1-2):203-14.
Database: /shared/MiSeq/ABC_Exon23.fasta
7,546 sequences; 4,346,452 total letters
Query= BarcodeS14--S14_Cluster0_Phase0_NumReads236
Length=1006
Score E
Sequences producing significant alignments: (Bits) Value
C_08_94 510 7e-145
C_08_90 510 7e-145
C_08_75 510 7e-145
C_08_74 510 7e-145
C_08_73 510 7e-145
C_08_69 510 7e-145
C_08_68 510 7e-145
C_08_67 510 7e-145
C_08_55N 510 7e-145
C_08_53 510 7e-145
C_08_52N 510 7e-145
C_08_45 510 7e-145
C_08_37 510 7e-145
C_08_28 510 7e-145
C_08_25 510 7e-145
C_08_18 510 7e-145
C_08_17 510 7e-145
C_08_05 510 7e-145
C_08_02_08 510 7e-145
<allele lists>
Lambda K H
1.33 0.621 1.12
Gapped
Lambda K H
1.28 0.460 0.850
Effective search space used: 4101954802
Database: /shared/MiSeq/ABC_Exon23.fasta
Posted date: Jul 22, 2014 5:00 PM
Number of letters in database: 4,346,452
Number of sequences in database: 7,546
I can get only top matching sequence so I can give option like below,
blastn -db .DB/ABC_Exon23.fasta.db -query ./test.fasta -out ./test.out -num_descriptions 1 -num_alignments 1-word_size 10
But then results will be
C_08_94 510 7e-145
Actually I want to get below
C_08_02_08 510 7e-145
So I want to order by ASC of results and get an top one.
PS. Top one has to be top score, order priory 1. Score 2. Name.
Anyone has good idea?
Thanks!