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= SPP00000048_1.0 # Protein # Selenophosphate synthetase 2 (SPS2) # Drosophila melanogaster # Complete (370 letters) Database: genome.fa 697 sequences; 71,881,135 total letters Searching.....................................................done Score E Sequences producing significant alignments: (bits) Value scaffold_36 133 4e-31 scaffold_84 32 1.7 scaffold_175 30 6.4 scaffold_70 30 8.3 scaffold_78 30 8.3 >scaffold_36 Length = 507435 Score = 133 bits (335), Expect = 4e-31 Identities = 82/222 (36%), Positives = 124/222 (55%), Gaps = 6/222 (2%) Frame = -3 Query: 153 IRQLKINPWCII------GGIATSVCRSEEIILPSNAQPGDVLVLTKPLGGQMAMDAHLW 206 IR + + P CI GG+ + ++ I+P+ A+ GD LVLTKPLG Q+A++ + W Sbjct: 409588 IRLILLGPKCINLFRFIRGGVGIASNHKDDFIMPNGAEAGDYLVLTKPLGVQLAVNTYQW 409409 Query: 207 QLNQTEKYKKLLSECSDADIKETFEIAVKSMTYLNKNAALLMHKYQAHCATDITGFGLLG 266 K+ + ++++I F+ A KSM LN AA L K+ +TD+TGFG+LG Sbjct: 409408 MTTIPSKWDAIKGITNESEIN*AFDHATKSMCTLNAAAASLFKKHNIKASTDVTGFGILG 409229 Query: 267 HANNLAQFQKEKVLFQINKLPIIKNVLKFSTLVGQSTKFRSGRSVETSGGLLICLPADAA 326 HA LA+ + KV F I+KLPII+N+ K + +F G + ETSGGLLI Sbjct: 409228 HARYLARA*RRKV*FVIDKLPIIQNLYKLDK-KARDFRFLEGLAAETSGGLLIA--TKDP 409058 Query: 327 DKFCRDFEEATNGEQKSFQIGHVTAANESDAVLCEDVEFIEV 368 D F +D++ + + IG V + E DAV+ ++V+ EV Sbjct: 409057 DSFLKDYKHLIG--EWGWVIGRVEKSQERDAVISKNVKIQEV 408938 Score = 100 bits (248), Expect = 5e-21 Identities = 55/139 (39%), Positives = 81/139 (58%) Frame = -2 Query: 4 PEKHGLEPDFQLTKFTTHTGUSCKIPQKVLEKYLRGTEIENKNNDGYLIGSGMDCAVIPL 63 PE +GL F +T FT G CK+P L KYL +G + D +V P+ Sbjct: 410570 PENYGLPKGFAMTNFTELKGCGCKVP*T*LLKYLG*I------GEGDIGKETPDVSVYPV 410409 Query: 64 KRHKDYLLIQTVDFFYPMVNDPELLGRIALANVLSDVYAVGVTQFDTVEMIVSTSTSFSE 123 KR D+ I T+D+FYP+V +P GRIA NVLSD+YA+GVTQ DT+ M++ T ++ Sbjct: 410408 KRT*DFHTISTIDYFYPLVANPYA*GRIACCNVLSDLYAMGVTQVDTMLMVLGVCTKMTD 410229 Query: 124 KERDVVIGLVMKGFQNSLK 142 +E+ + L+++GF + K Sbjct: 410228 EEK*ISTSLMLEGFNDCAK 410172 >scaffold_84 Length = 354111 Score = 32.0 bits (71), Expect = 1.7 Identities = 35/133 (26%), Positives = 56/133 (42%), Gaps = 6/133 (4%) Frame = +2 Query: 177 IILPSNAQPGDVLVLTKPLGGQMAMDAHLWQLNQTEKYKKLLSECSDADIKETFEIAVKS 236 ++ P N D+ +L K + Q E K+ L S A+ E FE +KS Sbjct: 44228 LVNPRNCVCIDIFILMKSANSNLPNSIQPDQNEINEYIKEYLRYSSFANTLECFEAEIKS 44407 Query: 237 MTYLNKNAALLMHKYQA---HCATDITGFGLLGHANNLA---QFQKEKVLFQINKLPIIK 290 NK L+ ++ + + G L + N L + KEK +++ NK P IK Sbjct: 44408 K*VSNKVFICLILPFRCLINNSSRSK*GMRLQDYFNCLRVTIRKLKEKSIWRRNKRPSIK 44587 Query: 291 NVLKFSTLVGQST 303 N+ K L +S+ Sbjct: 44588 NINKSYRLEDRSS 44626 >scaffold_175 Length = 139747 Score = 30.0 bits (66), Expect = 6.4 Identities = 17/47 (36%), Positives = 26/47 (55%), Gaps = 1/47 (2%) Frame = -1 Query: 271 LAQFQKEKVLFQINKLPIIKNVLKFSTL-VGQSTKFRSGRSVETSGG 316 L QF ++ VL +INK+ IIKN F L + + +F+ G + G Sbjct: 28144 LCQFYEKTVL*EINKIVIIKNERLFRMLEINKKERFKKGFRINNKYG 28004 >scaffold_70 Length = 388433 Score = 29.6 bits (65), Expect = 8.3 Identities = 18/41 (43%), Positives = 22/41 (53%) Frame = -1 Query: 220 ECSDADIKETFEIAVKSMTYLNKNAALLMHKYQAHCATDIT 260 +C + D + I VK LNK+A LL YQA A DIT Sbjct: 21944 KCRNQDCGDNGIIEVKIDKTLNKDALLLKKDYQAAKAIDIT 21822 >scaffold_78 Length = 358839 Score = 29.6 bits (65), Expect = 8.3 Identities = 18/41 (43%), Positives = 22/41 (53%) Frame = -3 Query: 220 ECSDADIKETFEIAVKSMTYLNKNAALLMHKYQAHCATDIT 260 +C + D + I VK LNK+A LL YQA A DIT Sbjct: 21196 KCRNQDCGDNGIIEVKIDKTLNKDALLLKKDYQAAKAIDIT 21074 Database: genome.fa Posted date: Feb 10, 2010 10:01 AM Number of letters in database: 71,881,135 Number of sequences in database: 697 Lambda K H 0.319 0.135 0.394 Gapped Lambda K H 0.267 0.0410 0.140 Matrix: BLOSUM62 Gap Penalties: Existence: 11, Extension: 1 Number of Hits to DB: 24,967,756 Number of Sequences: 697 Number of extensions: 336167 Number of successful extensions: 1145 Number of sequences better than 10.0: 5 Number of HSP's better than 10.0 without gapping: 23 Number of HSP's successfully gapped in prelim test: 10 Number of HSP's that attempted gapping in prelim test: 1084 Number of HSP's gapped (non-prelim): 148 length of query: 370 length of database: 23,960,378 effective HSP length: 107 effective length of query: 263 effective length of database: 23,885,799 effective search space: 6281965137 effective search space used: 6281965137 frameshift window, decay const: 40, 0.1 T: 13 A: 40 X1: 16 ( 7.4 bits) X2: 38 (14.6 bits) X3: 64 (24.7 bits) S1: 41 (21.8 bits) S2: 64 (29.3 bits)