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

High Error rate on short read simulation #1161

Closed
Hakuna-Mataha opened this issue Feb 14, 2024 · 5 comments
Closed

High Error rate on short read simulation #1161

Hakuna-Mataha opened this issue Feb 14, 2024 · 5 comments
Labels

Comments

@Hakuna-Mataha
Copy link

Hello,
I am doing schoolwork that requires me to reproduce Figure 1b in the minimap2 paper. I am using this human genome
https://www.ncbi.nlm.nih.gov/datasets/genome/GCF_000001405.26/ with simulation by mason2 and analysis by paf tools. I am using papers parameters.

I have attached results from paftools mapeval:
As you see error rates are 50 percent on all brackets. Is something wrong with my method?
Thank you very much

quality_new.txt

@lh3
Copy link
Owner

lh3 commented Feb 14, 2024

What are the exact command lines and input files?

@lh3 lh3 added the question label Feb 14, 2024
@Hakuna-Mataha
Copy link
Author

Here is the code.
mason_simulator --seed 2415 --illumina-prob-mismatch-scale 2.5 --illumina-read-length 150 -ir GCF_000001405.26_GRCh38_genomic.fna -n 10000000 -o tmp-l.fq -or tmp-r.fq -oa tmp.sam
~/k8/k8 ./paftools.js mason2fq tmp.sam | seqtk seq -1 > human_mason_1.fq
~/k8/k8 ./paftools.js mason2fq tmp.sam | seqtk seq -2 > human_mason_2.fq
minimap2 -ax sr GCF_000001405.26_GRCh38_genomic.fna human_mason_1.fq human_mason_2.fq > test.sam
~/k8/k8 paftools.js mapeval test.sam > quality.txt

@Hakuna-Mataha
Copy link
Author

Hi,
Are there any updates on this?
Thank you very much

@lh3
Copy link
Owner

lh3 commented Mar 11, 2024

Could you try older k8, not k8-1.0?

@lh3 lh3 closed this as completed in ba60aac Mar 12, 2024
@lh3 lh3 added bug and removed question labels Mar 12, 2024
@lh3
Copy link
Owner

lh3 commented Mar 12, 2024

Please try paftools.js from github HEAD. Let me know if you still have this problem. Thanks.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Projects
None yet
Development

No branches or pull requests

2 participants