packageVersion('GenomicFeatures')
## [1] ‘1.9.2’
packageVersion('rtracklayer')
## [1] ‘1.17.0’

Ok, without further ado:

snp135 <- makeFeatureDbFromUCSC('hg19','snp135Common','snp135Common')
## Download the snp135Common table ... Error in scan(file, what, nmax, sep,
dec, quote, skip, nlines, na.strings, :
## line 4129688 did not have 26 elements

This is rather unfortunate since it offers a very easy way to mask common
SNPs that intersect an interrogated di- or tri-nucleotide sequence.
Also, I need to get more of a feel for how a FeatureDb is supposed to
operate "in the wild" and this is an interesting, common use case.

Thoughts?

-- 
*A model is a lie that helps you see the truth.*
*
*
Howard Skipper<http://cancerres.aacrjournals.org/content/31/9/1173.full.pdf>

	[[alternative HTML version deleted]]

