<4D F736F F F696E74202D F90B695A8947A97F189F090CD8AEE91625F3189F196DA5F8E9197BF2E >

Size: px

Start display at page:

Download "<4D F736F F F696E74202D F90B695A8947A97F189F090CD8AEE91625F3189F196DA5F8E9197BF2E >"

ようたやなぎしま
5 years ago
Views:

1 2 物配列 = 塩基配列およびアミノ酸配列塩基配列 = DNA の塩基 (G A T C) の並び順どのようにして, 塩基配列 (GATCの並び順) を読むのか? 塩基配列の決定法 = DNA シークエンシング 30 3 ジデオキシ法別名 : サンガー法 Sanger et al., 1977 DNAポリメラーゼを使って相補鎖を合成する反応をう特定のヌクレオチドの位置で反応が停するようにしておく 4

2 5 デオキシヌクレオチド (dntp) ジデオキシヌクレオチド (ddntp) DNA ポリメラーゼの伸反応が進む伸反応が進まないヌクレオチド A: アデニン C: シトシン G: グアニン T: チミン蛍光標識試薬やサイクルシークエンスキャピラリー電気泳動などの技術によって精度や迅速性が飛躍的に向上したー陰極 7 4 キャピラリー電気泳動により分離する鋳型 DNA プライマー DNA ポリメラーゼ + 陽極 11 つのチューブで反応を行う 2ddNTP が取り込まれたところで反応が止まる 5 レーザー光を当てて蛍光色素を励起する 6 検出器で蛍光を測定する 3 未反応の蛍光試薬を取り除く 7 データがコンピューターに送られるこの DNA の塩基配列は A A C G T A T G 7

3 9 DNA シーケンサーからの出力ファイルの例ベースコール DNAシークエンサーから出力されたクロマトグラム ( 波形データ ) をもとに, 塩基配列を決定する作業. ATCTTAC GATC の文字列にしてしまうとシークエンスの読め具合がわからなくなってしまうので, 注意が必要

National Center for Biotechnology Information

なデータベースツールがまとめられている AP009356 と 13 14 All

これらのデータベースに対して横断検索ができる 15 16 例えば replication

4 National Center for Biotechnology Information 通称 :NCBI 国の国衛研究所 (NIH) の下の国医学図書館が運営するWebサイト GenbankやPubMed BLASTなど有なデータベースツールがまとめられている AP と All Databases データベースの統合検索システム主なデータベースは,PubMed 塩基配列データベースアミノ酸配列データベースゲノムデータベース 3D 分構造データベースなど All Databasesからはこれらのデータベースに対して横断検索ができる例えば replication protein phytoplasma としてみるファイトプラズマ Candidatus Phytoplasma 属細菌健全アジサイ罹病アジサイ植物の篩部細胞に寄する植物病原細菌感染植物ではがくや花弁が葉化する

5 FASTA フォーマット 17 データベースカタログ GenBank フォーマット 18 相同性検索 ( ホモロジー検索 ) 相同性検索 ( ホモロジー検索 ) とは? 相同性検索は, 配列の類似性から類縁の遺伝タンパク質を検索する法で, 進化系統分類の解析, 機能解析などを的とした配列解析の最も基本的な法のつである. SSEARCH 質問配列 (Query) MIGMMIT 質問配列と類似した ( 相同な ) 配列を, データベース上から探索する FASTA BLAST MMIGPIT MYLGPST MDIGIT MIIGMIT MMIQPMMDG 20

MIGMMIT MMIGPIT アラインメント ( 並置 ) アラインメント 2 つの配列を要素ごとに対応づけて並べる操作進化の過程で生じ得る配列要素の挿入欠失をギャップ (-) で対応づけるグローバルアラインメント配列全体の類似性を考慮 a = M-IGMMIT b = MMIGP-IT 二つのアミノ酸配列を整列化させるにはどのように並べればよいか?

C) = 8 3-7 3 9 s(a, A) + s(f, E) + s(d, E) + gap + s(c, C) = 0 3-7 3-8 9 完全に一致するアミノ酸や, 類似アミノ酸には高い点数を与えたい各アミノ酸の点数はどのように求めればよいか?

6 MIGMMIT MMIGPIT アラインメント ( 並置 ) アラインメント 2 つの配列を要素ごとに対応づけて並べる操作進化の過程で生じ得る配列要素の挿入欠失をギャップ (-) で対応づけるグローバルアラインメント配列全体の類似性を考慮 a = M-IGMMIT b = MMIGP-IT 二つのアミノ酸配列を整列化させるにはどのように並べればよいか? ローカルアラインメント局所的な類似性を考慮 a = MIGMMIT--- b = ---MMIGPIT 21 AFDC AEEC AFDGC AEE-C アラインメントスコアの計算配列の類似度 = アラインメントのスコアアラインメントのスコアの計算対応する各要素の類似度スコアの和ギャップの挿入にはペナルティを与える s(a, A) + s(f, E) + s(d, E) + s(c, C) = s(a, A) + s(f, E) + s(d, E) + gap + s(c, C) = 完全に一致するアミノ酸や, 類似アミノ酸には高い点数を与えたい各アミノ酸の点数はどのように求めればよいか? 22 BLOSUM スコア ( Henikoff らの方法 ) BLOSUM: BLOcks amino acid Substitution Matrix 同一ファミリータンパク質のギャップなしでアラインメントされた領域 ( ブロック ) に対しアミノ BLOSUM50マトリックス酸の置換の頻度を調べて作成良く似た配列の寄与が優勢になりすぎないように, 例えば62% 一致のパターンを一まとめにしてBLOSUM62を作るのに用いる. アラインメントのアルゴリズム Needleman-Wunschのアルゴリズム 2つの配列の最適なグローバルアラインメントを, ダイナミックプログラミング ( 動的計画法 ) により求める Smith-Waterman のアルゴリズム 2つの配列の部分配列間の致を探索する最もいスコアをもつ致箇所をすアラインメントを求めるダイナミックプログラミング ( 動的計画法 ) 23 24

FASTA と BLAST 動的計画法による検索法 (SSERACH) は mnに例した時間を要する (m, nは配列のさ ) 配列データベースに登録されている配列の数は膨時間がかかりすぎてしまう BLAST 検索他の法にべて速であり, ホモロジー検索の法として最もよく利されている質問配列を固定の断 ( ワード ) に区切るまずはワード単位で類似する断

AGPVFG GPVFGI 致する部分を検索 MSGPVFGIP ワードの切り出し Default の設定ではアミノ酸の場合は 6 字, 塩基配列は 28 字.

7 FASTA と BLAST 動的計画法による検索法 (SSERACH) は mnに例した時間を要する (m, nは配列のさ ) 配列データベースに登録されている配列の数は膨時間がかかりすぎてしまう BLAST 検索他の法にべて速であり, ホモロジー検索の法として最もよく利されている質問配列を固定の断 ( ワード ) に区切るまずはワード単位で類似する断をデータベース上から検索類似度が最になるまで両向にアラインメントを伸ばす最後にこれらの局所的なアラインメントを結合する 26 FASTA 最初に致する配列断を速に検索して絞り込む Lipman and Pearson (1985) BLAST 最初に局所的に類似の部分配列を速に検索して絞り込む Altschul (1990) 25 MAGPVFGIPSCSF MAGPVF AGPVFG GPVFGI 致する部分を検索 MSGPVFGIP ワードの切り出し Default の設定ではアミノ酸の場合は 6 字, 塩基配列は 28 字. 致したワードを中にして両向にアラインメントを伸ばしていく ( 類似度が下がってきたらアラインメントを終了する ) NCBI のトップページの右にあるリンクから BLAST 検索のページへ 27 BLASTP 検索 (protein blast) 2 貼り付ける >sample1 MNRVFLFGKLSFTPNRLQTKNGTLGATFSMECLDS SGFNNAKSFIRVTAWGKVASFIVAQNPGVMLFVEG RLTTYKITNSENKNTYALQVTADKIFHPDEKTTNE EPIKSTVVDSPFMNPKASVTEAEFEQAFPHQDETD FNNITPIFENDVQLEEESDD プログラム質問配列 (query) 検索対象 Protein BLAST アミノ酸配列アミノ酸配列データベース blastx 塩基配列アミノ酸配列データベース Nucleotide BLAST 塩基配列塩基配列データベース 1 配列をコピーする ( > の行は入れても入れなくてもよい ) 3 データベースを選ぶ (nr) 4 BLAST を押す tblastn アミノ酸配列塩基配列データベース tblastx 塩基配列塩基配列データベース nr : 冗長性をなくした (non-redundant) アミノ酸データベース 28

ヒットした配列の GenBank フォーマットデータへのリンク 29 E-value E-value は,

アラインメント Query : 質問配列中段 : 一致するアミノ酸あるいは +( 類似アミノ酸 ) Sbjct : Blast

8 ヒットした配列の GenBank フォーマットデータへのリンク 29 E-value E-value は, 現在のデータベースにおいて, 全く偶然に同じスコアに E-value 全長は 166 アミノ酸相同性 (identity) 相同性 (similarity) ギャップなる配列の数の期待値であり,E valueが小さいほど偶然には起こり得ないことを示している. スコア BLAST 検索の際にE-valueのしきい値を設定することで, その値よりも小さいE-valueの検索結果しか出力されなくなる. アラインメント Query : 質問配列中段 : 一致するアミノ酸あるいは +( 類似アミノ酸 ) Sbjct : Blast 検索の結果, ヒットした配列全長ではないので注意 ( 本当は 166 番目に E がある ) 30 検索結果の表示件数 blastx 塩基配列を入力 E-valueのしきい値 BLAST 検索時のWordサイズマトリックスの種類を選ぶギャップのスコア設定 6 通りの reading frame のすべてについて翻訳し, アミノ酸配列データベースに対して検索してくれる E-value 計算時の設定冗長配列を取り除く場合はチェック冗長配列を取り除く場合の設定小文字を無視する場合の設定塩基配列を決定したが, 何がコードされているかわからないとき non-coding 領域に, タンパク質がコードされていないかどうか, 調べたいときなど 32

>sample2 ATGAAATTAAGAATCTGCGAACTTGTTATTAATAAAACTTTAATTACTAAAACTAAA ATAGAAACTATTTTAGAAACTAAAAAAAAAGCCATTCAAAATTATGCCTATATTTTG CATGATAAAGATATTTATCAAAATGATAAAGAGGCTCAATTGAATGGTAAAAAAGTA

>sample3 TTGAAGAGGACTTGGAACTTCGAT 1 配列をコピーする ( > の行は入れても入れなくてもよい ) 2 貼り付ける 3 データベースを選ぶ (nr/nt) 4 BLAST を押す 33 と表示され, 短い配列用の設定で検索される 34 tblastn アミノ酸配列を入力データベース上の塩基配列を,6 通りの reading frame のすべてについて翻訳し,

LASENIHFYAINSDLQHLAFKTNVKNKLLIQDHTNKGFGAGG DPAKGASLAISFQEQFNTLTDGYDFCILVAGFGKGTGTGATP VFSKILKTKKILNVAIVTYPSLNEGLTVRNKATKGLEILNKA TDSYMLFCNEKCTNGIYQLANTEIVSAIKNLIELITIPLQQN

9 >sample2 ATGAAATTAAGAATCTGCGAACTTGTTATTAATAAAACTTTAATTACTAAAACTAAA ATAGAAACTATTTTAGAAACTAAAAAAAAAGCCATTCAAAATTATGCCTATATTTTG CATGATAAAGATATTTATCAAAATGATAAAGAGGCTCAATTGAATGGTAAAAAAGTA GGAGATATAAAAGCTCCTCATTGGCATATATATTTAAGATTTAATTATTCACATGAT ACAAAAAATATCGCTCAATGGTTTAATACTGAGGATAATTTTGTTTCCAAAATAAAA GGTAGATTTAGTGATGCCTTAATGTATATGATTCATGCTAATAGGTC blastx 検索 blastn (nucleotide blast) >sample3 TTGAAGAGGACTTGGAACTTCGAT 1 配列をコピーする ( > の行は入れても入れなくてもよい ) 2 貼り付ける 3 データベースを選ぶ (nr/nt) 4 BLAST を押す 33 と表示され, 短い配列用の設定で検索される 34 tblastn アミノ酸配列を入力データベース上の塩基配列を,6 通りの reading frame のすべてについて翻訳し, このアミノ酸配列データに対して検索してくれる EST 配列やドラフトゲノムなど, アノテーション情報が整備されていないデータから相同な配列を探したいときに便利 tblastx 塩基配列を入力 BLAST 検索 (GenomeNet) 35 >sample5 2 貼り付ける MDENETQFNKLNQVKNKLKIGVFGIGGAGNNIVDASLYHYPN LASENIHFYAINSDLQHLAFKTNVKNKLLIQDHTNKGFGAGG DPAKGASLAISFQEQFNTLTDGYDFCILVAGFGKGTGTGATP VFSKILKTKKILNVAIVTYPSLNEGLTVRNKATKGLEILNKA TDSYMLFCNEKCTNGIYQLANTEIVSAIKNLIELITIPLQQN IDFEDVRAFFQTKKTNQDQQLFTVTHPFSFSFDSKDSIEQFA KQFKNFEKVSYFDHSIVGAKKVLLKANINQKIVKLNFKQIQD IIWTKIDNYQLEIRLGVDFVTTIPNIQIFILSEHKNPVSLPI DNKSTENNQNKLKLLDELKELGMKYVKHQNQIY 1 配列をコピーする ( > の行は入れても入れなくてもよい ) 3Favorite organisms を選択 4 mge mpn uur と入力 mge: Mycoplasma genitalium mpn: Mycoplasma pneumoniae uur: Ureaplasma parvum 6 通りの reading frame のすべてについて翻訳 5 Compute を押すデータベース上の塩基配列も,6 通りの reading frame のすべてについて翻訳し, このアミノ酸配列データに対して検索質問配列, データベースとも, アノテーション情報が整備されていない場合に有効 Ureaplasma は,ftsZ を持っていないことがわかる

スタートすべてのプログラムアクセサリコマンドプロンプト C: Users iu> > > blastp -help 37 38

83 1995 Mycoplasma genitalium 0.58 1995 Mycoplasma pneumoniae 0.

10 スタートすべてのプログラムアクセサリコマンドプロンプト C: Users iu> > > blastp -help Windows の場合は, このファイルをダウンロードします Haemophilus influenzae Mycoplasma genitalium Mycoplasma pneumoniae Bacillus subtilis Escherichia coli Ureaplasma parvum

11 test1.seq test2.seq test3.seq Mgenitalium.faa Mpneumoniae.faa Ureaplasma.faa parse-blast7.pl の7つのファイルをダウンロードして作成したblastフォルダにれてくださいデータベースの準備 > cd C: Users iu Desktop blast C: Users iu Desktop blast> > more Mgenitalium.faa > dir more コマンドについて 2009/03/11 19:52 <DIR>. 2009/03/11 19:52 <DIR> /04/21 23:34 222,447 Mgenitalium.faa 2005/04/21 23:33 307,006 Mpneumoniae.faa 43 指定したファイルの内容を表示します. 次ページを見るには [Space] キー, 1 行ずつ見るには [Enter] キー, 終了するには [Q] キー押します. 44

12 データベースの準備 46 > makeblastdb -in Mgenitalium.faa -dbtype prot > more test1.seq >gi ref NP_ uracil-dna-glycosylase [Escherichia coli str. K-12 substr. MG1655] MANELTWHDVLAEEKQQPYFLNTLQTVASERQSGVTIYPPQKDVFNAFRFTELG DVKVVILGQDPYHGPGQAHGLAFSVRPGIAIPPSLLNMYKELENTIPGFTRPNH GYLESWARQGVLLLNTVLTVRAGQAHSHASLGWETFTDKVISLINQHREGVVFL LWGSHAQKKGAIIDKQRHHVLKAPHPSPLSAHRGFFGCNHFVLANQWLEQRGET PIDWMPVLPAESE 45 > blastp db Mgenitalium.faa query test1.seq > blastp db Mgenitalium.faa query test1.seq out result1.txt > more result1.txt -out : 出力ファイル指定 ( 上矢印 ) を押すと, 過去に入力したコマンドが出てきます 47 > blastp db Mgenitalium.faa query test1.seq > result1.txt 48

13 質問配列の名前検索対象として用いたデータベースアラインメント BLASTP [Oct ] 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: Query= hsa:7100 TLR5; toll-like receptor 5 (A) (858 letters) Database: nr-aa: Non-redundant protein sequence database Release ,952,394 sequences; 634,153,439 total letters スコア Searching...done Score E Sequences producing significant alignments: (bits) Value Top 100Top 50Top 20Top 10Top 5 Select operationclustalwmafftprrndraw alignmentsearch common motifs(pfam)search common motifs(prosite) gp:ab060695_1 [AB060695] Toll-like receptor 5 [Homo sapiens] sp:tlr5_human [O60602] Toll-like receptor 5 precursor (Toll/inte gp:ax590493_1 [AX590493] Sequence 5 from Patent WO [Hom gpu:ab208697_1 [AB208697] Toll-like receptor 5 [Sus scrofa] sp:tlr5_mouse [Q9JLF7] Toll-like receptor 5 precursor.>prf: tr:q8cb40_mouse [Q8CB40] Mus musculus adult female vagina cdna, tr:q5gda9_chick [Q5GDA9] Toll-like receptor tr:q5gr02_chick [Q5GR02] Toll-like receptor 5 precursor.>gpu:aj tr:q5u5b1_xenla [Q5U5B1] LOC protein.>gpu:bc084773_1 [BC gp:cq870716_1 [CQ870716] Sequence 9 from Patent EP [unid prf: a membrane-toll-like receptor - Oncorhynchus mykiss ( e-165 tr:q5h720_fugru [Q5H720] TLR5.>gpu:AC156437_1 [AC156437] TLR5 [T e-154 gp:ax590495_1 [AX590495] Sequence 7 from Patent WO [syn e-121 tr:q7zt81_oncmy [Q7ZT81] Toll-like receptor5.>gp:ab062504_1 [AB e-99 >gp:ab060695_1 [AB060695] Toll-like receptor 5 [Homo sapiens] Top Length = 858 Score = 1666 bits (4315), Expect = 0.0 Identities = 827/844 (97%), Positives = 827/844 (97%) Query: 15 AGPVFGIPSCSFDGRIAFYRFCNLTQVPQVLNTTERLLLSFNYIRTVTASSFPFXXXXXX 74 AGPVFGIPSCSFDGRIAFYRFCNLTQVPQVLNTTERLLLSFNYIRTVTASSFPF Sbjct: 15 AGPVFGIPSCSFDGRIAFYRFCNLTQVPQVLNTTERLLLSFNYIRTVTASSFPFLEQLQL 74 Query: 75 XXXGSQYTPLTIDKEAFRNLPNLRILDLGSSKIYFLHPDAFQGLFHLFELRLYFCGLSDA 134 GSQYTPLTIDKEAFRNLPNLRILDLGSSKIYFLHPDAFQGLFHLFELRLYFCGLSDA Sbjct: 75 LELGSQYTPLTIDKEAFRNLPNLRILDLGSSKIYFLHPDAFQGLFHLFELRLYFCGLSDA 134 E value E value 設定 > blastp db Mgenitalium.faa query test1.seq out result1.txt evalue 1e-10 > more result1.txt 1 l 50 BLASTX 大量 Query のホモロジー検索法 > more test2.seq > blastx db Mgenitalium.faa query test2.seq evalue 1e-10 out result2.txt > more result2.txt >gi ref NP_ phosphofructokinase II [Escherichia coli K12] MVRIYTLTLAPSLDSATITPQIYPEGKLRCTAPVFEPGGGGINVARAIAHLGGSATAIFPAGGATGEHLV SLLADENVPVATVEAKDWTRQNLHVHVEASGEQYRFVMPGAALNEDEFRQLEEQVLEIESGAILVISGSL PPGVKLEKLTQLISAAQKQGIRCIVDSSGEALSAALAIGNIELVKPNQKELSALVNRELTQPDDVRKAAQ EIVNSGKAKRVVVSLGPQGALGVDSENCIQVVPPPVKSQSTVGAGDSMVGAMTLKLAENASLEEMVRFGV AAGSAATLNQGTRLCSHDDTQKIYAYLSR >gi ref NP_ phosphoglyceromutase 2 [Escherichia coli K12] MLQVYLVRHGETQWNAERRIQGQSDSPLTAKGEQQAMQVATRAKELGITHIISSDLGRTRRTAEIIAQAC GCDIIFDSRLRELNMGVLEKRHIDSLTEEEENWRRQLVNGTVDGRIPEGESMQELSDRVNAALESCRDLP QGSRPLLVSHGIALGCLVSTILGLPAWAERRLRLRNCSISRVDYQESLWLASGWVVETAGDISHLDAPAL DELQR 51 >gi ref NP_ glucosephosphate isomerase [Escherichia coli K12] MKNINPTQTAAWQALQKHFDEMKDVTIADLFAKDGDRFSKFSATFDDQMLVDYSKNRITEETLAKLQDLA KECDLAGAIKSMFSGEKINRTENRAVLHVALRNRSNTPILVDGKDVMPEVNAVLEKMKTFSEAIISGEWK GYTGKAITDVVNIGIGGSDLGPYMVTEALRPYKNHLNMHFVSNVDGTHIAEVLKKVNPETTLFLVASKTF TTQETMTNAHSARDWFLKAAGDEKHVAKHFAALSTNAKAVGEFGIDTANMFEFWDWVGGRYSLWSAIGLS IVLSIGFDNFVELLSGAHAMDKHFSTTPAEKNLPVLLALIGIWYNNFFGAETEAILPYDQYMHRFAAYFQ QGNMESNGKYVDRNGNVVDYQTGPIIWGEPGTNGQHAFYQLIHQGTKMVPCDFIAPAITHNPLSDHHQKL LSNFFAQTEALAFGKSREVVEQEYRDQGKDPATLDYVVPFKVFEGNRPTNSILLREITPFSLGALIALYE HKIFTQGVILNIFTFDQWGVELGKQLANRILPELKDDKEISSHDSSTNGLINRYKAWRG 52

14 大量 Query のホモロジー検索法 > more test3.seq > blastp db Mgenitalium.faa query test3.seq evalue 1e-10 out result3.txt > more result3.txt 生物生物 > more Mpneumoniae.faa > blastp db Mgenitalium.faa query Mpneumoniae.faa evalue 1e-10 out result4.txt > more result4.txt 55 56

BLASTP 2.2.5 [Nov-16-2002] Reference: Altschul, Stephen F., Thomas L. Madden, Alejandro A. Schaffer, Jinghui Zhang, Zheng Zhang, Webb Miller, and David J.

15 BLASTP [Nov ] 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: Query= gi ref NP_ glucosephosphate isomerase [Escherichia coli K12] (549 letters) Database: yeast.aa 6298 sequences; 2,974,038 total letters > more parse-blast7.pl Score E Sequences producing significant alignments: (bits) Value ref NP_ Glucose-6-phosphate isomerase; Pgi1p ref NP_ Ygr130cp ref NP_ spindle pole body component; Stu2p ref NP_ (putative) involved in cell wall biogenesis; Ec ref NP_ Ylr419wp >ref NP_ Glucose-6-phosphate isomerase; Pgi1p Length = 554 parse-blast.pl #! /usr/local/bin/perl use strict; use warnings; use Getopt::Std; Score = 641 bits (1654), Expect = 0.0 Identities = 326/549 (59%), Positives = 401/549 (73%), Gaps = 16/549 (2%) Query: 7 TQTAAWQALQKHFDEM-KDVTIADLFAKDGDRFSKFSATFDD----QMLVDYSKNRITEE 61 T+ AW LQK ++ K +++ F KD RF K + TF + ++L DYSKN + +E Sbjct: 13 TELPAWSKLQKIYESQGKTLSVKQEFQKDAKRFEKLNKTFTNYDGSKILFDYSKNLVNDE 72 Query: 62 TLAKLQDLAKECDLAGAIKSMFSGEKINRTENRAVLHVALRNRSNTPILVDGKDVMPEVN 121 +A L +LAKE ++ G +MF GE IN TE+RAV HVALRNR+N P+ VDG +V PEV+ Sbjct: 73 IIAALIELAKEANVTGLRDAMFKGEHINSTEDRAVYHVALRNRANKPMYVDGVNVAPEVD my $mode = 0; my $name = "";.. 58 > perl parse-blast7.pl -i result4.txt -o list1.txt スタート list1.txt をExcel 上にすべてのプログラムドラッグ & ドロップ Microsoft Office Microsoft Office Excel 質問配列の情報 BLAST 検索でヒットした配列の情報 ( ヒットしなかった場合は空欄 ) スコア,E-value,Identity Ureaplasma.faaには,Ureaplasma parvumのゲノムにコードされる全タンパク質がmulti-fastaフォーマットで記述してあります Mpneumoniae.faa をデータベース, Ureaplasma.faa を質問配列にいてBLAST 検索をい,Ureaplasmaがコードするタンパク質と相同なものがM. pneumoniaeゲノム上にもあるかどうか, 調べてください (E-valueの閾値は,1e-3に設定してください) M. genitalium ゲノム上には, これらと相同なタンパク質がコードされていない parse-blast7.plを使って, ヒットしたアミノ酸配列のリストを作成してください作成したエクセルファイルを提出してください 60

<4D F736F F F696E74202D F90B695A8947A97F189F090CD8AEE91625F3189F196DA5F8E9197BF2E >

<4D F736F F F696E74202D F90B695A8947A97F189F090CD8AEE91625F3189F196DA5F8E9197BF2E > 1 2 kiso1 3 4 5 6 7 8 9 10 11 12 AP009356 と入力 National Center for Biotechnology Information http://www.ncbi.nlm.nih.gov/ 13 All Databases データベースの統合検索システム主なデータベースは,PubMed ヌクレオチドシークエンスデータベースタンパク質シークエンスデータベース