Command line: [exonerate -m p2g --showtargetgff -q selenoh_zf.aa.fa.canviat -t genomeselenoh_JXSJ01000003.1.fa] Hostname: [ws101100u] C4 Alignment: ------------ Query: SPP00000634_2.0 # Protein # Selenoprotein H (SELENOH) # Zebrafish Target: JXSJ01000003.1:subseq(9649972,98963) Miichthys miiuy scaffold3, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 315 Query range: 4 -> 127 Target range: 50869 -> 48962 5 : GlyLysSerAlaArgLysArgLysAlaAspSerAspGluLysGluLysLeuAspAspAla : 24 |||::: !!!.!! !||||||||||||!!: !! !||| |||||||||... GlyArgArgGlyThrLysArgLysAlaGluValValGlu---GluLysLeuSerValGlu 50869 : GGTCGTCGTGGGACGAAGCGTAAAGCGGAAGTCGTTGAG---GAGAAACTGTCcgtggag : 50815 25 : LysLysGluLysLeuGluAspLysAspGluGluThrGlyLeuArgValValIleGluHis : 44 ::: ... ||| ||||||:!::!! ! !!! !|||||||||||||||||| GluAspArgGlyGluGluLysLysAspAspGlnGluArgGlnArgValValIleGluHis 50814 : gaggacagaggagaggagaagaaggacgACCAGGAGAGACAAAGGGTGGTCATTGAACAC : 50755 45 : {Cy} >>>> Target Intron 1 >>>> {s}LysSerUnkArgValTyrGlyArg : 53 {||} 764 bp {|}|||||| ||||||||||||||| {Cy}++ ++{s}LysSer***ArgValTyrGlyArg 50754 : {TG}gt.........................ag{T}AAGAGCTGACGAGTGTATGGGCGT : 49964 54 : AsnAlaValValValArgGluAlaLeuAlaAspSerHisProGluLeuLysValMetIle : 73 ||||||! !!.!|||!:!...|||||| !! !:!!||||||! !|||! |||:!::!! AsnAlaGluAlaValLysSerAlaLeuLeuAlaAlaHisProGlyLeuThrValValLeu 49963 : AATGCCGAGGCGGTGAAATCTGCCCTCCTGGCTGCCCACCCCGGACTGACTGTGGTCCTC : 49904 74 : AsnProHisAsnProArgArgAsnSerPheGluIleThrLeuMetAspGlyGluArg{A} : 93 ||||||.!.!!.||||||||||||||||||||||||||||||:!!||||||! !:!{!} AsnProGluLysProArgArgAsnSerPheGluIleThrLeuLeuAspGlyGlyLys{G} 49903 : AACCCTGAAAAACCCCGCCGGAACAGCTTCGAGATCACTCTGCTGGATGGAGGCAAA{G} : 49844 94 : >>>> Target Intron 2 >>>> {la}AspValLeuTrpSerGlyIleLysLys : 102 777 bp { !} ! !||||||!!!|||||||||||| ++ ++{lu}ThrSerLeuTrpThrGlyIleLysLys 49843 : gt.........................ag{AA}ACTTCTCTGTGGACTGGAATAAAGAAG : 49040 103 : GlyProProArgLysLeuLysPheProGluProAlaGluValValThrAlaLeuLysGln : 122 |||||||||||||||||||||||||||:!!|||! !! ||||||.!!||||||:!!:!! GlyProProArgLysLeuLysPheProGlnProAspValValValAlaAlaLeuGlnGlu 49039 : GGTCCGCCTCGCAAGCTGAAATTTCCTCAACCTGATGTTGTAGTTGCTGCTCTACAGGAG : 48980 123 : AlaLeuGluLysGlu : 127 !.!|||:!!!!.||| ValLeuLysAsnGlu 48979 : GTTCTGAAGAATGAG : 48963 vulgar: SPP00000634_2.0 4 127 . JXSJ01000003.1:subseq(9649972,98963) 50869 48962 - 315 M 13 39 G 1 0 M 26 78 S 0 2 5 0 2 I 0 760 3 0 2 S 1 1 M 47 141 S 0 1 5 0 2 I 0 773 3 0 2 S 1 2 M 34 102 # --- 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 # JXSJ01000003.1:subseq(9649972,98963) exonerate:protein2genome:local gene 48963 50869 315 - . gene_id 1 ; sequence SPP00000634_2.0 ; gene_orientation + JXSJ01000003.1:subseq(9649972,98963) exonerate:protein2genome:local cds 50751 50869 . - . JXSJ01000003.1:subseq(9649972,98963) exonerate:protein2genome:local exon 50751 50869 . - . insertions 0 ; deletions 1 JXSJ01000003.1:subseq(9649972,98963) exonerate:protein2genome:local splice5 50749 50750 . - . intron_id 1 ; splice_site "GT" JXSJ01000003.1:subseq(9649972,98963) exonerate:protein2genome:local intron 49987 50750 . - . intron_id 1 JXSJ01000003.1:subseq(9649972,98963) exonerate:protein2genome:local splice3 49987 49988 . - . intron_id 0 ; splice_site "AG" JXSJ01000003.1:subseq(9649972,98963) exonerate:protein2genome:local cds 49844 49986 . - . JXSJ01000003.1:subseq(9649972,98963) exonerate:protein2genome:local exon 49844 49986 . - . insertions 0 ; deletions 0 JXSJ01000003.1:subseq(9649972,98963) exonerate:protein2genome:local splice5 49842 49843 . - . intron_id 2 ; splice_site "GT" JXSJ01000003.1:subseq(9649972,98963) exonerate:protein2genome:local intron 49067 49843 . - . intron_id 2 JXSJ01000003.1:subseq(9649972,98963) exonerate:protein2genome:local splice3 49067 49068 . - . intron_id 1 ; splice_site "AG" JXSJ01000003.1:subseq(9649972,98963) exonerate:protein2genome:local cds 48963 49066 . - . JXSJ01000003.1:subseq(9649972,98963) exonerate:protein2genome:local exon 48963 49066 . - . insertions 0 ; deletions 0 JXSJ01000003.1:subseq(9649972,98963) exonerate:protein2genome:local similarity 48963 50869 315 - . alignment_id 1 ; Query SPP00000634_2.0 ; Align 50870 5 39 ; Align 50831 19 78 ; Align 49986 46 141 ; Align 49065 94 102 # --- END OF GFF DUMP --- # -- completed exonerate analysis