Skip to content

minigraph --call output #134

Description

@panosioannidis

I'm running minigraph -cxasm --call out.gfa sample.fa > sample.bed to find the coverage of sample in the bubbles of the graph. The output is not hard to understand, but I am a bit confused with the coordinates.

For one particular bubble I get the following 3 different types of output (for 3 different samples):

Sample Th_026

CHM13#1#chr1	3690	3752	>s1	>s3	>s2:62:+:Th_026#1#ptg000057l:2751:2893

Sample Th_009

CHM13#1#chr1	3690	3752	>s1	>s3	.

Sample Th_137

CHM13#1#chr1	3690	3752	>s1	>s3	*:0:+:Th_137#1#ptg000037l:5644:5724

My interpretation is this:
For Th_026, there's an insertion. The insertion is on contig ptg000057l of Th_026, from position 2751 to 2893.

For Th_009, there's nothing (no insertion, no deletion); it contains the reference sequence.

For Th_137, there's a deletion (denoted by the *). The deletion is on contig ptg0000037l, but I don't understand what the coordinates mean; how can it be that the deletion goes from 5644 to 5724 on that contig of Th_137? In contrast, I would understand it if it had a range on the CHM13 genome (i.e. the reference used for building the graph). It would mean that this sequence from CHM13 is missing from sample Th_137...

What am I missing?

Activity

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

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