Command line: [exonerate --exhaustive yes -m p2g --showtargetgff -q /export/home/u114420/treball/fastas/SEPHS2.x.fa -t SEPHS2.extraction.KV884875.1.fa] Hostname: [sitdoc] C4 Alignment: ------------ Query: SEPHS2 Protein # Selenophosphate synthetase 2 (SEPHS2) # Zebrafish Target: KV884875.1:subseq(624694,100000) Monopterus albus unplaced genomic scaffold scaffold187.1, whole genome shotgun sequence Model: protein2genome:local Raw score: 1506 Query range: 7 -> 435 Target range: 36362 -> 51219 8 : SerGlyAspSerGlyGly >>>> Target Intron 1 >>>> AlaGlyGlyTyrA : 18 |||||| !!!!!||| ! 9178 bp !.!! !|||! !! SerGlyHisThrGlyGln++ ++ValValGlySerG 36363 : TCAGGCCACACAGGGCAGgt.........................agGTTGTGGGATCCG : 45571 19 : laValPheTyrProAlaGlyTyrGlnAlaLeuAsnProGluGluHisGlyLeuAspArgG : 38 !!.!! !!:!||| !!||||||:!! !!!!!!!.||||||||||||||||||!!:|||| luAlaCysPheProProGlyTyrLysProPheLysProGluGluHisGlyLeuGluArgG 45572 : AAGCGTGTTTCCCCCCTGGATACAAGCCTTTCAAACCCGAGGAGCACGGCCTGGAGCGCG : 45631 39 : lyPheArgLeuThrAlaPheSerAspMetLysGlyUnkGlyCysLysValProGlnGluT : 58 ||||||||||||||:!!|||||||||||||||||| |||||||||||||||||||||. lyPheArgLeuThrSerPheSerAspMetLysGly***GlyCysLysValProGlnGluA 45632 : GGTTCCGTCTCACATCCTTTTCGGACATGAAGGGATGAGGCTGCAAGGTCCCGCAGGAGG : 45691 59 : hrLeuLeuLysLeuLeuGlnGlyLeuGluProAspArgProProGlyGluAspGlyGlyL : 78 !!||||||||||||||| !||||||:!!||||||!:! !.!! !: laLeuLeuLysLeuLeuAlaGlyLeuGlnProAspGln---------------AlaAspM 45692 : CTCTGCTCAAACTCCTGGCGGGACTGCAGCCGGACCAG---------------GCCGACA : 45736 79 : euGlyThrGlyValGlyAspGluThrAlaAspPheGlyLeuValSerAlaAlaGlnGlyP : 98 !! ! ! !!.!||||||:!!! !:!!!!:||| ! ..! ! !! |||| etThrGlyLysAlaGlyAspGlnIleSerGluPheHis------GlnGlnLeuProGlyP 45737 : TGACAGGGAAGGCCGGGGACCAAATATCGGAGTTCCAC------CAACAGTTGCCCGGAC : 45790 99 : roArgLeu{G} >>>> Target Intron 2 >>>> {ly}IleGlyMetAspSerC : 107 ||||||||{|} 996 bp {||}:!!|||||||||!!!| roArgLeu{G}++ ++{ly}ValGlyMetAspCysC 45791 : CCCGTCTC{G}gt.........................ag{GT}GTAGGAATGGACTGCT : 46813 108 : ysValIleProLeuArgHisGlyGlyLeuSerLeuValGlnThrThrAspPhePheTyrP : 127 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| ysValIleProLeuArgHisGlyGlyLeuSerLeuValGlnThrThrAspPhePheTyrP 46814 : GTGTGATTCCCCTGAGACACGGAGGGCTCTCACTGGTCCAGACCACAGACTTCTTCTACC : 46873 128 : roLeuValGluAspProTyrMetMet >>>> Target Intron 3 >>>> GlyAr : 137 |||||||||||||||||||||||||| 606 bp ||||| roLeuValGluAspProTyrMetMet++ ++GlyAr 46874 : CTCTGGTGGAGGATCCCTACATGATGgt.........................agGGCAG : 47509 138 : gIleAlaCysAlaAsnValLeuSerAspLeuTyrAlaMetGlyIleThrGluCysAspAs : 157 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| gIleAlaCysAlaAsnValLeuSerAspLeuTyrAlaMetGlyIleThrGluCysAspAs 47510 : GATAGCGTGCGCCAACGTCCTCAGTGATCTCTACGCCATGGGCATCACCGAGTGTGACAA : 47569 158 : nMetLeuMetLeuLeuSerValSerGlnLysMetAsnGluLys >>>> Target Int : 172 ||||||||||||||||||||||||||||!:!||||||!!:||| 213 b nMetLeuMetLeuLeuSerValSerGlnArgMetAsnAspLys++ 47570 : CATGCTAATGCTGCTGAGTGTCAGCCAGAGGATGAATGACAAGgt............... : 47616 173 : ron 4 >>>> GluArgAspLeuValLeuProLeuMetMetLysGlyPheArgAspAla : 187 p !!:|||!!:! !|||:!!|||||||||!!:::!||||||||||||||| ++AspArgGluArgValMetProLeuMetIleArgGlyPheArgAspAla 47617 : ..........agGACCGTGAGCGGGTGATGCCACTTATGATTCGAGGCTTTCGAGACGCA : 47872 188 : AlaGluGluGlyGlyThrSerValThrGlyGlyGlnThrValIleAsnProTrpIleIle : 207 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| AlaGluGluGlyGlyThrSerValThrGlyGlyGlnThrValIleAsnProTrpIleIle 47873 : GCAGAGGAGGGTGGGACGTCCGTGACGGGTGGCCAGACGGTTATCAACCCCTGGATCATC : 47932 208 : IleGlyGlyValAlaSerValValCysGlnProAsnAspPheIle{Le} >>>> Targ : 223 ||||||||||||||||||||||||||||||||||||!!:||||||{:!} IleGlyGlyValAlaSerValValCysGlnProAsnGluPheIle{Me}++ 47933 : ATAGGAGGAGTAGCATCAGTCGTGTGCCAGCCCAACGAGTTCATC{AT}gt......... : 47984 224 : et Intron 5 >>>> {u}ProAspGlyAlaValProGlyAspValLeuValLeuThr : 236 598 bp {!}||||||||||||||||||||||||||||||||||||||| ++{t}ProAspGlyAlaValProGlyAspValLeuValLeuThr 47985 : ................ag{G}CCGGATGGTGCTGTACCTGGTGACGTCCTGGTCCTGACT : 48617 237 : LysProLeuGlyThrGlnValAlaValAsnAlaHisGlnTrpLeuAspIle >>>> Ta : 254 |||||||||||||||||||||||||||||||||||||||||||||||| LysProLeuGlyThrGlnValAlaValAsnAlaHisGlnTrpLeuAspGln++ 48618 : AAACCTCTGGGGACACAGGTCGCTGTAAATGCTCACCAGTGGCTTGATCAGgt....... : 48673 255 : rget Intron 6 >>>> ProGluLysTrpAsnLysIleLysLeuValIleSerArgG : 267 1025 bp ||||||!:!|||||||||||||||||||||:!!!!!!:!| ++ProGluArgTrpAsnLysIleLysLeuValValThrLysG 48674 : ..................agCCTGAAAGGTGGAACAAGATCAAACTAGTTGTCACCAAAG : 49735 268 : luGluValGluGlnAlaTyrGlnGluAlaMetLeuAsnMetAlaThrLeuAsnArgThr{ : 287 ||||||||:!!:!!||||||||||||||||||!!!::!|||||||||||||||||||||{ luGluValLysGluAlaTyrGlnGluAlaMetPheSerMetAlaThrLeuAsnArgThr{ 49736 : AGGAAGTTAAAGAGGCCTATCAGGAAGCCATGTTCTCCATGGCAACGCTGAATCGCACA{ : 49795 288 : A} >>>> Target Intron 7 >>>> {la}AlaAlaLeuMetHisLysPheAsnA : 296 |} 205 bp {||}|||||||||||||||||||||.!.| A}++ ++{la}AlaAlaLeuMetHisLysPheGlnA 49796 : G}gt.........................ag{CG}GCGGCCCTAATGCACAAGTTCCAGG : 50027 297 : laHisAlaAlaThrAspIleThrGlyPheGlyIleIleGlyHisAlaArgAsnLeuAlaL : 316 |||||||||||||||||:!:||||||||||||:!::!!|||||||||!.!|||||||||! laHisAlaAlaThrAspValThrGlyPheGlyLeuLeuGlyHisAlaAsnAsnLeuAlaT 50028 : CCCACGCTGCAACAGATGTGACAGGCTTTGGGTTGTTAGGGCACGCCAACAACCTGGCAA : 50087 317 : ysGlnGlnLysAsnAspValAlaPheValIleHisAsnLeuProIleIleSerLysMetA : 336 !||||||:!!|||!!:|||||||||||||||||||||||||||||||||:!!||||||| hrGlnGlnGlnAsnGluValAlaPheValIleHisAsnLeuProIleIleAlaLysMetA 50088 : CACAGCAACAGAATGAGGTGGCCTTTGTCATCCACAACCTACCAATCATAGCCAAGATGG : 50147 337 : laAlaIleSerLysAlaGlyGlyAsnLeuPheGlyLeuLeuGlnGlyThrSerSerGluT : 356 ||||||||||||||||| !!|||||| !|||.!!||||||||||||! !|||||||||| laAlaIleSerLysAlaCysGlyAsnAlaPheSerLeuLeuGlnGlyArgSerSerGluT 50148 : CCGCCATCAGTAAAGCCTGTGGAAACGCGTTCAGCCTGCTTCAGGGCAGATCATCAGAGA : 50207 357 : hrSer{G} >>>> Target Intron 8 >>>> {ly}GlyLeuLeuIleCysLeuP : 365 |||||{|} 791 bp {||}|||||||||:!:||||||| hrSer{G}++ ++{ly}GlyLeuLeuValCysLeuP 50208 : CTTCT{G}gt.........................ag{GT}GGGCTGCTGGTGTGTTTGC : 51025 366 : roArgGluGlnAlaAlaArgPheCysAlaGluMetLysSerSerArgMetGlyLeuLeuG : 385 |||||||||||||||||!:!||||||:!!||||||||||||..!!!! |||:!! !| roArgGluGlnAlaAlaLysPheCysSerGluMetLysSerGlnSerSerGlyValGlyG 51026 : CCAGAGAACAGGCAGCCAAGTTTTGCTCCGAAATGAAGAGCCAGAGCTCTGGCGTTGGAG : 51085 386 : lyAlaGlyGlnAspGlyGlyValGlyAspGlyGlnGlnAlaTrpIleIleGlyIleValG : 405 || !||| !! |||||| |||||||||||||||||||||| lyGlnGlyThrAlaGlyGly------------------AlaTrpIleIleGlyIleValG 51086 : GCCAGGGAACAGCGGGCGGA------------------GCATGGATAATTGGAATCGTAG : 51127 406 : luLysGlyAsnArgCysAlaArgIleIleAspLysProArgIleIleGluValProTyrA : 425 ||||||||:!!||| !|||||||||||||||||||||||||||||||||||| !! | luLysGlyAspArgHisAlaArgIleIleAspLysProArgIleIleGluValAlaProA 51128 : AGAAAGGCGACCGCCACGCTCGCATCATTGACAAGCCTCGTATTATTGAGGTTGCACCAC : 51187 426 : rgGlySerValValSerValGlnGluGlySer : 435 |||||||| !!.!:!! !||||||..!||| rgGlySerGlnAlaAlaAsnGlnGluAsnSer 51188 : GAGGAAGTCAGGCAGCGAATCAGGAGAACAGC : 51219 vulgar: SEPHS2 7 435 . KV884875.1:subseq(624694,100000) 36362 51219 + 1506 M 6 18 5 0 2 I 0 9174 3 0 2 M 57 171 G 5 0 M 15 45 G 2 0 M 8 24 S 0 1 5 0 2 I 0 992 3 0 2 S 1 2 M 34 102 5 0 2 I 0 602 3 0 2 M 36 108 5 0 2 I 0 209 3 0 2 M 51 153 S 0 2 5 0 2 I 0 594 3 0 2 S 1 1 M 30 90 5 0 2 I 0 1021 3 0 2 M 33 99 S 0 1 5 0 2 I 0 201 3 0 2 S 1 2 M 70 210 S 0 1 5 0 2 I 0 787 3 0 2 S 1 2 M 33 99 G 6 0 M 38 114 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2017-11-23 ##type DNA # # # seqname source feature start end score strand frame attributes # KV884875.1:subseq(624694,100000) exonerate:protein2genome:local gene 36363 51219 1506 + . gene_id 1 ; sequence SEPHS2 ; gene_orientation + KV884875.1:subseq(624694,100000) exonerate:protein2genome:local cds 36363 36380 . + . KV884875.1:subseq(624694,100000) exonerate:protein2genome:local exon 36363 36380 . + . insertions 0 ; deletions 0 KV884875.1:subseq(624694,100000) exonerate:protein2genome:local splice5 36381 36382 . + . intron_id 1 ; splice_site "GT" KV884875.1:subseq(624694,100000) exonerate:protein2genome:local intron 36381 45558 . + . intron_id 1 KV884875.1:subseq(624694,100000) exonerate:protein2genome:local splice3 45557 45558 . + . intron_id 0 ; splice_site "AG" KV884875.1:subseq(624694,100000) exonerate:protein2genome:local cds 45559 45799 . + . KV884875.1:subseq(624694,100000) exonerate:protein2genome:local exon 45559 45799 . + . insertions 0 ; deletions 7 KV884875.1:subseq(624694,100000) exonerate:protein2genome:local splice5 45800 45801 . + . intron_id 2 ; splice_site "GT" KV884875.1:subseq(624694,100000) exonerate:protein2genome:local intron 45800 46795 . + . intron_id 2 KV884875.1:subseq(624694,100000) exonerate:protein2genome:local splice3 46794 46795 . + . intron_id 1 ; splice_site "AG" KV884875.1:subseq(624694,100000) exonerate:protein2genome:local cds 46796 46899 . + . KV884875.1:subseq(624694,100000) exonerate:protein2genome:local exon 46796 46899 . + . insertions 0 ; deletions 0 KV884875.1:subseq(624694,100000) exonerate:protein2genome:local splice5 46900 46901 . + . intron_id 3 ; splice_site "GT" KV884875.1:subseq(624694,100000) exonerate:protein2genome:local intron 46900 47505 . + . intron_id 3 KV884875.1:subseq(624694,100000) exonerate:protein2genome:local splice3 47504 47505 . + . intron_id 2 ; splice_site "AG" KV884875.1:subseq(624694,100000) exonerate:protein2genome:local cds 47506 47613 . + . KV884875.1:subseq(624694,100000) exonerate:protein2genome:local exon 47506 47613 . + . insertions 0 ; deletions 0 KV884875.1:subseq(624694,100000) exonerate:protein2genome:local splice5 47614 47615 . + . intron_id 4 ; splice_site "GT" KV884875.1:subseq(624694,100000) exonerate:protein2genome:local intron 47614 47826 . + . intron_id 4 KV884875.1:subseq(624694,100000) exonerate:protein2genome:local splice3 47825 47826 . + . intron_id 3 ; splice_site "AG" KV884875.1:subseq(624694,100000) exonerate:protein2genome:local cds 47827 47981 . + . KV884875.1:subseq(624694,100000) exonerate:protein2genome:local exon 47827 47981 . + . insertions 0 ; deletions 0 KV884875.1:subseq(624694,100000) exonerate:protein2genome:local splice5 47982 47983 . + . intron_id 5 ; splice_site "GT" KV884875.1:subseq(624694,100000) exonerate:protein2genome:local intron 47982 48579 . + . intron_id 5 KV884875.1:subseq(624694,100000) exonerate:protein2genome:local splice3 48578 48579 . + . intron_id 4 ; splice_site "AG" KV884875.1:subseq(624694,100000) exonerate:protein2genome:local cds 48580 48670 . + . KV884875.1:subseq(624694,100000) exonerate:protein2genome:local exon 48580 48670 . + . insertions 0 ; deletions 0 KV884875.1:subseq(624694,100000) exonerate:protein2genome:local splice5 48671 48672 . + . intron_id 6 ; splice_site "GT" KV884875.1:subseq(624694,100000) exonerate:protein2genome:local intron 48671 49695 . + . intron_id 6 KV884875.1:subseq(624694,100000) exonerate:protein2genome:local splice3 49694 49695 . + . intron_id 5 ; splice_site "AG" KV884875.1:subseq(624694,100000) exonerate:protein2genome:local cds 49696 49795 . + . KV884875.1:subseq(624694,100000) exonerate:protein2genome:local exon 49696 49795 . + . insertions 0 ; deletions 0 KV884875.1:subseq(624694,100000) exonerate:protein2genome:local splice5 49796 49797 . + . intron_id 7 ; splice_site "GT" KV884875.1:subseq(624694,100000) exonerate:protein2genome:local intron 49796 50000 . + . intron_id 7 KV884875.1:subseq(624694,100000) exonerate:protein2genome:local splice3 49999 50000 . + . intron_id 6 ; splice_site "AG" KV884875.1:subseq(624694,100000) exonerate:protein2genome:local cds 50001 50213 . + . KV884875.1:subseq(624694,100000) exonerate:protein2genome:local exon 50001 50213 . + . insertions 0 ; deletions 0 KV884875.1:subseq(624694,100000) exonerate:protein2genome:local splice5 50214 50215 . + . intron_id 8 ; splice_site "GT" KV884875.1:subseq(624694,100000) exonerate:protein2genome:local intron 50214 51004 . + . intron_id 8 KV884875.1:subseq(624694,100000) exonerate:protein2genome:local splice3 51003 51004 . + . intron_id 7 ; splice_site "AG" KV884875.1:subseq(624694,100000) exonerate:protein2genome:local cds 51005 51219 . + . KV884875.1:subseq(624694,100000) exonerate:protein2genome:local exon 51005 51219 . + . insertions 0 ; deletions 6 KV884875.1:subseq(624694,100000) exonerate:protein2genome:local similarity 36363 51219 1506 + . alignment_id 1 ; Query SEPHS2 ; Align 36363 8 18 ; Align 45559 14 171 ; Align 45730 76 45 ; Align 45775 93 24 ; Align 46798 102 102 ; Align 47506 136 108 ; Align 47827 172 153 ; Align 48581 224 90 ; Align 49696 254 99 ; Align 50003 288 210 ; Align 51007 359 99 ; Align 51106 398 114 # --- END OF GFF DUMP --- # -- completed exonerate analysis