Displaying 1 result from an estimated 1 matches for "60bp".
Did you mean:
60b
2011 Apr 15
1
Whole genome searching of 100bp "D" sequence
...,
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 the alignments which match real genes
#(add this onto s...