Skip to content

Could not recreate the sequence for read #33

@Louis-MG

Description

@Louis-MG

Hello, I had this error occurred with a getdamage command, and I do not really understand how to solve it, besides deleting the line from the sam file.

Here is the output of metaDMG:

	-> metaDMG version: v0.4-123-g2cf6a01 (htslib: 1.17-28-ge13611a9) build(May  7 2025 09:54:35)
	-> ../../tools/metaDMG-cpp/metaDMG-cpp getdamage --threads 50 -o CarpeDeam/metaDMG/LDCbagSt6_S11_L001 /mnt/projects_tn03/AD_lab_old_dna/DATA/post_reu/CarpeDeam/alignments/LDCbagSt6_S11_L001.sorted_by_name.bam 
	-> ./metaDMG-cpp refName: (null) min_length: 35 print_length: 5 run_mode: 0 out_prefix: CarpeDeam/metaDMG/LDCbagSt6_S11_L001 nthreads: 50 ignore_errors: 0
skipping: M03353:64:000000000-CRD57:1:1101:1882:14211 unmapped, this msg is printed: 2 times more
skipping: M03353:64:000000000-CRD57:1:1101:1890:15773 unmapped, this msg is printed: 1 times more
skipping: M03353:64:000000000-CRD57:1:1101:1890:16478 unmapped, this msg is printed: 0 times more
Could not recreate the sequence for read: M03353:64:000000000-CRD57:1:2105:10530:2364 pp.first->s: CMGIIIIIIIIIIIIIIIMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMTMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMGMMMMMMMMMMMMMMMMMMMMMMMGMMMMMMMMMCMMMMMMMMMMMTMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMGMMMMMMAMMMMMMMMMMMMMAMCMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMGMMTMMCMMTMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMCMMMMMMMMGMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMMM strlen():484

Here is the guilty in the sam file:

M03353:64:000000000-CRD57:1:2105:10530:2364	16	NZ_CP024422.1	2000169	0	3M15I215M44D247M	*	0	0	ATCAACGGCATCAAGCTGATCGCCCGCCGCGTCGACGGCGTTTCCGGCAAGGAACTCGGCCCGCTGGTGGATGAGATGAAGGCCAGGCTGGGTTCGGGCGCCGTCGTGGTGCTGGCCGAGGCCGAGGGCAAGGCCACCGTCGCCGCCGCCGTCACCCCGGACCTGACCGCGCTCGTCAGCGCGGTGGACCTGGTCCAGGCCGCGACCGCGGCACTGGGCGGCAAGGACGGCGGTGGCCGACCAGGCCGTGGCCGCCGTCGAAACCCTGATCGGAGAGAAGGCATGACCGCGCTCTGGATTGCCCATGTGACCGTCACCGACCCCGAGGCTTACGGCGAATACGCCAGGCTCGCCGGTCCGGCCATCGCCGCCCATGGCGGGGTCTTCCTGGCGCGCGGCTCGCGCTATGTGCAGCTTGAAGGCGCCGACCGCGCCCGCAATGTCGTCGCCCGCTTTCCCAGCGTCGAGGATGCGGTGGCC	--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------==$$=====$==$=$==$=$=-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------	AS:i:-235	XS:i:-235	XN:i:0	XM:i:17	XO:i:2	XG:i:59	NM:i:76	MD:Z:0C1G80T49G23G9C11T31G6^CGGCCGCCCCGACCGTGCCCAGGGCGGCGCCCCCAGCCTGGCGC0A13A1C73G2T2C2T107C8G30	YT:Z:UU

EDIT: removing this sequence doe snot help, at least one other sequence has a similar problem.

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions