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]]