ITADN

combined quantification output from psirc

#201Closedjpfry327 创建于 2025-05-01
bug
J
jpfry327commented
### Description of the bug Hi, I've run the pipeline with default parameters, and I believe the main execution has worked. However, I have a question about the output in the `4_quantification` folder. When I look at the combined output from `psirc`, I get an empty `circular.tsv` file (in `4_quantification/psirc/combined`). From what I understand, this file is created by finding the circular transcript abundances from the individual samples. In my case, I have individual sample abundances at, for example, `4_quantification/psirc/samples/siCONTROL_REP1/psirc/siCONTROL_REP1`. When I inspect the `abundance.tsv` file at this location, I see abundances for transcripts such as "chr1:2191529-2192844:-", which I'm assuming are the circRNA transcripts. So the pipeline has quantified circRNAs for my individual samples, but it looks like they have not been combined. Looking through the code, it looks like you use `awk` to extract transcripts starting with "circ_". But it looks like none of the transcripts start with "circ_", but are of the form "chr1:2191529-2192844:-"? Anyway, regardless of the empty output at `4_quantification/psirc/combined`, can I just manually combine all the files at `4_quantification/psirc/samples/*/psirc/*/abundance.tsv`? Or does the combination step do any additional processing? ### Command used and terminal output ```console nextflow run nf-core/circrna -r dev \ -profile nu_genomics \ -resume \ --outdir nextflow_dir \ --input samplesheet.csv \ --gtf ~/ylab/jfp0082/ref/gencode.v42.annotation.gtf \ --fasta ~/ylab/jfp0082/ref/hg38.fa \ --bowtie2 "/home/jfp0082/ylab/jfp0082/ref/bowtie2" \ --igenomes_ignore \ --genome null ``` ### Relevant files _No response_ ### System information _No response_
关闭于 2025-07-21 1 条评论