Command line: [exonerate -m p2g --showtargetgff -q sephs1_zf.aa.fa.canviat -t genomesephs1_JXSJ01000007.1.fa] Hostname: [ws116078U] C4 Alignment: ------------ Query: SPP00000629_2.0 # Protein # Selenophosphate synthetase (SEPHS) # Zebrafish Target: JXSJ01000007.1:subseq(6430332,103396) Miichthys miiuy scaffold7, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 1933 Query range: 0 -> 392 Target range: 53396 -> 49999 1 : MetSerValArgGluSerPheAsnProGluSerTyrGluLeuAspLysAsnPheArgLeu : 20 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| MetSerValArgGluSerPheAsnProGluSerTyrGluLeuAspLysAsnPheArgLeu 53396 : ATGTCCGTGAGGGAATCCTTTAACCCTGAAAGCTATGAGCTGGACAAGAACTTCAGGCTC : 53339 21 : ThrArgPheThrGluLeuLysGlyThrGlyCysLys >>>> Target Intron 1 > : 33 |||||||||.!!|||||||||||||||||||||||| 111 bp ThrArgPheAlaGluLeuLysGlyThrGlyCysLys++ 53338 : ACACGCTTTGCAGAGCTCAAGGGCACAGGCTGCAAGgt...................... : 53298 34 : >>> ValProGlnAspValLeuGlnLysLeuLeuGluSerLeuGlnGluAsnHisTyrG : 51 |||||||||||||||||||||||||||||||||!!!||||||||||||||||||| ++ValProGlnAspValLeuGlnLysLeuLeuGluThrLeuGlnGluAsnHisTyrG 53297 : ...agGTGCCCCAAGATGTATTACAGAAACTGCTAGAAACCCTACAGGAGAATCACTATC : 53135 52 : lnGluAspGluGlnPheLeuGlyAlaValMetProArgLeu{G} >>>> Target In : 65 |||||||||||||||||||||||||||||||||||||||||{|} 312 lnGluAspGluGlnPheLeuGlyAlaValMetProArgLeu{G}++ 53134 : AAGAGGATGAACAGTTCCTGGGAGCAGTTATGCCTCGATTA{g}gt.............. : 53090 66 : tron 2 >>>> {ly}IleGlyMetAspThrCysValIleProLeuArgHisGlyGlyL : 80 bp {||}||||||||||||||||||||||||||||||||||||||||||| ++{ly}IleGlyMetAspThrCysValIleProLeuArgHisGlyGlyL 53089 : ...........ag{gc}aTAGGCATGGACACCTGTGTGATCCCCCTCAGACATGGAGGCC : 52736 81 : euSerLeuValGlnThrThrAspTyrIleTyrProIleValAspAspProTyrMetMet : 100 ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| euSerLeuValGlnThrThrAspTyrIleTyrProIleValAspAspProTyrMetMet+ 52735 : TTTCTCTGGTCCAGACGACAGATTACATCTACCCCATTGTGGATGACCCGTACATGATGg : 52676 101 : >>>> Target Intron 3 >>>> GlyArgIleAlaCysAlaAsnValLeuSerAs : 110 91 bp |||||||||||||||||||||||||||||||| + ++GlyArgIleAlaCysAlaAsnValLeuSerAs 52675 : t.........................agGGAAGAATTGCTTGTGCCAATGTTCTAAGTGA : 52555 111 : pLeuTyrAlaMetGlyValThrGluCysAspAsnMetLeuMetLeuLeuGlyValSerAs : 130 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| pLeuTyrAlaMetGlyValThrGluCysAspAsnMetLeuMetLeuLeuGlyValSerAs 52554 : CCTATATGCCATGGGAGTGACAGAGTGTGACAACATGTTGATGCTTCTGGGAGTcagcaa : 52495 131 : nLysLeuThrGluLys >>>> Target Intron 4 >>>> GluArgAspLysVal : 140 ||||::::::|||||| 199 bp ||||||||||||||| nLysMetSerGluLys++ ++GluArgAspLysVal 52494 : taaaatgtcagagaaggt.........................agGAAAGAGACAAAGTT : 52266 141 : MetProLeuValIleGlnGlyPheLysAspAlaSerGluGluAlaGlyThrSerValThr : 160 |||||||||:!!|||||||||||||||||||||||||||||||||||||||||||||||| MetProLeuIleIleGlnGlyPheLysAspAlaSerGluGluAlaGlyThrSerValThr 52265 : ATGCCACTGATCATCCAGGGATTCAAGGATGCATCGGAGGAGGCAGGCACATCTGTAACA : 52206 161 : GlyGlyGlnThrValIleAsnProTrpIleValLeuGlyGlyValAlaThrThrValCys : 180 |||||||||||||||:!!|||||||||:!:|||:!!|||||||||||||||||||||||| GlyGlyGlnThrValLeuAsnProTrpValValMetGlyGlyValAlaThrThrValCys 52205 : GGAGGACAGACGGTGCTCAACCCATGGGTGGTGATGGGAGGAGTTGCGACAACAGTCTGC : 52146 181 : GlnProAsnGluPheIle{Me} >>>> Target Intron 5 >>>> {t}ProAsp : 189 ||||||||||||||||||{||} 96 bp {|}|||||| GlnProAsnGluPheIle{Me}++ ++{t}ProAsp 52145 : CAGCCCAACGAATTTATC{AT}gt.........................ag{G}CCTGAC : 52023 190 : AsnAlaValProGlyAspValLeuValLeuThrLysProLeuGlyThrGlnValAlaVal : 209 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| AsnAlaValProGlyAspValLeuValLeuThrLysProLeuGlyThrGlnValAlaVal 52022 : AATGCAGTACCAGGAGACGTGTTAGTGTTGACCAAGCCGCTTGGAACACAAGTTGCAGTT : 51963 210 : AlaValHisGlnTrpLeuAspIle >>>> Target Intron 6 >>>> ProGluL : 220 |||||||||||||||||||||||| 84 bp ||||||| AlaValHisGlnTrpLeuAspIle++ ++ProGluL 51962 : GCAGTACACCAGTGGCTAGATATTgt.........................agCCTGAGA : 51846 221 : ysTrpAsnLysIleLysLeuValValThrGlnGluAspValGluLeuAlaTyrGlnGluA : 240 |||||||||||||||||||||||||||||||||||||||||||||||||||||!!.|||| ysTrpAsnLysIleLysLeuValValThrGlnGluAspValGluLeuAlaTyrHisGluA 51845 : AGTGGAACAAGATCAAATTGGTGGTAACCCAGGAAGATGTAGAGCTGGCATACCATGAAG : 51786 241 : laMetLeuAsnMetAlaArgLeuAsnArgThr{A} >>>> Target Intron 7 >> : 251 |||||:!!||||||||||||||||||||||||{|} 589 bp laMetMetAsnMetAlaArgLeuAsnArgThr{A}++ 51785 : CCATGATGAACATGGCTCGGCTCAACAGGACA{G}gt....................... : 51750 252 : >> {la}AlaGlyLeuMetHisThrPheAsnAlaHisAlaAlaThrAspIleThrGlyP : 269 {||}|||||||||||||||||||||||||||||||||||||||||||||||||||| ++{la}AlaGlyLeuMetHisThrPheAsnAlaHisAlaAlaThrAspIleThrGlyP 51749 : ..ag{CG}GCCGGTCTCATGCACACCTTCAATGCTCATGCAGCCACAGACATCACAGGGT : 51110 270 : heGlyIleLeuGlyHisAlaGlnAsnLeuAlaArgGlnGlnArgThrGluValSerPheV : 289 |||||||||||||||||||||||!..||||||||||||||||||!:!||||||||||||| heGlyIleLeuGlyHisAlaGlnThrLeuAlaArgGlnGlnArgSerGluValSerPheV 51109 : TTGGCATCCTTGGTCACGCTCAAACGTTGGCACGGCAACAGCGAAGCGAGGTGTCGTTTG : 51050 290 : alIleHisAsnLeuProValLeuAlaLysMetAlaAlaValSerLysAlaCysGlyAsnM : 309 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| alIleHisAsnLeuProValLeuAlaLysMetAlaAlaValSerLysAlaCysGlyAsnM 51049 : TCATCCACAACCTTCCGGTCCTTGCCAAGATGGCCGCAGTGTCCAAGGCCTGTGGGAACA : 50990 310 : etPheGlyLeuMetHisGlyThrCysProGluThrSer{G} >>>> Target Intro : 322 ||||||||||||||||||||||||||||||||||||||{|} 739 bp etPheGlyLeuMetHisGlyThrCysProGluThrSer{G}++ 50989 : TGTTTGGACTGATGCACGGCACGTGCCCTGAAACAtca{g}gt................. : 50948 323 : n 8 >>>> {ly}GlyLeuLeuIleCysLeuProArgGluGlnAlaAlaArgPheCysA : 338 {||}|||||||||||||||||||||||||||||||||||||||||||||| ++{ly}GlyLeuLeuIleCysLeuProArgGluGlnAlaAlaArgPheCysA 50947 : ........ag{GA}GGCCTGCTTATCTGTCTGCCTCGGGAGCAGGCCGCTCGCTTCTGTG : 50164 339 : laGluIleLysSerProLysTyrGlyGluGlyHisGlnAlaTrpIleIleGlyIleValG : 358 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| laGluIleLysSerProLysTyrGlyGluGlyHisGlnAlaTrpIleIleGlyIleValG 50163 : CCGAGATCAAATCCCCAAAATACGGCGAGGGCCACCAAGCGTGGATCATTGGGATTGTAG : 50104 359 : luLysGlyAsnArgThrAlaArgIleIleAspLysProArgIleIleGluValAlaProG : 378 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| luLysGlyAsnArgThrAlaArgIleIleAspLysProArgIleIleGluValAlaProG 50103 : AGAAAGGAAACCGCACTGCTCGCATCATTGACAAACCACGGATCATAGAGGTGGCACCAC : 50044 379 : lnValAlaThrGlnAsnValAsnThrThrProGlyAlaThrSer : 392 ||!.!|||||||||||||||||| !!|||||||||||||||||| lnAlaAlaThrGlnAsnValAsnProThrProGlyAlaThrSer 50043 : AAGCGGCCACGCAGAACGTCAACCCAACTCCTGGTGCAACTTCT : 50000 vulgar: SPP00000629_2.0 0 392 . JXSJ01000007.1:subseq(6430332,103396) 53396 49999 - 1933 M 32 96 5 0 2 I 0 107 3 0 2 M 32 96 S 0 1 5 0 2 I 0 308 3 0 2 S 1 2 M 34 102 5 0 2 I 0 87 3 0 2 M 36 108 5 0 2 I 0 195 3 0 2 M 51 153 S 0 2 5 0 2 I 0 92 3 0 2 S 1 1 M 30 90 5 0 2 I 0 80 3 0 2 M 33 99 S 0 1 5 0 2 I 0 585 3 0 2 S 1 2 M 70 210 S 0 1 5 0 2 I 0 735 3 0 2 S 1 2 M 70 210 # --- 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 # JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local gene 50000 53396 1933 - . gene_id 1 ; sequence SPP00000629_2.0 ; gene_orientation + JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local cds 53301 53396 . - . JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local exon 53301 53396 . - . insertions 0 ; deletions 0 JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local splice5 53299 53300 . - . intron_id 1 ; splice_site "GT" JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local intron 53190 53300 . - . intron_id 1 JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local splice3 53190 53191 . - . intron_id 0 ; splice_site "AG" JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local cds 53093 53189 . - . JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local exon 53093 53189 . - . insertions 0 ; deletions 0 JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local splice5 53091 53092 . - . intron_id 2 ; splice_site "gt" JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local intron 52781 53092 . - . intron_id 2 JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local splice3 52781 52782 . - . intron_id 1 ; splice_site "ag" JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local cds 52677 52780 . - . JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local exon 52677 52780 . - . insertions 0 ; deletions 0 JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local splice5 52675 52676 . - . intron_id 3 ; splice_site "GT" JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local intron 52586 52676 . - . intron_id 3 JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local splice3 52586 52587 . - . intron_id 2 ; splice_site "AG" JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local cds 52478 52585 . - . JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local exon 52478 52585 . - . insertions 0 ; deletions 0 JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local splice5 52476 52477 . - . intron_id 4 ; splice_site "GT" JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local intron 52279 52477 . - . intron_id 4 JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local splice3 52279 52280 . - . intron_id 3 ; splice_site "AG" JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local cds 52124 52278 . - . JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local exon 52124 52278 . - . insertions 0 ; deletions 0 JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local splice5 52122 52123 . - . intron_id 5 ; splice_site "GT" JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local intron 52028 52123 . - . intron_id 5 JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local splice3 52028 52029 . - . intron_id 4 ; splice_site "AG" JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local cds 51937 52027 . - . JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local exon 51937 52027 . - . insertions 0 ; deletions 0 JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local splice5 51935 51936 . - . intron_id 6 ; splice_site "GT" JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local intron 51853 51936 . - . intron_id 6 JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local splice3 51853 51854 . - . intron_id 5 ; splice_site "AG" JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local cds 51753 51852 . - . JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local exon 51753 51852 . - . insertions 0 ; deletions 0 JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local splice5 51751 51752 . - . intron_id 7 ; splice_site "GT" JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local intron 51164 51752 . - . intron_id 7 JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local splice3 51164 51165 . - . intron_id 6 ; splice_site "AG" JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local cds 50951 51163 . - . JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local exon 50951 51163 . - . insertions 0 ; deletions 0 JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local splice5 50949 50950 . - . intron_id 8 ; splice_site "gt" JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local intron 50212 50950 . - . intron_id 8 JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local splice3 50212 50213 . - . intron_id 7 ; splice_site "AG" JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local cds 50000 50211 . - . JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local exon 50000 50211 . - . insertions 0 ; deletions 0 JXSJ01000007.1:subseq(6430332,103396) exonerate:protein2genome:local similarity 50000 53396 1933 - . alignment_id 1 ; Query SPP00000629_2.0 ; Align 53397 1 96 ; Align 53190 33 96 ; Align 52779 66 102 ; Align 52586 100 108 ; Align 52279 136 153 ; Align 52027 188 90 ; Align 51853 218 99 ; Align 51162 252 210 ; Align 50210 323 210 # --- END OF GFF DUMP --- # -- completed exonerate analysis