This is a test version of Biostars. For the public version, visit https://www.biostars.org.
MIRA: Assembly from fasta file with identical names

I have a fasta file that includes sequences with identical names. I need to assemble contigs from these sequences using MIRA. However, Mira seems to have a problem with identical sequence names. Is there a way to solve this problem?

mira assembly fasta
• 2,301 views • link
written by Tamandua • Tamandua

Are these sequences (with identical names) identical? if so, you can use fastx_collapser from fastx-toolkit or rmdup from seqkit. If not, append an unique string to identical ids.

• 0 views • link
written by cpad0112 cpad0112

No, they are not identical. I used 

blastn -db reference.fasta -query input.fasta -evalue 1e-60 -outfmt '6 qseqid qseq'     | awk 'BEGIN{FS="\t"; OFS="\n"}{gsub(/-/, "", $2); print ">"$1,$2}'     > output.fastq

to get the best matching hits. Maybe there is a better way of doing this? And how can I rename identical ids?

• 0 views • link
written by Tamandua • Tamandua

example fasta with identical headers

$ cat test.fa
>gene1
ATGCGGG
>gene1
TAGCTGT

Rename identical fasta ids:

$ seqkit rename test.fa

output: 

>gene1
ATGCGGG
>gene1_2 gene1
TAGCTGT

download seqkit from here. Output can be redirected to either fa or fa.gz with -o option.

• 0 views • link
written by cpad0112 cpad0112

Just rename all sequences - just google rename fasta header. There are dozens of ways, see two at this post. With awk:

awk '/^>/ { printf("%s_%s\n",$0,i++);next;} { print $0;}' teste.fas > out.fas

or with R:

library(Biostrings)
fa = read.DNAStringSet(...)
names(fa) = make.unique(names(fa))
write.XStringSet(fa, ...)
• 0 views • link
written by h.mon h.mon

Ah, thanks a lot! That's it!

• 0 views • link
written by Tamandua • Tamandua

0 answers

No answers yet.

Log in to answer this question.

More posts like this