24小时热门版块排行榜    

查看: 946  |  回复: 5

fang8232335

银虫 (初入文坛)

[交流] 全基因组测序 已有4人参与

新手,自己工科出身,但现在老板让弄分子生物学,测一个细菌的全基因组。现在拼接好的序列已经有了,想问下,我现在想查下这个全基因组序列里有没有我想要的基因或蛋白质序列要怎么操作哇(比如我知道这个菌有农药的降解效果,降解酶的蛋白质序列已经有了,现在想拿我测的基因组序列跟这个进行比对,看是不是有这个序列,以及确定这个序列的位置),谢谢大家了!自己用NCBI上blast怎么也用不好
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

growlywolf

金虫 (小有名气)


小木虫: 金币+0.5, 给个红包,谢谢回帖
用基因组序列建立本地blast数据库,以蛋白序列为query做本地tblast

[ 发自小木虫客户端 ]
2楼2015-05-10 10:54:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fang8232335

银虫 (初入文坛)

引用回帖:
2楼: Originally posted by growlywolf at 2015-05-10 10:54:24
用基因组序列建立本地blast数据库,以蛋白序列为query做本地tblast

我后来又找到了那个基因的序列,然后用BioEdit 进行本地blast,得到下面的数据,这是代表有还是没有哇,感觉都是片段匹配,是什么意思哇?
BLASTN 2.2.10 [Oct-19-2004]


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= gi|148712|gb|M29593.1|FVBOPD Flavobacterium sp. parathion
hydrolase gene, complete cds
         (1693 letters)

Database: D:\BioEdit\database\X1-01.fas
           1 sequences; 2,450,153 total letters



                                                                 Score    E
Sequences producing significant alignments:                      (bits) Value

X1-01.seq                                                              34   0.21

>X1-01.seq
          Length = 2450153

Score = 34.2 bits (17), Expect = 0.21
Identities = 17/17 (100%)
Strand = Plus / Minus

                              
Query: 466    tctgctcggcggcctgg 482
              |||||||||||||||||
Sbjct: 894223 tctgctcggcggcctgg 894207



Score = 34.2 bits (17), Expect = 0.21
Identities = 17/17 (100%)
Strand = Plus / Plus

                              
Query: 986    gccagcttggccaccgg 1002
              |||||||||||||||||
Sbjct: 659980 gccagcttggccaccgg 659996



Score = 32.2 bits (16), Expect = 0.83
Identities = 16/16 (100%)
Strand = Plus / Minus

                              
Query: 83      ggcttgccagcagccg 98
               ||||||||||||||||
Sbjct: 1669234 ggcttgccagcagccg 1669219



Score = 32.2 bits (16), Expect = 0.83
Identities = 16/16 (100%)
Strand = Plus / Plus

                              
Query: 449     gcggccgccgcaggaa 464
               ||||||||||||||||
Sbjct: 1562756 gcggccgccgcaggaa 1562771



Score = 32.2 bits (16), Expect = 0.83
Identities = 16/16 (100%)
Strand = Plus / Minus

                              
Query: 860    cagttcttcctgcgtg 875
              ||||||||||||||||
Sbjct: 699001 cagttcttcctgcgtg 698986



Score = 32.2 bits (16), Expect = 0.83
Identities = 16/16 (100%)
Strand = Plus / Plus

                              
Query: 486    ggtgcgcgagcgtggc 501
              ||||||||||||||||
Sbjct: 393049 ggtgcgcgagcgtggc 393064



Score = 32.2 bits (16), Expect = 0.83
Identities = 19/20 (95%)
Strand = Plus / Minus

                                 
Query: 467   ctgctcggcggcctggctgg 486
             ||||||||||||||| ||||
Sbjct: 25911 ctgctcggcggcctgactgg 25892


  Database: D:\BioEdit\database\X1-01.fas
    Posted date:  May 10, 2015  1:46 PM
  Number of letters in database: 2,450,153
  Number of sequences in database:  1
  
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 Hits to DB: 3969
Number of Sequences: 1
Number of extensions: 3969
Number of successful extensions: 7
Number of sequences better than  1.0: 1
Number of HSP's better than  1.0 without gapping: 1
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): 7
length of query: 1693
length of database: 2,450,153
effective HSP length: 16
effective length of query: 1677
effective length of database: 2,450,137
effective search space: 4108879749
effective search space used: 4108879749
T: 0
A: 0
X1: 11 (21.8 bits)
X2: 15 (29.7 bits)
S1: 12 (24.3 bits)
S2: 16 (32.2 bits)
3楼2015-05-10 16:19:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nanhu1984

银虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
最好选择BlastP或者BlastX,这样获得的结果更多,灵敏度更高
https://www.microbialgenomic.cn/
4楼2015-05-12 11:24:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

王志云006

新虫 (初入文坛)

本帖内容被屏蔽

5楼2015-05-12 20:19:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ruociwang

铁虫 (初入文坛)

MicrobeTrakr
不论付出多少努力,一定要把上海生命科学院的大门给撞开了!
6楼2015-06-02 12:08:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fang8232335 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 关于filecode,很负责任的告诉大家 +6 爱看书的可乐 2026-08-08 7/350 2026-08-08 22:13 by a_niu
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 qTvzQBRAHjCi 2026-08-07 4/200 2026-08-08 17:22 by oEVWOejN9taj
[基金申请] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 qTvzQBRAHjCi 2026-08-07 5/250 2026-08-08 17:02 by oEVWOejN9taj
[考研] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +3 qTvzQBRAHjCi 2026-08-07 5/250 2026-08-08 16:47 by oEVWOejN9taj
[考博] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:27 by oEVWOejN9taj
[基金申请] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:19 by oEVWOejN9taj
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 Vi50GxzrFcSG 2026-08-07 4/200 2026-08-08 16:02 by oEVWOejN9taj
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 5/250 2026-08-08 15:59 by oEVWOejN9taj
[论文投稿] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 WQyTGMbfH7Cx 2026-08-07 4/200 2026-08-08 15:02 by oEVWOejN9taj
[基金申请] 这样的filecode谁见过 +5 布布和一二 2026-08-08 8/400 2026-08-08 14:54 by archvillain
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 HEQlVqMTIA7d 2026-08-07 5/250 2026-08-08 14:27 by oEVWOejN9taj
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 KXLV3nuBVBY7 2026-08-07 5/250 2026-08-08 14:07 by oEVWOejN9taj
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 Vi50GxzrFcSG 2026-08-07 6/300 2026-08-08 10:19 by 3OOjAIS77qg2
[基金申请] 基金中了 +14 laoda193707 2026-08-06 14/700 2026-08-08 00:23 by 实验小白ha
[基金申请] 听说今天filecode变了 +24 布布和一二 2026-08-06 47/2350 2026-08-07 16:02 by zhiyanjiang
[基金申请] filecode变化情况 +6 布布和一二 2026-08-07 22/1100 2026-08-07 14:45 by 且听虎啸
[基金申请] 求各位大神看下 100+6 hpkpkpkp 2026-08-05 33/1650 2026-08-06 14:49 by zhiyanjiang
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 纯娱乐,不喜欢勿喷 +7 Tide man 2026-08-04 10/500 2026-08-04 15:10 by loufangrui
[基金申请] 什么时候能放榜呀? +3 Jacob678 2026-08-03 3/150 2026-08-03 16:14 by gltch
信息提示
请填处理意见