Error: The specified genome ./data/BSgenome.Hsapiens.UCSC.hg19 does not have the extension of a fasta file (fa,fasta,fna)
1
0
Entering edit mode
ezgi_can00 • 0
@ezgi_can00-9811
Last seen 5.8 years ago

I am very new in rna-seq analysis in R so I am sorry that I have confused somethings in the analysis.

I am using quasr package for hg19 genome. I have generated sample file and download BSgenome.Hsapiens.UCSC.hg19 for genome file. When I run 

>proj <- qAlign(sampleFile, genome=genomeFile, maxHits=1, splicedAlignment=FALSE, alignmentsDir=results, clObj=cl, cacheDir=results)

command, I get an exactly error like this:

Error: The specified genome ./data/BSgenome.Hsapiens.UCSC.hg19 does not have the extension of a fasta file (fa,fasta,fna)

and then I have download chrom.Fa file and used this commands,

tar -zxvf chromFa.tar.gz
and
cat chr*.fa > hg19.fa

then I have run qAlign again and still get same error. I have confused what is wrong.

 

quasr error • 862 views
ADD COMMENT
0
Entering edit mode
@hotz-hans-rudolf-3951
Last seen 13 months ago
Switzerland

Hi

You don't say how you have definded "genomeFile". Nevertheless, try this:

proj <- qAlign(sampleFile, genome="BSgenome.Hsapiens.UCSC.hg19", maxHits=1, splicedAlignment=FALSE, alignmentsDir=results, clObj=cl, cacheDir=results)

 

Hope this helps, Hans-Rudolf

 

PS: and for future references: please include the output of  sessionInfo()

ADD COMMENT

Login before adding your answer.

Traffic: 200 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6