RecGraph icon indicating copy to clipboard operation
RecGraph copied to clipboard

How to interprete the output information

Open hangsuUNC opened this issue 1 year ago • 7 comments

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

hangsuUNC avatar Aug 01 '24 15:08 hangsuUNC

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

gdv avatar Aug 01 '24 16:08 gdv

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!

hangsuUNC avatar Aug 01 '24 17:08 hangsuUNC

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

gdv avatar Aug 02 '24 10:08 gdv

That works! Thanks so much!

hangsuUNC avatar Aug 02 '24 13:08 hangsuUNC

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

gdv avatar Aug 02 '24 14:08 gdv

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

hangsuUNC avatar Aug 02 '24 14:08 hangsuUNC

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

dcmonti avatar Aug 03 '24 07:08 dcmonti