makeblastdb注意事項
本地blast,需要自行構建數據庫。那麼就要用到makeblastdb。這個命令很簡單。
USAGE
makeblastdb [-h] [-help] [-in input_file] [-input_type type]
-dbtype molecule_type [-title database_title] [-parse_seqids]
[-hash_index] [-mask_data mask_data_files] [-mask_id mask_algo_ids]
[-mask_desc mask_algo_descriptions] [-gi_mask]
[-gi_mask_name gi_based_mask_names] [-out database_name]
[-max_file_sz number_of_bytes] [-logfile File_Name] [-taxid TaxID]
[-taxid_map TaxIDMapFile] [-version]
DESCRIPTION
Application to create BLAST databases, version 2.4.0+
Use '-help' to print detailed descriptions of command line arguments
建庫命令如下:
makeblastdb -in IWGSC_v1.0_blastdb.fasta -dbtype nucl -title CS_v1.0_full -parse_seqids -out CS_v1.0_full
makeblastdb時一定要加上 -parse_seqids ,否則調取序列會有問題。調取序列的命令是:
blastdbcmd -db CS_v1.0_full -entry chr1A -range 2-100 #因未寫輸出文件,默認輸出到屏幕
>chr1A
TAAACCCTAAACCCTAAACCCTAAACCCTAAACCCCTAACCCTAAACCCTAACCAAAACCCTAAACCCTAAACCCCTAAA
CCCTAAACCCTAACCTAAA
如果不加 -parse_seqids,會出現 OID not found的錯誤。
當要調取的序列較大時可以將序列輸出☞文件,見如下命令。然後
blastdbcmd -db /data2/Fshare/IWGSC_v1.0_Formatdb/CS_v1.0_merge -entry chr1A_part1 -range 100000-110000 -out result.fasta
scp mollab@192.168.3.253:/home/mollab/result.fasta ~/Desktop/ #此條命令需要在本機執行,不是在服務器端執行,並且需要輸入密碼,文件保存在你自己電腦的桌面上