Edit:
Never mind. It was a simple database error. Fixed it, but don't know how to delete post...
I am trying to map around 300 cDNA sequences to a genome using GMAP. I encountered some problems installing GMAP, but eventually I think it worked out. And I also managed to build an index of my genome. But when I run GMAP like this:
The program never progresses through:
Any idea what could be the problem?
Edit:
It seems that the problem must be with some of the options I used for output, because when I run it like this:
the program progresses further, but gives me an error:
Never mind. It was a simple database error. Fixed it, but don't know how to delete post...
I am trying to map around 300 cDNA sequences to a genome using GMAP. I encountered some problems installing GMAP, but eventually I think it worked out. And I also managed to build an index of my genome. But when I run GMAP like this:
Code:
gmap -d My_genome -2 ~/cDNAs.fasta --format=3
Code:
Checking compiler assumptions for popcnt: 000041A7 __builtin_clz=17 __builtin_ctz=0 _mm_popcnt_u32=7 __builtin_popcount=7 Checking compiler assumptions for SSE2: 000041A7 10D63AF1 xor=10D67B56 Checking compiler assumptions for SSE4.1: -89 -15 max=241 => compiler zero extends Finished checking compiler assumptions
Edit:
It seems that the problem must be with some of the options I used for output, because when I run it like this:
Code:
gmap -d My_genome ~/cDNAs.fasta
Code:
Checking compiler assumptions for popcnt: 000041A7 __builtin_clz=17 __builtin_ctz=0 _mm_popcnt_u32=7 __builtin_popcount=7 Checking compiler assumptions for SSE2: 000041A7 10D63AF1 xor=10D67B56 Checking compiler assumptions for SSE4.1: -89 -15 max=241 => compiler zero extends Finished checking compiler assumptions Note: >1 sequence detected, so index files are being memory mapped. GMAP can run slowly at first while the computer starts to accumulate pages from the hard disk into its cache. To copy index files into RAM instead of memory mapping, use -B 3, -B 4, or -B 5, if you have enough RAM. For more speed, also try multiple threads (-t <int>), if you have multiple processors or cores. Pre-loading compressed genome (oligos).....,......,...,...,...,...,...done (58,449,588 bytes, 14270 pages, 0.09 sec) Pre-loading compressed genome (bits).....,......,...,...,...,...,...done (58,449,600 bytes, 14270 pages, 0.10 sec) Cannot find genomic index files in either current or old format. Looking for files containing ref
Comment