13 votes
Accepted

Using the t-SNE algorithm on microarray data + an error bonus

Converting your data.frame to a matrix (and then removing the data.frame) will often free up ...
Devon Ryan's user avatar
  • 19.5k
13 votes
Accepted

Volcano plot in R

This plot is clearly done using core R functions. There are smoother alternatives how to make a pretty volcano plot (like ggplot with example here), but if you ...
Kamil S Jaron's user avatar
6 votes

Interpreting Intergrative Genomic Viewer (IGV)

You've apparently colored your alignments by read strand. In this case, red indicates "+ (watson) strand" and blue indicates "- (strand) strand". This strand association is determined by the ...
Devon Ryan's user avatar
  • 19.5k
6 votes
Accepted

Visualization tools for 3C/Hi-C long-range interaction data?

I would also recommend two very recent Hi-C visualization frameworks (with some public data available in both): HiGlass and JuiceBox.
golobor's user avatar
  • 176
6 votes
Accepted

full visualisation of draft genomes alignment

Found a solution, using D-Genies, worked great. Some examples from their website: Thanks to @user172818.
aechchiki's user avatar
  • 2,656
6 votes
Accepted

What is the name of this type of figure?

This kind of visualization showing the banding patterns of chromosomes are called "chromosome ideograms". You can use for example IdeoViz in R to generate them.
PPK's user avatar
  • 876
5 votes

How to enlarge a section of y-axis?

I would strongly discourage you from making discontinuous axis, it's going to be very confusing for a reader. The facet plot you proposed seems like a good solution to me. Alternatively you can use ...
Kamil S Jaron's user avatar
5 votes
Accepted

3D PCA group labelling

You can easily color 3D pca plots in R based on the code given below: ...
Ammar Sabir Cheema's user avatar
5 votes

Visualisation of long read RNA-Seq splicing

DNASTAR's software is for purchase, but high quality. GenVision Pro does genomic visualization, including Sashimi plots. Edit: not sure why this answer is being downvoted, unless it's because the ...
Rosalind Was Robbed's user avatar
5 votes

A good tool for gene locus visualization

I think in R, Gviz is the best choice. I have used it for human sequences but it will probably also be suitable for Prokaryotes. Another option is MS powerpoint, if you only need the arrows...
benn's user avatar
  • 3,571
5 votes
Accepted

What do the read colors in IGV mean?

Could it be that these are reads with their mates found on other chromosomes? That's the part of the IGV documentation which seems to make sense in this context: for paired end reads that are coded ...
Ram RS's user avatar
  • 1,993
5 votes

Need a alternative or more complex version of venn diagram in python for matching dna sequences

To get around limitations in using Venn diagrams for set overlaps, when there are more than three sets, back around 2013 or so I created something I called an Eulergrid plot (example at the bottom of ...
Alex Reynolds's user avatar
5 votes
Accepted

Create identicon of a DNA sequence

Sequenticon is a Python library for generating identicons for DNA sequences. For instance, the sequence ATGGTGCA gets converted into the following icon: A web ...
Peter's user avatar
  • 2,624
4 votes
Accepted

Highlight CpG sites with Texshade

The best solution I could find was to use the \shaderegion function on all the sequences so it ended up looking something like this: ...
quantik's user avatar
  • 255
4 votes
Accepted

creating graph of distribution of blast hits on the query sequence

You probably want to include query start (qstart) and query end (qend) in your blast output. Something like this: ...
benn's user avatar
  • 3,571
4 votes
Accepted

Customizing bigWig file

There's no equivalent to the wiggle header in bigWig (or bigBed) files, which is why UCSC uses the file name. This is actually the reason for the track line stuff that you linked to, since you can ...
Devon Ryan's user avatar
  • 19.5k
4 votes
Accepted

Extracting genes from corrplot and adding labels based on high and low corelation

Here is a solution using the pheatmap library to cluster and visualise the correlation matrix, then extract the groups from the cluster dendrograms: ...
user3479780's user avatar
4 votes

Visualize similar interactions in two homologous networks

An easy way to visualize this is to make a tab-delimited table that contains edge information for both species and both networks, color the edges depending on the species, and make the thickness and ...
conchoecia's user avatar
  • 3,141
4 votes
Accepted

Occupancy of TFs with the target genes

I would ignore peak calling for this and instead compute enrichment of ChIP/input for the genome (e.g., with deepTools or presumably homer) and then plot it for the genes of interest individually (e.g....
Devon Ryan's user avatar
  • 19.5k
4 votes
Accepted

How can I export a full alignment from IGV as an image?

Running an IGV batch script helps increase the number of reads per image, compared with File | Save image, but it does not solve the issue by itself. To run a ...
Timur Shtatland's user avatar
3 votes

Seurat Merged objects tSNE - How to paint on original IDs?

Below a few lines of code that accompany BC Wang's answer. After using MergeSeurat the sample name will be added to meta data under orig.ident. this can then be used to color the tSNE either using ...
Mack123456's user avatar
3 votes

Volcano plot in R

Without a minimal reproducible example (MWE), I can't reproduce the plot but I would suggest using existing volcanoplot functions such as the ...
Tom Kelly's user avatar
  • 873
3 votes
Accepted

Visualisation of long read RNA-Seq splicing

IGB gives the intended result without much hassle. It doesn't have the bells and whistles of IGV, but presents a clean and intuitive view of the individual reads. In IGB, each read is presented as a ...
Scott Gigante's user avatar
3 votes

Visualisation of long read RNA-Seq splicing

The Sashmi Plot feature built into IGV. It gives a nice summary of the spliced transcripts and the coverage of each exon. For example:
GWW's user avatar
  • 752
3 votes

Visualization tools for 3C/Hi-C long-range interaction data?

If you have more than a dozen or even hundreds of regions that you want to compare you could check a companion tool of HiGlass called HiPiler. HiPiler lets you arrange, cluster, and pile up to ...
F Lekschas's user avatar
3 votes

A good tool for gene locus visualization

gggenes is a ggplot2 based solution and looks very appealing. Plus, it has a nice introduction with examples. ...
Alex's user avatar
  • 41
3 votes
Accepted

A good tool for gene locus visualization

I ended up using Easyfig, which is made for the task and requires little tinkering around, once fed the genbank files. (Granted, being so simple, it's hard to make the image more advanced, but it ...
Mgall's user avatar
  • 101
3 votes

A good tool for gene locus visualization

GFF2PS was developed for exactly this sort of thing. It takes a GFF file as input, and will produce a pretty PS output showing the annotations. It was used to create the poster for Celera's human ...
terdon's user avatar
  • 9,442
3 votes
Accepted

How to plot p-values in a circular barplot?

The adjustcolor function from grDevices can be used to generate colours with transparency. What you want is colours based on the ...
Tom Kelly's user avatar
  • 873
3 votes
Accepted

Plotting coverage of annotation over collection of region

If one assumes that repeat elements of a given type (e.g., LINEs) don't overlap each other, then the following will work: Split your BED file by repeat element, such that you have a LINE.bed, SINE....
Devon Ryan's user avatar
  • 19.5k

Only top scored, non community-wiki answers of a minimum length are eligible