Models for newer versions of Guppy with sup basecalls

[hasindu2008]

trafficstars

HI NanoSim developers,

I am wondering if you have a pre-trained model for simulating genomic reads from the human genomes close to ones basecalled from newer versions of Guppy (e.g. Guppy 6) under the super accuracy mode (R9.4.1 chemistry). Or will the default -c human option already suffice to emulate this kind of data?

[kmnip]

The existing models in NanoSim are very old. If you have an ONT dataset with the right basecaller+chemistry combination, then you can subsample to 1 million reads and run the characterization stage. It should be relatively straightforward, e.g.

# subsample 1 million reads
seqtk sample reads.fq.gz 1000000 > reads_sub1M.fq.gz

# NanoSim characterization stage
read_analysis.py genome -i reads_sub1M.fq.gz ...

[hasindu2008]

Can we use NA12878 data to train the human genome and then use the trained model as a generic human model? To explain further, say I create a custom human genome reference file by manually incorporating some variants, can the model trained with NA12878 be used on such a reference? Or do I need to train on data samples exactly with those variants?

[saberhq]

Can we use NA12878 data to train the human genome and then use the trained model as a generic human model? To explain further, say I create a custom human genome reference file by manually incorporating some variants, can the model trained with NA12878 be used on such a reference? Or do I need to train on data samples exactly with those variants?

I am not sure if I got your question correctly. But if you are interested in simulating nanopore reads having the characteristics of the latest basecaller and chemistry, I would suggest you train NanoSim using such reads and then use the simulator to generate reads given a reference genome/transcriptome.

Btw, if you are interested in newer basecaller, dorado, we provide a trained model for H. sapiens NA24385 - AshkenazimTrio - Son (hg002) which is sequenced using Kit v14 (R10 chemistry) and basecalled by dorado. You may find the trained model on 1 Million subsampled reads on the GitHub page (available along with the other models at pre-trained_models folder). If you are interested in the trained model based on the whole dataset, you can get it through Zenodo - DOI: 10.5281/zenodo.10064740. The model is trained using NanoSim v3.0.2 with scikit-learn v0.23.2 and python v3.7.10.

If you have any issues using the pre-trained models, check the dependencies section for some information and tips.