import logging
import os
try:
from urllib.parse import unquote as unquote
except ImportError:
from urllib import unquote as unquote
from flask import Flask, render_template, request, jsonify, Response, session
from svviz import export
logging.getLogger('werkzeug').setLevel(logging.ERROR)
dataHub = None
# Initialize the Flask application
app = Flask(__name__,
static_folder=os.path.join(os.path.dirname(__file__), "static"),
template_folder=os.path.join(os.path.dirname(__file__), "templates")
)
app.secret_key = '\xfd{H\xe5<\x95\xf9\xe3\x96.5\xd1\x01O<!\xd5\xa2\xa0\x9fR"\xa1\xa8'
def getRandomPort():
import socket
sock = socket.socket(socket.AF_INET, socket.SOCK_STREAM)
sock.bind(("", 0))
port = sock.getsockname()[1]
sock.close()
return port
def checkPortIsClosed(port):
import socket;
sock = socket.socket(socket.AF_INET, socket.SOCK_STREAM)
result = sock.connect_ex(('127.0.0.1', port))
if result == 0:
return False
else:
return True
@app.route('/')
def index():
if not "last_format" in session:
session["last_format"] = "svg"
session.permanent = True
try:
variantDescription = str(dataHub.variant).replace("::", " ").replace("-", "–")
return render_template('index.html',
samples=list(dataHub.samples.keys()),
annotations=dataHub.annotationSets,
results_table=dataHub.getCounts(),
insertSizeDistributions=[sample.name for sample in dataHub if sample.insertSizePlot],
dotplots=dataHub.dotplots,
variantDescription=variantDescription)
except Exception as e:
logging.error("ERROR:{}".format(e))
raise
@app.route('/<path:path>')
def static_proxy(path):
# send_static_file will guess the correct MIME type
return app.send_static_file(path)
@app.route('/_export', methods=["POST"])
def do_export():
format = request.form.get("format", "svg").lower()
session["last_format"] = format
filename = "export.{}".format(format)
svg = dataHub.trackCompositor.render()
converter = export.getExportConverter(dataHub.args, format)
if format == "svg":
mimetype = "image/svg+xml"
data = svg
elif format == "png":
mimetype = "image/png"
data = export.convertSVG(svg, "png", converter)
elif format == "pdf":
mimetype = "application/pdf"
data = export.convertSVG(svg, "pdf", converter)
else:
raise Exception("unknown format")
response = Response(data,
mimetype=mimetype,
headers={"Content-Disposition": "attachment;filename={}".format(filename)})
return response
@app.route('/_haspdfexport')
def _hasPDFExport():
if export.getExportConverter(dataHub.args, "pdf"):
return jsonify({"haspdfexport":True})
return jsonify({"haspdfexport":False})
@app.route('/_haspngexport')
def _hasPNGExport():
if export.getExportConverter(dataHub.args, "png"):
return jsonify({"haspngexport":True})
return jsonify({"haspngexport":False})
def _getsvg(track):
track.render()
svgText = track.svg.asString("web")
return svgText
@app.route('/_disp')
def display():
req = request.args.get('req', 0)
if req == "progress":
return jsonify(result="done")
if req in ["alt", "ref", "amb"]:
allele = req
results = []
for name, sample in dataHub.samples.items():
# svg = open("{}.{}.svg".format(req, name)).read()
track = sample.tracks[allele]
track.render()
svg = track.svg.asString("web")
results.append({"name":name, "svg":svg})
for annotation in dataHub.alleleTracks[allele]:
track = dataHub.alleleTracks[allele][annotation]
track.render(spacing=5)
annoSVG = track.svg.asString("web")
results.append({"name":annotation, "svg":annoSVG})
return jsonify(results=results)
if req == "counts":
return jsonify(result=dataHub.getCounts())
return jsonify(result="unknown request: {}".format(req))
@app.route('/_info')
def info():
from . import alignment
readid = unquote(request.args.get('readid', 0))
alnSet = dataHub.getAlignmentSetByName(readid)
if alnSet:
if dataHub.args.verbose > 3:
print("\n")
print(alnSet.parentCollection.name)
for setName, moreAlnSet in alnSet.parentCollection.sets.items():
desc = [setName]
for aln in moreAlnSet.getAlignments():
desc.append("{}:{}-{}{} ({}[{}])".format(aln.regionID, aln.start, aln.end, aln.strand, aln.score, aln.score2))
print("\t".join(desc))
reads = alnSet.getAlignments()
result = []
for read in reads:
html = "{}<br/>".format(alignment.getBlastRepresentation(read).replace("\n", "<br/>"))
html = html.replace(" ", ".")
result.append(html)
result.append("<br/>Regions={}".format(",".join(map(str, [aln.regionID for aln in alnSet.getAlignments()]))))
result.append("<br/>Total length={}".format(len(alnSet)))
result.append(" Reason={}".format(alnSet.parentCollection.why))
if len(alnSet.getAlignments()) > 1:
result += " Alignment Scores: ref={} alt={}".format(
alnSet.parentCollection.sets["ref"].evidences["alignmentScore"],
alnSet.parentCollection.sets["alt"].evidences["alignmentScore"])
result += " Insert Size Scores: ref={} alt={}".format(
alnSet.parentCollection.sets["ref"].evidences["insertSizeScore"],
alnSet.parentCollection.sets["alt"].evidences["insertSizeScore"])
result += " Lengths: ref={} alt={}".format(
len(alnSet.parentCollection.sets["ref"]),
len(alnSet.parentCollection.sets["alt"]))
result += " mapq={}".format(",".join(str(read.mapq) for read in reads))
# result.append(" Log odds={:.3g}".format(float(READ_INFO[readid].prob)))
result = "".join(result)
result = "<div style='font-family:Courier;'>" + result + "</div>"
result = jsonify(result=result)
return result
else:
logging.debug("NOT FOUND:{}".format(readid))
@app.route('/_isizes/<name>')
def displayIsizes(name):
if not dataHub.samples[name].insertSizePlot:
return None
return Response(dataHub.samples[name].insertSizePlot, mimetype="image/png")
@app.route("/_dotplots/<name>")
def get_dotplot(name):
if name in dataHub.dotplots:
return Response(dataHub.dotplots[name], mimetype="image/png")
return None
def run(port=None):
import webbrowser, threading
if port is None:
port = getRandomPort()
# load()
url = "http://127.0.0.1:{}/".format(port)
logging.info("Starting browser at {}".format(url))
# webbrowser.open_new(url)
threading.Timer(1.25, lambda: webbrowser.open(url) ).start()
app.run(
port=port#,
# debug=True
)
if __name__ == '__main__':
pass
