Can Exonerate return only the Cigar string?
1
1
Entering edit mode
7 weeks ago
helloword • 0

Hello. Is there any possibility that Exonerate, the sequence alignment tool, returns only cigar string? I mean how could we print only the string without printing unnecessary stuff like below (or without parsing output)?

Command line:

[exonerate --model affine:global --exhaustive yes tube.fasta calm.fasta --showcigar yes --showalignment no]
Hostname: [beep]

(process:30145): WARNING **: 01:19:52.429: Exhaustively generating suboptimal alignments will be VERY SLOW
cigar: Query 0 100 + Target 0 100 + 266  M 13 D 3 M 13 I 2 M 5 I 1 M 66
vulgar: Query 0 100 + Target 0 100 + 266 M 13 13 G 0 3 M 13 13 G 2 0 M 5 5 G 1 0 M 66 66
-- completed exonerate analysis

Wanted output:

M 13 D 3 M 13 I 2 M 5 I 1 M 66
Exonerate • 312 views
ADD COMMENT
0
Entering edit mode

Please do not delete posts that have received feedback. If your post was answered and the answer worked, accept it. If it did not work, provide feedback on what went wrong.

ADD REPLY
1
Entering edit mode
7 weeks ago
seidel 8.3k

Have you tried the roll your own (--ryo) output format (see the manual)? There are detailed options.

--showalignment no --showvulgar no --ryo "%C\n"

 M 5 I 2 M 3 D 2 M 10 I 1 M 15 I 1 M 10 I 1 M 18 D 1 M 16 I 3 M 7 I 1 M 15 D 2 M 5 I 1 M 10 I 1 M 38 D 4 M 47 D 1 M 2 D 1 M 7 D 2 M 22 I 2 M 4 I 2 M 18 I 2 M 2 I 1 M 24 I 1 M 4 I 1 M 21 I 2 M 11 I 1 M 25 D 1 M 9 D 1 M 12 I 1 M 20
ADD COMMENT

Login before adding your answer.

Traffic: 1620 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6