[genoPlotR-help] genoPlotR question
Lionel Guy
lionel.guy at icm.uu.se
Mon Aug 27 10:05:11 CEST 2012
Hi Wei,
Thanks for your interest in genoPlotR.
It's perfectly possible to plot a segment in reverse orientation with genoPlotR. However, the DNA sequence itself is not read by genoPlotR and thus won't be reverse complemented.
The argument to use is xlims, putting first the 3'-end of the segment. From the help page:
'xlims' allow the user to plot subsegments of a 'dna_seg'. 'xlims'
consists of a list composed of as many numeric vectors as there
are segments. Each of these numeric vectors give pairs of left and
right borders, and gives the direction. For example, c(1,2,6,4)
will plot two subsegments, segment 1 to 2 which is plotted left to
right and segment 4 to 6, plotted right to left. '-Inf' and 'Inf'
values are accepted. 'NULL' values will result in plotting the
whole segment.
From the examples:
data("three_genes")
## xlims
## Just returning a segment
plot_gene_map(dna_segs, comparisons,
xlims=list(NULL, NULL, c(Inf,-Inf)),
dna_seg_scale=TRUE)
Note the third element of the xlims list: c(Inf,-Inf) which will plot the whole third segment reversed.
HTH,
Lionel
On 27 Aug 2012, at 4:25 , Wei Lu wrote:
> Dear Lionel,
>
> Thanks for your great work on this R package. We are now using this tool to help our comparative genomic analysis. I've got a quick question to ask you. Can I do a reverse and complement of dan_seg in genoPlotR? For example, I read a genbank file and plot a 5kb segment of the genome and I would like to reverse and complement this region for the comparison. Can we do this in genoPlotR?
>
> Thanks,
>
> Dr. Wei LU
> Research Associate
> University of Western Australia | The Marshall Centre for Infectious
> Diseases Research and Training
> Reception Room 1.2, M Block, QEII Medical Centre | Monash Avenue |
> Nedlands 6009 | Western Australia, Australia.
> T: 08-92876933 | F: 08-92876213 | E: wei.lu at uwa.edu.au
>
--
Lionel Guy
Molecular Evolution, Uppsala University, Uppsala, Sweden
postal address: Box 596, SE-751 24 Uppsala; visiting address: BMC C6:305d, Husargatan 3, SE-752 37 Uppsala
phone: +46 18 471 6129, mobile +46 73 976 0618
lionel.guy at icm.uu.se
More information about the genoPlotR-help
mailing list