Can Exonerate return only the Cigar string?
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
• 1,363 views
•
link
updated
by
Ram
•
written
by
helloword •
•
written
by
helloword •
1 answer
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
Log in to answer this question.
More posts like this
-
Invalid CIGAR after using bam clipOverlap
written by MboiTui •
-
How to properly mock a (Pysam) read
written by ManuelDB
-
How can Biostrings calculate consensus sequence without counting gaps?
written by MatthewP
-
How to convert amino acid sequence to numbers
written by nrizzolo1796 •
-
Exonerate - issues when running it
written by helloword •
-
Adding information tags with types, positions and counts of matches and mismatches in every read in…
written by lechu •
-
New Column based on another data frame
written by j.lunger18 •
-
combine columns from multiple files
written by bk11 •
-
needleman-wunsch algorithm on python
written by elisheva
-
Deciphering Cigar
written by Leszek
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.