Command line: [exonerate --exhaustive yes -m p2g --showtargetgff -q /export/home/u114420/treball/fastas/SELENOP2.x.fa -t SELENOP2.extraction.KV884708.1.fa] Hostname: [sitdoc] C4 Alignment: ------------ Query: SELENOP2 # Protein # Selenoprotein P (SELENOP) # Zebrafish Target: KV884708.1:subseq(991377,100000) Monopterus albus unplaced genomic scaffold scaffold20.1, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 641 Query range: 0 -> 348 Target range: 50563 -> 19856 1 : MetTrpLysAlaLeuSerLeuThrLeuAlaLeuCysLeuLeuValGlyCysSerAlaGlu : 20 |||||| !..!||||||||| ||| |||||||||||| !||| !!!!!|||||| MetTrpAlaCysLeuSerLeuLeuLeuProLeuCysLeuLeuHisGlyGlyArgAlaGlu 50563 : ATGTGGGCGTGCCTCAgccttctcctccctctctgcctgctCCATGGGGGCAGAGCAGAG : 50506 21 : SerGluThrGluGlyAlaArgCysLysLeuProProGluTrpLysValGlyAspValGlu : 40 |||||| !! ||| !!||||||:!!|||||||||! |||!!.:!!|||!!:|||||| SerGluGlyValGlyProArgCysGlnLeuProProAlaTrpAsnIleGlyGluValGlu 50505 : AGTGAGGGGGTTGGGCCTCGCTGTCAACTGCCGCCAGCCTGGAACATAGGGGAGGTGGAG : 50446 41 : ProMetLysAsnAlaLeuGlyGlnValThrValValAlaTyrLeuGlnAlaSerUnkLeu : 60 |||||||||....!!:!!||||||||||||||||||||| !|||||||||||| ||| ProMetLysGlyThrMetGlyGlnValThrValValAlaLeuLeuGlnAlaSerCysLeu 50445 : CCCATGAAGGGGACAATGGGCCAGGTGACGGTGGTGGCTCTTTTGCAGGCCAGCTGCCTG : 50386 61 : PheCysLeuGluGlnAlaSer{Ly} >>>> Target Intron 1 >>>> {s}Leu : 69 |||||||||! !|||||||||{!:} 377 bp {!}:!! PheCysLeuValGlnAlaSer{Ar}++ ++{g}Met 50385 : TTCTGCCTGGTGCAGGCTTCC{AG}gt.........................ag{A}ATG : 49982 70 : AsnAspLeuLeuLeuLysLeuGluLysGlnGlyTyrProAsnIleAlaTyrMetValVal : 89 :!!..!|||! !! !! !:!!|||!..|||||| !|||:!:||||||||||||:!! AspSerLeuHisGlnMetMetGluSerGlnGlyLeuLysAsnValAlaTyrMetValIle 49981 : GACAGCCTGCACCAGATGATGGAGAGTCAGGGCCTGAAGAATGTGGCTTACATGGTGATT : 49922 90 : AsnAsnArgGluGluArgSerGlnArgLeuHisHisLeuLeuGlnGluArgLeuLeu<-> : 109 |||:!!!:!! !|||!:!:!!|||! !||||||! !:!!||| !:!!|||:!!! ! AsnHisGlnGlyGluGlnAlaGlnLeuLeuHisProMetLeuAlaGlnArgMetSerGlu 49921 : AACCACCAGGGGGAGCAAGCACAGCTCCTGCACCCTATGCTGGCACAGAGAATGTCAGAG : 49862 110 : AsnIleThrLeuTyrAlaGlnAspLeuSerGlnProAspAlaTrpGlnAlaValAsnAla : 128 ||||||||||||||| !||||||! !..!||| !|||!.!|||:!!.!!:!!..!!.! AsnIleThrLeuTyrLysGlnAspGlnGlnGlnValAspValTrpLysThrLeuGlyGly 49861 : AACATCACACTGTACAAACAGGACCAGCAACAGGTTGATGTCTGGAAGACACTGGGTGGA : 49802 129 : GluLysAspAspIleLeuValTyrAsp{Ar} >>>> Target Intron 2 >>>> : 138 :!!|||||||||.!!|||:!!||||||{||} 159 bp GlnLysAspAspPheLeuIleTyrAsp{Ar}++ ++ 49801 : CAAAAAGATGACTTTCTCATTTATGAC{AG}gt.........................ag : 49613 139 : {g}CysGlyArgLeuThrTyrHisLeuSerLeuProTyrThrIleLeuIleHisProHis : 157 {|}|||||||||||||||:!!|||:!!||||||||||||:!!|||:!! !!!. !||| {g}CysGlyArgLeuThrHisHisIleSerLeuProTyrSerIleIleGlyGlnGlyHis 49612 : {A}TGTGGCCGTCTCACCCATCACATTTCTCTTCCATACTCCATCATTGGACAAGGCCAT : 49556 158 : ValGluGluAlaIleLysHisThrTyrCysAspArgIleCysGlyGluCysSer{L} > : 176 ||||||! ||||||!:! !!|||||||||:!!||||||||||||!!:|||!!!{!} ValGluGlyAlaIleArgAspThrTyrCysAsnArgIleCysGlyAspCysThr{G}-+ 49555 : GTTGAGGGCGCCATCAGAGATACTTACTGTAATCGCATATGTGGGGACTGCACA{C}at. : 49496 177 : >>> Target Intron 3 >>>> {eu}GluSerSerAlaGlnLeuGluGluCysLys : 186 184 bp { !}|||!!!!.!!.!:!! !!!:! ! ++{ln}GluCysLysGlyLysAlaAspValGlnPro 49495 : ........................ag{AA}GAGTGCAAAGGGAAAGCAGATGTACAGCCT : 49285 187 : LysAlaThrGluGluValAsnLysProValGluGluGluProArgGlnAspHisGlyHis : 206 ! ||| !! ! !! ! ! !.!. ! !! ! !!!!. !!||||||||| AspAlaAspGlyThrProAlaIleGluHisAsnThrGlyHisGlyHisHisHisGlyHis 49284 : GATGCGGATGGTACCCCAGCTATAGAGCACAACACTGGACATGGTCACCATCATGGACAT : 49225 207 : HisGluGlnGlyHisHisGluHisGlnGlyGluAlaGluArgHisArgHisGlyHisHis : 226 |||.!. ! !|||! |||!!.|||!!: !... !! !!.!||||||||| ! HisHisGlyHisGlyHisGlyHisHisGlyAspAsnSerGlyPheHisHisGlyHisGly 49224 : CACCATGGTCATGGCCATGGCCACCATGGGGATAATAGCGGTTTTCATCATGGTCATGGC : 49165 227 : HisProHisHisHisHisHisHisHisArg{G} >>>> Target Intron 4 >>>> : 237 ||| !|||:!!||| !||||||||| !!{|} 78 bp HisAspHisAsnHisGlyHisHisHisGly{G}++ 49164 : CATGATCATAACCATGGCCATCACCATGGA{G}gt......................... : 49129 238 : {ly}Gln<->GlnGlnValAspValAspGlnGlnValLeuSerGlnValAspPheGly : 254 {||}||| ||||||!.!! :!! !||||||:!!|||..!:!!!.! !! !.! ++{ly}GlnMetGlnGlnAlaValLeuThrGlnGlnMetLeuGlnGluAlaHisArgAla 49128 : ag{GC}CAGATGCAGCAAGCCGTGCTCACTCAACAGATGTTACAGGAGGCCCACAGAGCC : 49000 255 : GlnValAlaValGluThrProMetMetLysArgProUnkAlaLysHisSerArgUnkLys : 274 ! ||| ! ! !.!! !! !! ! !||| ! :!!||||||||| !! ! ProValArgPro***AlaSerGluArgAlaArg***LysSerLysHisSer***GlnTrp 48999 : CCTGTAAGACCTTGAGCATCAGAGAGGGCAAGGTGAAAGTCAAAGCACAGCTGACAGTGG : 48940 275 : ValGlnTyrSerUnkGlnGlnGlyAlaAspSerProValAlaSerUnkCysUnkHisUnk : 294 ..! ! !||| .!.:!!!.! |||||| ! !||| ||| ||| ThrAlaGlySerAspAsnGluAla------SerProLysIleSer***Cys***His*** 48939 : ACAGCAGGCTCTGACAATGAAGCC------TCTCCTAAGATCAGCTGATGCTGACACTGA : 48886 295 : ArgGlnLeuPheGlyGlyGluGlyAsnGlyArgValAlaGlyLeuUnkHisCysAspGlu : 314 |||::!|||||||||! !! !|||!:!! !!:! !!.!|||||| ||||||:!!||| ArgArgLeuPheGlyAspAlaGlySerGluGlnProValGlyLeu***HisCysAsnGlu 48885 : CGCAGGCTGTTTGGCGATGCAGGGAGTGAACAGCCGGTCGGTCTCTGACACTGTAATGAG : 48826 315 : ProLeuProAlaSerUnkProUnkGlnGlyLeuLys{G} >>>> Target Intron : 327 !!|||||||||||| ! ! !!.||||||! {|} 28865 bp AlaLeuProAlaSer***Gln***HisGlyLeuIle{G}++ 48825 : GCGCTGCCCGCCTCCTGACAGTGACACGGACTGATT{G}gt................... : 48784 328 : 5 >>>> {lu}GlnAspAsnHisIleLysGluThrUnkGlnUnkArgProAlaProPro : 343 {||}::: ...|||::: ||| |||... ... ||| ++{lu}GluProGlyHisLeuProGluArgLeuGlnAlaSerGlyThrHisPro 48783 : ......ag{aa}gagcctggtcacctgcctgaacgactacaggccagtggcactcacccc : 19874 344 : AlaGluUnkGluLeu : 348 |||::: AsnTyrTyrGluVal 19873 : aattattatgaagtg : 19857 vulgar: SELENOP2 0 348 . KV884708.1:subseq(991377,100000) 50563 19856 - 641 M 67 201 S 0 2 5 0 2 I 0 373 3 0 2 S 1 1 M 40 120 G 0 3 M 29 87 S 0 2 5 0 2 I 0 155 3 0 2 S 1 1 M 37 111 S 0 1 5 0 2 I 0 180 3 0 2 S 1 2 M 60 180 S 0 1 5 0 2 I 0 74 3 0 2 S 1 2 M 1 3 G 0 3 M 44 132 G 2 0 M 42 126 S 0 1 5 0 2 I 0 28861 3 0 2 S 1 2 M 21 63 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2017-11-22 ##type DNA # # # seqname source feature start end score strand frame attributes # KV884708.1:subseq(991377,100000) exonerate:protein2genome:local gene 19857 50563 641 - . gene_id 1 ; sequence SELENOP2 ; gene_orientation + KV884708.1:subseq(991377,100000) exonerate:protein2genome:local cds 50361 50563 . - . KV884708.1:subseq(991377,100000) exonerate:protein2genome:local exon 50361 50563 . - . insertions 0 ; deletions 0 KV884708.1:subseq(991377,100000) exonerate:protein2genome:local splice5 50359 50360 . - . intron_id 1 ; splice_site "GT" KV884708.1:subseq(991377,100000) exonerate:protein2genome:local intron 49984 50360 . - . intron_id 1 KV884708.1:subseq(991377,100000) exonerate:protein2genome:local splice3 49984 49985 . - . intron_id 0 ; splice_site "AG" KV884708.1:subseq(991377,100000) exonerate:protein2genome:local cds 49771 49983 . - . KV884708.1:subseq(991377,100000) exonerate:protein2genome:local exon 49771 49983 . - . insertions 3 ; deletions 0 KV884708.1:subseq(991377,100000) exonerate:protein2genome:local splice5 49769 49770 . - . intron_id 2 ; splice_site "GT" KV884708.1:subseq(991377,100000) exonerate:protein2genome:local intron 49612 49770 . - . intron_id 2 KV884708.1:subseq(991377,100000) exonerate:protein2genome:local splice3 49612 49613 . - . intron_id 1 ; splice_site "AG" KV884708.1:subseq(991377,100000) exonerate:protein2genome:local cds 49499 49611 . - . KV884708.1:subseq(991377,100000) exonerate:protein2genome:local exon 49499 49611 . - . insertions 0 ; deletions 0 KV884708.1:subseq(991377,100000) exonerate:protein2genome:local splice5 49497 49498 . - . intron_id 3 ; splice_site "AT" KV884708.1:subseq(991377,100000) exonerate:protein2genome:local intron 49315 49498 . - . intron_id 3 KV884708.1:subseq(991377,100000) exonerate:protein2genome:local splice3 49315 49316 . - . intron_id 2 ; splice_site "AG" KV884708.1:subseq(991377,100000) exonerate:protein2genome:local cds 49132 49314 . - . KV884708.1:subseq(991377,100000) exonerate:protein2genome:local exon 49132 49314 . - . insertions 0 ; deletions 0 KV884708.1:subseq(991377,100000) exonerate:protein2genome:local splice5 49130 49131 . - . intron_id 4 ; splice_site "GT" KV884708.1:subseq(991377,100000) exonerate:protein2genome:local intron 49054 49131 . - . intron_id 4 KV884708.1:subseq(991377,100000) exonerate:protein2genome:local splice3 49054 49055 . - . intron_id 3 ; splice_site "AG" KV884708.1:subseq(991377,100000) exonerate:protein2genome:local cds 48787 49053 . - . KV884708.1:subseq(991377,100000) exonerate:protein2genome:local exon 48787 49053 . - . insertions 3 ; deletions 2 KV884708.1:subseq(991377,100000) exonerate:protein2genome:local splice5 48785 48786 . - . intron_id 5 ; splice_site "GT" KV884708.1:subseq(991377,100000) exonerate:protein2genome:local intron 19922 48786 . - . intron_id 5 KV884708.1:subseq(991377,100000) exonerate:protein2genome:local splice3 19922 19923 . - . intron_id 4 ; splice_site "ag" KV884708.1:subseq(991377,100000) exonerate:protein2genome:local cds 19857 19921 . - . KV884708.1:subseq(991377,100000) exonerate:protein2genome:local exon 19857 19921 . - . insertions 0 ; deletions 0 KV884708.1:subseq(991377,100000) exonerate:protein2genome:local similarity 19857 50563 641 - . alignment_id 1 ; Query SELENOP2 ; Align 50564 1 201 ; Align 49983 69 120 ; Align 49860 109 87 ; Align 49611 139 111 ; Align 49313 177 180 ; Align 49052 238 3 ; Align 49046 239 132 ; Align 48914 285 126 ; Align 19920 328 63 # --- END OF GFF DUMP --- # -- completed exonerate analysis