-
Notifications
You must be signed in to change notification settings - Fork 5
Open
Description
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():484Here 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:UUEDIT: removing this sequence doe snot help, at least one other sequence has a similar problem.
Reactions are currently unavailable
Metadata
Metadata
Assignees
Labels
No labels