TBLASTN 2.2.13 [Nov-27-2005] Reference: Altschul, Stephen F., Thomas L. Madden, Alejandro A. Schaffer, Jinghui Zhang, Zheng Zhang, Webb Miller, and David J. Lipman (1997), "Gapped BLAST and PSI-BLAST: a new generation of protein database search programs", Nucleic Acids Res. 25:3389-3402. Query= SPP00000001_1.0 # Protein # Iodothyronine deiodinase 1 (DI1) # Homo sapiens # Complete (249 letters) Database: /users/rg/didac/GENOMES/A.rara/genome.fa 3231 sequences; 1,450,095 total letters Searching..................................................done Score E Sequences producing significant alignments: (bits) Value gb|ADNL01000617.1| Astrammina rara contig00691, whole genome sho... 27 0.50 gb|ADNL01001134.1| Astrammina rara contig01331, whole genome sho... 23 5.5 gb|ADNL01000769.1| Astrammina rara contig00880, whole genome sho... 23 7.2 gb|ADNL01000200.1| Astrammina rara contig00210, whole genome sho... 23 7.2 gb|ADNL01002808.1| Astrammina rara contig03299, whole genome sho... 22 9.4 gb|ADNL01001283.1| Astrammina rara contig01512, whole genome sho... 22 9.4 >gb|ADNL01000617.1| Astrammina rara contig00691, whole genome shotgun sequence Length = 644 Score = 26.6 bits (57), Expect = 0.50 Identities = 12/33 (36%), Positives = 16/33 (48%), Gaps = 4/33 (12%) Frame = -2 Query: 65 FFSTQYFWFVLKVRWQRL----EDTTELGGLAP 93 FF +QYFW+ + + W L T G L P Sbjct: 232 FFQSQYFWYCVVLEWTYLSIFIHQRTLFGDLIP 134 >gb|ADNL01001134.1| Astrammina rara contig01331, whole genome shotgun sequence Length = 241 Score = 23.1 bits (48), Expect = 5.5 Identities = 13/46 (28%), Positives = 23/46 (50%), Gaps = 5/46 (10%) Frame = -3 Query: 38 KRNILAMGEKT-----GMTRNPHFSHDNWIPTFFSTQYFWFVLKVR 78 K N+L+ + T T F+H + T+F+T F +V+ +R Sbjct: 179 KDNLLSRADHTTHVTLNFTTIATFNHATFFLTYFTTFSFTYVISLR 42 >gb|ADNL01000769.1| Astrammina rara contig00880, whole genome shotgun sequence Length = 406 Score = 22.7 bits (47), Expect = 7.2 Identities = 8/17 (47%), Positives = 10/17 (58%) Frame = +1 Query: 105 CNIWEFMQGNRPLVLNF 121 CNIW F+ P +L F Sbjct: 124 CNIWYFLNNXTPALL*F 174 >gb|ADNL01000200.1| Astrammina rara contig00210, whole genome shotgun sequence Length = 1151 Score = 22.7 bits (47), Expect = 7.2 Identities = 10/26 (38%), Positives = 12/26 (46%), Gaps = 5/26 (19%) Frame = -1 Query: 89 GGLAPNCPVVRLSGQRCN-----IWE 109 G +CP V G RCN +WE Sbjct: 113 GSCGYHCPAVVFRGSRCNWSMNGLWE 36 >gb|ADNL01002808.1| Astrammina rara contig03299, whole genome shotgun sequence Length = 743 Score = 22.3 bits (46), Expect = 9.4 Identities = 10/24 (41%), Positives = 12/24 (50%) Frame = +3 Query: 190 RSPQCPVVVDTMQNQSSQLYAALP 213 R PQ PVVV + +LY P Sbjct: 657 RPPQSPVVVGRAWQRREELYGESP 728 >gb|ADNL01001283.1| Astrammina rara contig01512, whole genome shotgun sequence Length = 522 Score = 22.3 bits (46), Expect = 9.4 Identities = 16/45 (35%), Positives = 18/45 (40%), Gaps = 4/45 (8%) Frame = -2 Query: 72 WFVLKVRWQRLEDTTELGGLAPNCPVVRLSGQRCN----IWEFMQ 112 W K+RWQRL T G P C +CN W F Q Sbjct: 143 WQARKLRWQRL---TVSGVAYPTCKC------KCNQDEDEWHFAQ 36 Database: /users/rg/didac/GENOMES/A.rara/genome.fa Posted date: Sep 29, 2011 8:09 PM Number of letters in database: 1,450,095 Number of sequences in database: 3231 Lambda K H 0.324 0.138 0.448 Gapped Lambda K H 0.267 0.0410 0.140 Matrix: BLOSUM62 Gap Penalties: Existence: 11, Extension: 1 Number of Hits to DB: 435,840 Number of Sequences: 3231 Number of extensions: 5582 Number of successful extensions: 29 Number of sequences better than 10.0: 6 Number of HSP's better than 10.0 without gapping: 29 Number of HSP's successfully gapped in prelim test: 0 Number of HSP's that attempted gapping in prelim test: 0 Number of HSP's gapped (non-prelim): 29 length of query: 249 length of database: 483,365 effective HSP length: 72 effective length of query: 177 effective length of database: 250,733 effective search space: 44379741 effective search space used: 44379741 frameshift window, decay const: 40, 0.1 T: 13 A: 40 X1: 15 ( 7.0 bits) X2: 38 (14.6 bits) X3: 64 (24.7 bits) S1: 40 (21.5 bits) S2: 46 (22.3 bits)