Command line: [exonerate -m p2g --showtargetgff -q ../SelU3-human.fa -t SelU3-NIVO01017783.1.subseq] Hostname: [sitdoc] C4 Alignment: ------------ Query: SPP00000028_1.0 # Protein # Selenoprotein U3 (SelU3) # Homo sapiens # Complete Target: NIVO01017783.1:subseq(851831,1164089) Ammotragus lervia isolate NAZ scaffold4661_len2021309_cov71, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 866 Query range: 0 -> 193 Target range: 52396 -> 50002 1 : MetSerThrValAspLeuAlaArgValGlyAlaCysIleLeuLysHisAlaValThrGly : 20 ||||||||||||||||||||||||||||||||||||:!!||||||||||||||||||||| MetSerThrValAspLeuAlaArgValGlyAlaCysValLeuLysHisAlaValThrGly 52396 : ATGAGTACGGTGGACCTTGCCCGCGTTGGCGCGTGTGTCCTGAAGCACGCGGTCACCGGG : 52339 21 : Glu >>>> Target Intron 1 >>>> AlaValGluLeuArgSerLeuTrpArgG : 31 ||| 199 bp ||||||||||||||||||||||||!:!| Glu++ ++AlaValGluLeuArgSerLeuTrpGlnG 52338 : GAGgt.........................agGCCGTAGAGCTGCGGAGCCTGTGGCAGG : 52107 32 : luHisAlaCysValValAlaGlyLeuArgArgPheGlyCysValValCysArgTrpIleA : 51 ||!!.||||||||||||||||||||||||||||||||||||:!!|||||||||||||||| luGlnAlaCysValValAlaGlyLeuArgArgPheGlyCysMetValCysArgTrpIleA 52106 : AGCAGGCGTGCGTGGTGGCCGGCTTGCGGCGCTTCGGCTGCATGGTGTGTCGCTGGATCG : 52047 52 : laGlnAspLeuSerSerLeuAlaGlyLeuLeuAspGlnHisGlyValArgLeuValGlyV : 71 ||!::|||||||||!:!||| !||||||||||||||||||||||||||||||||||||| laArgAspLeuSerAsnLeuLysGlyLeuLeuAspGlnHisGlyValArgLeuValGlyV 52046 : CCCGCGACCTCAGCAACCTCAAGGGGCTCCTGGACCAGCACGGCGTGCGCCTGGTGGGCG : 51987 72 : alGlyProGluAlaLeuGlyLeuGlnGluPheLeuAspGlyAspTyrPheAlaGly{G} : 90 |||||||||||||||||||||||||||||||||||||||||! !||||||||||||{|} alGlyProGluAlaLeuGlyLeuGlnGluPheLeuAspGlyGlyTyrPheAlaGly{G}+ 51986 : TGGGGCCAGAGGCCCTGGGCCTGCAGGAGTTCCTGGATGGTGGCTACTTTGCAGGA{G}g : 51929 91 : >>>> Target Intron 2 >>>> {lu}LeuTyrLeuAspGluSerLysGlnLeuT : 100 700 bp {||}||||||||||||||||||||||||!!!| + ++{lu}LeuTyrLeuAspGluSerLysGlnPheT 51928 : t.........................ag{AG}CTCTACCTGGACGAGAGCAAGCAGTTTT : 51200 101 : yrLysGluLeuGlyPheLys{Ar} >>>> Target Intron 3 >>>> {g}TyrA : 109 ||||||||||||||||||||{||} 543 bp {|}|||| yrLysGluLeuGlyPheLys{Ar}++ ++{g}TyrA 51199 : ACAAGGAGCTGGGCTTCAAG{CG}gt.........................ag{G}TACA : 50630 110 : snSerLeuSerIleLeuProAlaAlaLeuGlyLysProValArgAspValAlaAlaLys : 129 ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| snSerLeuSerIleLeuProAlaAlaLeuGlyLysProValArgAspValAlaAlaLys+ 50629 : ACAGCTTGAGCATCCTGCCGGCTGCCCTGGGAAAGCCTGTTCGTGATGTGGCCGCCAAGg : 50570 130 : >>>> Target Intron 4 >>>> AlaLysAlaValGlyIleGlnGlyAsnLeuSe : 139 80 bp |||||||||||||||||||||||||||||||| + ++AlaLysAlaValGlyIleGlnGlyAsnLeuSe 50569 : t.........................agGCCAAGGCTGTCGGCATCCAGGGGAATTTGTC : 50460 140 : rGlyAspLeuLeuGlnSerGlyGlyLeuLeuValValSerLys{G} >>>> Target : 154 |||||||||||||||||||||||||||||||||||||:!!|||{|} 29 rGlyAspLeuLeuGlnSerGlyGlyLeuLeuValValAlaLys{G}++ 50459 : TGGGGATCTGCTGCAGAGCGGAGGGCTGCTGGTGGTTGCCAAA{G}gt............ : 50412 155 : Intron 5 >>>> {ly}GlyAspLysValLeuLeuHisPheValGlnLysSerProGl : 168 3 bp {||}||||||||||||||||||||||||||||||||||||||||| ++{ly}GlyAspLysValLeuLeuHisPheValGlnLysSerProGl 50411 : .............ag{GT}GGTGATAAAGTGCTGCTGCACTTTGTCCAGAAGTCCCCGGG : 50080 169 : yAspTyrValProLysGluHisIleLeuGlnValLeuGlyIleSerAlaGluValCysAl : 188 |||||||!.!||| ||| !|||||||||!.!||||||||||||||||||||| !! ! yAspTyrAlaProProGluSerIleLeuGlnAlaLeuGlyIleSerAlaGluValGlyPr 50079 : TGACTACGCCCCGCCTGAGAGCATCCTGCAGGCCTTGGGCATCTCTGCAGAGGTCGGCCC : 50020 189 : aSerAspProProGln : 193 !|||!!:||||||||| oSerGluProProGln 50019 : CAGCGAACCGCCCCAG : 50003 vulgar: SPP00000028_1.0 0 193 . NIVO01017783.1:subseq(851831,1164089) 52396 50002 - 866 M 21 63 5 0 2 I 0 195 3 0 2 M 68 204 S 0 1 5 0 2 I 0 696 3 0 2 S 1 2 M 16 48 S 0 2 5 0 2 I 0 539 3 0 2 S 1 1 M 21 63 5 0 2 I 0 76 3 0 2 M 25 75 S 0 1 5 0 2 I 0 289 3 0 2 S 1 2 M 39 117 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2017-11-21 ##type DNA # # # seqname source feature start end score strand frame attributes # NIVO01017783.1:subseq(851831,1164089) exonerate:protein2genome:local gene 50003 52396 866 - . gene_id 1 ; sequence SPP00000028_1.0 ; gene_orientation + NIVO01017783.1:subseq(851831,1164089) exonerate:protein2genome:local cds 52334 52396 . - . NIVO01017783.1:subseq(851831,1164089) exonerate:protein2genome:local exon 52334 52396 . - . insertions 0 ; deletions 0 NIVO01017783.1:subseq(851831,1164089) exonerate:protein2genome:local splice5 52332 52333 . - . intron_id 1 ; splice_site "GT" NIVO01017783.1:subseq(851831,1164089) exonerate:protein2genome:local intron 52135 52333 . - . intron_id 1 NIVO01017783.1:subseq(851831,1164089) exonerate:protein2genome:local splice3 52135 52136 . - . intron_id 0 ; splice_site "AG" NIVO01017783.1:subseq(851831,1164089) exonerate:protein2genome:local cds 51930 52134 . - . NIVO01017783.1:subseq(851831,1164089) exonerate:protein2genome:local exon 51930 52134 . - . insertions 0 ; deletions 0 NIVO01017783.1:subseq(851831,1164089) exonerate:protein2genome:local splice5 51928 51929 . - . intron_id 2 ; splice_site "GT" NIVO01017783.1:subseq(851831,1164089) exonerate:protein2genome:local intron 51230 51929 . - . intron_id 2 NIVO01017783.1:subseq(851831,1164089) exonerate:protein2genome:local splice3 51230 51231 . - . intron_id 1 ; splice_site "AG" NIVO01017783.1:subseq(851831,1164089) exonerate:protein2genome:local cds 51178 51229 . - . NIVO01017783.1:subseq(851831,1164089) exonerate:protein2genome:local exon 51178 51229 . - . insertions 0 ; deletions 0 NIVO01017783.1:subseq(851831,1164089) exonerate:protein2genome:local splice5 51176 51177 . - . intron_id 3 ; splice_site "GT" NIVO01017783.1:subseq(851831,1164089) exonerate:protein2genome:local intron 50635 51177 . - . intron_id 3 NIVO01017783.1:subseq(851831,1164089) exonerate:protein2genome:local splice3 50635 50636 . - . intron_id 2 ; splice_site "AG" NIVO01017783.1:subseq(851831,1164089) exonerate:protein2genome:local cds 50571 50634 . - . NIVO01017783.1:subseq(851831,1164089) exonerate:protein2genome:local exon 50571 50634 . - . insertions 0 ; deletions 0 NIVO01017783.1:subseq(851831,1164089) exonerate:protein2genome:local splice5 50569 50570 . - . intron_id 4 ; splice_site "GT" NIVO01017783.1:subseq(851831,1164089) exonerate:protein2genome:local intron 50491 50570 . - . intron_id 4 NIVO01017783.1:subseq(851831,1164089) exonerate:protein2genome:local splice3 50491 50492 . - . intron_id 3 ; splice_site "AG" NIVO01017783.1:subseq(851831,1164089) exonerate:protein2genome:local cds 50415 50490 . - . NIVO01017783.1:subseq(851831,1164089) exonerate:protein2genome:local exon 50415 50490 . - . insertions 0 ; deletions 0 NIVO01017783.1:subseq(851831,1164089) exonerate:protein2genome:local splice5 50413 50414 . - . intron_id 5 ; splice_site "GT" NIVO01017783.1:subseq(851831,1164089) exonerate:protein2genome:local intron 50122 50414 . - . intron_id 5 NIVO01017783.1:subseq(851831,1164089) exonerate:protein2genome:local splice3 50122 50123 . - . intron_id 4 ; splice_site "AG" NIVO01017783.1:subseq(851831,1164089) exonerate:protein2genome:local cds 50003 50121 . - . NIVO01017783.1:subseq(851831,1164089) exonerate:protein2genome:local exon 50003 50121 . - . insertions 0 ; deletions 0 NIVO01017783.1:subseq(851831,1164089) exonerate:protein2genome:local similarity 50003 52396 866 - . alignment_id 1 ; Query SPP00000028_1.0 ; Align 52397 1 63 ; Align 52135 22 204 ; Align 51228 91 48 ; Align 50634 108 63 ; Align 50491 129 75 ; Align 50120 155 117 # --- END OF GFF DUMP --- # -- completed exonerate analysis