buttery-eel
buttery-eel copied to clipboard
RNA m6a called reads seem not to align?
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?