xqtl-protocol
xqtl-protocol copied to clipboard
RNA_calling.ipynb rsem_call error when using BAM files from STAR with wasp filter
Running alignment with the STAR wasp option creates *toTranscriptome.out_wasp_qc.bam files. When running RSEM, the code expect the files to not include "_wasp_qc" in the name, causing an error.
@grennfp is this fixed now? It looks like a file naming issue that can be addressed with a bit of SoS/python syntax?