gfapy icon indicating copy to clipboard operation
gfapy copied to clipboard

Missing coverage information after merge_linear_paths()

Open jamiemcg opened this issue 2 years ago • 0 comments

Is it possible to keep coverage/count information when using merge_linear_paths?

E.g., this is snippet of my starting graph which has FC tags:

S	35	TATTTTTATTAGAGACCTTGATTCCGATAATATTCATGTTGTATTAGCCGATTTTGGTCTAGCCAAAGCAAGTATGGGTTCAGTTGACAGATCCTATGCTGGCACTCCTTTATACATGTCTTTAGAACTCGCAACCGGAGGAAAGTACAGTTTCAATACCGATATCTATTCATTGGGAGTGGCTATCTATCAAATCATGACCAAAGACACAACCACAGCTGTCAGTCACATGTTTATTTCGAATAGAGAGGAAGAAGCAATTGATAAATTAACAAAGGAAATGGAAAAATCAGGAAATTATTCAACAACATTGATTAATATTTGTTTATCTATGCTGGAGAAGGATGCCAGTAAGAGACCAAGTGCCAAGCAAATCCTTTCCAATACTTTCTTGGAAGATATTCTCATCACACCAACAGAATAAACAATTTCTTGTGTTTTACCGCTGAAGTTAACCGTTAGAATTATAATAATTTTATTTTTCACAACAAAAAAAACAAGCCAACAACAATAGAGCGGGTGCAAACTCGGGTGTATCTACCGTTTTTCGTTTCTGAAAACGACTTTCGTTTTTGTTCGTGTTTGAATGAAGATTTTTGGTGTTTCTTTCAAATTCATACAAGATCGTTTTTCGTTTTTGTAAAAAGTTAGAAAAATACGAAAACGTGGCAAGTTGCACCTGCTCTAATGTGGGCAACAAACTTGAGGTTACCATTACCTTCTTGAATGATCATCAGAAGATTATTCATCTTCGCTTACCATCAGAGATCTGATTATTAACCGCTCTTGACTGTCAACTGAGATTTGGCAGGTGAAATTCAGTTTTCACATTCGCCAAAAGAAAATGAAGATAAATGATGATGAATTGAATAATTCTTGAATTTCTTCCTACGGATAGCTATGCATCATTCTCACATATCAATAGGCAGTGTCAATCTATTCTTGTTTCGGAAAAACCAAACTTTTGGAAACAGAAAAAGAGTGAATATGATAGAACTATTGCCAACCAAGAACGATACCTCGACCATGTTATTGATACAATGACAAGAAATGCTTCTAAGAACAATTATTGTTTATTTGCTTGTTTTGTCGGGATGTTGTTTTCGTTTCTTGTGTTTTTACCAGCATTTATGATCATGGTGATTATGGCATCTAGTTCCCTTTCAGGAAATATATTGAATTGTGAGAGGAATGTTTCTTCTATTCAAGTTTTGAATAGAACTATTGAAAGGGTTTCTTATCCAGGATGTGTAGAATGTAACACAATTTGGAATACTAGTTGCAATAATTTCAATGTTATGCCATGTTGTCAACCAAATAATAGTTATTCTTCAATTTTCTCATATGAATGCAATGTAACTGAATGCATTAAGGAAAATTATCAAGATGGACCTCTCAGATGTCGAGTAATTAGTTACAACAACAAAAAAGGGCCATCCACATATTCACTAGTCATTTACAAACCAGAAAGCAATTATTCAGGGTCTATAGGACTATTATCTGGTTAGCTACCAATACTTATGATTTGAAAGTTTAGTCAGATTACTGAAATACTGACTTTATTTCACTAACAACACGTAGGTCGGCAATGTGATTCCTTCTTCATAATTGAGTATTAATCGTACGGTATTTGCAAGTAAATAATAACTAATAATGGTCAATTAACAAAATATTAAATAAATTATTTAAAAAAAATGAAAAGAAGCAGAGAGGAGGAACCTGCTGTTACCAATACATCATGTACTGATCAAAGTAATACAAGTTCCGATCCTCTCTTGGAAACGAATCTAAAAAAACAAGTGACTGGCAATATTT	ll:f:44.5	FC:i:80768	kl:f:44.5
S	11508	ACTGTCAACTGAGATTTGGCAGGTGAAATTCAGTTTTCACATTCGCCAAAAGAAAATGAAGATAAATGATGATGAATTGAATAATTCTTGAATTTCTTCCTACGGATAGCTATGCATCATTCTCACATATCAATAGGCAGTGTCAATCTATTCTTGTTTCGGAAAAACCAAACTTTTGGAAACAGAAAAAGAGTGAATATGATAGAACTATTGCCAACCAAGAACGATACCTCGACCATGTTATTGATACAATGACAAGAAATGCTTCTAAGAACAATTATTGTTTATTTGCTTGTTTTGTCGGGATGTTGTTTTCGTTTCTTGTGTTTTTACCAGCATTTATGATCATGGTGATTATGGCATCTAGTTCCCTTTCAGGAAATATATTGAATTGTGAGAGGAATGTTTCTTCTATTCAAGTTTTGAATAGAACTATTGAAAGGGTTTCTTATCCAGGATGTGTAGAATGTAACACAATTTGGAATACTAGTTGCAATAATTTCAATGTTATGCCATGTTGTCAACCAAATAATAGTTATTCTTCAATTTTCTCATATGAATGCAATGTAACTGAATGCATTAAGGAAAATTATCAAGATGGACCTCTCAGATGTCGAGTAATTAGTTACAACAACAAAAAAGGGCCATCCACATATTCACTAGTCATTTACAAACCAGAAAGCAATTATTCAGGGTCTATAGGACTATTATCTGGTTAGCTACCAATACTTATGATTTGAAAGTTTAGTCAGATTACTGAAATACTGACTTTATTTCACTAACAACACGTAGGTCGGCAATGTGATTCCTTCTTCATAATTGAGTATTAATCGTACGGTATTTGCAAGTAAATAATAACTAATAATGGTCAATTAACAAAATATTAAATAAATTATTTAAAAAAAATGAAAAGAAGCAGAGAGGAGGAACCTGCTGTTACCAATACATCATGTACTGATCAAAGTAATACAAGTTCCGATCCTCTCTTGGAAACGAATCTAAAAAAACAAGTGACTGGCAATATTTCTAGTGAATCATCGTATTTTAGTCTTGGTTGTAATTATATGGATGGTGTTAATGGTGTTGAAAAAGATTATTTAAAAGCATTTGAATTATTTTCTAAAGGTGTTGAGGAGGGATGCTTGCGATCAATACATAAAATTGGTTACTTTTATAGCAAAGGATTGGGAGTTGAACAAGATTATGATAAAGCAGTCCAGTTCTATATGAAATCTGCTGAAGGGGGTTTTGTGGCCTCATATGCTAAGGGAAAGGCGTTCAACAGGACCATGTTAAAGCATTTGAATGGACTTTAAAAGCTGCTGAAAATGATAATACTGATGCACAATTTAACATTGGAAGATGCTATAGTAAAGGGACAGGAGTTGAGCAAGACCATGTTAAAGCATTTGAATGGTATTTAAAAGCTGCTGAAAAGGGTGATACAGATGCACAGTTTGTGATTGGA	ll:f:43	FC:i:63124	kl:f:43

But when I run graph.merge_linear_paths() the FC tag is lost:

S	35_11508	TATTTTTATTAGAGACCTTGATTCCGATAATATTCATGTTGTATTAGCCGATTTTGGTCTAGCCAAAGCAAGTATGGGTTCAGTTGACAGATCCTATGCTGGCACTCCTTTATACATGTCTTTAGAACTCGCAACCGGAGGAAAGTACAGTTTCAATACCGATATCTATTCATTGGGAGTGGCTATCTATCAAATCATGACCAAAGACACAACCACAGCTGTCAGTCACATGTTTATTTCGAATAGAGAGGAAGAAGCAATTGATAAATTAACAAAGGAAATGGAAAAATCAGGAAATTATTCAACAACATTGATTAATATTTGTTTATCTATGCTGGAGAAGGATGCCAGTAAGAGACCAAGTGCCAAGCAAATCCTTTCCAATACTTTCTTGGAAGATATTCTCATCACACCAACAGAATAAACAATTTCTTGTGTTTTACCGCTGAAGTTAACCGTTAGAATTATAATAATTTTATTTTTCACAACAAAAAAAACAAGCCAACAACAATAGAGCGGGTGCAAACTCGGGTGTATCTACCGTTTTTCGTTTCTGAAAACGACTTTCGTTTTTGTTCGTGTTTGAATGAAGATTTTTGGTGTTTCTTTCAAATTCATACAAGATCGTTTTTCGTTTTTGTAAAAAGTTAGAAAAATACGAAAACGTGGCAAGTTGCACCTGCTCTAATGTGGGCAACAAACTTGAGGTTACCATTACCTTCTTGAATGATCATCAGAAGATTATTCATCTTCGCTTACCATCAGAGATCTGATTATTAACCGCTCTTGACTGTCAACTGAGATTTGGCAGGTGAAATTCAGTTTTCACATTCGCCAAAAGAAAATGAAGATAAATGATGATGAATTGAATAATTCTTGAATTTCTTCCTACGGATAGCTATGCATCATTCTCACATATCAATAGGCAGTGTCAATCTATTCTTGTTTCGGAAAAACCAAACTTTTGGAAACAGAAAAAGAGTGAATATGATAGAACTATTGCCAACCAAGAACGATACCTCGACCATGTTATTGATACAATGACAAGAAATGCTTCTAAGAACAATTATTGTTTATTTGCTTGTTTTGTCGGGATGTTGTTTTCGTTTCTTGTGTTTTTACCAGCATTTATGATCATGGTGATTATGGCATCTAGTTCCCTTTCAGGAAATATATTGAATTGTGAGAGGAATGTTTCTTCTATTCAAGTTTTGAATAGAACTATTGAAAGGGTTTCTTATCCAGGATGTGTAGAATGTAACACAATTTGGAATACTAGTTGCAATAATTTCAATGTTATGCCATGTTGTCAACCAAATAATAGTTATTCTTCAATTTTCTCATATGAATGCAATGTAACTGAATGCATTAAGGAAAATTATCAAGATGGACCTCTCAGATGTCGAGTAATTAGTTACAACAACAAAAAAGGGCCATCCACATATTCACTAGTCATTTACAAACCAGAAAGCAATTATTCAGGGTCTATAGGACTATTATCTGGTTAGCTACCAATACTTATGATTTGAAAGTTTAGTCAGATTACTGAAATACTGACTTTATTTCACTAACAACACGTAGGTCGGCAATGTGATTCCTTCTTCATAATTGAGTATTAATCGTACGGTATTTGCAAGTAAATAATAACTAATAATGGTCAATTAACAAAATATTAAATAAATTATTTAAAAAAAATGAAAAGAAGCAGAGAGGAGGAACCTGCTGTTACCAATACATCATGTACTGATCAAAGTAATACAAGTTCCGATCCTCTCTTGGAAACGAATCTAAAAAAACAAGTGACTGGCAATATTTCTAGTGAATCATCGTATTTTAGTCTTGGTTGTAATTATATGGATGGTGTTAATGGTGTTGAAAAAGATTATTTAAAAGCATTTGAATTATTTTCTAAAGGTGTTGAGGAGGGATGCTTGCGATCAATACATAAAATTGGTTACTTTTATAGCAAAGGATTGGGAGTTGAACAAGATTATGATAAAGCAGTCCAGTTCTATATGAAATCTGCTGAAGGGGGTTTTGTGGCCTCATATGCTAAGGGAAAGGCGTTCAACAGGACCATGTTAAAGCATTTGAATGGACTTTAAAAGCTGCTGAAAATGATAATACTGATGCACAATTTAACATTGGAAGATGCTATAGTAAAGGGACAGGAGTTGAGCAAGACCATGTTAAAGCATTTGAATGGTATTTAAAAGCTGCTGAAAAGGGTGATACAGATGCACAGTTTGTGATTGGA	ll:f:44.5	kl:f:44.5	LN:i:2257

jamiemcg avatar May 15 '23 11:05 jamiemcg