24小时热门版块排行榜    

查看: 949  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 小木虫上这么多卖论文的,真有人买论文么?感觉没必要啊 +4 Tide man 2026-08-10 4/200 2026-08-11 07:45 by emeraldine
[基金申请] 帮忙看看fileCode +4 wwncly 2026-08-10 8/400 2026-08-10 22:07 by demonstreets
[基金申请] FileCode能看出啥? +7 要乐观耀哥 2026-08-10 15/750 2026-08-10 19:29 by 要乐观耀哥
[基金申请] 确定了,国自然21号放榜 +6 布布和一二 2026-08-10 7/350 2026-08-10 19:15 by 2000zf36392
[基金申请] 关于代码变化问题,想知道的进来 +17 且听虎啸 2026-08-07 24/1200 2026-08-10 18:35 by zhangduo2008
[基金申请] 静等基金结果 +5 gjjjzhong 2026-08-10 16/800 2026-08-10 17:19 by Tide man
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +14 Lanmanbaby 2026-08-09 22/1100 2026-08-10 14:57 by wwncly
[基金申请] 好奇怪的filecode +4 布布和一二 2026-08-08 5/250 2026-08-10 14:20 by 冰心玉壶晴
[基金申请] 这样的filecode谁见过 +11 布布和一二 2026-08-08 22/1100 2026-08-10 11:10 by wmfsnow
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
[基金申请] 关于豆爷回答的JTJC与%2F数量 +5 yang182083 2026-08-06 7/350 2026-08-08 18:29 by zhanghaozhu
[基金申请] 固定端突然变了,今天 +6 archvillain 2026-08-06 10/500 2026-08-07 16:03 by 医学老男孩
[基金申请] 听说今天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 且听虎啸
[基金申请] 大家散了吧,后缀研究没有意义,别浪费时间了,过好目前的每一天,不要焦虑 +5 Tide man 2026-08-06 7/350 2026-08-07 13:11 by 医学老男孩
[基金申请] filecode +14 等待解的谜 2026-08-06 19/950 2026-08-07 12:20 by wlwhappy
[基金申请] filecode +8 布布和一二 2026-08-06 11/550 2026-08-06 20:41 by tangpu318
[基金申请] 影响面上的因素 +8 布布和一二 2026-08-05 11/550 2026-08-06 10:41 by 宝贝虫子
[基金申请] 好消息?这个有何含义??? +8 Tide man 2026-08-05 10/500 2026-08-05 16:14 by xmuxiaoyu
信息提示
请填处理意见