Command line: [exonerate -m p2g --showtargetgff -q ./SPP00000358_2.0.fa -t ./SPP00000358_2.0/fastasubseq/KV389679.1_fastasubseq.fa]
Hostname: [ws101092U]
C4 Alignment:
------------
Query: SPP00000358_2.0 # Protein # Selenoprotein H (SELENOH) # Marmoset
Target: KV389679.1:subseq(2553092,100146) Cebus capucinus imitator isolate Cc_AM_T3 unplaced genomic scaffold Scaffold253, whole genome shotgun sequence:[revcomp]
Model: protein2genome:local
Raw score: 581
Query range: 0 -> 122
Target range: 50354 -> 49749
1 : MetAlaProArgGlyArgLysArgLysAlaGluAlaAlaMetValAlaThrThrGluLys : 20
|||||| !!||||||||||||||||||||||||||||||:!!||||||||||||||||||
MetAlaSerArgGlyArgLysArgLysAlaGluAlaAlaValValAlaThrThrGluLys
50354 : ATGGCTTCCCGCGGGAGGAAGCGGAAGGCTGAGGCCGCGGTGGTCGCCACGACCGAGAAG : 50297
21 : ArgGluLysLeuAlaAsnGlyGlyGluGlyValGluGluAlaThrIleValIleGluHis : 40
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
ArgGluLysLeuAlaAsnGlyGlyGluGlyValGluGluAlaThrIleValIleGluHis
50296 : CGAGAGAAGCTGGCGAACGGCGGGGAGGGAGTGGAGGAGGCGACCATTGTTATTGAGCAT : 50237
41 : {Cy} >>>> Target Intron 1 >>>> {s}ThrSerUnkArgValTyrGlyArg : 49
{||} 85 bp {|}|||||| |||||||||||||||
{Cy}++ ++{s}ThrSer***ArgValTyrGlyArg
50236 : {TG}gt.........................ag{C}ACTAGCTGACGCGTTTATGGGCGC : 50125
50 : AsnAlaAlaAlaLeuSerGlnAlaLeuArgLeuGluAlaProGluLeuProValLysVal : 69
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
AsnAlaAlaAlaLeuSerGlnAlaLeuArgLeuGluAlaProGluLeuProValLysVal
50124 : AACGCCGCGGCCCTGAGCCAGGCGCTGCGCCTGGAAGCCCCAGAACTTCCAGTAAAGGTG : 50065
70 : AsnProAlaLysProArgArgGlySerPheGluValThrLeuLeuArgProAspGlySer : 89
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
AsnProAlaLysProArgArgGlySerPheGluValThrLeuLeuArgProAspGlySer
50064 : AACCCAGCCAAGCCCCGGAGGGGCAGCTTCGAGGTGACGCTGTTGCGCCCGGACGGCAGC : 50005
90 : {S} >>>> Target Intron 2 >>>> {er}AlaGluLeuTrpThrGlyIleLys : 98
{|} 154 bp {||}||||||||||||||||||||||||
{S}++ ++{er}AlaGluLeuTrpThrGlyIleLys
50004 : {A}gt.........................ag{GT}GCGGAGCTCTGGACTGGGATTAAG : 49824
99 : LysGlyProProArgLysLeuLysPheProGluProGlnGluValValGluGluLeuLys : 118
||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
LysGlyProProArgLysLeuLysPheProGluProGlnGluValValGluGluLeuLys
49823 : AAGGGGCCCCCGCGCAAACTCAAATTCCCTGAGCCTCAAGAGGTGGTCGAAGAGTTGAAG : 49764
119 : LysTyrLeuSer : 122
||||||||||||
LysTyrLeuSer
49763 : AAGTACCTGTCG : 49750
vulgar: SPP00000358_2.0 0 122 . KV389679.1:subseq(2553092,100146) 50354 49749 - 581 M 40 120 S 0 2 5 0 2 I 0 81 3 0 2 S 1 1 M 48 144 S 0 1 5 0 2 I 0 150 3 0 2 S 1 2 M 32 96
# --- 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
#
KV389679.1:subseq(2553092,100146) exonerate:protein2genome:local gene 49750 50354 581 - . gene_id 1 ; sequence SPP00000358_2.0 ; gene_orientation +
KV389679.1:subseq(2553092,100146) exonerate:protein2genome:local cds 50233 50354 . - .
KV389679.1:subseq(2553092,100146) exonerate:protein2genome:local exon 50233 50354 . - . insertions 0 ; deletions 0
KV389679.1:subseq(2553092,100146) exonerate:protein2genome:local splice5 50231 50232 . - . intron_id 1 ; splice_site "GT"
KV389679.1:subseq(2553092,100146) exonerate:protein2genome:local intron 50148 50232 . - . intron_id 1
KV389679.1:subseq(2553092,100146) exonerate:protein2genome:local splice3 50148 50149 . - . intron_id 0 ; splice_site "AG"
KV389679.1:subseq(2553092,100146) exonerate:protein2genome:local cds 50002 50147 . - .
KV389679.1:subseq(2553092,100146) exonerate:protein2genome:local exon 50002 50147 . - . insertions 0 ; deletions 0
KV389679.1:subseq(2553092,100146) exonerate:protein2genome:local splice5 50000 50001 . - . intron_id 2 ; splice_site "GT"
KV389679.1:subseq(2553092,100146) exonerate:protein2genome:local intron 49848 50001 . - . intron_id 2
KV389679.1:subseq(2553092,100146) exonerate:protein2genome:local splice3 49848 49849 . - . intron_id 1 ; splice_site "AG"
KV389679.1:subseq(2553092,100146) exonerate:protein2genome:local cds 49750 49847 . - .
KV389679.1:subseq(2553092,100146) exonerate:protein2genome:local exon 49750 49847 . - . insertions 0 ; deletions 0
KV389679.1:subseq(2553092,100146) exonerate:protein2genome:local similarity 49750 50354 581 - . alignment_id 1 ; Query SPP00000358_2.0 ; Align 50355 1 120 ; Align 50147 42 144 ; Align 49846 91 96
# --- END OF GFF DUMP ---
#
-- completed exonerate analysis