Hey everyone, I am trying to identify circular RNAs in my seq data using CIRCquant. I have read and understood the usage and all the files required I have generated. However, my reads are single end reads and the usage and manual for CIRIquant do not provide any option for single end reads input even though in the published paper they mentioned that the program can take in signle end reads. can any one help ? thanks and see usage below for paired end reads
CIRIquant -t 20 --config ./chr1.yml -1 ./fastq1 -2 ./fastq2 --config ./chr1.yml --no-gene -o ./GCF -p GCFloc
•
written
by
jaqx008
1 answer
Hi jaqx, you can still use CIRI2 for identification of circRNAs from single end data, CIRIquant is designed for quantification of circRNAs based on alignment result of paired-end reads, and does not support single end reads yet.
By the way, our manuscript is still under review, I think what you mentioned is another tool CIRI-full, which is designed to detect internal splicing patterns of circRNAs.
Log in to answer this question.
-
CIRIquant: ValueError: file has no sequences defined (mode='rb') - is it SAM/BAM format? Consider o…
written by Atul K. •
-
How to calculate insert size and Stdev of insert file for Configure file of MaSuRCA assembly ?
written by Sony Nguyen •
-
mapping reads using kallisto - rna seq analysis
written by bioinformatics •
-
Running STAR-SEQR with star outputs
written by avelarbio46 •
-
CWL - glob an array of output files based on a single input string
written by drkennetz
-
Tophat Usage difficulty?
written by jaqx008
-
HiSAT2 edit distance option
written by gokberk
-
Pindel Config File Error
written by pwiner •
-
How to print stats from Python filtering script
written by arabhorsegal18 •
-
We Need Advice On Choosing The Appropriate Parameters For Smallrna Sequencing
written by sethugunja •
If this is the published paper (https://genomemedicine.biomedcentral.com/articles/10.1186/s13073-019-0614-1), where do they mention that the program can take in single end reads?
Given that the method is based on (reverse) overlap detection of paired-end reads
feeding it single-end reads does not really make sense to me...
Edit: citation
this is the publication and they mentioned it here
Hmm, have you tried to only use the -1 option?
Yes it doesnt work. at this point I just assume its not possible. Thanks for trying to help.
Just judging from the first couple lines of code at the CIRIquant github repo, the tool will terminate if it doesn't receive -1 and -2, so yea, not without modifying that (which might be easy enough, though!).
Yes I agree. I tested with a paired end it ran just fine.