RetroSeq : Fatal error : unrecognised symbol

Hi everybody!

I don't know if you use RetroSeq to detect transposable element but I have an issue with it :

** FATAL ERROR **: Unrecognised symbol 'x' (ascii:120) file:[17846.allrefs.fasta] seq:[Endogenous_retrovirus] pos:[204]
exiting ...

I checked the file (which is the file with the TE sequences) and saw that in fact there are some x, k, etc. in it. How can I deal with it?

Thank you for your help!

Lou

te ascii

• 2,221 views • link

updated 3.1 years ago by Ram

• written 9.1 years ago by FryHandiz •

1 answer

try to replace the non-ATGCN bases to N ?

 sed '/^[^>]/s/[^ATGCNatgcn]/N/g' in.fasta > out.fasta

• 0 views • link

written 9.1 years ago by Pierre Lindenbaum

Huhu... I'm stupid! I didn't think about that, I try it! Thank you.

• 0 views • link

written 9.1 years ago by FryHandiz •

More posts like this

AMRFinderPlus input fasta file for clinical isolates
written 1.0 years ago by Adriana •

I recently ran NCBI's AMRFinderPlus to detect antimicrobial resistance on two different datasets of clinical isolates of *Mycobacterium tuberculosis*. Since the input file requires is …
VEP doesn't count CADD
written 1.5 years ago by balanannie •

I've been stuck with counting CADD scores for a while now. I'm working on server with VEP, it annotates my vcf, and uses CADD.pm plugin …
GFF Error: overlapping duplicate dispersed_repeat feature in stringtie
written 2.7 years ago by Aki •

Hi. I got following error when I use stringtie. with repeatmasker annotation gff file and RNA-seq bam files which is already sorted with samtools. GFF …
STARsolo
written 3.4 years ago by Nodilan •

Good morning everybody, I need you expertise : In fact I'm trying to write a script using python object oriented programming to perform STARsolo to …
No CR tag in the sam file mapped by STARsolo (10x v3.1 data)
written 3.6 years ago by shaoqian ma •

I mapped the 10x data (v3.1) using STARsolo but the strange thing is that output sam file does not have the tags including CR CY …
HISAT2- (ERR): hisat2-align exited with value 134
written 5.9 years ago by Chvatil

Hello I'm using **hisat2 version 2-2.1.0** I use hisat2 in order to map fastq reads into my assembly.fa file. here is the command I use …
How to remove reads with only Ns from bam/sam/cram file?
written 6.1 years ago by erik.stricker •

When using my current sam file (which I generated from a cram file, which in turn was generated from a bam file) in my analysis …
repeatmasker empty output
written 7.3 years ago by Chironex •

hi everybody, finally I run Repeatmasker but after all, it generates a table like this: file name: OB_100DEC.fa sequences: 100 total length: 209466439 bp (184235452 …
SNP plotting from VCF file in (non database) organism
written 11.8 years ago by justinfanyi •

Hello everyone, I have called SNP and annotated;results in vcf format, I have been trying to work out visualization of this SNP using various tools …
Exonerate ** FATAL ERROR **: Unrecognised symbol 'J' (ascii:74)
written 12.2 years ago by jomaco

Hi, I am receiving the following error when running exonerate-2.2.0: ** FATAL ERROR **: Unrecognised symbol 'J' (ascii:74) Here are a couple of examples of …