Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Problem with running poxomis #24

Open
ookil opened this issue Oct 25, 2018 · 2 comments
Open

Problem with running poxomis #24

ookil opened this issue Oct 25, 2018 · 2 comments

Comments

@ookil
Copy link

ookil commented Oct 25, 2018

I ran into some problems and not sure what to do now.
I installed pomoxis how it's written in guide but I get an error with racon. I then activated venv and wanted to check if mini_assemble will work but it's not working.

Here's a screenshot.

zrzut ekranu z 2018-10-25 12-22-50
Thank you for any solutions ;)

@cjw85
Copy link
Member

cjw85 commented Oct 25, 2018

The screenshot indicates that cmake is not installed on your system. This is required to compile the third-party racon code. You should install this using your system package manager. We will update the README to make this clear. Sorry for the inconvenience.

@ookil
Copy link
Author

ookil commented Oct 25, 2018

Thank you. But I got another problem, hope it's okay to post it here

(pomoxis) umian@umach:~/pomoxis$ mini_assemble -i /media/umian/Toshiba/candida_pacbio.fasta -r /media/umian/Toshiba/candida_reference1.fasta -o /media/umian/Toshiba/Ref_nanopore -p ref_assembly
Copying FASTA input to workspace: /media/umian/Toshiba/candida_pacbio.fasta > /media/umian/Toshiba/Ref_nanopore/ref_assembly.fasta
Skipped adapter trimming.
Skipped pre-assembly correction.
Using supplied reference to perform reference-guided consensus.
Running racon read shuffle 1...
Running round 1 consensus...
[M::mm_idx_gen::0.3780.80] collected minimizers
[M::mm_idx_gen::0.585
0.87] sorted minimizers
[M::main::0.5850.87] loaded/built the index for 14 target sequence(s)
[M::mm_mapopt_update::0.613
0.88] mid_occ = 14
[M::mm_idx_stat] kmer size: 15; skip: 10; is_hpc: 0; #seq: 14
[M::mm_idx_stat::0.6380.88] distinct minimizers: 2204001 (96.66% are singletons); average occurrences: 1.048; average spacing: 5.343
[M::worker_pipeline::28.674
0.31] mapped 19 sequences
[M::main] Version: 2.10-r761
[M::main] CMD: minimap2 -t1 /media/umian/Toshiba/candida_reference1.fasta ref_assembly.fasta
[M::main] Real time: 28.740 sec; CPU: 9.016 sec
[racon::Polisher::initialize] loaded target sequences
[bioparser::FastaParser] error: invalid file format!

Im not sure which file has invalid format. Both are in fasta which should be okay comparing with docs. But a little background:
1.Reference file was from NCBI so I changed it to fasta --> cp candida_reference.fna candida_reference.fasta (At first I put it in .fna format but got an error early on)
2. Thought maybe I should change it from multiline sequences to one so I did it with cat candida_reference.fasta | awk '{if (substr($0,1,1)==">"){if (p){print "\n";} print $0} else printf("%s",$0);p++;}END{print "\n"}' > candida_reference1.fasta as recommended on some forum
3. As a last hope did the same to input fasta

Well in each case I recived the same error. Any idea what's wrong? Or should I seek help on racon sites?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants