当前位置:网站首页>Solve blast database error: error pre fetching sequence data
Solve blast database error: error pre fetching sequence data
2022-07-28 14:55:00 【javafalcon】
1、 background :
download Uniprot From the protein sequence in the database to the local , Use this file to create BLAST Search database for , And then use blastp Command to search the database .
2、 error
according to blast Instructions in the manual , To speed up the search process , Use mask file . from Uniprot Downloaded sequence file sequences.fasta, The content is as follows :
>sid1
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
>sid2
xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx
.................................................................
1) establish mask file :
segmasker -in sequences.fasta -infmt fasta -outfmt maskingfo_asnl_bin \
-out refsequences.asnb -parse_seqids
2) Create database
makeblastdb -in sequences.fasta -input_type fasta -dbtype prot -out uniprot \
-title "Uniprot" -parse_seqids
3) Sequence alignment
blastp -db uniprot -query input.fasta -out output.txt -outfmt 7
The result is wrong :BLAST Database error: Error pre-fetching sequences data
3、 resolvent
Creating mask Files and databases are not used parse_seqids.
segmasker -in sequences.fasta -infmt fasta -outfmt maskingfo_asnl_bin \
-out refsequences.asnb
makeblastdb -in sequences.fasta -input_type fasta -dbtype prot -out uniprot -title "Uniprot"
Why? ? I don't know why , Looking up many articles on the Internet can't solve , I guess there is something wrong with creating the database , Try to remove unnecessary parameters one by one , It was found that parse_seqids After the parameter , Can use balstp Get the comparison results correctly .
边栏推荐
- Qtableview in QT sets three methods of paging display [easy to understand]
- 八、picker用法 下拉框选择效果
- 8、 Picker usage drop-down box selection effect
- 如何只降3D相机不降UI相机的分辨率
- SSL socket cross platform solution libevent OpenSSL
- Installing redis in Linux
- Afnetworking crash course
- SwiftUI 4.0 的全新导航系统
- VTK vtkcontourwidget extracts regions of interest
- Namespace conflict problem
猜你喜欢

10、 Timestamp

It's so hot that solar power can't take off? Hello, head

21、 TF coordinate transformation (I): coordinate MSG message

企鹅一面:为什么不建议使用SELECT * ?

VTK vtkcontourwidget extracts regions of interest

使用Weka与Excel进行简单的数据分析

Creating, deleting and viewing Anaconda virtual environment

【七夕】七夕孤寡小青蛙究极版?七夕节最终章!

2022 safety officer-a certificate operation certificate examination question bank simulated examination platform operation

爆肝整理JVM十大模块知识点总结,不信你还不懂
随机推荐
面试官:ThreadLocal使用场景有哪些?内存泄露问题如何避免?
How to effectively conduct the review meeting (Part 1)?
看了就会的 Rainbond 入门教程
19、 ROS parameter name setting
58 sub station Anju, broker marketing management platform login interface encryption reverse
8、 C scope rules
卡方分布和伽马函数(Chi-Square Distribution)
&0xffffffff(0x08)
SwiftUI 布局 —— 尺寸( 上 )
On July 29, apachecon | apachepulsar's exploration and practice in vivo will be broadcast soon
Qt development tips
Namespace conflict problem
OKR and grad
如何只降3D相机不降UI相机的分辨率
6、 C language circular statement
Core Data 是如何在 SQLite 中保存数据的
C language exercises
如何让照片中的人物笑起来?HMS Core视频编辑服务一键微笑功能,让人物笑容更自然
14、 ROS meta function package
Creating, deleting and viewing Anaconda virtual environment