nanoporetech / megalodon

Megalodon is a research command line tool to extract high accuracy modified base and sequence variant calls from raw nanopore reads by anchoring the information rich basecalling neural network output to a reference genome/transriptome.
Other
197 stars 30 forks source link

RNA variation models #168

Open franrodalg opened 3 years ago

franrodalg commented 3 years ago

Hi @marcus1487 !

I'm aware that the models for modified RNA bases are still in preparation, but is there any way of calling variants from direct RNA data? There doesn't seem to be associated calibration files at the moment:

ERROR: No default sequence variant calibration file found for guppy config: rna_r9.4.1_70bps_hac.cfg

Cheers, Fran

marcus1487 commented 3 years ago

This is not currently implemented. There is a bit of logic that needs to be filled in to enable variant calling on RNA/3’ to 5’ signal direction reads. This is missing at the moment. This feature is of interest though, so we will look into adding this feature in a future megalodon release.

franrodalg commented 3 years ago

Thanks, Marcus. It would indeed be really really useful if this could be incorporated into future releases.