Command line: [exonerate -m p2g --showtargetgff -q ../Sel.human/DI1.fa -t gi|397456552|gb|JH738948.1|subseq.fa -E] Hostname: [ws101062] ** (process:3183): WARNING **: Exhaustively generating suboptimal alignments will be VERY SLOW C4 Alignment: ------------ Query: SPP00000002_2.0 # Protein # Iodothyronine deiodinase 1 (DI1) # Human (Homo sapiens) Target: gi|397456552|gb|JH738948.1|:subseq(11286,101000) Alligator mississippiensis unplaced genomic scaffold scaffold-13070, whole genome shotgun sequence Model: protein2genome:local Raw score: 392 Query range: 25 -> 246 Target range: 7156 -> 50482 26 : GlyLysValLeuLeuIleLeuPheProAspArgValLysArgAsnIleLeuAlaMetGly : 45 !.!|||:!!!!!|||:!::!! ! ! !!|||..!!:!!!!||| !||| !:!! ! AlaLysLeuPheLeuValValValAspHisArgThrArgSerAsnAspLeuLysLeuGln 7157 : GCTAAACTATTCTTAGTGGTGGTAGATCATAGAACAAGGAGCAATGATCTCAAGCTGCAG : 7214 46 : Glu{L} >>>> Target Intron 1 >>>> {ys}Thr{G} >>>> Target I : 49 :!!{ } 437 bp { !}|||{|} 1301 Gln{V}+- ++{al}Thr{G}++ 7215 : CAA{G}gg.........................ag{TG}ACA{G}gt............. : 7666 50 : ntron 2 >>>> {ly}MetThrArgAsnProHisPheSerHis{A} >>>> Targe : 59 7 bp {||}! ! !!:!||| !||||||..!...{ } 5 ++{ly}LysPheGlnAsnPheHisPheGlnArg{A}++ 7667 : ............ag{GA}AAGTTTCAAAACTTTCATTTTCAAAGA{A}gt.......... : 20713 60 : t Intron 3 >>>> {sp}AsnTrpIleProThrPhePheSerThrGlnTyrPheTrp : 72 105 bp { }||| :!!||| !|||!:!!!! !.!.! !.!.!:: ++{rg}AsnIleValProLeuPheTyrCysLeuAspCysLeuTyr 20714 : ...............ag{GA}AATATAGTTCCTCTATTTTATTGCTTAGACTGCCTATAC : 25854 73 : {Ph} >>>> Target Intron 4 >>>> {e}ValLeuLysValArgTrpGlnArg : 81 {!:} 23260 bp {:} !! ! ! !! !||| ||| {Tr}++ ++{p}Arg***AlaGlyProTrpGlyArg 25855 : {TG}gt.........................ag{G}CGGTGAGCCGGTCCCTGGGGCCGT : 49141 82 : LeuGluAspThr{T} >>>> Target Intron 5 >>>> {hr}GluLeuGlyGly : 90 ! ! !!{.} 838 bp {!!}.!.:!!|||!.! AlaProArgPro{A}++ -+{la}HisValGlyAla 49142 : GCCCCGAGGCCA{G}gt.........................gg{CG}CACGTGGGCGCG : 50006 91 : LeuAlaProAsnCysProValValArgLeuSerGlyGln<->ArgCysAsnIleTrpGlu : 109 ! !||||||||| ! !|||:!!||||||..!|||||| ||| ..!:!! !!!: ArgAlaProAsnProGluValIleArgLeuAspGlyGlnLysArgLeuArgLeuLeuAsp 50007 : CGGGCGCCCAACCCCGAGGTCATCCGGCTGGACGGGCAGAAGCGGCTGCGCCTCCTGGAC : 50066 110 : PheMetGlnGlyAsnArgProLeuValLeuAsnPheGlySerCysThrUnkProSerPhe : 129 ||| !::|||!!.||||||||||||||||||||||||||||||||| ||| !!||| PheAlaArgGlyLysArgProLeuValLeuAsnPheGlySerCysThr***ProProPhe 50067 : TTCGCCCGCGGCAAGCGGCCGCTCGTGCTCAACTTCGGCAGCTGCACCTGACCGCCCTTC : 50126 130 : MetPheLysPheAspGlnPheLysArgLeuIleGluAspPheSerSerIleAlaAspPhe : 149 ||| !:::.!. ! |||:!!|||||| !! ! !!||| !..!|||||||||||| MetAlaArgLeuArgAlaPheGlnArgLeuAlaAlaHisPheLeuAspIleAlaAspPhe 50127 : ATGGCCCGCCTGCGCGCCTTCCAGCGCCTGGCCGCGCACTTCCTGGACATCGCCGACTTC : 50186 150 : LeuValIleTyrIleGluGluAlaHisAlaSerAspGlyTrpAlaPheLysAsn<->Asn : 168 |||:!!:!:|||||||||||||||||| !!||||||||||||!.! !..!:!! ! LeuLeuValTyrIleGluGluAlaHisProSerAspGlyTrpValSerSerAspAlaAla 50187 : CTGCTGGTGTACATCGAGGAGGCGCACCCCTCCGACGGCTGGGTCAGCTCGGACGCCGCC : 50246 169 : MetAspIleArgAsnHisGlnAsnLeuGlnAspArgLeuGlnAlaAlaHisLeuLeuLeu : 188 ! |||! !!!.|||||| !|||||||||||||||!:!||||||!!.|||:!!! ! TyrLysIleProLysHisGlnCysLeuGlnAspArgLeuArgAlaAlaGlnLeuMetArg 50247 : TACAAGATCCCCAAGCACCAGTGCCTGCAGGACCGGCTGCGGGCGGCGCAGCTCATGCGG : 50306 189 : AlaArgSerProGlnCysProValValValAspThrMetGlnAsnGlnSerSerGlnLeu : 208 ! ! !!:!!||| !||||||:!!!.!||||||||||||.!.||| |||||| ! GluGlyAlaProGlyCysProLeuAlaValAspThrMetAspAsnAlaSerSerAlaAla 50307 : GAAGGCGCCCCCGGGTGCCCGCTGGCGGTGGACACCATGGACAACGCCTCCAGCGCCGCC : 50366 209 : TyrAlaAlaLeuProGluArgLeuTyrIleIleGlnGluGlyArgIleLeuTyrLysGly : 228 |||!.!|||! ! !|||||||||||||||:!!||||||! !|||:!!:!!|||:!!||| TyrGlyAlaTyrPheGluArgLeuTyrIleValGlnGluGluArgValMetTyrGlnGly 50367 : TACGGAGCCTACTTCGAGCGGCTGTACATCGTCCAGGAGGAGCGCGTCATGTACCAGGGC : 50426 229 : LysSerGlyProTrpAsnTyrAsnProGluGluValArgAlaValLeuGluLys : 246 !!|||||| !..!|||!!. !..!|||:!!|||::! !|||!!::!! GlyArgGlyProGluGlyTyrLysIleSerGluLeuArgSerTrpLeuAspGln 50427 : GGCCGCGGCCCCGAGGGCTACAAGATCTCGGAGCTGCGCAGCTGGCTCGACCAG : 50482 vulgar: SPP00000002_2.0 25 246 . gi|397456552|gb|JH738948.1|:subseq(11286,101000) 7156 50482 + 392 M 21 63 S 0 1 5 0 2 I 0 433 3 0 2 S 1 2 M 1 3 S 0 1 5 0 2 I 0 13013 3 0 2 S 1 2 M 9 27 S 0 1 5 0 2 I 0 5101 3 0 2 S 1 2 M 13 39 S 0 2 5 0 2 I 0 23256 3 0 2 S 1 1 M 12 36 S 0 1 5 0 2 I 0 834 3 0 2 S 1 2 M 17 51 G 0 3 M 64 192 G 0 3 M 79 237 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2013-11-15 ##type DNA # # # seqname source feature start end score strand frame attributes # gi|397456552|gb|JH738948.1|:subseq(11286,101000) exonerate:protein2genome:local gene 7157 50482 392 + . gene_id 1 ; sequence SPP00000002_2.0 ; gene_orientation + gi|397456552|gb|JH738948.1|:subseq(11286,101000) exonerate:protein2genome:local cds 7157 7220 . + . gi|397456552|gb|JH738948.1|:subseq(11286,101000) exonerate:protein2genome:local exon 7157 7220 . + . insertions 0 ; deletions 0 gi|397456552|gb|JH738948.1|:subseq(11286,101000) exonerate:protein2genome:local splice5 7221 7222 . + . intron_id 1 ; splice_site "GG" gi|397456552|gb|JH738948.1|:subseq(11286,101000) exonerate:protein2genome:local intron 7221 7657 . + . intron_id 1 gi|397456552|gb|JH738948.1|:subseq(11286,101000) exonerate:protein2genome:local splice3 7656 7657 . + . intron_id 0 ; splice_site "AG" gi|397456552|gb|JH738948.1|:subseq(11286,101000) exonerate:protein2genome:local cds 7658 7663 . + . gi|397456552|gb|JH738948.1|:subseq(11286,101000) exonerate:protein2genome:local exon 7658 7663 . + . insertions 0 ; deletions 0 gi|397456552|gb|JH738948.1|:subseq(11286,101000) exonerate:protein2genome:local splice5 7664 7665 . + . intron_id 2 ; splice_site "GT" gi|397456552|gb|JH738948.1|:subseq(11286,101000) exonerate:protein2genome:local intron 7664 20680 . + . intron_id 2 gi|397456552|gb|JH738948.1|:subseq(11286,101000) exonerate:protein2genome:local splice3 20679 20680 . + . intron_id 1 ; splice_site "AG" gi|397456552|gb|JH738948.1|:subseq(11286,101000) exonerate:protein2genome:local cds 20681 20710 . + . gi|397456552|gb|JH738948.1|:subseq(11286,101000) exonerate:protein2genome:local exon 20681 20710 . + . insertions 0 ; deletions 0 gi|397456552|gb|JH738948.1|:subseq(11286,101000) exonerate:protein2genome:local splice5 20711 20712 . + . intron_id 3 ; splice_site "GT" gi|397456552|gb|JH738948.1|:subseq(11286,101000) exonerate:protein2genome:local intron 20711 25815 . + . intron_id 3 gi|397456552|gb|JH738948.1|:subseq(11286,101000) exonerate:protein2genome:local splice3 25814 25815 . + . intron_id 2 ; splice_site "AG" gi|397456552|gb|JH738948.1|:subseq(11286,101000) exonerate:protein2genome:local cds 25816 25858 . + . gi|397456552|gb|JH738948.1|:subseq(11286,101000) exonerate:protein2genome:local exon 25816 25858 . + . insertions 0 ; deletions 0 gi|397456552|gb|JH738948.1|:subseq(11286,101000) exonerate:protein2genome:local splice5 25859 25860 . + . intron_id 4 ; splice_site "GT" gi|397456552|gb|JH738948.1|:subseq(11286,101000) exonerate:protein2genome:local intron 25859 49118 . + . intron_id 4 gi|397456552|gb|JH738948.1|:subseq(11286,101000) exonerate:protein2genome:local splice3 49117 49118 . + . intron_id 3 ; splice_site "AG" gi|397456552|gb|JH738948.1|:subseq(11286,101000) exonerate:protein2genome:local cds 49119 49156 . + . gi|397456552|gb|JH738948.1|:subseq(11286,101000) exonerate:protein2genome:local exon 49119 49156 . + . insertions 0 ; deletions 0 gi|397456552|gb|JH738948.1|:subseq(11286,101000) exonerate:protein2genome:local splice5 49157 49158 . + . intron_id 5 ; splice_site "GT" gi|397456552|gb|JH738948.1|:subseq(11286,101000) exonerate:protein2genome:local intron 49157 49994 . + . intron_id 5 gi|397456552|gb|JH738948.1|:subseq(11286,101000) exonerate:protein2genome:local splice3 49993 49994 . + . intron_id 4 ; splice_site "GG" gi|397456552|gb|JH738948.1|:subseq(11286,101000) exonerate:protein2genome:local cds 49995 50482 . + . gi|397456552|gb|JH738948.1|:subseq(11286,101000) exonerate:protein2genome:local exon 49995 50482 . + . insertions 6 ; deletions 0 gi|397456552|gb|JH738948.1|:subseq(11286,101000) exonerate:protein2genome:local similarity 7157 50482 392 + . alignment_id 1 ; Query SPP00000002_2.0 ; Align 7157 26 63 ; Align 7660 48 3 ; Align 20683 50 27 ; Align 25818 60 39 ; Align 49120 74 36 ; Align 49997 87 51 ; Align 50051 104 192 ; Align 50246 168 237 # --- END OF GFF DUMP --- # -- completed exonerate analysis