BLASTP 2.2.26+ Reference: Stephen F. Altschul, 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 compositional score matrix adjustment: Stephen F. Altschul, John C. Wootton, E. Michael Gertz, Richa Agarwala, Aleksandr Morgulis, Alejandro A. Schaffer, and Yi-Kuo Yu (2005) "Protein database searches using compositionally adjusted substitution matrices", FEBS J. 272:5101-5109. RID: GNZ3KF5H016 Database: All non-redundant GenBank CDS translations+PDB+SwissProt+PIR+PRF excluding environmental samples from WGS projects 16,839,230 sequences; 5,784,003,470 total letters Query= contig_36:subseq(2370043,20236) [translate(1)] L.Major Length=113 Score E Sequences producing significant alignments: (Bits) Value gb|ABB82550.1| selenoprotein 1 [Leishmania major] 159 3e-49 gb|EGC99887.1| putative sigma-54 modulation protein [Burkhold... 33.5 4.5 ref|YP_001578715.1| sigma 54 modulation protein/ribosomal pro... 33.1 5.6 ALIGNMENTS >gb|ABB82550.1| selenoprotein 1 [Leishmania major] Length=117 Score = 159 bits (401), Expect = 3e-49, Method: Compositional matrix adjust. Identities = 102/116 (88%), Positives = 104/116 (90%), Gaps = 3/116 (3%) Query 1 MRVLRHTGTCSDIAVDATESSAATHHKQADSQ---QRNEYPFSTILLIMDLXENWRTALY 57 MRVLRHTGTCSDIAVDATESSAATHHKQADSQ + NEYPFSTILLIMDL ENWR ALY Sbjct 1 MRVLRHTGTCSDIAVDATESSAATHHKQADSQATTKENEYPFSTILLIMDLUENWRIALY 60 Query 58 FVLAVRRSSHSPSSLAPAVLHQLPPLVLLLAALCFYSSTAAAAHTAIESHERWSWT 113 FVLAV RSSHSPSS AP VL QLPPLV LLAALCFYSSTAAAAHTAIE+ ERWSWT Sbjct 61 FVLAVSRSSHSPSSPAPTVLRQLPPLVPLLAALCFYSSTAAAAHTAIETRERWSWT 116 >gb|EGC99887.1| putative sigma-54 modulation protein [Burkholderia sp. TJI49] Length=119 Score = 33.5 bits (75), Expect = 4.5, Method: Compositional matrix adjust. Identities = 29/99 (29%), Positives = 44/99 (44%), Gaps = 24/99 (24%) Query 2 RVLRHTGTCSDIAVDATESSAATHHKQADSQQRNEYP-----------------FSTILL 44 RVLRH SD +D T + +HK+ D QQR E ++ I L Sbjct 26 RVLRH----SDQVIDGTVILSVDNHKEKDKQQRAEINLHLKGKDIFVESSNGNLYAAIDL 81 Query 45 IMDLXENWRTALYFVLAVRRSSHSPSSLAPAVLH-QLPP 82 ++D + R + + ++ +H P L PAV +LPP Sbjct 82 LIDKLD--RQVVKHMERLQTHAHDPIKLQPAVEQIELPP 118 >ref|YP_001578715.1| sigma 54 modulation protein/ribosomal protein S30EA [Burkholderia multivorans ATCC 17616] ref|YP_001947161.1| putative sigma-54 modulation protein [Burkholderia multivorans ATCC 17616] ref|ZP_03573039.1| ribosomal subunit interface protein [Burkholderia multivorans CGD2M] ref|ZP_03579861.1| ribosomal subunit interface protein [Burkholderia multivorans CGD2] ref|ZP_03584182.1| ribosomal subunit interface protein [Burkholderia multivorans CGD1] gb|ABX14218.1| sigma 54 modulation protein/ribosomal protein S30EA [Burkholderia multivorans ATCC 17616] dbj|BAG44625.1| putative sigma54 modulation protein [Burkholderia multivorans ATCC 17616] gb|EEE01032.1| ribosomal subunit interface protein [Burkholderia multivorans CGD1] gb|EEE05939.1| ribosomal subunit interface protein [Burkholderia multivorans CGD2] gb|EEE12639.1| ribosomal subunit interface protein [Burkholderia multivorans CGD2M] Length=119 Score = 33.1 bits (74), Expect = 5.6, Method: Compositional matrix adjust. Identities = 29/99 (29%), Positives = 44/99 (44%), Gaps = 24/99 (24%) Query 2 RVLRHTGTCSDIAVDATESSAATHHKQADSQQR-----------------NEYPFSTILL 44 RVLRH SD +D T + +HK+ D QQR N ++ I L Sbjct 26 RVLRH----SDQVIDGTVILSVDNHKEKDKQQRAEINLHLKGKDIFVESSNGNLYAAIDL 81 Query 45 IMDLXENWRTALYFVLAVRRSSHSPSSLAPAVLH-QLPP 82 ++D + R + + ++ +H P L PAV +LPP Sbjct 82 LIDKLD--RQVVKHMERLQTHAHDPIKLQPAVDQIELPP 118 Database: All non-redundant GenBank CDS translations+PDB+SwissProt+PIR+PRF excluding environmental samples from WGS projects Posted date: Jan 6, 2012 10:12 AM Number of letters in database: 1,489,036,174 Number of sequences in database: 16,839,230 Lambda K H 0.326 0.129 0.404 Gapped Lambda K H 0.267 0.0410 0.140 Matrix: BLOSUM62 Gap Penalties: Existence: 11, Extension: 1 Number of Sequences: 16839230 Number of Hits to DB: 114572459 Number of extensions: 3497482 Number of successful extensions: 10888 Number of sequences better than 100: 2 Number of HSP's better than 100 without gapping: 0 Number of HSP's gapped: 10888 Number of HSP's successfully gapped: 2 Length of query: 113 Length of database: 5784003470 Length adjustment: 81 Effective length of query: 32 Effective length of database: 4420025840 Effective search space: 141440826880 Effective search space used: 141440826880 T: 11 A: 40 X1: 15 (7.0 bits) X2: 38 (14.6 bits) X3: 64 (24.7 bits) S1: 40 (20.0 bits) S2: 67 (30.4 bits)