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

How to interprete the output information #8

Open
hangsuUNC opened this issue Aug 1, 2024 · 7 comments
Open

How to interprete the output information #8

hangsuUNC opened this issue Aug 1, 2024 · 7 comments

Comments

@hangsuUNC
Copy link

Hi,

I tried the program using the files in the example folder. The output is listed as follows. Is this something expected? How to interprete this output?
image

Thanks,

Hang

@gdv
Copy link
Member

gdv commented Aug 1, 2024

May you please provide us how you have installed recgraph and the complete command line that you have used to run the program?

@hangsuUNC
Copy link
Author

Installation command:

git clone https://github.com/AlgoLab/RecGraph.git
cd RecGraph
cargo build --release

commandline:
cargo run --release ./example/reads.fa ./example/graph.gfa

Thanks!

@gdv
Copy link
Member

gdv commented Aug 2, 2024

You need to use the -m option to determine the type of alignment you want to compute.
For example

target/release/recgraph -m 9 ./example/reads.fa ./example/graph.gfa

if you want to allow recombinations. See https://github.com/AlgoLab/RecGraph#alignment-modes

Best regards

@hangsuUNC
Copy link
Author

That works! Thanks so much!

@gdv
Copy link
Member

gdv commented Aug 2, 2024

We are going to release the new version soon, which is way faster.

@hangsuUNC hangsuUNC reopened this Aug 2, 2024
@hangsuUNC
Copy link
Author

When I use the api.rs provided in this crate, I found the same output issues when using the functions as follows:

        let score_matrix = create_score_matrix_f32(Some(0), Some(1), None);
        let alignment = align_global_no_gap(&read_seq_upper, &g, Some(h_tuple), Some(score_matrix), Some(0.1));

        let gaf_output = alignment.to_string();
        gaf_file.write_all(gaf_output.as_bytes()).expect("Unable to write to file");

Is there any parameters that I should set here?

Thanks,

Hang

@dcmonti
Copy link
Collaborator

dcmonti commented Aug 3, 2024

You can try by setting the last parameter of align_global_no_gap() to Some(1.0), it should work

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

3 participants