Using Rsubread buildindex with GRCh37.p13.genome.fa.gz gives me an error

0

Entering edit mode

efernandez • 0

@efernandez-22025

Last seen 2.5 years ago

Argentina

Hi I am triying to build the human index using

ftp.ebi.ac.uk/pub/databases/gencode/Gencode_human/release_19/GRCh37.p13.genome.fa.gz

I am using Rsubread 2.4.3 an it gives me the following error

//================================= Running ==================================\ || || || Check the integrity of provided reference sequences ... || ERROR: A fasta file cannot have a line longer than 1000 bytes. You need to split a very long line into many lines. || || \============================================================================//

Any advise?

Thank you very much in advance.

Rsubread • 1.0k views

ADD COMMENT • link updated 2.6 years ago by Gordon Smyth 51k • written 2.6 years ago by efernandez • 0

0

Entering edit mode

Is there a reason why you're using such an old version of the human genome and such an old version of Rsubread? The current version of the genome is available from https://www.gencodegenes.org/human/ and the current version of Rsubread is here: https://bioconductor.org/packages/release/bioc/html/Rsubread.html

ADD REPLY • link 2.6 years ago Gordon Smyth 51k

Login before adding your answer.