Alineamiento de B.Bigemina contra selenoproteínas de Drosophila


SLA/LP

	TBLASTN 2.2.17 [Aug-26-2007]


	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.

	Reference for composition-based statistics:
	Schaffer, Alejandro A., L. Aravind, Thomas L. Madden,
	Sergei Shavirin, John L. Spouge, Yuri I. Wolf,  
	Eugene V. Koonin, and Stephen F. Altschul (2001), 
	"Improving the accuracy of PSI-BLAST protein database searches with 
	composition-based statistics and other refinements",  Nucleic Acids Res. 29:2994-3005.

	Query= gi|45550694|ref|NP_649556.2| CG1427 CG1427-PA, isoform A
	[Drosophila melanog
	         (448 letters)

	Database: B_bigemina.contigs.fa 
	           4396 sequences; 19,733,690 total letters

	Searching..................................................done



	                                                                 Score    E
	Sequences producing significant alignments:                      (bits) Value

	Contig3744                                                             33   0.31 
	Contig3897                                                             33   0.33 
	Contig3683.0                                                           32   0.60 
	Contig2932                                                             32   0.79 
	Contig854                                                              32   0.81 
	Contig1362                                                             31   0.88 
	Contig846                                                              31   0.97 

	>Contig3744
	          Length = 5928

	 Score = 32.7 bits (73), Expect = 0.31,   Method: Composition-based stats.
	 Identities = 31/90 (34%), Positives = 39/90 (43%), Gaps = 9/90 (10%)
	 Frame = +3

	Query: 113  LLARLTNALI--LDLIRGIGLPSC-----AGCFL--VPMCTGMTLTLCLQSLRKRRPGAR 163
	            L+ARLT  L   L L  G  LP C     A C +  + +C G +L LCL +     PGAR
	Sbjct: 2316 LVARLTLRLCASLALCLGARLPLCLVASLARCLVASLALCPGASLALCLVACLALWPGAR 2495

	Query: 164  YVLWSRIDQKSCFKAITATGLVPVVIPCLI 193
	              L+       C  A     LV     CL+
	Sbjct: 2496 LALFLVARLARCLVARLTRCLVASSTRCLV 2585


	>Contig3897
	          Length = 12828

	 Score = 32.7 bits (73), Expect = 0.33,   Method: Composition-based stats.
	 Identities = 25/68 (36%), Positives = 29/68 (42%), Gaps = 7/68 (10%)
	 Frame = -2

	Query: 123  LDLIRGIGLPSCAGCFL-------VPMCTGMTLTLCLQSLRKRRPGARYVLWSRIDQKSC 175
	            L L  G  LP C G  L       + +C G  L LCL +   R  GAR  L S +D    
	Sbjct: 5594 LALCLGARLPLCLGACLALCLVASLALCLGARLALCLGARSTRCLGARLALSSIVDSSLN 5415

	Query: 176  FKAITATG 183
	              AI  TG
	Sbjct: 5414 LNAIQFTG 5391


	>Contig3683.0
	          Length = 4137

	 Score = 32.0 bits (71), Expect = 0.60,   Method: Composition-based stats.
	 Identities = 25/69 (36%), Positives = 31/69 (44%)
	 Frame = -3

	Query: 131  LPSCAGCFLVPMCTGMTLTLCLQSLRKRRPGARYVLWSRIDQKSCFKAITATGLVPVVIP 190
	            LP C  C  +P C G TLTLCL        GAR  L        C +A  A  LV  +  
	Sbjct: 3607 LPLCL-CASLPRCLGATLTLCL--------GARSTLCLGASLPLCLRARLAPCLVARLAL 3455

	Query: 191  CLIKGESLN 199
	            CL+   +L+
	Sbjct: 3454 CLVASLALS 3428


	>Contig2932
	          Length = 11681

	 Score = 31.6 bits (70), Expect = 0.79,   Method: Composition-based stats.
	 Identities = 34/106 (32%), Positives = 42/106 (39%), Gaps = 5/106 (4%)
	 Frame = -1

	Query: 132  PSCAGCFLVPMCTGMTLTLCLQSLRKRRPGARYVLWSRIDQKSCFKAITATGLVPVVIPC 191
	            P C G     +C G TLT CL +      GA   L        C  A +A  LV  +  C
	Sbjct: 6803 PLCLGAS-SALCLGATLTRCLVASLALCLGASSALCLVARSPRCLSARSALSLVASLPLC 6627

	Query: 192  LIKGESLNTNVDLFREKIKS----LGVDSILCL-YTTTSCFAPRNS 232
	            L    +L     L R  + S    LG  S LCL  T T C   R++
	Sbjct: 6626 LDASLALCLGATLTRCLVASLALCLGASSALCLGATLTLCLGARST 6489


	>Contig854
	          Length = 2089

	 Score = 31.6 bits (70), Expect = 0.81,   Method: Composition-based stats.
	 Identities = 31/84 (36%), Positives = 37/84 (44%), Gaps = 2/84 (2%)
	 Frame = -3

	Query: 112 TLLARLTNALI--LDLIRGIGLPSCAGCFLVPMCTGMTLTLCLQSLRKRRPGARYVLWSR 169
	           +LLARLT  L+  L L     L  C G  L P+C G   TLCL +      GAR  L   
	Sbjct: 782 SLLARLTRCLVASLALCLCARLALCLGASL-PLCLGAWSTLCLGASLALCLGARLALCLG 606

	Query: 170 IDQKSCFKAITATGLVPVVIPCLI 193
	                C  A  A  LV  +  CL+
	Sbjct: 605 ASLPRCLVARLALCLVARLPRCLV 534


	>Contig1362
	          Length = 9418

	 Score = 31.2 bits (69), Expect = 0.88,   Method: Composition-based stats.
	 Identities = 21/67 (31%), Positives = 30/67 (44%)
	 Frame = -3

	Query: 132  PSCAGCFLVPMCTGMTLTLCLQSLRKRRPGARYVLWSRIDQKSCFKAITATGLVPVVIPC 191
	            P+   C  +P+C G +LTLCL +      GAR  L        C  A  +  LV  +  C
	Sbjct: 6365 PNDGSCASLPLCLGASLTLCLGASLALCLGARLALCLGASLPLCLVACLSR*LVASLARC 6186

	Query: 192  LIKGESL 198
	            L+   +L
	Sbjct: 6185 LVASLAL 6165


	>Contig846
	          Length = 1960

	 Score = 31.2 bits (69), Expect = 0.97,   Method: Composition-based stats.
	 Identities = 23/56 (41%), Positives = 28/56 (50%), Gaps = 2/56 (3%)
	 Frame = -1

	Query: 113 LLARLTNALILDLIRGIG--LPSCAGCFLVPMCTGMTLTLCLQSLRKRRPGARYVL 166
	           L+A L   L+  L R +G  LP C G  L P C G +L LCL +   R  GA   L
	Sbjct: 934 LIASLPLCLVASLPRCLGACLPRCLGARL-PRCLGASLALCLSASLPRCLGASLAL 770


	  Database: B_bigemina.contigs.fa
	    Posted date:  Feb 28, 2008 11:56 AM
	  Number of letters in database: 19,733,690
	  Number of sequences in database:  4396
	  
	Lambda     K      H
	   0.320    0.136    0.403 

	Gapped
	Lambda     K      H
	   0.267   0.0410    0.140 


	Matrix: BLOSUM62
	Gap Penalties: Existence: 11, Extension: 1
	Number of Sequences: 4396
	Number of Hits to DB: 10,559,300
	Number of extensions: 167152
	Number of successful extensions: 1271
	Number of sequences better than  1.0: 30
	Number of HSP's gapped: 1178
	Number of HSP's successfully gapped: 31
	Length of query: 448
	Length of database: 6,577,896
	Length adjustment: 100
	Effective length of query: 348
	Effective length of database: 6,138,296
	Effective search space: 2136127008
	Effective search space used: 2136127008
	Neighboring words threshold: 13
	Window for multiple hits: 40
	X1: 16 ( 7.4 bits)
	X2: 38 (14.6 bits)
	X3: 64 (24.7 bits)
	S1: 41 (21.8 bits)
	S2: 43 (21.2 bits)

		

SPS1

	TBLASTN 2.2.17 [Aug-26-2007]


	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.

	Reference for composition-based statistics:
	Schaffer, Alejandro A., L. Aravind, Thomas L. Madden,
	Sergei Shavirin, John L. Spouge, Yuri I. Wolf,  
	Eugene V. Koonin, and Stephen F. Altschul (2001), 
	"Improving the accuracy of PSI-BLAST protein database searches with 
	composition-based statistics and other refinements",  Nucleic Acids Res. 29:2994-3005.

	Query= gi|2411487|gb|AAB88790.1| selenophosphate synthetase
	[Drosophila melanogaster]
	         (398 letters)

	Database: B_bigemina.contigs.fa 
	           4396 sequences; 19,733,690 total letters

	Searching..................................................done

	 ***** No hits found ******


	  Database: B_bigemina.contigs.fa
	    Posted date:  Feb 28, 2008 11:56 AM
	  Number of letters in database: 19,733,690
	  Number of sequences in database:  4396
	  
	Lambda     K      H
	   0.321    0.136    0.409 

	Gapped
	Lambda     K      H
	   0.267   0.0410    0.140 


	Matrix: BLOSUM62
	Gap Penalties: Existence: 11, Extension: 1
	Number of Sequences: 4396
	Number of Hits to DB: 9,204,411
	Number of extensions: 136460
	Number of successful extensions: 669
	Number of sequences better than  1.0: 0
	Number of HSP's gapped: 669
	Number of HSP's successfully gapped: 0
	Length of query: 398
	Length of database: 6,577,896
	Length adjustment: 99
	Effective length of query: 299
	Effective length of database: 6,142,692
	Effective search space: 1836664908
	Effective search space used: 1836664908
	Neighboring words threshold: 13
	Window for multiple hits: 40
	X1: 16 ( 7.4 bits)
	X2: 38 (14.6 bits)
	X3: 64 (24.7 bits)
	S1: 41 (21.9 bits)
	S2: 43 (21.2 bits)

	

SPS2

	TBLASTN 2.2.17 [Aug-26-2007]


	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.

	Reference for composition-based statistics:
	Schaffer, Alejandro A., L. Aravind, Thomas L. Madden,
	Sergei Shavirin, John L. Spouge, Yuri I. Wolf,  
	Eugene V. Koonin, and Stephen F. Altschul (2001), 
	"Improving the accuracy of PSI-BLAST protein database searches with 
	composition-based statistics and other refinements",  Nucleic Acids Res. 29:2994-3005.

	Query= gi|28574119|ref|NP_477478.4| Selenophosphate synthetase
	2[Drosophila melan]
	         (370 letters)

	Database: B_bigemina.contigs.fa 
	           4396 sequences; 19,733,690 total letters

	Searching..................................................done

	 ***** No hits found ******


	  Database: B_bigemina.contigs.fa
	    Posted date:  Feb 28, 2008 11:56 AM
	  Number of letters in database: 19,733,690
	  Number of sequences in database:  4396
	  
	Lambda     K      H
	   0.320    0.136    0.395 

	Gapped
	Lambda     K      H
	   0.267   0.0410    0.140 


	Matrix: BLOSUM62
	Gap Penalties: Existence: 11, Extension: 1
	Number of Sequences: 4396
	Number of Hits to DB: 7,626,367
	Number of extensions: 106203
	Number of successful extensions: 392
	Number of sequences better than  1.0: 0
	Number of HSP's gapped: 392
	Number of HSP's successfully gapped: 0
	Length of query: 370
	Length of database: 6,577,896
	Length adjustment: 98
	Effective length of query: 272
	Effective length of database: 6,147,088
	Effective search space: 1672007936
	Effective search space used: 1672007936
	Neighboring words threshold: 13
	Window for multiple hits: 40
	X1: 16 ( 7.4 bits)
	X2: 38 (14.6 bits)
	X3: 64 (24.7 bits)
	S1: 41 (21.8 bits)
	S2: 43 (21.2 bits)

	

SelH

	TBLASTN 2.2.17 [Aug-26-2007]


	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.

	Reference for composition-based statistics:
	Schaffer, Alejandro A., L. Aravind, Thomas L. Madden,
	Sergei Shavirin, John L. Spouge, Yuri I. Wolf,  
	Eugene V. Koonin, and Stephen F. Altschul (2001), 
	"Improving the accuracy of PSI-BLAST protein database searches with 
	composition-based statistics and other refinements",  Nucleic Acids Res. 29:2994-3005.

	Query= gi|28571177|ref|NP_572903.3| BthD selenoprotein CG11177-PA
	[Drosophila mr]
	         (249 letters)

	Database: B_bigemina.contigs.fa 
	           4396 sequences; 19,733,690 total letters

	Searching..................................................done



	                                                                 Score    E
	Sequences producing significant alignments:                      (bits) Value

	Contig3887                                                             31   0.48 

	>Contig3887
	          Length = 25891

	 Score = 31.2 bits (69), Expect = 0.48,   Method: Composition-based stats.
	 Identities = 23/98 (23%), Positives = 44/98 (44%)
	 Frame = +2

	Query: 149   KKSESTEEAQENEAPTSTSTSRKSKKEQKSEEEPTQVDSKEAKQSKELVKTKRQPKAQKK 208
	             +KS    E   ++A  ++  + +   E+  E       +KEA++SKE  K     + Q  
	Sbjct: 14789 RKSAVIVEVAADDASNASQPAAEGTAEEAGERAANAEQAKEAQESKEHDKADSAKEEQAA 14968

	Query: 209   QAKASESQEEVAEDKPPSSQKRKRTTRSSTDEATAGAK 246
	             + ++ ES +     +  S  +  R  + S D+A AG +
	Sbjct: 14969 EPESKESSDAEQTTESGSEDQGNRENKKSDDKAEAGKE 15082


	  Database: B_bigemina.contigs.fa
	    Posted date:  Mar 4, 2008 12:01 PM
	  Number of letters in database: 19,733,690
	  Number of sequences in database:  4396
	  
	Lambda     K      H
	   0.305    0.121    0.320 

	Gapped
	Lambda     K      H
	   0.267   0.0410    0.140 


	Matrix: BLOSUM62
	Gap Penalties: Existence: 11, Extension: 1
	Number of Sequences: 4396
	Number of Hits to DB: 4,625,062
	Number of extensions: 54524
	Number of successful extensions: 397
	Number of sequences better than  1.0: 30
	Number of HSP's gapped: 389
	Number of HSP's successfully gapped: 47
	Length of query: 249
	Length of database: 6,577,896
	Length adjustment: 94
	Effective length of query: 155
	Effective length of database: 6,164,672
	Effective search space: 955524160
	Effective search space used: 955524160
	Neighboring words threshold: 13
	Window for multiple hits: 40
	X1: 16 ( 7.0 bits)
	X2: 38 (14.6 bits)
	X3: 64 (24.7 bits)
	S1: 43 (21.9 bits)
	S2: 41 (20.4 bits)

	

SelK

	TBLASTN 2.2.17 [Aug-26-2007]


	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.

	Reference for composition-based statistics:
	Schaffer, Alejandro A., L. Aravind, Thomas L. Madden,
	Sergei Shavirin, John L. Spouge, Yuri I. Wolf,  
	Eugene V. Koonin, and Stephen F. Altschul (2001), 
	"Improving the accuracy of PSI-BLAST protein database searches with 
	composition-based statistics and other refinements",  Nucleic Acids Res. 29:2994-3005.

	Query= gi|28571165|ref|NP_572763.3| Selenoprotein G CG1844-PA
	[Drosophila melanogast
	         (110 letters)

	Database: B_bigemina.contigs.fa 
	           4396 sequences; 19,733,690 total letters

	Searching..................................................done



	                                                                 Score    E
	Sequences producing significant alignments:                      (bits) Value

	Contig4132                                                             28   1.00 

	>Contig4132
	          Length = 50512

	 Score = 27.7 bits (60), Expect = 1.00,   Method: Composition-based stats.
	 Identities = 15/30 (50%), Positives = 16/30 (53%)
	 Frame = +1

	Query: 59    GGGGGWGGGGGGGGGGGGGRPGSGSGGLRP 88
	             GGG G   G G  GG   G PGSGS   +P
	Sbjct: 13156 GGGKGPASGAGQQGGTSKGSPGSGSKATQP 13245


	  Database: B_bigemina.contigs.fa
	    Posted date:  Feb 28, 2008 11:56 AM
	  Number of letters in database: 19,733,690
	  Number of sequences in database:  4396
	  
	Lambda     K      H
	   0.322    0.149    0.534 

	Gapped
	Lambda     K      H
	   0.267   0.0410    0.140 


	Matrix: BLOSUM62
	Gap Penalties: Existence: 11, Extension: 1
	Number of Sequences: 4396
	Number of Hits to DB: 3,676,195
	Number of extensions: 75439
	Number of successful extensions: 8125
	Number of sequences better than  1.0: 209
	Number of HSP's gapped: 4572
	Number of HSP's successfully gapped: 967
	Length of query: 110
	Length of database: 6,577,896
	Length adjustment: 82
	Effective length of query: 28
	Effective length of database: 6,217,424
	Effective search space: 174087872
	Effective search space used: 174087872
	Neighboring words threshold: 13
	Window for multiple hits: 40
	X1: 16 ( 7.4 bits)
	X2: 38 (14.6 bits)
	X3: 64 (24.7 bits)
	S1: 41 (21.8 bits)
	S2: 38 (19.2 bits)

	

SBP2

	TBLASTN 2.2.17 [Aug-26-2007]


	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.

	Reference for composition-based statistics:
	Schaffer, Alejandro A., L. Aravind, Thomas L. Madden,
	Sergei Shavirin, John L. Spouge, Yuri I. Wolf,  
	Eugene V. Koonin, and Stephen F. Altschul (2001), 
	"Improving the accuracy of PSI-BLAST protein database searches with 
	composition-based statistics and other refinements",  Nucleic Acids Res. 29:2994-3005.

	Query= gi|21355989|ref|NP_648204.1| CG7066 CG7066-PA [Drosophila
	melanogaster]
	         (313 letters)

	Database: B_bigemina.contigs.fa 
	           4396 sequences; 19,733,690 total letters

	Searching..................................................done

	 ***** No hits found ******


	  Database: B_bigemina.contigs.fa
	    Posted date:  Feb 28, 2008 11:56 AM
	  Number of letters in database: 19,733,690
	  Number of sequences in database:  4396
	  
	Lambda     K      H
	   0.319    0.132    0.377 

	Gapped
	Lambda     K      H
	   0.267   0.0410    0.140 


	Matrix: BLOSUM62
	Gap Penalties: Existence: 11, Extension: 1
	Number of Sequences: 4396
	Number of Hits to DB: 6,676,898
	Number of extensions: 104494
	Number of successful extensions: 561
	Number of sequences better than  1.0: 2
	Number of HSP's gapped: 560
	Number of HSP's successfully gapped: 2
	Length of query: 313
	Length of database: 6,577,896
	Length adjustment: 96
	Effective length of query: 217
	Effective length of database: 6,155,880
	Effective search space: 1335825960
	Effective search space used: 1335825960
	Neighboring words threshold: 13
	Window for multiple hits: 40
	X1: 16 ( 7.4 bits)
	X2: 38 (14.6 bits)
	X3: 64 (24.7 bits)
	S1: 41 (21.8 bits)
	S2: 42 (20.8 bits)

	

Secp43

	TBLASTN 2.2.17 [Aug-26-2007]


	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.

	Reference for composition-based statistics:
	Schaffer, Alejandro A., L. Aravind, Thomas L. Madden,
	Sergei Shavirin, John L. Spouge, Yuri I. Wolf,  
	Eugene V. Koonin, and Stephen F. Altschul (2001), 
	"Improving the accuracy of PSI-BLAST protein database searches with 
	composition-based statistics and other refinements",  Nucleic Acids Res. 29:2994-3005.

	Query= gi|28574707|ref|NP_608837.2| CG15440 CG15440-PA [Drosophila
	melanogaster]
	         (336 letters)

	Database: B_bigemina.contigs.fa 
	           4396 sequences; 19,733,690 total letters

	Searching..................................................done



	                                                                 Score    E
	Sequences producing significant alignments:                      (bits) Value

	Contig4184.3                                                           67   2e-11
	Contig3945                                                             48   5e-06
	Contig3934                                                             48   5e-06
	Contig4104.1                                                           38   0.005
	Contig3713                                                             38   0.007
	Contig4166.0                                                           34   0.11 
	Contig4158                                                             32   0.57 

	>Contig4184.3
	          Length = 137674

	 Score = 66.6 bits (161), Expect = 2e-11,   Method: Composition-based stats.
	 Identities = 49/158 (31%), Positives = 78/158 (49%)
	 Frame = -1

	Query: 8     LWMGSLESYMTENFIIAAFRKMGEDPTTVRLMRNKYTGEPAGYCFVNFISDDHALDAMHK 67
	             L++G L   +TE  +   F  +G    ++R+ R+  T +  GY +VN+ +   A  A+  
	Sbjct: 30733 LYVGDLLPEVTEAMLYEVFNGVGP-VASIRVCRDSVTRKSLGYAYVNYYNVQDAETALEC 30557

	Query: 68    LNGKPIPGTNPIVRFRLNSASNSYKLPGNEREFSVWVGDLSSDVDDYQLYKVFSSKFTSI 127
	             LN   I G    + +     S      GN     ++V +L   +D   LY  FS  F +I
	Sbjct: 30556 LNYIDIKGHPARIMWSDRDPSLRKSGTGN-----IFVKNLDKAIDTKALYDTFS-HFGTI 30395

	Query: 128   KTAKVILDSLGFSKGYGFVRFGIEDEQKSALYDMNGYI 165
	              + KV +D+LG SKGYGFV +  E+  K A+  +NG +
	Sbjct: 30394 LSCKVAMDALGNSKGYGFVHYTTEESAKEAIEKVNGML 30281



	 Score = 37.0 bits (84), Expect = 0.012,   Method: Composition-based stats.
	 Identities = 31/118 (26%), Positives = 54/118 (45%), Gaps = 1/118 (0%)
	 Frame = +2

	Query: 49     GYCFVNFISDDHALDAMHKLNGKPIPGTNPIVRFRLNSASNSYKLPGNEREFSVWVGDLS 108
	              GYCF+ F + + A D   KL+G    G    +R  +   + S     +  +  V+V ++ 
	Sbjct: 104231 GYCFLEFSTAELA-DMCFKLDGINYKGQTLKIRRPIEYGTTS-----SSDDTKVFVQNIP 104392

	Query: 109    SDVDDYQLYKVFSSKFTSIKTAKVILD-SLGFSKGYGFVRFGIEDEQKSALYDMNGYI 165
	                + D  + K    K   + +A ++ D   G +KGYGF  F      K A+  +NG++
	Sbjct: 104393 PTMSDADI-KALLEKHGKLTSANLVKDLKTGQNKGYGFFEFDDSRAAKMAVCHLNGHV 104563


	>Contig3945
	          Length = 13847

	 Score = 48.1 bits (113), Expect = 5e-06,   Method: Composition-based stats.
	 Identities = 41/128 (32%), Positives = 66/128 (51%), Gaps = 6/128 (4%)
	 Frame = +3

	Query: 53    VNFISDDHALDAMHKLNGKPIPGTNPIVRFRLNSASNS-YKLPG----NEREFSVWVGDL 107
	             +N I+D+ A     ++  +   GT    R  L  A+ S +  P      E +F ++ GDL
	Sbjct: 13401 INRITDEEAA----RMGSRMAVGTKGNGRIHLRKAAGSVWNDPTLDDWPENDFRIFCGDL 13568

	Query: 108   SSDVDDYQLYKVFSSKFTSIKTAKVILD-SLGFSKGYGFVRFGIEDEQKSALYDMNGYIG 166
	              +DV D  L   F  K+ + + A+VI D + G SKGYGFV     ++  +AL +MN    
	Sbjct: 13569 GNDVTDDVLANAFK-KYPTFQRARVIRDRNSGKSKGYGFVSMMSAEDMLAALNEMNNKF- 13742

	Query: 167   LGTKPIKI 174
	             +G +PI++
	Sbjct: 13743 VGNRPIRV 13766


	>Contig3934
	          Length = 22845

	 Score = 48.1 bits (113), Expect = 5e-06,   Method: Composition-based stats.
	 Identities = 30/79 (37%), Positives = 47/79 (59%), Gaps = 1/79 (1%)
	 Frame = -3

	Query: 97    EREFSVWVGDLSSDVDDYQLYKVFSSKFTSIKTAKVILD-SLGFSKGYGFVRFGIEDEQK 155
	             E +F ++ GDL +DV D  L   F  K+ + + A+VI D + G SKGYGFV     ++  
	Sbjct: 21145 ENDFRIFCGDLGNDVTDDVLANAFK-KYPTFQRARVIRDRNSGKSKGYGFVSMMSAEDML 20969

	Query: 156   SALYDMNGYIGLGTKPIKI 174
	             +AL +MN    +G +PI++
	Sbjct: 20968 AALNEMNNKF-VGNRPIRV 20915


	>Contig4104.1
	          Length = 37990

	 Score = 38.1 bits (87), Expect = 0.005,   Method: Composition-based stats.
	 Identities = 29/101 (28%), Positives = 55/101 (54%), Gaps = 3/101 (2%)
	 Frame = -3

	Query: 98   REFSVWVGDLSSDVDDYQLYKVFSSKFTSIKTAKVILDSL-GFSKGYGFVRFGIEDEQKS 156
	            +E +++VG++ + VD+  L++ F  +   +K   +  D + G  +GY FV F  ED+   
	Sbjct: 8885 QEATLYVGNVDTQVDEELLWEFFV-QVGPVKHVHIPRDKVTGHHQGYAFVEFESEDDADY 8709

	Query: 157  ALYDMNGYIGLGTKPIKICNAVPKPKS--ELGGAVGEGNTN 195
	            A+  +N ++ L  KP++ CN   + K   E+G  +  GN +
	Sbjct: 8708 AIRILN-FVKLYNKPLR-CNKASRDKQNYEIGANLFIGNLD 8592


	>Contig3713
	          Length = 2421

	 Score = 37.7 bits (86), Expect = 0.007,   Method: Composition-based stats.
	 Identities = 26/86 (30%), Positives = 43/86 (50%), Gaps = 11/86 (12%)
	 Frame = -2

	Query: 8   LWMGSLESYMTENFIIAAFRKMGEDPTT-VRLMRNKYTGEPAGYCFVNFISDDHALDAMH 66
	           L++G++++ + E  +   F ++G  P   V + R+K TG   GY FV F S+DHA   + 
	Sbjct: 593 LYVGNVDTQVDEELLWEFFVQVG--PVKHVHIPRDKVTGHHQGYAFVEFESEDHADYPLR 420

	Query: 67  KLN----------GKPIPGTNPIVRF 82
	            LN          G  +P T+   +F
	Sbjct: 419 FLNFLKRVKKPFRGNKVPRTSQTSKF 342


	>Contig4166.0
	          Length = 90045

	 Score = 33.9 bits (76), Expect = 0.11,   Method: Composition-based stats.
	 Identities = 21/88 (23%), Positives = 44/88 (50%), Gaps = 1/88 (1%)
	 Frame = -1

	Query: 84    LNSASNSYKLPGNEREFSVWVGDLSSDVDDYQLYKVFSSKFTSIKTAKVILD-SLGFSKG 142
	              ++ S + + P      ++++  + ++   Y L + FS +F  I ++++  D S G  KG
	Sbjct: 30189 FDNPSGASQAPAGPPGANLFIFHIPNEWTQYDLVRTFS-QFGKILSSRIASDKSTGRHKG 30013

	Query: 143   YGFVRFGIEDEQKSALYDMNGYIGLGTK 170
	             Y FV +   D    A+  +NG+  +G +
	Sbjct: 30012 YAFVSYDSPDSAAQAIQHLNGFTVMGKR 29929



	 Score = 32.7 bits (73), Expect = 0.25,   Method: Composition-based stats.
	 Identities = 19/58 (32%), Positives = 32/58 (55%)
	 Frame = -1

	Query: 18    TENFIIAAFRKMGEDPTTVRLMRNKYTGEPAGYCFVNFISDDHALDAMHKLNGKPIPG 75
	             T+  ++  F + G+  ++ R+  +K TG   GY FV++ S D A  A+  LNG  + G
	Sbjct: 30105 TQYDLVRTFSQFGKILSS-RIASDKSTGRHKGYAFVSYDSPDSAAQAIQHLNGFTVMG 29935


	>Contig4158
	          Length = 87318

	 Score = 31.6 bits (70), Expect = 0.57,   Method: Composition-based stats.
	 Identities = 19/64 (29%), Positives = 36/64 (56%)
	 Frame = +3

	Query: 101   SVWVGDLSSDVDDYQLYKVFSSKFTSIKTAKVILDSLGFSKGYGFVRFGIEDEQKSALYD 160
	             +++V  L ++ DD      F  +F ++ +A+V  DS G +KG+GF+ +       +A+  
	Sbjct: 67644 NLFVFHLPAEWDDADFVLNFQ-QFGTVLSARVQRDSAGRNKGFGFISYDNPLSALNAIKS 67820

	Query: 161   MNGY 164
	             MNG+
	Sbjct: 67821 MNGF 67832


	  Database: B_bigemina.contigs.fa
	    Posted date:  Feb 28, 2008 11:56 AM
	  Number of letters in database: 19,733,690
	  Number of sequences in database:  4396
	  
	Lambda     K      H
	   0.314    0.131    0.390 

	Gapped
	Lambda     K      H
	   0.267   0.0410    0.140 


	Matrix: BLOSUM62
	Gap Penalties: Existence: 11, Extension: 1
	Number of Sequences: 4396
	Number of Hits to DB: 7,395,256
	Number of extensions: 101067
	Number of successful extensions: 451
	Number of sequences better than  1.0: 9
	Number of HSP's gapped: 447
	Number of HSP's successfully gapped: 11
	Length of query: 336
	Length of database: 6,577,896
	Length adjustment: 97
	Effective length of query: 239
	Effective length of database: 6,151,484
	Effective search space: 1470204676
	Effective search space used: 1470204676
	Neighboring words threshold: 13
	Window for multiple hits: 40
	X1: 16 ( 7.3 bits)
	X2: 38 (14.6 bits)
	X3: 64 (24.7 bits)
	S1: 42 (22.0 bits)
	S2: 42 (20.8 bits)

	

PSTK

	TBLASTN 2.2.17 [Aug-26-2007]


	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.

	Reference for composition-based statistics:
	Schaffer, Alejandro A., L. Aravind, Thomas L. Madden,
	Sergei Shavirin, John L. Spouge, Yuri I. Wolf,  
	Eugene V. Koonin, and Stephen F. Altschul (2001), 
	"Improving the accuracy of PSI-BLAST protein database searches with 
	composition-based statistics and other refinements",  Nucleic Acids Res. 29:2994-3005.

	Query= gi|19920352|ref|NP_608329.1| CG12788 CG12788-PA [Drosophila
	melanogaster]
	         (292 letters)

	Database: B_bigemina.contigs.fa 
	           4396 sequences; 19,733,690 total letters

	Searching..................................................done



	                                                                 Score    E
	Sequences producing significant alignments:                      (bits) Value

	Contig4120                                                             31   0.70 

	>Contig4120
	          Length = 42572

	 Score = 30.8 bits (68), Expect = 0.70,   Method: Composition-based stats.
	 Identities = 23/70 (32%), Positives = 35/70 (50%)
	 Frame = +2

	Query: 130   MASSLDSCLQANSLRSDATRVPVDVVRQMNERLEVPDTSEAWERNSLTLNGLDMDTTGSA 189
	             + S+LDS    +++     +   D+V+    RLE P    A++R S  L+G+D D T  A
	Sbjct: 31088 LQSALDSSNLPSTIVLLEGKSDGDIVKAAMRRLEYP-LQVAFKRLSFALSGVDADDTSRA 31264

	Query: 190   LLAFIASLLD 199
	             L      LLD
	Sbjct: 31265 LTQRSRDLLD 31294


	  Database: B_bigemina.contigs.fa
	    Posted date:  Mar 4, 2008 11:40 AM
	  Number of letters in database: 19,733,690
	  Number of sequences in database:  4396
	  
	Lambda     K      H
	   0.322    0.138    0.407 

	Gapped
	Lambda     K      H
	   0.267   0.0410    0.140 


	Matrix: BLOSUM62
	Gap Penalties: Existence: 11, Extension: 1
	Number of Sequences: 4396
	Number of Hits to DB: 6,722,194
	Number of extensions: 111013
	Number of successful extensions: 613
	Number of sequences better than  1.0: 1
	Number of HSP's gapped: 613
	Number of HSP's successfully gapped: 1
	Length of query: 292
	Length of database: 6,577,896
	Length adjustment: 96
	Effective length of query: 196
	Effective length of database: 6,155,880
	Effective search space: 1206552480
	Effective search space used: 1206552480
	Neighboring words threshold: 13
	Window for multiple hits: 40
	X1: 16 ( 7.4 bits)
	X2: 38 (14.6 bits)
	X3: 64 (24.7 bits)
	S1: 41 (21.9 bits)
	S2: 42 (20.8 bits)

	

tRNASec

	BLASTN 2.2.17 [Aug-26-2007]


	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= Predicted RNA/mRNA Sequence (introns spliced out)  Hide
	         (90 letters)

	Database: B_bigemina.contigs.fa 
	           4396 sequences; 19,733,690 total letters

	Searching..................................................done

	 ***** No hits found ******


	  Database: B_bigemina.contigs.fa
	    Posted date:  Feb 28, 2008 11:56 AM
	  Number of letters in database: 19,733,690
	  Number of sequences in database:  4396
	  
	Lambda     K      H
	    1.37    0.711     1.31 

	Gapped
	Lambda     K      H
	    1.37    0.711     1.31 


	Matrix: blastn matrix:1 -3
	Gap Penalties: Existence: 5, Extension: 2
	Number of Sequences: 4396
	Number of Hits to DB: 13,786
	Number of extensions: 591
	Number of successful extensions: 171
	Number of sequences better than  1.0: 0
	Number of HSP's gapped: 171
	Number of HSP's successfully gapped: 0
	Length of query: 90
	Length of database: 19,733,690
	Length adjustment: 15
	Effective length of query: 75
	Effective length of database: 19,667,750
	Effective search space: 1475081250
	Effective search space used: 1475081250
	X1: 11 (21.8 bits)
	X2: 15 (29.7 bits)
	X3: 25 (49.6 bits)
	S1: 12 (24.3 bits)
	S2: 16 (32.2 bits)