Mercurial > repos > drosofff > msp_sr_readmap_and_size_histograms
annotate readmap.xml @ 8:4e7f59d5ee18 draft
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
| author | drosofff |
|---|---|
| date | Tue, 02 Feb 2016 11:50:13 -0500 |
| parents | 9784e5046f04 |
| children | 44a0b0fa52a9 |
| rev | line source |
|---|---|
|
7
9784e5046f04
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit 450d1b3d954cc7d96e53652d468811546951f995
drosofff
parents:
6
diff
changeset
|
1 <tool id="Readmap" name="Generate readmap and histograms from alignment files" version="1.0.5"> |
| 0 | 2 <description>from sRbowtie aligment</description> |
| 3 <requirements> | |
| 4 <requirement type="package" version="0.12.7">bowtie</requirement> | |
| 5 <requirement type="package" version="0.7.7">pysam</requirement> | |
|
6
0f7cf7fb3465
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit fe40dec87779c1fcfbd03330e653aa886f4a2cda
drosofff
parents:
5
diff
changeset
|
6 <requirement type="package" version="3.1.2">R</requirement> |
| 0 | 7 <requirement type="package" version="2.14">biocbasics</requirement> |
| 1 | 8 <requirement type="package" version="1.9">numpy</requirement> |
| 0 | 9 </requirements> |
| 10 <command interpreter="python"> | |
| 11 readmap.py | |
| 12 #if $refGenomeSource.genomeSource == "history": | |
| 13 --reference_fasta ## sys.argv[2] | |
| 14 $refGenomeSource.ownFile ## index source | |
| 15 #else: | |
| 16 #silent reference= filter( lambda x: str( x[0] ) == str( $refGenomeSource.series[0].input.dbkey ), $__app__.tool_data_tables[ 'bowtie_indexes' ].get_fields() )[0][-1] | |
| 17 --reference_bowtie_index | |
| 18 $reference | |
| 19 #end if | |
| 20 --rcode | |
| 21 $plotCode | |
| 22 --output_readmap | |
| 23 $readmap_dataframe | |
| 24 --output_size_distribution | |
| 25 $size_distribution_dataframe | |
| 26 --minquery | |
| 27 $minquery | |
| 28 --maxquery | |
| 29 $maxquery | |
| 30 --input | |
| 31 #for $i in $refGenomeSource.series | |
| 32 $i.input | |
| 33 #end for | |
| 34 --ext | |
| 35 #for $i in $refGenomeSource.series | |
| 36 $i.input.ext | |
| 37 #end for | |
| 38 --label | |
| 39 #for $i in $refGenomeSource.series | |
| 40 "$i.input.name" | |
| 41 #end for | |
| 42 --normalization_factor | |
| 43 #for $i in $refGenomeSource.series | |
| 44 $i.norm | |
| 45 #end for | |
| 46 #if $gff: | |
| 47 --gff | |
| 48 $gff | |
| 49 #end if | |
| 50 | |
| 51 </command> | |
| 52 <inputs> | |
| 53 <conditional name="refGenomeSource"> | |
| 54 <param name="genomeSource" type="select" label="Will you select a reference genome from your history or use a built-in index?" help="Built-ins were indexed using default options"> | |
| 55 <option value="indexed">Use a built-in index</option> | |
| 56 <option value="history">Use one from the history</option> | |
| 57 </param> | |
| 58 <when value="indexed"> | |
| 59 <repeat name="series" title="Add alignment files"> | |
|
4
dffa22efc6a8
planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
3
diff
changeset
|
60 <param name="input" type="data" label="Select multiple alignments to parse" format="tabular,sam,bam"> |
| 0 | 61 <validator type="dataset_metadata_in_data_table" table_name="bowtie_indexes" metadata_name="dbkey" metadata_column="0" message="database not set for this bowtie output. Select the database(=genome used for matching) manually, or select a reference fasta from your history."/> |
| 62 </param> | |
| 63 <param name="norm" type="float" value="1" label="Indicate a normalization factor to compare multiple aligments"/> | |
| 64 </repeat> | |
| 65 </when> | |
| 66 <when value="history"> | |
| 67 <param name="ownFile" type="data" format="fasta" label="Select a fasta file, that served as the reference index for the alignments" /> | |
| 68 <repeat name="series" title="Add alignment files"> | |
|
4
dffa22efc6a8
planemo upload for repository https://bitbucket.org/drosofff/gedtools/
drosofff
parents:
3
diff
changeset
|
69 <param name="input" type="data" label="Select multiple alignments to parse" format="tabular,sam,bam"/> |
| 0 | 70 <param name="norm" type="float" value="1" label="Indicate a normalization factor to compare multiple aligments"/> |
| 71 </repeat> | |
| 72 </when> | |
| 73 </conditional> | |
| 74 <param name="gff" type="data" format="gff3" optional="true" label="Optional: select a GFF to investigate regions of interest" help="GFF must match genome build"/> | |
| 75 <!-- <validator type="dataset_metadata_in_data_table" table_name="bowtie_indexes" metadata_name="dbkey" metadata_column="0" message="GFF database and alignment file databse do not match!"/> --> | |
| 76 <param name="minquery" type="integer" size="3" value="18" label="Min size of query small RNAs" help="'18' = 18 nucleotides"/> | |
| 77 <param name="maxquery" type="integer" size="3" value="28" label="Max size of query small RNAs" help="'28' = 28 nucleotides"/> | |
| 78 <param name="title" type="text" size="15" value= "Readmaps and size distributions" label="Main Titles"/> | |
| 79 <param name="xlabel" type="text" size="15" value="Coordinates/read size" label="x axis label"/> | |
|
8
4e7f59d5ee18
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
drosofff
parents:
7
diff
changeset
|
80 <param name="yrange" type="integer" size="6" value="0" label="y axis range" help="leave at 0 for autoscaling"/> |
| 0 | 81 <param name="ylabel" type="text" size="15" value="Number of reads" label="y axis label"/> |
| 82 <param name="rows_per_page" type="text" size="9" value="8" label="How many items to display per page?"> | |
| 83 <validator type="in_range" min="6" max="20" message="Select between 6 and 20 rows, as the readability will suffer otherwise."/> | |
| 84 </param> | |
| 85 </inputs> | |
| 86 <configfiles> | |
| 87 <configfile name="plotCode"> | |
| 88 ## Setup R error handling to go to stderr | |
| 89 options( show.error.messages=F, | |
| 90 error = function () { cat( geterrmessage(), file=stderr() ); q( "no", 1, F ) } ) | |
| 91 library(RColorBrewer) | |
| 92 library(lattice) | |
| 93 library(latticeExtra) | |
| 94 library(grid) | |
| 95 library(gridExtra) | |
| 96 | |
| 97 ## data frames implementation | |
| 98 | |
| 99 rm=read.delim("${readmap_dataframe}", header=T, row.names=NULL) | |
| 100 n_samples=length(unique(rm\$sample)) | |
| 101 genes=unique(levels(rm\$gene)) | |
| 102 per_gene_readmap=lapply(genes, function(x) subset(rm, gene==x)) ####### ? | |
| 103 n_genes=length(per_gene_readmap) | |
| 104 | |
| 105 size=read.delim("${size_distribution_dataframe}", header=T, row.names=NULL) | |
| 106 per_gene_size=lapply(genes, function(x) subset(size, gene==x)) ###### ? | |
| 107 | |
| 108 ## end of data frames implementation | |
| 109 | |
| 110 ## functions | |
| 111 | |
| 112 plot_readmap=function(df, ...) { | |
| 113 combineLimits(xyplot(count~coord|factor(sample, levels=unique(sample))+reorder(gene, count, function(x) -sum(abs(x))), | |
| 114 data=df, | |
| 115 type='h', | |
| 116 scales= list(relation="free", x=list(rot=0, cex=0.7, axs="i", tck=0.5), y=list(tick.number=4, rot=90, cex=0.7)), | |
| 117 xlab=NULL, main=NULL, ylab=NULL, | |
| 118 as.table=T, | |
| 119 origin = 0, | |
| 120 horizontal=FALSE, | |
| 121 group=polarity, | |
| 122 col=c("red","blue"), | |
| 123 par.strip.text = list(cex=0.7), | |
| 124 ...)) | |
| 125 } | |
| 126 | |
| 127 plot_size_distribution= function(df, ...) { | |
| 128 smR.prepanel=function(x,y,...){; yscale=c(-max(abs(y)), max(abs(y)));list(ylim=yscale);} | |
| 129 bc= barchart(count~as.factor(size)|factor(sample, levels=unique(sample))+gene, data = df, origin = 0, | |
| 130 horizontal=FALSE, | |
|
8
4e7f59d5ee18
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
drosofff
parents:
7
diff
changeset
|
131 group=polarity, |
|
4e7f59d5ee18
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit e27d18d58ae095e7fad4b08b04370857a1d37964-dirty
drosofff
parents:
7
diff
changeset
|
132 stack=TRUE, |
| 0 | 133 col=c('red', 'blue'), |
| 134 cex=0.75, | |
| 135 scales=list(y=list(tick.number=4, rot=90, relation="free", cex=0.7), x=list(cex=0.7) ), | |
| 136 prepanel=smR.prepanel, | |
| 137 xlab = NULL, | |
| 138 ylab = NULL, | |
| 139 main = NULL, | |
| 140 as.table=TRUE, | |
| 141 newpage = T, | |
| 142 par.strip.text = list(cex=0.7), | |
| 143 ...) | |
| 144 combineLimits(bc) | |
| 145 } | |
| 146 | |
| 147 ## end of functions | |
| 148 | |
| 149 ## function parameters' | |
| 150 | |
| 151 par.settings.readmap=list(layout.heights=list(top.padding=0, bottom.padding=-2.5), strip.background = list(col=c("lightblue","lightgreen")) ) | |
| 152 par.settings.size=list(layout.heights=list(top.padding=-1, bottom.padding=-2.5), strip.background = list(col=c("lightblue","lightgreen")) ) | |
| 153 par.settings.combination.readmap=list(layout.heights=list(top.padding=0, bottom.padding=-3), strip.background=list(col=c("lightblue","lightgreen")) ) | |
| 154 par.settings.combination.size=list(layout.heights=list(top.padding=-2, bottom.padding=-0.5), strip.background=list(col=c("lightblue", "lightgreen")) ) | |
| 155 | |
| 156 ## end of function parameters' | |
| 157 | |
| 158 ## GRAPHS | |
| 159 | |
| 160 if (n_genes > 7) {page_height_simple = 11.69; page_height_combi=11.69; rows_per_page=${rows_per_page}; extrarow=0 } else { | |
|
7
9784e5046f04
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit 450d1b3d954cc7d96e53652d468811546951f995
drosofff
parents:
6
diff
changeset
|
161 rows_per_page= 8; page_height_simple = 11.69; page_height_combi=11.69; extrarow=0 } |
|
9784e5046f04
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit 450d1b3d954cc7d96e53652d468811546951f995
drosofff
parents:
6
diff
changeset
|
162 ## rows_per_page= 8; page_height_simple = 11.69/7*n_genes; page_height_combi=11.69/9*(n_genes*2); extrarow=0 } |
|
9784e5046f04
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit 450d1b3d954cc7d96e53652d468811546951f995
drosofff
parents:
6
diff
changeset
|
163 ## rows_per_page= n_genes; page_height_simple = 11.69/n_genes/4; page_height_combi=11.69/(n_genes*2); extrarow=1 } |
|
9784e5046f04
planemo upload for repository https://github.com/ARTbio/tools-artbio/tree/master/tools/msp_sr_readmap_and_size_histograms commit 450d1b3d954cc7d96e53652d468811546951f995
drosofff
parents:
6
diff
changeset
|
164 if (n_samples > 4) {page_width = 8.2677*n_samples/4} else {page_width = 8.2677*n_samples/3} # to test |
| 0 | 165 |
| 166 pdf(file="${readmap_PDF}", paper="special", height=page_height_simple, width=page_width) | |
| 167 for (i in seq(1,n_genes,rows_per_page)) { | |
| 168 start=i | |
| 169 end=i+rows_per_page-1 | |
| 170 if (end>n_genes) {end=n_genes} | |
| 171 readmap_plot.list=lapply(per_gene_readmap[start:end], function(x) plot_readmap(x, par.settings=par.settings.readmap)) | |
| 172 args.list=c(readmap_plot.list, list(nrow=rows_per_page, ncol=1, | |
| 173 main=textGrob("Read Maps (nucleotide coordinates)", gp=gpar(cex=1), just="top"), | |
| 174 left=textGrob("${ylabel}", gp=gpar(cex=1), vjust=1, rot=90) | |
| 175 #sub=textGrob("readmap coordinates", gp=gpar(cex=.75), just="bottom") | |
| 176 ) | |
| 177 ) | |
| 178 do.call(grid.arrange, args.list) | |
| 179 } | |
| 180 devname=dev.off() | |
| 181 | |
| 182 | |
| 183 pdf(file="${size_PDF}", paper="special", height=page_height_simple, width=page_width) | |
| 184 for (i in seq(1,n_genes,rows_per_page)) { | |
| 185 start=i | |
| 186 end=i+rows_per_page-1 | |
| 187 if (end>n_genes) {end=n_genes} | |
| 188 plot.list=lapply(per_gene_size[start:end], function(x) plot_size_distribution(x, par.settings=par.settings.size) ) | |
| 189 args.list=c(plot.list, list(nrow=rows_per_page, ncol=1, | |
| 190 main=textGrob("Size distributions (in nucleotides)", gp=gpar(cex=1), just="top"), | |
| 191 left=textGrob("${ylabel}", gp=gpar(cex=1), vjust=1, rot=90) | |
| 192 #sub="readsize in nucleotides" | |
| 193 ) | |
| 194 ) | |
| 195 do.call(grid.arrange, args.list) | |
| 196 } | |
| 197 devname=dev.off() | |
| 198 | |
| 199 pdf(file="${combi_PDF}", paper="special", height=page_height_combi, width=page_width) | |
| 200 for (i in seq(1,n_genes,rows_per_page/2)) { | |
| 201 start=i | |
| 202 end=i+rows_per_page/2-1 | |
| 203 if (end>n_genes) {end=n_genes} | |
| 204 read_plot.list=lapply(per_gene_readmap[start:end], function(x) plot_readmap(x, par.settings=par.settings.combination.readmap)) | |
| 205 size_plot.list=lapply(per_gene_size[start:end], function(x) plot_size_distribution(x, strip=FALSE, par.settings=par.settings.combination.size)) | |
| 206 plot.list=rbind(read_plot.list, size_plot.list ) | |
| 207 args.list=c(plot.list, list(nrow=rows_per_page + extrarow, ncol=1, | |
| 208 main=textGrob("${title}", gp=gpar(cex=1), just="top"), | |
| 209 left=textGrob("${ylabel}", gp=gpar(cex=1), vjust=1, rot=90), | |
| 210 sub=textGrob("${xlabel}", gp=gpar(cex=1), just="bottom") | |
| 211 ) | |
| 212 ) | |
| 213 do.call(grid.arrange, args.list) | |
| 214 } | |
| 215 devname=dev.off() | |
| 216 | |
| 217 | |
| 218 </configfile> | |
| 219 </configfiles> | |
| 220 | |
| 221 <outputs> | |
| 222 <data format="tabular" name="readmap_dataframe" label="Readmap dataframe"/> | |
| 223 <data format="tabular" name="size_distribution_dataframe" label="Size distribution dataframe"/> | |
| 224 <data format="pdf" name="readmap_PDF" label="Readmaps"/> | |
| 225 <data format="pdf" name="size_PDF" label="Size distribution"/> | |
| 226 <data format="pdf" name="combi_PDF" label="Size distribution and Readmaps"/> | |
| 227 </outputs> | |
| 228 <help> | |
| 229 | |
| 230 **What it does** | |
| 231 | |
| 232 Takes one or more alignment files (BAM, SAM or tabular bowtie output) as input and produces a "Readmap", | |
| 233 where by default for each "chromosome" the position of the read is recorded on the x-axis, and the y-axis indicates | |
| 234 the number of reads per position. Reads that map in sense are on the top, reads that map antisense are on the bottom. | |
| 235 | |
| 236 | |
| 237 .. class:: warningmark | |
| 238 | |
| 239 '''TIP''' The input data can be produced using the sRbowtie tool. | |
| 240 | |
| 241 ---- | |
| 242 | |
| 243 '''Example''' | |
| 244 | |
| 245 Query sequence:: | |
| 246 For a SAM file as the following: | |
| 247 | |
| 248 5 16 2L_79 24393 255 17M * 0 0 CCTTCATCTTTTTTTTT IIIIIIIIIIIIIIIII XA:i:0 MD:Z:17 NM:i:0 | |
| 249 | |
| 250 11 0 2R_1 12675 255 21M * 0 0 AAAAAAAACGCGTCCTTGTGC IIIIIIIIIIIIIIIIIIIII XA:i:0 MD:Z:21 NM:i:0 | |
| 251 | |
| 252 2 16 2L_5 669 255 23M * 0 0 TGTTGCTGCATTTCTTTTTTTTT IIIIIIIIIIIIIIIIIIIIIII XA:i:0 MD:Z:23 NM:i:0 | |
| 253 | |
| 254 produce a plot like this: | |
| 255 | |
| 256 ---- | |
| 257 | |
| 258 .. image:: static/images/readmap.png | |
| 259 :height: 800 | |
| 260 :width: 500 | |
| 261 | |
| 262 </help> | |
| 1 | 263 <tests> |
| 0 | 264 <test> |
| 265 <param name="genomeSource" value="history" /> | |
| 266 <param name="ownFile" value ="transposons.fasta" ftype="fasta" /> | |
| 1 | 267 <param name="series_0|input" value="sample1.srbowtie_out" ftype="tabular"/> |
| 268 <param name="series_0|norm" value="1" /> | |
| 269 <param name="series_1|input" value="sample2.srbowtie_out" ftype="tabular"/> | |
| 270 <param name="series_1|norm" value="1" /> | |
| 271 <param name="series_2|input" value="sample3.srbowtie_out" ftype="tabular"/> | |
| 272 <param name="series_2|norm" value="1" /> | |
| 0 | 273 <param name="minquery" value="20" /> |
| 274 <param name="maxquery" value="30" /> | |
| 275 <param name="title" value="Readmaps and size distributions" /> | |
| 276 <param name="xlabel" value="Coordinates/read size" /> | |
| 277 <param name="ylabel" value="Number of reads" /> | |
| 278 <param name="rows_per_page" value="8" /> | |
| 1 | 279 <output name="readmap_dataframe" ftype="tabular" file="Readmap_dataframe.tab" /> |
| 280 <output name="size_distribution_dataframe" ftype="tabular" file="Size_distribution_dataframe.tab" /> | |
| 281 <output name="readmap_PDF" ftype="pdf" file="Readmaps.pdf" /> | |
| 282 <output name="size_PDF" ftype="pdf" file="Size_distribution.pdf" /> | |
| 283 <output name="combi_PDF" ftype="pdf" file="Size_distribution_and_Readmaps.pdf" /> | |
| 0 | 284 </test> |
| 1 | 285 </tests> |
| 0 | 286 </tool> |