From 877f6aa573f176222633d5ef9b118aba5e827688 Mon Sep 17 00:00:00 2001 From: Hasindu Gamaarachchi Date: Mon, 19 Feb 2024 22:59:12 +1100 Subject: [PATCH] rna004 test --- test/test_extensive.sh | 12 ++++++++++-- 1 file changed, 10 insertions(+), 2 deletions(-) diff --git a/test/test_extensive.sh b/test/test_extensive.sh index 7dd32b9..d26b0cf 100755 --- a/test/test_extensive.sh +++ b/test/test_extensive.sh @@ -63,7 +63,7 @@ echo "R10 DNA ideal-time" eel -i new.blow5 --config dna_r10.4.1_e8.2_400bps_sup.cfg --device cuda:all -o new.fastq &>> a.log || die "eel failed" identitydna.sh $REF_HG38_IDX new.fastq > a.acc 2>> a.log || die "identitydna failed" cat a.acc -CHECK_ACC 0.89 a.acc +CHECK_ACC 0.91 a.acc REMOVE_TMP echo "R10 DNA" @@ -71,7 +71,15 @@ echo "R10 DNA" eel -i new.blow5 --config dna_r10.4.1_e8.2_400bps_sup.cfg --device cuda:all -o new.fastq &>> a.log || die "eel failed" identitydna.sh $REF_HG38_IDX new.fastq > a.acc 2>> a.log || die "identitydna failed" cat a.acc -CHECK_ACC 0.85 a.acc +CHECK_ACC 0.90 a.acc samtools sort -o new.bam new.sam || die "samtools failed" samtools index new.bam || die "samtools failed" REMOVE_TMP + +echo "RNA004 RNA" +./squigulator -x rna004-prom $REF_GENCODE -o new.blow5 2> a.log || die "squigulator failed" +/install/buttery-eel-0.4.2+dorado7.2.13/scripts/eel -i new.blow5 --config rna_rp4_130bps_sup.cfg --device cuda:all -o new.fastq &>> a.log|| die "eel failed" +identityrna.sh $REF_GENCODE new.fastq > a.acc 2>> a.log || die "identitydna failed" +cat a.acc +CHECK_ACC 0.77 a.acc +REMOVE_TMP \ No newline at end of file