Command line: [exonerate --exhaustive yes -m p2g --showtargetgff -q ../DI1-human.fa -t DI1-NIVO01071179.1.subseq] Hostname: [sitdoc] C4 Alignment: ------------ Query: SPP00000001_1.0 # Protein # Iodothyronine deiodinase 1 (DI1) # Homo sapiens # Complete Target: NIVO01071179.1:subseq(592045,106696) Ammotragus lervia isolate NAZ scaffold350633_len810912_cov72, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 647 Query range: 109 -> 249 Target range: 56696 -> 49999 110 : PheMetGlnGlyAsnArgProLeuValLeuAsnPheGlySerCysThrUnkProSerPhe : 129 |||:!: !||||||||||||||||||||||||||||||||||||||| ||||||||| PheValLeuGlyAsnArgProLeuValLeuAsnPheGlySerCysThr***ProSerPhe 56696 : TTTGTTTTAGGCAACAGACCCTTGGTGCTGAATTTTGGAAGTTGCACCTgaccttcattt : 56639 130 : MetPheLysPheAspGlnPheLysArgLeuIleGluAspPheSerSerIleAlaAspPhe : 149 :::|||||||||||||||||||||||||||||||||||||||.!!|||:!!||||||||| IlePheLysPheAspGlnPheLysArgLeuIleGluAspPheGlySerValAlaAspPhe 56638 : attttcaaatttgacCAATTCAAGAGACTTATTGAAGATTTTGGTTCCGTAGCAGATTTT : 56579 150 : LeuValIleTyrIleGluGluAlaHisAlaSer{A} >>>> Target Intron 1 > : 161 |||:!!|||||||||||||||||||||||||||{|} 4048 bp LeuIleIleTyrIleGluGluAlaHisAlaSer{A}++ 56578 : CTCATCATTTACATCGAAGAAGCACATGCATCA{G}gt...................... : 56540 162 : >>> {sp}GlyTrpAlaPheLysAsnAsnMetAspIleArgAsnHisGlnAsnLeuGln : 178 {||}|||||||||||||||||||||:!!||||||!:!||||||!:!||||||||| ++{sp}GlyTrpAlaPheLysAsnAsnValAspIleLysAsnHisArgAsnLeuGln 56539 : ...ag{AT}GGCTGGGCTTTTAAGAACAACGTGGACATCAAGAATCACCGGAATCTCCAG : 52444 179 : AspArgLeuGlnAlaAlaHisLeuLeuLeuAlaArgSerProGlnCysProValValVal : 198 |||||||||!:!||||||||||||||||||! !||||||||||||||||||||||||||| AspArgLeuArgAlaAlaHisLeuLeuLeuAspArgSerProGlnCysProValValVal 52443 : GACCGCTTGCGGGCAGCCCACCTGCTGCTGGACAGGAGCCCCCAGTGCCCCGTGGTGGTA : 52384 199 : AspThrMetGlnAsnGlnSerSerGlnLeuTyrAlaAlaLeuProGluArgLeuTyrIle : 218 |||||||||:!!||||||||||||..!! !|||||||||||||||||||||||||||:!: AspThrMetLysAsnGlnSerSerSerCysTyrAlaAlaLeuProGluArgLeuTyrVal 52383 : GACACAATGAAGAACCAGAGCAGCTCGTGCTACGCAGCACTGCCAGAGCGGCTCTACGTG : 52324 219 : IleGlnGluGlyArgIleLeuTyrLys >>>> Target Intron 2 >>>> GlyL : 229 :!!!:!||||||||||||||||||||| 2229 bp |||| LeuArgGluGlyArgIleLeuTyrLys++ ++GlyL 52323 : CTCCGGGAGGGCAGGATCCTGTACAAGgt.........................agGGTA : 50062 230 : ysSerGlyProTrpAsnTyrAsnProGluGluValArgAlaValLeuGluLysLeuHisS : 249 || !!|||||||||||||||:!!||||||||||||||||||||||||||||||||||||| ysProGlyProTrpAsnTyrHisProGluGluValArgAlaValLeuGluLysLeuHisS 50061 : AACCTGGTCCTTGGAACTACCATCCAGAGGAAGTTCGTGCTGTTCTGGAAAAGCTCCACA : 50002 250 : er : 249 || er 50001 : GT : 50000 vulgar: SPP00000001_1.0 109 249 . NIVO01071179.1:subseq(592045,106696) 56696 49999 - 647 M 51 153 S 0 1 5 0 2 I 0 4044 3 0 2 S 1 2 M 66 198 5 0 2 I 0 2225 3 0 2 M 22 66 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2017-11-16 ##type DNA # # # seqname source feature start end score strand frame attributes # NIVO01071179.1:subseq(592045,106696) exonerate:protein2genome:local gene 50000 56696 647 - . gene_id 1 ; sequence SPP00000001_1.0 ; gene_orientation + NIVO01071179.1:subseq(592045,106696) exonerate:protein2genome:local cds 56543 56696 . - . NIVO01071179.1:subseq(592045,106696) exonerate:protein2genome:local exon 56543 56696 . - . insertions 0 ; deletions 0 NIVO01071179.1:subseq(592045,106696) exonerate:protein2genome:local splice5 56541 56542 . - . intron_id 1 ; splice_site "gt" NIVO01071179.1:subseq(592045,106696) exonerate:protein2genome:local intron 52495 56542 . - . intron_id 1 NIVO01071179.1:subseq(592045,106696) exonerate:protein2genome:local splice3 52495 52496 . - . intron_id 0 ; splice_site "AG" NIVO01071179.1:subseq(592045,106696) exonerate:protein2genome:local cds 52295 52494 . - . NIVO01071179.1:subseq(592045,106696) exonerate:protein2genome:local exon 52295 52494 . - . insertions 0 ; deletions 0 NIVO01071179.1:subseq(592045,106696) exonerate:protein2genome:local splice5 52293 52294 . - . intron_id 2 ; splice_site "GT" NIVO01071179.1:subseq(592045,106696) exonerate:protein2genome:local intron 50066 52294 . - . intron_id 2 NIVO01071179.1:subseq(592045,106696) exonerate:protein2genome:local splice3 50066 50067 . - . intron_id 1 ; splice_site "AG" NIVO01071179.1:subseq(592045,106696) exonerate:protein2genome:local cds 50000 50065 . - . NIVO01071179.1:subseq(592045,106696) exonerate:protein2genome:local exon 50000 50065 . - . insertions 0 ; deletions 0 NIVO01071179.1:subseq(592045,106696) exonerate:protein2genome:local similarity 50000 56696 647 - . alignment_id 1 ; Query SPP00000001_1.0 ; Align 56697 110 153 ; Align 52493 162 198 ; Align 50066 228 66 # --- END OF GFF DUMP --- # -- completed exonerate analysis