I try to map illumina results to genome by bwa.
However, I have stalled at the entrance.
I prepared concatenated human genome sequence derived from UCSC.
It was removed random, Un, hap files in advance.
Although I entered following command, I encountered error message.
$ bwa index -a bwtsw hg19.fa
$ [bwa_index] fail to open 'hg19.fa'. Abort!
When I used only chromosome 1 instead of whole chromosomes, program was not stopped.
I guess bwa on my PC is functional and I felt point is file size.
Nevertheless, size of hg19.fa is 2.9 GB.
In addition, everyone reported on the net can use hg19.fa.
Could you please help me?
However, I have stalled at the entrance.
I prepared concatenated human genome sequence derived from UCSC.
It was removed random, Un, hap files in advance.
Although I entered following command, I encountered error message.
$ bwa index -a bwtsw hg19.fa
$ [bwa_index] fail to open 'hg19.fa'. Abort!
When I used only chromosome 1 instead of whole chromosomes, program was not stopped.
I guess bwa on my PC is functional and I felt point is file size.
Nevertheless, size of hg19.fa is 2.9 GB.
In addition, everyone reported on the net can use hg19.fa.
Could you please help me?
Comment