empty gfa file

[rthapa1]

Hi,I am trying to create a graph file using multiple files. I used the following command but got an empty gfa output. Could you provide any suggestions? Thanks

minigraph -xggs -t16 reference.fa sample1.fa sample2.fa sample3.fa sample4.fa sample5.fa sample6.fa> out.gfa

[rthapa1]

If I use only two query samples at a time, it works and I got gfa output. But the bed file with structural variants looks weird.

[lh3]

The BED file looks ok.

What is the log message when you run for multiple input fasta files? Thanks.

[rthapa1]

Thank you, how could we know if there is deletion and translocation? Is the position in bed file corresponding to position in genome assembly?

When I run multiple fast files, it says segmentation fault.

[lh3]

BED shows the reference coordinate. Minigraph doesn't generate translocations. I may need your input files to understand the segfault.

[rthapa1]

Thank you. I have sent you the files at [email protected].

[lh3]

Thanks for the email. Which minigraph version are you using? I tried the minigraph but didn't see an error in the asan debugging mode.

[rthapa1]

It is minigraph 0.14-r415. There is no error but the output gfa file is empty.