Command line: [exonerate -m p2g --showtargetgff -q ./SPP00000360_2.0.fa -t ./SPP00000360_2.0/fastasubseq/KV389684.1_fastasubseq.fa]
Hostname: [ws101092U]
C4 Alignment:
------------
Query: SPP00000360_2.0 # Protein # Selenoprotein M (SELENOM) # Marmoset
Target: KV389684.1:subseq(1556918,100428) Cebus capucinus imitator isolate Cc_AM_T3 unplaced genomic scaffold Scaffold258, whole genome shotgun sequence
Model: protein2genome:local
Raw score: 678
Query range: 0 -> 141
Target range: 47951 -> 50428
1 : LeuProProLeuProLeuLeuLeuLeuLeuAlaAlaLeuAlaAlaProAlaThrAlaAla : 20
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
LeuProProLeuProLeuLeuLeuLeuLeuAlaAlaLeuAlaAlaProAlaThrAlaAla
47952 : TTGCCTccgctgccgctgctgctgcttctcGCGGCGCTTGCGGCTCCAGCCACAGCCGCC : 48009
21 : ThrThrTyrArgProAspTrpAsnArgLeuArgGlyLeuThrArgAlaArgValGlu > : 40
|||||||||||||||||||||||||||||||||||||||||||||||||||||||||
ThrThrTyrArgProAspTrpAsnArgLeuArgGlyLeuThrArgAlaArgValGlu++
48010 : ACTACCTATCGGCCGGACTGGAACCGTCTGAGAGGCCTGACCCGCGCCCGGGTAGAGgt. : 48071
41 : >>> Target Intron 1 >>>> ThrCysGlyGlyUnkGlnLeuAsnArgLeuLysG : 51
1403 bp |||||||||||| |||||||||||||||||||
++ThrCysGlyGly***GlnLeuAsnArgLeuLysG
48072 : ........................agACCTGTGGGGGATGACAGCTGAACCGCCTAAAGG : 49505
52 : lu >>>> Target Intron 2 >>>> ValLysAlaPheValThrGlnAspIlePr : 61
|| 216 bp |||||||||||||||||||||||||||||
lu++ ++ValLysAlaPheValThrGlnAspIlePr
49506 : AGgt.........................agGTGAAGGCTTTCGTCACTCAGGACATTCC : 49751
62 : oPhe{Ty} >>>> Target Intron 3 >>>> {r}HisAsnLeuValMetLysHi : 70
||||{||} 357 bp {|}||||||||||||||||||||
oPhe{Ty}++ ++{r}HisAsnLeuValMetLysHi
49752 : ATTC{TA}gt.........................ag{T}CACAACCTGGTGATGAAACA : 50135
71 : sLeuProGlyAlaAspProGluLeuValLeuLeuGlyArgArgTyrGluGluLeuGlu : 90
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
sLeuProGlyAlaAspProGluLeuValLeuLeuGlyArgArgTyrGluGluLeuGlu++
50136 : CCTCCCTGGGGCCGACCCCGAGCTCGTGCTGCTGGGCCGCCGCTACGAGGAACTGGAGgt : 50195
91 : >>>> Target Intron 4 >>>> ArgIleProLeuSerGluMetThrArgGluGlu : 100
78 bp |||||||||||||||||||||||||||||||||
++ArgIleProLeuSerGluMetThrArgGluGlu
50196 : .........................agCGCATTCCGCTCAGCGAAATGACCCGCGAAGAG : 50303
101 : IleAsnAlaLeuValGlnGluLeuGlyPheTyrArgLysAlaAlaProAspAlaGlnVal : 120
|||||||||||||||||||||||||||||||||||||||!.!||||||||||||||||||
IleAsnAlaLeuValGlnGluLeuGlyPheTyrArgLysValAlaProAspAlaGlnVal
50304 : ATCAACGCGCTGGTGCAGGAGCTCGGCTTCTACCGCAAGGTGGCGCCCGACGCGCAGGTG : 50363
121 : ProProGluTyrValTrpAlaProAlaArgProProValGluThrSerAspHisAlaAsp : 140
||||||||||||||||||||||||||||||||||||! !|||||||||||||||.!!|||
ProProGluTyrValTrpAlaProAlaArgProProGluGluThrSerAspHisThrAsp
50364 : CCCCCCGAGTACGTGTGGGCGCCCGCTAGGCCCCCGGAGGAAACTTCAGACCACACTGAC : 50423
141 : Leu : 141
|||
Leu
50424 : CTG : 50428
vulgar: SPP00000360_2.0 0 141 . KV389684.1:subseq(1556918,100428) 47951 50428 + 678 M 39 117 5 0 2 I 0 1399 3 0 2 M 12 36 5 0 2 I 0 212 3 0 2 M 11 33 S 0 2 5 0 2 I 0 353 3 0 2 S 1 1 M 26 78 5 0 2 I 0 74 3 0 2 M 52 156
# --- START OF GFF DUMP ---
#
#
##gff-version 2
##source-version exonerate:protein2genome:local 2.2.0
##date 2016-11-18
##type DNA
#
#
# seqname source feature start end score strand frame attributes
#
KV389684.1:subseq(1556918,100428) exonerate:protein2genome:local gene 47952 50428 678 + . gene_id 1 ; sequence SPP00000360_2.0 ; gene_orientation +
KV389684.1:subseq(1556918,100428) exonerate:protein2genome:local cds 47952 48068 . + .
KV389684.1:subseq(1556918,100428) exonerate:protein2genome:local exon 47952 48068 . + . insertions 0 ; deletions 0
KV389684.1:subseq(1556918,100428) exonerate:protein2genome:local splice5 48069 48070 . + . intron_id 1 ; splice_site "GT"
KV389684.1:subseq(1556918,100428) exonerate:protein2genome:local intron 48069 49471 . + . intron_id 1
KV389684.1:subseq(1556918,100428) exonerate:protein2genome:local splice3 49470 49471 . + . intron_id 0 ; splice_site "AG"
KV389684.1:subseq(1556918,100428) exonerate:protein2genome:local cds 49472 49507 . + .
KV389684.1:subseq(1556918,100428) exonerate:protein2genome:local exon 49472 49507 . + . insertions 0 ; deletions 0
KV389684.1:subseq(1556918,100428) exonerate:protein2genome:local splice5 49508 49509 . + . intron_id 2 ; splice_site "GT"
KV389684.1:subseq(1556918,100428) exonerate:protein2genome:local intron 49508 49723 . + . intron_id 2
KV389684.1:subseq(1556918,100428) exonerate:protein2genome:local splice3 49722 49723 . + . intron_id 1 ; splice_site "AG"
KV389684.1:subseq(1556918,100428) exonerate:protein2genome:local cds 49724 49758 . + .
KV389684.1:subseq(1556918,100428) exonerate:protein2genome:local exon 49724 49758 . + . insertions 0 ; deletions 0
KV389684.1:subseq(1556918,100428) exonerate:protein2genome:local splice5 49759 49760 . + . intron_id 3 ; splice_site "GT"
KV389684.1:subseq(1556918,100428) exonerate:protein2genome:local intron 49759 50115 . + . intron_id 3
KV389684.1:subseq(1556918,100428) exonerate:protein2genome:local splice3 50114 50115 . + . intron_id 2 ; splice_site "AG"
KV389684.1:subseq(1556918,100428) exonerate:protein2genome:local cds 50116 50194 . + .
KV389684.1:subseq(1556918,100428) exonerate:protein2genome:local exon 50116 50194 . + . insertions 0 ; deletions 0
KV389684.1:subseq(1556918,100428) exonerate:protein2genome:local splice5 50195 50196 . + . intron_id 4 ; splice_site "GT"
KV389684.1:subseq(1556918,100428) exonerate:protein2genome:local intron 50195 50272 . + . intron_id 4
KV389684.1:subseq(1556918,100428) exonerate:protein2genome:local splice3 50271 50272 . + . intron_id 3 ; splice_site "AG"
KV389684.1:subseq(1556918,100428) exonerate:protein2genome:local cds 50273 50428 . + .
KV389684.1:subseq(1556918,100428) exonerate:protein2genome:local exon 50273 50428 . + . insertions 0 ; deletions 0
KV389684.1:subseq(1556918,100428) exonerate:protein2genome:local similarity 47952 50428 678 + . alignment_id 1 ; Query SPP00000360_2.0 ; Align 47952 1 117 ; Align 49472 40 36 ; Align 49724 52 33 ; Align 50117 64 78 ; Align 50273 90 156
# --- END OF GFF DUMP ---
#
-- completed exonerate analysis