buttery-eel icon indicating copy to clipboard operation
buttery-eel copied to clipboard

RNA m6a called reads seem not to align?

Open hasindu2008 opened this issue 4 months ago • 0 comments

I tried this rna_rp4_130bps_modbases_m6a_drach_sup model for m6a calling and aligned the basecalled reads to the reference. Seems no alignments. See the commands below:

/install/buttery-eel-0.5.0+dorado7.4.12/scripts/eel --config rna_rp4_130bps_modbases_m6a_drach_sup.cfg --device cuda:all -i /data/hasindu/hasindu2008.git/nanodegrade/PNXRXX240011_reads_500k.blow5 -o tmp/PNXRXX240011_reads_500k.sam --call_mods
samtools fastq -TMM,ML tmp/PNXRXX240011_reads_500k.sam | minimap2 -ax map-ont /genome/gencode.v40.transcripts.fa -uf -t 8 --secondary=no -y - | samtools sort -@20 - > tmp/PNXRXX240011_reads_500_mapped.bam
samtools flagstats PNXRXX240011_reads_500_mapped.bam

Note that using rna_rp4_130bps_sup.cfg is working as intended. So something happens when modcalls are enabled?

hasindu2008 avatar Oct 09 '24 00:10 hasindu2008