Skip to content

Resquiggle TSV formatting issue: kmer base with "." as end_raw_idx includes additional tab #180

Open
@williamphu

Description

@williamphu

Hi,

I haven't found a similar issue, so apologies if this has already been addressed.

The output of f5c tsv resquiggle appears to add an additional tab at the end of the line in cases where there is no corresponding signal for a kmer, at least in the case of running f5c resquiggle using a direct RNA sequencing dataset.

Quick look at the code suggest maybe this is the incorrect section?

f5c/src/resquiggle.c

Lines 411 to 412 in a88df46

if(signal_end_point<0){
printf(".\t");

Here's an example of what I mean:
Screenshot 2024-10-08 175310

sed 's/\.\t\.\t/\.\t\./g' seemed to fix it.

Thanks.

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions