pore-c icon indicating copy to clipboard operation
pore-c copied to clipboard

multi-way contacts file

Open erbon7 opened this issue 4 years ago • 3 comments

Hi, I was wondering if the multi-way contacts derived from the alignment file are stored in a dedicated file and if yes which one ? The "alignments.parquet" file contains only filtered pairwise alignments if I understood correctly ? Thx

erbon7 avatar Sep 22 '20 09:09 erbon7

@erbon7 sorry for the delayed response. In version 0.3.0 there's *.contacts.parquet file with all pairwise contacts. You can use the read_id column to group these to see all the contacts for a given contcatemer/read. Is this what you meant?

eharr avatar Feb 11 '21 20:02 eharr

Thanks Eoghan for the answer. I just have another question. In the preprint there is a graph analysis procedure to select the best alignments. Is the code to do that part of the pore-c distribution ?

erbon7 avatar Feb 24 '21 11:02 erbon7

Yes - there the 2 functions linked here carry it out: https://github.com/nanoporetech/pore-c/blob/a865b7ea64d7a972fe878e8f9ea4aa7553def363/pore_c/analyses/alignments.py#L181-L230

eharr avatar Feb 24 '21 15:02 eharr