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

Cigar strings altered after extracting subgraph #10

Open
rlorigro opened this issue Apr 6, 2020 · 0 comments
Open

Cigar strings altered after extracting subgraph #10

rlorigro opened this issue Apr 6, 2020 · 0 comments

Comments

@rlorigro
Copy link

rlorigro commented Apr 6, 2020

Before using the command:

gfatools view -l 303726,303738 -r 20 Assembly-BothStrands.gfa > Assembly-BothStrands_303726-303738_r20.gfa

my GFA has the following elements:

S       5924062 GTTTCAAAAAAAAAAAAAGAGCATGGCTCTGG        RC:i:400
L       5209751 +       5924062 +       18M1D6M
L       5420387 +       5924062 +       18M1D6M

After running the command, they look like this:

S	5924062	GTTTCAAAAAAAAAAAAAGAGCATGGCTCTGG	LN:i:32	RC:i:400
L	5209751	+	5924062	+	24M	L1:i:33	L2:i:8
L	5420387	+	5924062	+	24M	L1:i:33	L2:i:8

These overlap Cigars are incorrect:

Cigar: 24M

5924062    1                                   GTTTCAAAAAAAAAAAAAGAGCATGGCTCTGG
                                               |||||||||||||||||||||||||
5420387    34 ACTGCATTCCAGCCTGGGTGATAGAGTGAGGCTGTTTCAAAAAAAAAAAAAAGAGCAT

Compared to a local alignment with Smith-Waterman:

Cigar:  5M1I19M

5924062    1                                   GTTTC-AAAAAAAAAAAAAGAGCATGGCTCTGG
                                               ||||| |||||||||||||||||||
5420387    34 ACTGCATTCCAGCCTGGGTGATAGAGTGAGGCTGTTTCAAAAAAAAAAAAAAGAGCAT
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

1 participant