xa6xa6 / metaOthello

Other
16 stars 7 forks source link

segfault during build #7

Open 1618gold opened 7 years ago

1618gold commented 7 years ago

_3_mytax.txt I'd like to construct a custom database with a couple million sequences. To start, I decided to try and make a custom database with just 3 sequences from greengenes to get going. This is the output I get from the build:

OpenFile Kmer_files/1.Kmer 17c73b0 OpenFile to write Kmers /srv/data1/greengenes/classifier_benchmarks/metaOthello/gg_test/_3_seqs/1.bintmp 17c73b0 OpenFile to read Kmers /srv/data1/greengenes/classifier_benchmarks/metaOthello/gg_test/_3_seqs/1.bintmp 17c93e0 OpenFile Kmer_files/2.Kmer 17c9620 OpenFile to write Kmers /srv/data1/greengenes/classifier_benchmarks/metaOthello/gg_test/_3_seqs/2.bintmp 17c9620 OpenFile to read Kmers /srv/data1/greengenes/classifier_benchmarks/metaOthello/gg_test/_3_seqs/2.bintmp 17d1590 OpenFile Kmer_files/3.Kmer 17d17d0 OpenFile to write Kmers /srv/data1/greengenes/classifier_benchmarks/metaOthello/gg_test/_3_seqs/3.bintmp 17d17d0 OpenFile to read Kmers /srv/data1/greengenes/classifier_benchmarks/metaOthello/gg_test/_3_seqs/3.bintmp 17d17d0 Segmentation fault (core dumped)

I constructed the kmer files as instructed (I used Jellyfish 2 with a kmer size of 31). This is an example command line for jellyfish that I ran:

~/jellyfish-2.2.6/bin/jellyfish count -o _1_gg.jf -m 31 -t 3 -s 10M -C _1_gg.fa

This is the build command that I ran:

build _3_mytax.txt Kmer_files/ .Kmer 31 6 mybacterial_index ./

Since its definitely possible I messed up the taxonomy info file, I've attached it here (_3_mytax.txt). Obviously extremely simple given my 3 sequences.

xa6xa6 commented 7 years ago

Thanks for your interest in MetaOthello. I am busy with other duties these days but will get back to you soon.