meth test

hasindu2008 · Jul 18, 2024 · 86e60c5 · 86e60c5
1 parent d7a3784
commit 86e60c5
Show file tree

Hide file tree

Showing 3 changed files with 22 additions and 0 deletions.
diff --git a/scripts/test.sh b/scripts/test.sh
@@ -127,6 +127,10 @@ diff -q test/ont_friendly.exp a.slow5 || die "diff failed"
 ex ./squigulator -x dna-r10-min -o a.slow5 -n 1 --seed 1 -t1 test/rnasequin_sequences_2.4.fa --digitisation 4096 --sample-rate 10000 --range 300 --offset-mean -1000 --offset-std 0 --median-before-mean 100 --median-before-std 0
 diff -q test/dev.exp a.slow5 || die "diff failed"
 
+#meth r9
+ex ./squigulator -x dna-r9-prom -o a.slow5 --seed 1 -t1 -n 2 -r 29000 test/nCoV-2019.reference.fasta --meth-freq test/mfreq.tsv
+diff -q test/r9_mfreq.exp a.slow5  || die "diff failed"
+
 # threads and batch size
 redundancy_check () {
     N=$(grep -v ^[@#] a.slow5 | cut -f ${1}  | sort | uniq -c | sort -nr -k1,1 | head -1 | awk '{print $1}')

diff --git a/test/mfreq.tsv b/test/mfreq.tsv
@@ -0,0 +1,5 @@
+MN908947.3	70	1.0
+MN908947.3	2093	0.5
+MN908947.3	20132	0.8
+MN908947.3	27999	0.2
+MN908947.3	10231	0.0
diff --git a/test/r9_mfreq.exp b/test/r9_mfreq.exp