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

unitigger segfaults #1

Open
imranfanaswala opened this issue May 14, 2015 · 2 comments
Open

unitigger segfaults #1

imranfanaswala opened this issue May 14, 2015 · 2 comments

Comments

@imranfanaswala
Copy link

Hello,

'tigger' segfaults on the attached FASTA file (rename the file to .png). Here are the commands:

nucmer -maxmatch -c 40 sample.fa sample.fa -p sample
show-coords -H -c -l -o -r sample.delta | nucmerAnnotate -ignore 180 | egrep 'BEGIN|END|CONTAIN|IDENTITY' > sample.coords
nucmer2ovl -ignore 180 -tab sample.coords | sort2 > sample.ovl
ovl2OVL sample.ovl > sample.OVL
toAmos -o sample.afg -s sample.fa
bank-transact -c -z -b sample.bnk -m sample.afg
bank-transact    -z -b sample.bnk -m sample.OVL
tigger -v 3 -g -b sample.bnk
![sample](https://cloud.githubusercontent.com/assets/422028/7638044/763d9414-fa72-11e4-8fbe-813602f61eca.png)


The stack trace shows:
#0 Graph::new_edge (this=0x87e1048, p_n1=0x87e8260, p_n2=0x0, p_element=0x87e8bf0) at Graph.cc:217

217 p_n2->add_iedge(e);

@imranfanaswala
Copy link
Author

sample

@imranfanaswala
Copy link
Author

For posterity, I leaving the solution here... Issue can be fixed by first calling 'dumpreads', therefore the correct calls should be:

toAmos -o sample.afg -s sample.fa
bank-transact -c -z -b sample.bnk -m sample.afg
dumpreads sample.bnk > sample_amos.fa
nucmer -maxmatch -c 150 sample_amos.fa sample_amos.fa -p sample
show-coords -H -c -l -o -r sample.delta | nucmerAnnotate -ignore 180 | egrep 'BEGIN|END|CONTAIN|IDENTITY' > sample.coords
nucmer2ovl -ignore 180 -tab sample.coords | sort2 > sample.ovl
ovl2OVL sample.ovl > sample.OVL
bank-transact -z -b sample.bnk -m sample.OVL
tigger -v 3 -g -b sample.bnk

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

1 participant