在如何在IGV上使用BLAT搜索非模式物種中，我們討論如何用Apache的CGI服務(wù)器，來響應(yīng)IGV發(fā)送的BLAT請(qǐng)求?？紤]到我們大部分的時(shí)候都是個(gè)人使用，并不需要Apache這種重量級(jí)的Web服務(wù)器，因此完全可以省去這一組件。

我之前學(xué)過Python和Django和Flask這兩種網(wǎng)頁應(yīng)用框架，其中Flask比較輕量，非常適合我們這種小型應(yīng)用。因此，我就用Flask編寫了一個(gè)Blat網(wǎng)頁應(yīng)用，用來響應(yīng)IGV的請(qǐng)求。

首先，安裝Flask（為了避免沖突，我用了虛擬環(huán)境)

mkdir myproject
cd myproject
python3 -m venv venv

# 啟動(dòng)
. venv/bin/activate
pip3 install flask

接著，從我的GitHub上下載腳本，見https://github.com/xuzhougeng/myscripts/blob/master/misc/myBlat.py

或者復(fù)制下面的代碼進(jìn)行修改，命名為myBlast.py

from flask import Flask
from flask import request
import tempfile
import subprocess
from os.path import basename

app = Flask(__name__)


"""
注意:
    1. 修改ref_dict中的ref_name1和ref_name2fa的位置
    2. 修改query_seq_db中的 ref1 和 ref2, 改成你的物種里的特有關(guān)鍵詞, 例如A.thaliana 我用的是thaliana

"""


ref_dict = {
    "ref_name1":  "/path/to/your/reference1.fa",
    "ref_name2":  "/path/to/your/reference2.fa"
}

def query_seq_db(db):
    #if "\\" in db:
    #    dbname = db.split("\\")[-1]
    #else:
    #    dbname = basename(db)
    #dbname = dbname.split(".")[0]
    #ref = ref_dict[dbname]
    #return ref
    db = db.lower()
    if "ref1" in db:
        return ref_dict["ref_name1"]
    elif "ref2" in db:
        return ref_dict["ref_name2"]
    else:
        return ref_dict["ref_name1"]
   
    

def blat(seq,db,type):
    # 需要改成你實(shí)際的blat路徑
    blat = "/opt/biosoft/ucsctools/blat/blat"
    # 需要改成你實(shí)際的參考位置
    ref = query_seq_db(db)

    fd1, fa_file = tempfile.mkstemp(suffix=".fa")
    fd2, psl_file = tempfile.mkstemp(suffix=".psl")

    with open(fa_file,"w") as f:
        f.write(f">query\n{seq}\n")
    process = subprocess.Popen([blat, ref, fa_file, psl_file],
            stdout=subprocess.PIPE,
            stderr=subprocess.PIPE)
    
    stdout,stderr = process.communicate()

    psl_ret = [line for line in open(psl_file, "r") ]
    psl_rec = [line.strip().split('\t') for line in psl_ret[5:] ]

    return psl_rec

@app.route("/myBlat", methods=['GET'])
def myBlat():

    seq = request.args.get('userSeq','')
    db = request.args.get('db','')
    output = request.args.get('output','json')
    type = request.args.get('type', 'DNA')

    psl_rec = blat(seq, db, type)
    data =  {
        "track": "blat",
        "genome": "ath",
        "fields": ["matches", "misMatches", "repMatches", "nCount", "qNumInsert", "qBaseInsert", \
                "tNumInsert", "tBaseInsert", "strand", "qName", "qSize", "qStart", "qEnd", "tName", \
                "tSize", "tEnd", "blockCount", "blockSizes", "qStarts", "tStarts"],
        "blat" : psl_rec
        }
    return data

```

最后運(yùn)行Flask

```nginx
flask run -p 5000
```

需要注意的是，我用的端口是5000，而不是80，因此還需要在IGV偏好中BLAT URL改成

```nginx
http://127.0.0.1:5000/myBlat?userSeq=$SEQUENCE&type=DNA&db=$DB&output=json
```