24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 867  |  回复: 5
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

fang8232335

银虫 (初入文坛)

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

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

» 猜你喜欢

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

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

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

nanhu1984

银虫 (正式写手)


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

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的回帖

王志云006

新虫 (初入文坛)

本帖内容被屏蔽

5楼2015-05-12 20:19:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 275求调剂 +10 Micky11223 2026-03-25 13/650 2026-03-27 22:42 by Micky11223
[考研] 265求调剂 +8 小木虫085600 2026-03-27 8/400 2026-03-27 22:16 by 无际的草原
[考研] 287求调剂 +10 land xuxu 2026-03-26 10/500 2026-03-27 15:33 by 帕尔马拉特
[考研] 安徽大学专硕生物与医药专业(086000)324分,英语已过四六级,六级521,求调剂 +4 美味可乐鸡翅 2026-03-26 4/200 2026-03-27 15:27 by 星空星月
[考研] 08开头275求调剂 +4 拉谁不重要 2026-03-26 4/200 2026-03-27 14:12 by Delta2012
[考研] 化学308分求调剂 +8 你好明天你好 2026-03-23 9/450 2026-03-27 14:01 by 杨光于青云
[考研] 348求调剂 +4 小懒虫不懒了 2026-03-27 5/250 2026-03-27 12:47 by 果果妈咪
[考研] 303求调剂 +7 安忆灵 2026-03-22 8/400 2026-03-27 11:46 by sanrepian
[考研] 312求调剂 +9 上岸吧ZJY 2026-03-22 13/650 2026-03-27 11:24 by sanrepian
[考研] 一志愿武汉理工,总分321,英一数二,求老师收留。 +5 nnnnnnn5 2026-03-25 5/250 2026-03-27 04:42 by wxiongid
[考研] 349求调剂 +5 杰斯塔里斯 2026-03-21 5/250 2026-03-27 00:31 by wxiongid
[考研] 一志愿 南京邮电大学 288分 材料考研 求调剂 +3 jl0720 2026-03-26 3/150 2026-03-26 13:39 by zzll406
[考研] 一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研经历 +6 cccchenso 2026-03-23 6/300 2026-03-25 22:25 by 544594351
[考研] 求调剂 +3 QiMing7 2026-03-25 3/150 2026-03-25 21:13 by 给你你注意休息
[考研] 网络空间安全0839招调剂 +4 w320357296 2026-03-25 6/300 2026-03-25 17:59 by 255671
[考研] 318求调剂 +5 plum李子 2026-03-21 8/400 2026-03-25 09:26 by aa331100
[考研] 311求调剂 +3 冬十三 2026-03-24 3/150 2026-03-24 21:31 by peike
[考研] 307求调剂 +3 余意卿 2026-03-21 6/300 2026-03-24 15:03 by 余意卿
[论文投稿] 急发核心期刊论文 +3 贤达问津 2026-03-23 5/250 2026-03-23 17:13 by 妹子不好惹
[考研] 求老师收我 +3 zzh16938784 2026-03-23 3/150 2026-03-23 12:56 by ztnimte
信息提示
请填处理意见