search

HadrienG / InSilicoSeq

:rocket: A sequencing simulator
https://insilicoseq.readthedocs.io
MIT License
176 stars 32 forks source link

Seems like the error free profile performs worse than the error profile #200

Closed Kumereng closed 3 years ago

Kumereng commented 3 years ago

When simulating from bacterial genomes, it appears that when no error profile is used (--model basic), it performs worse after classification on kraken2 than when a Miseq profile is used for instance.

HadrienG commented 3 years ago

--model basic is not an error-free error model. All bases have Q30 and all nucleotides have equal substitution rate. For error-free reads, you can use --model perfect.