Minzhe / Clirc

Clirc: a user-friendly bioinformatics software to identify RBP-bound circRNAs through analysis of CLIP-Seq data
1 stars 0 forks source link

Problem with extracting circRNA sequences #1

Open MSajek opened 3 years ago

MSajek commented 3 years ago

I try to use Clirc and I have a problem in the first step. When I try to extract circRNA sequences from my file with coordinates and reference genome that were used to detect circRNA i have the following error:

perl bin/Clirc_library.pl -coord /media/marcin/data/KasiasDataRean/RBP/edger/Clirc/DGEcircs.csv -genome /media/marcin/data/MartaAndAgaAnalysis/Genome/GRCh37_primary_assembly.fa -library /media/marcin/data/KasiasDataRean/RBP/edger/Clirc/lib.txt

I used CIRI2 do detect circRNAs, circMeta package for annotation and DGE analysis and GRCh37_primary_assembly human genome from gencode. I will be grateful for any help. Best, Marcin

MSajek commented 3 years ago

I found the reason. In the original genomic file chromosome names looks like:

chr1 1 chr10 10 chr11 11 chr12 12 chr13 13 chr14 14 chr15 15 chr16 16 chr17 17 chr18 18 chr19 19 chr2 2 chr20 20 chr21 21 chr22 22 chr3 3 chr4 4 chr5 5 chr6 6 chr7 7 chr8 8 chr9 9 chrM MT chrX X chrY Y When I change them to: chr1 chr10 .... everythink starts working