[genoPlotR-help] genoplotR help with graphic presentation or CDS and names
Artemis H
hatziiod at gmail.com
Wed Jan 25 19:18:29 CET 2017
Hello,
I've just tried to use genoplotR for the first time. I used
read_dna_seg_from_genbank to import CDS info of cluster genes and
read_comparison_from_blast to read blastn comparison files. I checked them
all with is.dna_seg and is.comparison and they all came out TRUE. Finally I
used plot_gene_map(dna_segs=*all my seg files* ), comparisons=*all my
comparison files*),override_color_schemes=TRUE,
global_color_scheme=c("e_value", "auto", "grey", "0")).
With this process I end up with a nice image which seems to show the blastn
conserved regions but instead of CDS boxes I get thin blue dispersed lines
along my stick cluster and I can't seem to find how to give each cluster a
name. I've attached the usual output.
I would like to ask first, how can I label the clusters, secondly how can I
make the CDS show for each cluster as a box or arrow, thirdly how does one
go about generating a newick2phylog file with the raw data?
I'm pasting an example of a dna_seg object in case that helps:
> geneA.seq
name start end strand length pid gene synonym product proteinid
feature
1 geneA 828 1001 1 57 NA geneA NA geneA CD352.1
CDS
2 geneB 1109 4090 1 993 NA geneB NA geneB CD353.1
CDS
3 geneT 4101 5903 1 600 NA geneT NA geneT CD354.1
CDS
4 geneC 5896 7140 1 414 NA geneC NA geneC CD355.1
CDS
5 geneI 7137 7874 1 245 NA geneI NA geneI CD356.1
CDS
6 geneP 7876 9924 1 682 NA geneP NA geneP CD357.1
CDS
7 geneR 9993 10679 1 228 NA geneR NA geneR CD358.1
CDS
8 geneK 10672 12015 1 447 NA geneK NA geneK CD359.1
CDS
9 geneF 12114 12791 1 225 NA geneF NA geneF CD360.1
CDS
10 geneE 12793 13521 1 242 NA geneE NA geneE CD361.1
CDS
11 geneG 13508 14152 1 214 NA geneG NA geneG CD362.1
CDS
gene_type col lty lwd pch cex
1 bars blue 1 1 8 1
2 bars blue 1 1 8 1
3 bars blue 1 1 8 1
4 bars blue 1 1 8 1
5 bars blue 1 1 8 1
6 bars blue 1 1 8 1
7 bars blue 1 1 8 1
8 bars blue 1 1 8 1
9 bars blue 1 1 8 1
10 bars blue 1 1 8 1
11 bars blue 1 1 8 1
Thank you in advance, any help and tips would be appreciated.
Diane
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.r-forge.r-project.org/pipermail/genoplotr-help/attachments/20170125/49e60a04/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Rplots_v1.pdf
Type: application/pdf
Size: 7909 bytes
Desc: not available
URL: <http://lists.r-forge.r-project.org/pipermail/genoplotr-help/attachments/20170125/49e60a04/attachment.pdf>
More information about the genoPlotR-help
mailing list