Enhort - Genomic Position Profiling

Info

Current File: HIV

Positions395

Background

Data

Change tracks Settings

Upload background sites

Results

Download results

Annotation group ( 7 / 23 )

Max effect size

Value distribution

Genetic (6)

1.81

Significant Tracks

P Value

Effect Size

Relative Integration Data vs. Background

Housekeeping genes

0.0

1.81

	In	Out
HIV	73 (18.48%)	322 (81.52%)
Background	508 (5.08%)	9492 (94.92%)

Housekeeping genes
Cell line Unknown (hg19).
The original data source is UCSC http://www.tau.ac.il/~elieis/HKG/HK_genes.txt.

RefSeq genes

0.0

1.38

	In	Out
HIV	315 (79.75%)	80 (20.25%)
Background	4699 (46.99%)	5301 (53.01%)

RefSeq genes
Cell line Unknown (hg19).
The original data source is UCSC https://genome.ucsc.edu/cgi-bin/hgTables?clade=mammal&org=Human&db=hg19&hgta_group=allTables&hgta_table=refGene.

3' UTR

1.27E-4

1.32

	In	Out
HIV	14 (3.54%)	381 (96.46%)
Background	127 (1.27%)	9873 (98.73%)

3' UTR
Cell line Unknown (hg19).
The original data source is UCSC https://genome.ucsc.edu/cgi-bin/hgTables?clade=mammal&org=Human&db=hg19&hgta_group=genes&hgta_track=refSeqComposite&hgta_table=refGene.

Coding Exons

2.54E-4

1.3

	In	Out
HIV	13 (3.29%)	382 (96.71%)
Background	119 (1.19%)	9881 (98.81%)

Coding Exons
Cell line Unknown (hg19).
The original data source is UCSC https://genome.ucsc.edu/cgi-bin/hgTables?clade=mammal&org=Human&db=hg19&hgta_group=genes&hgta_track=refSeqComposite&hgta_table=refGene.

Exons

2.33E-6

1.24

	In	Out
HIV	27 (6.84%)	368 (93.16%)
Background	276 (2.76%)	9724 (97.24%)

Exons
Cell line Unknown (hg19).
The original data source is UCSC https://genome.ucsc.edu/cgi-bin/hgTables?clade=mammal&org=Human&db=hg19&hgta_group=genes&hgta_track=refSeqComposite&hgta_table=refGene.

Introns

0.0

1.09

	In	Out
HIV	293 (74.18%)	102 (25.82%)
Background	4473 (44.73%)	5527 (55.27%)

Introns
Cell line Unknown (hg19).
The original data source is UCSC https://genome.ucsc.edu/cgi-bin/hgTables?clade=mammal&org=Human&db=hg19&hgta_group=genes&hgta_track=refSeqComposite&hgta_table=refGene.

Restriction Sites (1)

0.99

Significant Tracks

P Value

Effect Size

Relative Integration Data vs. Background

BamHI

1.99E-4

0.99

	In	Out
HIV	24 (6.08%)	371 (93.92%)
Background	1229 (12.29%)	8771 (87.71%)

Restriction enzyme
Cell line Unknown (hg19).
The original data source is Other

Covariates (0)

chr1

chr2

chr3

chr4

chr5

chr6

chr7

chr8

chr9

chr10

chr11

chr12

chr13

chr14

chr15

chr16

chr17

chr18

chr19

chr20

chr21

chr22

chrX

chrY

Each blue bar represents the relative count of integration sites in this part of the genome. Lighter bars show few, and darker bars show an increased count of sites.

Positional intersect of the most significant tracks

Copy this code and run it on your machine using R to generate a barplot of the results of your analysis. You can select a subset of tracks by clicking 'change tracks' on the left.

dat = data.frame(
c(1.27E-4, 0.0268, 0.505, 2.54E-4, 1.0, 0.0691, 0.0803, 0.141, 0.00566, 0.931, 2.33E-6, 0.426, 0.0, 0.0, 0.0, 0.0906, 0.0283, 1.99E-4, 0.758, 0.705),
c(1.32, 1.04, 0.55, 1.3, 0.0, 2.09, 1.04, 0.61, 0.54, 0.01, 1.24, 0.7, 1.81, 1.09, 1.38, 0.5, 0.42, 0.99, 0.02, 0.03),
c(14, 6, 1, 13, 395, 0, 4, 10, 49, 64, 27, 0, 73, 293, 315, 22, 52, 24, 217, 213),
c(127, 61, 49, 119, 10000, 83, 236, 400, 843, 1604, 276, 16, 508, 4473, 4699, 388, 980, 1229, 5415, 5489),
c(3.54, 1.52, 0.25, 3.29, 100.0, 0.0, 1.01, 2.53, 12.41, 16.2, 6.84, 0.0, 18.48, 74.18, 79.75, 5.57, 13.16, 6.08, 54.94, 53.92),
c(1.27, 0.61, 0.49, 1.19, 100.0, 0.83, 2.36, 4.0, 8.43, 16.04, 2.76, 0.16, 5.08, 44.73, 46.99, 3.88, 9.8, 12.29, 54.15, 54.89)
)

row.names(dat) <- c("3' UTR", "5' UTR", "Blacklisted regions", "Coding Exons", "Contigs", "CpG islands", "CpG islands with 1kb window", "CpG islands with 2kb window", "CpG islands with 5kb window", "DNAse cluster regions", "Exons", "FANTOM5 CAGE Peaks", "Housekeeping genes", "Introns", "RefSeq genes", "TSS with 2kb window", "TSS with 5kb window", "BamHI", "EcoRII-1", "EcoRII-2")
colnames(dat) <- c('p value','effect size','count positions measured inside','count positions expected inside ','percent positions measured inside','percent positions expected inside')

#### Barplot ####

#bottom, left, top and right margins
par(mar=c(15,5,5,5))
fontsize=1

perc <- as.matrix(rbind(dat$`percent positions measured inside`,dat$`percent positions expected inside`))
colnames(perc) <- row.names(dat)

mp <- barplot(perc,las=2,
              ylim=c(0,115),
              col = c('blue', 'lightgray'),
              #log='y',
              beside=TRUE,
              cex=fontsize,
              cex.axis=fontsize,
              cex.lab=fontsize,
              ylab = '% inside')

#### pvalues####
pval = 0.05 # significance level 5%
pvals <- dat$`p value` #pvalues

#Plot stars above each significant bar group:
stars <- lapply(pvals,function(x)(ifelse(x < pval, '*', '')))
text(mp[1,] + 0.5, apply(perc,2,max)+2, stars,col='red', cex=2)

#Or plot the p value for each two bars:
#text(mp[1,] + 0.5, apply(perc,2,max)+4, pvals,col='red')

#### Legend ####-
legend('topright', legend=c('% inside HIV', '% inside control', 'p value < 0.05'), col=c('blue','lightgray','red'), bty='n', pch = c(15, 15, 8))

#### Vertical percent values, delete or comment out both lines to remove from plot ####
percVals <- lapply(perc, function(x) paste(toString(x[1]),'%'))
text(mp, apply(perc, 1:2, function(x) max(x-5,20)), percVals,col='coral',srt=90)

### column number warning ####
if (length(colnames(perc)) > 20) {
    title(main = 'You selected too many tracks in Enhort, 
 try to run the analysis with fewer tracks to get a better plot.')
}

Download background positions as .bed-file Download results as .csv-file