search for: chr7

I have a dataframe in the general format: chr1 0.5 chr1 0 chr1 0.75 chr2 0 chr2 0 chr3 1 chr3 1 chr3 0.5 chr7 0.75 chr9 1 chr9 1 chr22 0.5 chr22 0.5 where the first column is the chromosome location and the second column is some value. What I'd like to do is have a histogram created for each chr location (i.e. a separate histogram for chr1, chr2, chr3, chr7, chr9, and chr22). I am just having a hard t...

...3 V4 V5 V6 V7 1 chr1 500 500 CHH 0 0.5 + 2 chr1 550 550 CHH 0 0.0 + 3 chr2 700 700 CHH 0 0.0 + 4 chr2 1000 1000 CHH 0 0.0 + 5 chr3 100 100 CHH 0 0.0 + 6 chr4 450 450 CG 0 0.0 + 7 chr5 450 450 CHH 0 0.0 + 8 chr5 50034 50034 CHG 0 0.0 + 9 chr7 50055 50055 CHG 0 0.0 + 10 chr10 50063 50063 CHH 0 0.0 + > dput(data) structure(list(V1 = structure(c(1L, 1L, 3L, 3L, 4L, 5L, 6L, 6L, 7L, 2L), .Label = c("chr1", "chr10", "chr2", "chr3", "chr4", "chr5", "chr7"), class =...

Hi, I was wondering I'm going about this in the correct way. I need to test if there are coding sequences or exons in hg19 which match a string of 100bp "D" i.e. [A,G or T]. However I'm getting a strange result. I get a hit on chr7, using the 100bp search however when I search with 60bp sequence of "D" I don't get any hits. library("BSgenome") library("Biostrings") library("BSgenome.Hsapiens.UCSC.hg19") library("biomaRt") library("GenomicFeatures") #extract...

Dear all, I have the following data frame: V1 V2 aaa chr1 aaa chr2 aaa NM aaa QC aaa chr10 att NM att chr7 What I want to do is to count the string (V1). But the condition of counting is: if the V2 of the string is "NM" or "QC" then the count is not increased. Hence the contigency table will look like this: #tag count aaa 3 att 1 Is there a compact way to achieve that...

...t;chr1", "chr10", "chr11", "chr12", "chr13", "chr14", "chr15", "chr16", "chr17", "chr18", "chr19", "chr2", "chr3", "chr4", "chr5", "chr6", "chr7", "chr8", "chr9", "chrX"), class = "factor"), start = c(5291000L, 10988025L, 11767950L, 11840900L, 12267450L, 12276675L), end = c(5291926L, 10988526L, 11768676L, 11841851L, 12268076L, 12277051L), peak = c(8L, 7L, 8L, 8L, 12L, 7L)), .Names = c("chrom...

...it out. I'm trying to make a fairly simple figure, but also trying to make it look "better" using ggplot2. I have data.frame that contains chrom length 1 chr1 249250621 2 chr2 243199373 3 chr3 198022430 4 chr4 191154276 5 chr5 180915260 6 chr6 171115067 7 chr7 159138663 8 chrX 155270560 9 chr8 146364022 10 chr9 141213431 11 chr10 135534747 12 chr11 135006516 13 chr12 133851895 14 chr13 115169878 15 chr14 107349540 16 chr15 102531392 17 chr16 90354753 18 chr17 81195210 19 chr18 78077248 20 chr20 63025520 21 chrY 59373566 22 chr19...

Hello, I have an input file: http://r.789695.n4.nabble.com/file/n3700031/testOut.txt testOut.txt where col 1 is chromosome, column2 is start of region, column 3 is end of region, column 4 and 5 is base position, column 6 is total reads, column 7 is methylation data, and column 8 is the strand. I would like a summary output file such as: