mappingQuality.ssp 4.03 KB
Newer Older
Peter van 't Hof's avatar
Peter van 't Hof committed
1
#import(nl.lumc.sasc.biopet.utils.summary.db.SummaryDb)
2
#import(nl.lumc.sasc.biopet.core.report.ReportPage)
Peter van 't Hof's avatar
Peter van 't Hof committed
3
4
#import(scala.concurrent.Await)
#import(scala.concurrent.duration.Duration)
5
6
#import(nl.lumc.sasc.biopet.pipelines.bammetrics.BammetricsReport)
#import(java.io.File)
Peter van 't Hof's avatar
Peter van 't Hof committed
7
8
#import(scalaz._)
#import(Scalaz._)
Peter van 't Hof's avatar
Peter van 't Hof committed
9
10
11
<%@ var summary: SummaryDb %>
<%@ var sampleId: Option[Int] = None %>
<%@ var libId: Option[Int] = None %>
12
<%@ var sampleLevel: Boolean = false %>
Peter van 't Hof's avatar
Peter van 't Hof committed
13
<%@ var rootPath: String %>
14
15
16
17
<%@ var outputDir: File %>
<%@ var showPlot: Boolean = false %>
<%@ var showTable: Boolean = true %>
<%@ var showIntro: Boolean = true%>
Peter van 't Hof's avatar
Peter van 't Hof committed
18
19
20
<%@ var runId: Int %>
<%@ var fields: List[String] = List("min", "max", "mean", "median", "modal")%>

21
#{
Peter van 't Hof's avatar
Peter van 't Hof committed
22
    val samples = Await.result(summary.getSamples(runId = Some(runId)), Duration.Inf)
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
}#

#if (showIntro)
    <br/>
    <div class="row">
    <div class="col-md-1"></div>
    <div class="col-md-6">
        <p>
        #if (sampleId.isDefined && libId.isDefined)
            This plot shows the mapping quality distribution for all libraries combined in sample <b>${sampleId}</b>.
        #elseif(sampleId.isDefined)
            This plot shows the mapping quality distribution for the libraries of sample <b>${sampleId}</b>.
        #else
            This plot shows the mapping quality distribution for each of the <b>${samples.size}</b> samples.
        #end
        </p>
    </div>
    </div>
#end

#if (showPlot)
Peter van 't Hof's avatar
Peter van 't Hof committed
44
    #{ BammetricsReport.mappingQualityPlot(outputDir, "mapping_quality", summary, !sampleLevel, sampleId = sampleId, libraryId = libId) }#
45
46
47
48
49
50
51
52
53
54
55
56

    <div class="panel-body">
    <img src="mapping_quality.png" class="img-responsive" />
    </div>
    <div class="panel-footer">
        #if (showTable)
            <button type="button" class="btn btn-info" data-toggle="collapse" data-target="#mapping_qualityTable">
                <i class="glyphicon glyphicon-eye-close"></i> Hide table</button>
        #else
            <button type="button" class="btn btn-info" data-toggle="collapse" data-target="#mapping_qualityTable">
            <i class="glyphicon glyphicon-eye-open"></i> Show table</button>
        #end
Peter van 't Hof's avatar
Peter van 't Hof committed
57
    <a href="mapping_quality.tsv"><button type="button" class="btn btn-info"><i class="glyphicon glyphicon-cloud-download"></i> TSV file</button></a>
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72

    </div>
#end

<div class="panel-body collapse #if (showTable)in#end" id="mapping_qualityTable">
<!-- Table -->
<table class="table">
    <thead><tr>
        <th data-sorted="true" data-sorted-direction="ascending">Sample</th>
        #if (!sampleLevel) <th>Library</th> #end
        #for (field <- fields)
            <th>${field.replaceAll("_", " ")}</th>
        #end
    </tr></thead>
    <tbody>
Peter van 't Hof's avatar
Peter van 't Hof committed
73
        #for (sample <- samples.sortBy(_.name))
74
            #{
Peter van 't Hof's avatar
Peter van 't Hof committed
75
76
77
78
                val libs: List[Option[Int]] = (libId, sampleLevel) match {
                case (_, true) => List(None)
                case (Some(_), _) => List(libId)
                case _ => Await.result(summary.getLibraries(sampleId = Some(sample.id), runId = Some(runId)), Duration.Inf).map(x => Some(x.id)).toList
79
80
                }
            }#
81
            <tr><td rowspan="${libs.size}"><a href="${rootPath}Samples/${sample.name}/index.html">${sample.name}</a></td>
82
            #for (libId <- libs)
83
                #{ val libName = libId.map(l => Await.result(summary.getLibraryName(l), Duration.Inf)) }#
84
                #if (libs.head != libId) <tr> #end
85
                    #if (!sampleLevel) <td><a href="${rootPath}Samples/${sample.name}/Libraries/${libName}/index.html">${libName}</a></td> #end
86
                    #{
Peter van 't Hof's avatar
Peter van 't Hof committed
87
                        val statsPaths = fields.map(x => x -> List("mapping_quality", "general", x)).toMap
Peter van 't Hof's avatar
Peter van 't Hof committed
88
                        val results = summary.getStatKeys(runId, "bammetrics".right, Some("bamstats".right), sample = Some(sample.id.left), library = libId.map(_.left), keyValues = statsPaths)
89
                    }#
Peter van 't Hof's avatar
Peter van 't Hof committed
90
91
                    #for (field <- fields)
                        <td>${results(field)}</td>
92
93
94
95
96
97
98
99
                    #end
                </tr>
            #end
        #end
    </tbody>
</table>

</div>