Nice tool! I have noticed some genes (e.g. HNRNPA1) gives an assertionerror, and I have not been able to figure out why.
code used:
./seqcover report --genes HNRNPA1 \
--fasta $fasta samples/*_recal.per-base.d4 \
--background seqcover_out/seqcover_p5.d4 \
-r my_genes_report.html \
--hg19
# Gives:
fatal.nim(49) sysFatal
Error: unhandled exception: transcript.nim(143, 14) `r_off - l_off == o_exon[1] - o_exon[0]` [AssertionError]
Nice tool! I have noticed some genes (e.g. HNRNPA1) gives an assertionerror, and I have not been able to figure out why.
code used: