Command line: [exonerate -m p2g --showtargetgff -q ./SPP00000028_1.0.fa -t ./SPP00000028_1.0/fastasubseq/KV389571.1_fastasubseq.fa]
Hostname: [ws101072u]
C4 Alignment:
------------
Query: SPP00000028_1.0 # Protein # Selenoprotein U3 (SelU3) # Homo sapiens # Complete
Target: KV389571.1:subseq(4650514,100083) Cebus capucinus imitator isolate Cc_AM_T3 unplaced genomic scaffold Scaffold145, whole genome shotgun sequence
Model: protein2genome:local
Raw score: 526
Query range: 20 -> 189
Target range: 48087 -> 50083
21 : GluAlaValGluLeuArgSerLeuTrpArgGluHisAlaCys--ValValAlaGlyLeuA : 40
:!!||||||! !||| !!||||||||||||||||||!.!|||## ||||||||||||
GlnAlaValAlaLeuTrpSerLeuTrpArgGluHisValCys##---ValAlaGlyLeuT
48088 : CAGGCTGTGGCGCTGTGGAGCCTGTGGCGGGAGCACGTGTGTGG---GTGGCCGGGCTGT : 48144
41 : rgArgPheGlyCysValValCysArgTrpIleAlaGlnAspLeuSerSerLeuAlaGlyL : 60
!!!.!||| !!|||! !||||||!.!||||||||| !|||||||||.!!|||!.!|||!
rpHisPheArgCysGlyValCysHisTrpIleAlaTrpAspLeuSerGlyLeuValGlyP
48145 : GGCACTTCAGGTGCGGGGTGTGCCACTGGATCGCCTGGGACCTCAGTGGCCTCGTGGGGT : 48204
61 : euLeuAspGlnHisGlyVal-ArgLeuValGlyValGlyProGluAlaLeuGlyLeuGln : 79
!!|||!!:||||||||||||# |||||||||:!!||||||||||||||||||||||||
heLeuGluGlnHisGlyVal#---LeuValGlyLeuGlyProGluAlaLeuGlyLeuGln
48205 : TCCTGGAGCAGCACGGTGTGC---CTGGTGGGCTTGGGGCCGGAGGCCCTGGGTCTGCAG : 48259
80 : GluPheLeuAspGlyAspTyrPheAlaGly{G} >>>> Target Intron 1 >>>> : 90
||||||||| ! |||..!||||||||||||{|} 435 bp
GluPheLeuLysGlySerTyrPheAlaGly{G}++
48260 : GAGTTCCTGAAGGGCAGCTACTTCGCGGGA{G}gt......................... : 48295
91 : {lu}LeuTyrLeuAspGluSerLysGlnLeuTyrLysGluLeuGlyPheLys{Ar} : 107
{||}||||||||||||||||||||||||||||||||||||||||||||||||{||}
++{lu}LeuTyrLeuAspGluSerLysGlnLeuTyrLysGluLeuGlyPheLys{Ar}++
48296 : ag{AG}CTCTACCTGGATGAGAGCAAGCAGCTTTATAAGGAGCTGGGCTTCAAG{CG}gt : 48780
108 : >>>> Target Intron 2 >>>> {g}TyrAsnSerLeuSerIleLeuProAlaAla : 117
579 bp {|}|||!:!||||||||||||||||||||||||
++{g}TyrSerSerLeuSerIleLeuProAlaAla
48781 : .........................ag{G}TACAGCAGCCTGAGCATCCTCCCTGCGGCC : 49387
118 : LeuGlyLysProValArgAspValAlaAlaLys >>>> Target Intron 3 >>>> : 129
||||||||||||||||||||||||:!!|||:!! 88 bp
LeuGlyLysProValArgAspValSerAlaGln++
49388 : CTGGGGAAGCCCGTGCGCGATGTGTCTGCCCAGgt......................... : 49425
130 : AlaLysAlaValGlyIleGlnGlyAsnLeuSer--GlyAspLeuLeuGlnSerGlyGl : 147
|||||||||:!!|||||||||||||||||||||## ||||||||||||||||||||
++AlaLysAlaIleGlyIleGlnGlyAsnLeuSer##---AspLeuLeuGlnSerGlyGl
49426 : agGCCAAGGCTATTGGCATCCAGGGGAACCTGTCTGG---GACCTGCTGCAGAGTGGAGG : 49564
148 : yLeuLeuValValSerLys{G} >>>> Target Intron 4 >>>> {ly}GlyAs : 156
|||||||||||||!!!|||{|} 388 bp {||}|||||
yLeuLeuValValCysLys{G}++ ++{ly}GlyAs
49565 : GCTGCTGGTGGTCTGCAAA{G}gt.........................ag{GT}GGTGA : 49979
157 : p<->LysValLeuLeuHisPheValGlnLysSerProGlyAspTyrValProLysGluHi : 175
| :!!|||||||||! !||||||||||||||| !!|||||||||||||||::!|||||
pThrGluValLeuLeuProPheValGlnLysSerAlaGlyAspTyrValProArgGluHi
49980 : CACAGAAGTGCTGCTGCCTTTCGTCCAGAAGTCTGCAGGCGACTATGTTCCCCGGGAGCA : 50039
176 : sIleLeuGlnValLeuGlyIleSerAlaGluValCysAlaSer : 189
|||||||||||||||||||!!:|||||||||||||||||||||
sIleLeuGlnValLeuGlyMetSerAlaGluValCysAlaSer
50040 : TATCTTGCAGGTCCTGGGCATGTCTGCAGAGGTGTGTGCCAGT : 50083
vulgar: SPP00000028_1.0 20 189 . KV389571.1:subseq(4650514,100083) 48087 50083 + 526 M 14 42 F 0 2 G 1 0 M 31 93 F 0 1 G 1 0 M 22 66 S 0 1 5 0 2 I 0 431 3 0 2 S 1 2 M 16 48 S 0 2 5 0 2 I 0 575 3 0 2 S 1 1 M 21 63 5 0 2 I 0 84 3 0 2 M 11 33 F 0 2 G 1 0 M 13 39 S 0 1 5 0 2 I 0 384 3 0 2 S 1 2 M 2 6 G 0 3 M 33 99
# --- START OF GFF DUMP ---
#
#
##gff-version 2
##source-version exonerate:protein2genome:local 2.2.0
##date 2016-11-16
##type DNA
#
#
# seqname source feature start end score strand frame attributes
#
KV389571.1:subseq(4650514,100083) exonerate:protein2genome:local gene 48088 50083 526 + . gene_id 1 ; sequence SPP00000028_1.0 ; gene_orientation +
KV389571.1:subseq(4650514,100083) exonerate:protein2genome:local cds 48088 48292 . + .
KV389571.1:subseq(4650514,100083) exonerate:protein2genome:local exon 48088 48292 . + . insertions 0 ; deletions 2 ; frameshifts 3
KV389571.1:subseq(4650514,100083) exonerate:protein2genome:local splice5 48293 48294 . + . intron_id 1 ; splice_site "GT"
KV389571.1:subseq(4650514,100083) exonerate:protein2genome:local intron 48293 48727 . + . intron_id 1
KV389571.1:subseq(4650514,100083) exonerate:protein2genome:local splice3 48726 48727 . + . intron_id 0 ; splice_site "AG"
KV389571.1:subseq(4650514,100083) exonerate:protein2genome:local cds 48728 48779 . + .
KV389571.1:subseq(4650514,100083) exonerate:protein2genome:local exon 48728 48779 . + . insertions 0 ; deletions 0
KV389571.1:subseq(4650514,100083) exonerate:protein2genome:local splice5 48780 48781 . + . intron_id 2 ; splice_site "GT"
KV389571.1:subseq(4650514,100083) exonerate:protein2genome:local intron 48780 49358 . + . intron_id 2
KV389571.1:subseq(4650514,100083) exonerate:protein2genome:local splice3 49357 49358 . + . intron_id 1 ; splice_site "AG"
KV389571.1:subseq(4650514,100083) exonerate:protein2genome:local cds 49359 49422 . + .
KV389571.1:subseq(4650514,100083) exonerate:protein2genome:local exon 49359 49422 . + . insertions 0 ; deletions 0
KV389571.1:subseq(4650514,100083) exonerate:protein2genome:local splice5 49423 49424 . + . intron_id 3 ; splice_site "GT"
KV389571.1:subseq(4650514,100083) exonerate:protein2genome:local intron 49423 49510 . + . intron_id 3
KV389571.1:subseq(4650514,100083) exonerate:protein2genome:local splice3 49509 49510 . + . intron_id 2 ; splice_site "AG"
KV389571.1:subseq(4650514,100083) exonerate:protein2genome:local cds 49511 49585 . + .
KV389571.1:subseq(4650514,100083) exonerate:protein2genome:local exon 49511 49585 . + . insertions 0 ; deletions 1 ; frameshifts 2
KV389571.1:subseq(4650514,100083) exonerate:protein2genome:local splice5 49586 49587 . + . intron_id 4 ; splice_site "GT"
KV389571.1:subseq(4650514,100083) exonerate:protein2genome:local intron 49586 49973 . + . intron_id 4
KV389571.1:subseq(4650514,100083) exonerate:protein2genome:local splice3 49972 49973 . + . intron_id 3 ; splice_site "AG"
KV389571.1:subseq(4650514,100083) exonerate:protein2genome:local cds 49974 50083 . + .
KV389571.1:subseq(4650514,100083) exonerate:protein2genome:local exon 49974 50083 . + . insertions 3 ; deletions 0
KV389571.1:subseq(4650514,100083) exonerate:protein2genome:local similarity 48088 50083 526 + . alignment_id 1 ; Query SPP00000028_1.0 ; Align 48088 21 42 ; Align 48132 36 93 ; Align 48226 68 66 ; Align 48730 91 48 ; Align 49360 108 63 ; Align 49511 129 33 ; Align 49546 141 39 ; Align 49976 155 6 ; Align 49985 157 99
# --- END OF GFF DUMP ---
#
-- completed exonerate analysis