8000 make MAXK=31 NUM_COLS=1 cortex_var · Issue #22 · iqbal-lab/cortex · GitHub
[go: up one dir, main page]
More Web Proxy on the site http://driver.im/
Skip to content

make MAXK=31 NUM_COLS=1 cortex_var #22

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.

Sign up for GitHub

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

Jump to bottom
Open
EmilyWesel opened this issue Jul 25, 2018 · 5 comments
Open

make MAXK=31 NUM_COLS=1 cortex_var #22

EmilyWesel opened this issue Jul 25, 2018 · 5 comments

Comments

@EmilyWesel
Copy link

As instructed in #0 of the README in example, I used the commands make clean and make MAXK=31 NUM_COLS=1 cortex_var.

I received the error:
collect2: ld returned 1 exit status
make: *** [cortex_var] Error 1

I'm on a Mac.
@iqbal-lab
Copy link
Owner

Hi there, can I ask which version of Cortex you are using? Where did you get it from?

@EmilyWesel
Copy link
Author

I'm using CORTEX_release_v1.0.5.21. I got it from https://sourceforge.net/projects/cortexassembler/?source=typ_redirect

@iqbal-lab
Copy link
Owner

Can you try using git clone with this repo, and see if that fixes it?
Am also going to make a singularity container (I guess I should do docker too) too make this easier for people

@EmilyWesel
Copy link
Author

Hi,

I used git clone and got the tool installed and running! Thank you for your help. However, when I run the tool, the output is only one file, bubbles_called_from_genome, with the following contents:

var_1_5p_flank length:276 INFO:KMER=31
ACTAGCAATGAATATGCTGAGTCACTGCATATACATACTCACAACTTTATTCTATAATGTAATATTCTATAAAGTAGCTGTATCAGTTTATACTTTAACCAGTAATGGACAAGATTTTCTATTACTTCCCATCTTTGTTATTACTTTTAGACTCTAACTTTTATCAGGCTCATGGATGTAAAAAGCATCTCAGGGTGGTTTTAATTTGCATTTATCTGCTCATCTATGAAGATGAGCTTCTTTTCATATAATTATGAGTCATTATTTTTGTTTTGC
var_1_branch_1 length:32 kmer:31
ATTCTTTTGTTTATGCATTTTGCTTGTTCTAT
var_1_branch_2 length:32 kmer:31
CTTCTTTTGTTTATGCATTTTGCTTGTTCTAT
var_1_3p_flank length:352 kmer:31
GTCTTATTTTTCCTGTTGATTTTTGGGAGTTCATATATATTCTAAATGTATATTTATTCACTTATATATATGTTGTAAATATTACAGTTTATGATTTGTCACCTTATGATATCTTCCAAATAGAGAAGCTTTATATTTTGATGTAGTCATATGTTCATTTTTCCTCCTTAATGTTTGTTTTTCTTGGTTCTATGACCTACCAAAAGTAACAAAAATTCTCATTTATTTTTAATCTAAATGTTTTAAGTATTTTCCTGGAATTCACCTTGAATTGATTTCTATTGGAGATAGGTATCCAATCTAATTTGCCTCATATGGATAACCACTTGTTCTATTACTGCTGTAACAAATT

The manual advertises three output files that detail the structural variants in the input fasta files. How can I get these three files as my output?

Thanks,
Emily

@iqbal-lab
Copy link
Owner

Hi Emily
Can you give the command line you ran and copy the text which describes these 3 files?
The output you show is the raw output of the bubble caller. Usually people use the run_calls wrapper script that takes you all the way to vcf files. Is vcf what you want?
I'm actually on vacation fir the next 2 weeks, will reply when I can

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Assignees
No one assigned
Labels
None yet
Projects
None yet
Milestone
No milestone
Development

No branches or pull requests
2 participants
@iqbal-lab @EmilyWesel
0