實(shí)驗(yàn)中經(jīng)常需要從基因組中提取蚊逢、序列用于設(shè)計(jì)引物，或者進(jìn)行基因功能研究等箫章，之前看到別人使用的一個(gè)腳本烙荷，自己拿來(lái)用一下，感覺(jué)確實(shí)很方便
在網(wǎng)站 http://hgdownload.cse.ucsc.edu/admin/exe/linux.x86_64/（網(wǎng)站內(nèi)還有很多其他腳本可以下載使用檬寂，有時(shí)間可以研究一下）
下載faSomeRecords腳本：faSomeRecords.txt
faSomeRecords
將txt后綴刪除后拷貝到指定的文件夾下即可運(yùn)行终抽，運(yùn)行命令行加載腳本：
chmod +x faSomeRecords
# 賦予文件可執(zhí)行權(quán)限，為L(zhǎng)inux系統(tǒng)下執(zhí)行
執(zhí)行如下命令桶至，找到自己需要的序列昼伴，一般使用CDs序列文件
./faSomeRecords genome.cds.fa ID.txt out.gene.fa # 執(zhí)行中 ./ 不能刪除
# 其中g(shù)enome.cds.fa是原始的fasta文件,包含自己需要的基因
# ID.txt 是需要查找并導(dǎo)出的基因ID，每行一個(gè)
# out.gene.fa 為輸出文件镣屹，包含對(duì)應(yīng)ID的序列信息圃郊。
CDs文件示意圖
ID.txt文件示意圖
out.gene.fa輸出文件示意圖
