-
Notifications
You must be signed in to change notification settings - Fork 33
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
Lachesis failed at the step of LoadDenovoCLMSFromSAMs. #16
Comments
Hi, |
Hi, I remembered the problem was result from un-sorted bam. Maybe it is a problem of SortSam.jar in picard (so I change it to SortSame). Or you had to sort the bam before the step of "samtools flagstat" (I had a new file called *consistent.bam). It's a long time since I solve this problem and I couldn't get the files. So I would give the relate code that was used to revised the file of PreprocessSAMs.pl blow, Hope it would help. By the way, you could sorted bam file before "samtools flagstat", then continue the following process if it failed without repeat the previously steps.
|
Hi, Thanks for your message. |
Hi,@lamz138138, |
@Epigenetics-Wang Hi, I remember it was caused by failing to sort the bam, that is the output bam isn't sorted by name as expected. And I think filterInconsistentPairReads.pl is my script, but I forget what is used for (it's a pity I can't find this scripts now). I had solved this problem by running the PrepocessSAMs.pl step by step, then figure out whether bam is sorted by name. Then if it still failed to work, try to find read that make "c2.tid == c1.mtid" is false (read name is different which may result from failing to sort bam), then write filterInconsistentPairReads.pl. Hope this useful. |
Hi,@lamz138138, |
Hi, everyone!
Lachesis is a grate software to link, order scaffolds. After test the software with the demo, I failed to run the script with my data. It failed at the step of LoadDenovoCLMSFromSAMs. According to the script of "ChromLinkMatrix.cc", it said "an internal inconsistency in the SAM file. Maybe the file is truncated or it has an incorrect header". Considering it may be caused by sort problem, I had re-run the sort step with picard by sorting the file with reads name. However, it was still failed, so I had following question:
Any suggestion would be grateful!
Error and log were list below:
\1. Error:
\2. Log:
The text was updated successfully, but these errors were encountered: