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

(Uncaught exception of type std::runtime_error: Illumina read is too long, increase fragment length) #2516

Open
mariaelf97 opened this issue Apr 11, 2024 · 1 comment

Comments

@mariaelf97
Copy link

mariaelf97 commented Apr 11, 2024

Hello, I am trying to simulate reads from a multi fasta file containing a config with length about 1000 bps but I get the following error.
Could you explain what the issue is?

mason_simulator -ir ~/GCA_014899445.1_amplicon_0.fasta -n 20000 -o read_1.fq -or read_2.fq --fragment-mean-size 450 --fragment-size-std-dev 50 --illumina-read-length 500 --fragment-min-size 400 --fragment-max-size 500

error:

WARNING: Tried 1000 times to sample fragment from contig of size 1019. Giving up.
(Uncaught exception of type std::runtime_error: Illumina read is too long, increase fragment length)

This is the fasta file for your reference
`

CP041800.1_amplicon_0
CCGATCCCGTGCGCAGAAACTGCAGCATCATCGCCGGCAACTCCCAGTACGCCCGCAGTGCGGCAATCTGCCACTCGCCGGTAACCGGTCGTAGGTCATAGCGTAGGAAGGCGGGAATGAACACCGTCACAGCCGAGTCCATCGCGACCTCGAGTTCGAGATCGCGCAGCACCACCGTGCCGGAGACGATATCCAGATCGCGATGGAACGTGATATCCCGCGGCCCGATGAAGGTGTCGTAGAAGCGGCCGATGGCCTCATGCCCCACCTGCGGCTGCGAACCCACCGGGTCTTCGACCCGCGCGTCACCGGTGAACAACCCGACCCAGCCGGCGCGGTCGTGCGCGGCGGCCGCTTGCGGCGAGCGCTCCACCGCCGCCAACAGTTCATCCCGGTTCGGCGGTGCCATCAGGAGCTGCAAACCAACTCGACGCTGGCGGTGCGCATCTCCTCCAGCGCGGCGACGGTGGTATCGGCCGACACACCCGCTGTCAGGTCCACCAGCACCCTGGTGGCCAAGCCATTGCGTACCGCGTCCTCGGCCGTCTGGCGCACACAATGATCGGTGGCAATACCGACCACATCGACCTCATCGACGCCGCGTTGCCGCAGCCAATTCAGCAGTGGCGTGCCGTTCTCGTCGACTCCTTCGAAGCCGCTGTACGCTCCGGTGTAGGCACCCTTGTAGAACACCGCCTCGATTGCCGACGTGTCCAGACTGGGATGGAAGTCCGCGCCGGGAGTACCGCTGACGCAATGCGGTGGCCACGACGAGGAATAGTCCGGTGTGCCGGAGAAGTGGTCACCCGGGTCGATGTGGAAGTCCTTGGTTGCCACGACGTGATGGTAGTCCGCCGCTTCGGCCAGGTAGTCGCTGATGGCGCGGGCCAGCGCGGCGCCACCGGTTACCGCCAGCGAGCCACCCTCGCAGAAGTCGTTCTGCACGTCGACGATGATCAACGCCCGCATACGTCCACCATACGTTCGGGCGACTGCCCGGGCAGTTTGCCTACCGACGC

`

@smehringer
Copy link
Member

Hi @mariaelf97,

thanks for reaching out. The error message basically tells you whats going wrong.

You are trying to simulate paired-end sequencing data with reads of length 500 but your fragment size is <= 500.

See this:
image
(Source https://images.app.goo.gl/Zm4aDLFdZGdhaSC98)

Or this:
image
(Source https://images.app.goo.gl/qZrQ8njVfYhiuvcb7)

Not sure myself if mason_simulator is including Adaptor length but you should still get the idea.

You need to increase your fragment size or lower the read length.

Best,
Svenja

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

No branches or pull requests

2 participants