Command line: [exonerate -m p2g --showtargetgff -q ../SelP-human.fa -t SelP-NIVO01048714.1.subseq] Hostname: [sitdoc] C4 Alignment: ------------ Query: SPP00000020_1.0 # Protein # Selenoprotein P (SelP) # Homo sapiens # Complete Target: NIVO01048714.1:subseq(460434,105985) Ammotragus lervia isolate NAZ scaffold4518_len1342136_cov72, whole genome shotgun sequence Model: protein2genome:local Raw score: 1305 Query range: 0 -> 381 Target range: 49999 -> 55985 1 : MetTrpArgSerLeuGlyLeuAlaLeuAlaLeuCysLeuLeuProSerGlyGlyThrGlu : 20 |||||||||.!!|||||||||||||||||||||||||||||| !!!! !!||||||||| MetTrpArgGlyLeuGlyLeuAlaLeuAlaLeuCysLeuLeuLeuThrArgGlyThrGlu 50000 : ATGTGGAGAGGCCTGGGGCTGGCTCtggctctctgcctcctcctAACCAGAGGAACAGAG : 50057 21 : SerGlnAspGlnSerSerLeuCysLysGlnProProAlaTrpSerIleArgAspGlnAsp : 40 ||||||! !|||||||||! !|||||||||||||||||||||||||||!:!||||||||| SerGlnGlyGlnSerSerTyrCysLysGlnProProAlaTrpSerIleLysAspGlnAsp 50058 : AGCCAGGGCCAAAGCTCCTATTGTAAGCAACCTCCAGCCTGGAGCATAAAGGATCAAGAC : 50117 41 : ProMetLeuAsnSerAsnGlySerValThrValValAlaLeuLeuGlnAlaSerUnkTyr : 60 ||||||||||||||| !!|||||||||||||||||||||||||||||||||||| ||| ProMetLeuAsnSerTyrGlySerValThrValValAlaLeuLeuGlnAlaSer***Tyr 50118 : CCTATGCTGAACTCCTATGGTTCAGTGACCGTGGTTGCTCTTCTTCAAGCCAGCTGATAC : 50177 61 : LeuCysIleLeuGlnAlaSer{Ly} >>>> Target Intron 1 >>>> {s}Leu : 69 |||||||||||||||||||||{!:} 1001 bp {!}||| LeuCysIleLeuGlnAlaSer{Ar}++ ++{g}Leu 50178 : CTGTGCATTCTGCAGGCCTCT{AG}gt.........................ag{a}tTG : 51205 70 : GluAspLeuArgValLysLeuLysLysGluGlyTyrSerAsnIleSerTyrIleValVal : 89 |||||||||||||||||||||:!!|||||||||||||||||||||||||||:::|||||| GluAspLeuArgValLysLeuGluLysGluGlyTyrSerAsnIleSerTyrValValVal 51206 : GAGGACCTGCGAGTAAAACTGGAGAAAGAAGGATATTCTAAcatctcttatgttgttgtt : 51265 90 : AsnHisGlnGlyIleSerSerArgLeuLysTyrThrHisLeuLysAsnLysValSerGlu : 109 |||||||||! !|||||||||||||||||||||..!|||||||||||||||||||||||| AsnHisGlnGluIleSerSerArgLeuLysTyrValHisLeuLysAsnLysValSerGlu 51266 : AATCATCAAGAAATCTCTTCTCGATTAAAATATGTGCATCTTAAAAATAAGGTTTCAGAA : 51325 110 : HisIleProValTyrGlnGlnGluGluAsnGlnThrAspValTrpThrLeuLeuAsnGly : 129 ||||||||||||||||||||||||!!::!!||| !!|||||||||||||||||||||||| HisIleProValTyrGlnGlnGluAspAspGlnProAspValTrpThrLeuLeuAsnGly 51326 : CATATTCCTGTTTACCAACAAGAAGACGACCAACCAGATGTCTGGACTCtcttaaatgga : 51385 130 : SerLysAspAspPheLeuIleTyrAsp{Ar} >>>> Target Intron 2 >>>> : 139 :::||||||||||||||||||||||||{||} 1982 bp AsnLysAspAspPheLeuIleTyrAsp{Ar}++ ++ 51386 : aataaagatgACTTCCTCATATATGAC{AG}gt.........................ag : 53397 140 : {g}CysGlyArgLeuValTyrHisLeuGlyLeuProPheSerPheLeuThrPheProTyr : 158 {|}|||||||||||||||||||||||||||||||||!:!||||||||||||||| !!||| {g}CysGlyArgLeuValTyrHisLeuGlyLeuProTyrSerPheLeuThrPheThrTyr 53398 : {A}TGTGGCCGCCTTGTATATCATCTTGGTTTGCCTTATTCCTTCCTAACTTTCACATAT : 53454 159 : ValGluGluAlaIleLysIleAlaTyrCysGluLysLysCysGlyAsnCysSerLeuThr : 178 ||||||!!::!!||||||! !!.!||||||||| ! |||||||||||||||||||||! ! ValGluAspSerIleLysThrValTyrCysGluAspLysCysGlyAsnCysSerLeuLys 53455 : GTAGAAGATTCCATTAAGACTGTTTACTGTGAAGATAAATGTGGAAACTGCTCTCTCAAG : 53514 179 : >>>> Target Intron 3 >>>> ThrLeuLysAspGluAspPheCysLysArgV : 189 1773 bp ||||||:::||||||||||||||||||!.!| ++ ++ThrLeuGluAspGluAspPheCysLysAsnV 53515 : gt.........................agacactggaagatgaagacTTCTGTAAAAATG : 55320 190 : alSerLeuAlaThrValAspLysThrValGluThrProSerProHisTyrHisHisGluH : 209 |||||||||||||| !!!:||||||!.!|||.!! !!..! |||:::||||||...| alSerLeuAlaThrLysGluLysThrAlaGluAlaSerGlnArgHisHisHisHisHisH 55321 : TATCTCTGGCTACCAAGGAGAAAACAGCTGAGGCTTCACAGcgacatcaccatcaccatc : 55380 210 : is{H} >>>> Target Intron 4 >>>> {is}HisAsnHisGlyHisGlnHisL : 218 ||{|} 72 bp {||} :::||||||||||||! !! is{H}++ --{is}ProHisHisGlyHisGlnLeuH 55381 : ac{c}gt.........................cc{at}ccccatcaCGGGCATCAGCTTC : 55479 219 : euGlySerSerGluLeuSerGluAsnGlnGlnProGlyAlaProAsnAlaProThrHisP : 238 !! !!:!:!!.!.|||||||||::!! !.!.|||! !.!!|||:!!.!!||| !:!!| isGluAsnAlaHisLeuSerGluSerProAsnProAspThrProAspThrProGluAsnP 55480 : ATGAAAATGCTCATCTTTCAGAGTCTCCAAACCCAGACACACCAGATACCCCTGAGAATC : 55539 239 : roAlaProProGlyLeu<->HisHisHisHisLysHisLysGlyGlnHisArgGlnGlyH : 257 || !!||| !!|||||| ||||||||||||:::|||||||||!!.!!.|||||||||| roProProSerGlyLeuHisHisHisHisHisArgHisLysGlyHisGlnArgGlnGlyH 55540 : CCCCTCCTTCAGGtcttcatcaccaccaccataggCACAAGGGTCACCAAAGACAGGGTC : 55599 258 : isProGluAsnArgAspMetPro<->AlaSerGluAspLeuGln<->AspLeu<->GlnL : 274 || !!!!:||| !!|||! ||| !.!||||||:!!|||||| ..!||| |||| isSerAspAsnCysAspThrProValGlySerGluAsnLeuGlnLeuSerLeuProGlnL 55600 : ACTCAGATAACTGTGATACACCAGTAGGAAGTGAAAATTTACAACTTTCTCTTCCACAAA : 55659 275 : ysLysLeuCysArgLysArgCysIleAsnGlnLeuLeuCysLysLeuProThrAspSerG : 294 ||||||||!! |||||||||||||||||||||||||||!! :!!!!!|||! !||||||| ysLysLeu***ArgLysArgCysIleAsnGlnLeuLeu***GlnPheProLysAspSerG 55660 : AGAAGCTCTGACGAAAGAGATGCATAAATCAGTTACTCTGACAGTTTCCCAAAGATTCAG : 55719 295 : luLeuAlaProArgSerUnkCysCysHisCysArgHisLeuIlePheGluLysThrGlyS : 314 ||! !||| !!!!||| |||||||||||||||||||||:!!|||||||||||||||| luSerAlaLeuSerSerCysCysCysHisCysArgHisLeuValPheGluLysThrGlyS 55720 : AATCTGCTTTGAGTAGCTGCTGTTGCCACTGTCGACATCTGGTATTTGAAAAAACAGGGT : 55779 315 : erAlaIleThrUnkGlnCysLysGluAsnLeuProSerLeuCysSerUnkGlnGlyLeuA : 334 ||||||||||| ||||||! !|||||||||||||||||||||||| ||||||||| erAlaIleThr***GlnCysThrGluAsnLeuProSerLeuCysSer***GlnGlyLeuL 55780 : CTGCAATCACCTGACAGTGTACAGAAAACCTCCCCTCTTTATGTAGCTGACAGGGCCTTT : 55839 335 : rgAlaGluGluAsnIleThrGluSerCysGlnUnkArgLeuProProAlaAlaUnkGlnI : 354 !||||||||||||:!!! !||||||!! ||| |||||||||||||||||| ||| euAlaGluGluAsnValIleGluSer***Gln***ArgLeuProProAlaAla***GlnA 55840 : TGGCAGAGGAGAACGTCATTGAATCTTGACAGTGACGTTTGCCTCCAGCTGCCTGACAGG : 55899 355 : leSer<->GlnGlnLeuIleProThrGluAlaSerAlaSerUnkArgUnkLysAsnGlnA : 373 !:!! |||||||||! !|||||||||||||||.!!!.! !!! ||||||:!!| laAlaGlyGlnGlnLeuAsnProThrGluAlaSerThrLys***Ser***LysAsnLysA 55900 : CAGCAGGTCAGCAGCTCAATCCCACAGAAGCCAGCACCAAGTGAAGCTGAAAAAATAAGG : 55959 374 : laLysLysUnkGluUnkProSerAsn : 381 |||||! ! :!! ||||||||| laLysMet***Lys***ProSerAsn 55960 : CCAAAATGTGAAAATGACCTTcaaat : 55985 vulgar: SPP00000020_1.0 0 381 . NIVO01048714.1:subseq(460434,105985) 49999 55985 + 1305 M 67 201 S 0 2 5 0 2 I 0 997 3 0 2 S 1 1 M 70 210 S 0 2 5 0 2 I 0 1978 3 0 2 S 1 1 M 39 117 5 0 2 I 0 1769 3 0 2 M 31 93 S 0 1 5 0 2 I 0 68 3 0 2 S 1 2 M 33 99 G 0 3 M 21 63 G 0 3 M 6 18 G 0 3 M 2 6 G 0 3 M 83 249 G 0 3 M 26 78 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2017-11-28 ##type DNA # # # seqname source feature start end score strand frame attributes # NIVO01048714.1:subseq(460434,105985) exonerate:protein2genome:local gene 50000 55985 1305 + . gene_id 1 ; sequence SPP00000020_1.0 ; gene_orientation + NIVO01048714.1:subseq(460434,105985) exonerate:protein2genome:local cds 50000 50202 . + . NIVO01048714.1:subseq(460434,105985) exonerate:protein2genome:local exon 50000 50202 . + . insertions 0 ; deletions 0 NIVO01048714.1:subseq(460434,105985) exonerate:protein2genome:local splice5 50203 50204 . + . intron_id 1 ; splice_site "GT" NIVO01048714.1:subseq(460434,105985) exonerate:protein2genome:local intron 50203 51203 . + . intron_id 1 NIVO01048714.1:subseq(460434,105985) exonerate:protein2genome:local splice3 51202 51203 . + . intron_id 0 ; splice_site "ag" NIVO01048714.1:subseq(460434,105985) exonerate:protein2genome:local cds 51204 51416 . + . NIVO01048714.1:subseq(460434,105985) exonerate:protein2genome:local exon 51204 51416 . + . insertions 0 ; deletions 0 NIVO01048714.1:subseq(460434,105985) exonerate:protein2genome:local splice5 51417 51418 . + . intron_id 2 ; splice_site "GT" NIVO01048714.1:subseq(460434,105985) exonerate:protein2genome:local intron 51417 53398 . + . intron_id 2 NIVO01048714.1:subseq(460434,105985) exonerate:protein2genome:local splice3 53397 53398 . + . intron_id 1 ; splice_site "AG" NIVO01048714.1:subseq(460434,105985) exonerate:protein2genome:local cds 53399 53516 . + . NIVO01048714.1:subseq(460434,105985) exonerate:protein2genome:local exon 53399 53516 . + . insertions 0 ; deletions 0 NIVO01048714.1:subseq(460434,105985) exonerate:protein2genome:local splice5 53517 53518 . + . intron_id 3 ; splice_site "gt" NIVO01048714.1:subseq(460434,105985) exonerate:protein2genome:local intron 53517 55289 . + . intron_id 3 NIVO01048714.1:subseq(460434,105985) exonerate:protein2genome:local splice3 55288 55289 . + . intron_id 2 ; splice_site "ag" NIVO01048714.1:subseq(460434,105985) exonerate:protein2genome:local cds 55290 55383 . + . NIVO01048714.1:subseq(460434,105985) exonerate:protein2genome:local exon 55290 55383 . + . insertions 0 ; deletions 0 NIVO01048714.1:subseq(460434,105985) exonerate:protein2genome:local splice5 55384 55385 . + . intron_id 4 ; splice_site "gt" NIVO01048714.1:subseq(460434,105985) exonerate:protein2genome:local intron 55384 55455 . + . intron_id 4 NIVO01048714.1:subseq(460434,105985) exonerate:protein2genome:local splice3 55454 55455 . + . intron_id 3 ; splice_site "cc" NIVO01048714.1:subseq(460434,105985) exonerate:protein2genome:local cds 55456 55985 . + . NIVO01048714.1:subseq(460434,105985) exonerate:protein2genome:local exon 55456 55985 . + . insertions 15 ; deletions 0 NIVO01048714.1:subseq(460434,105985) exonerate:protein2genome:local similarity 50000 55985 1305 + . alignment_id 1 ; Query SPP00000020_1.0 ; Align 50000 1 201 ; Align 51205 69 210 ; Align 53400 140 117 ; Align 55290 179 93 ; Align 55458 211 99 ; Align 55560 244 63 ; Align 55626 265 18 ; Align 55647 271 6 ; Align 55656 273 249 ; Align 55908 356 78 # --- END OF GFF DUMP --- # -- completed exonerate analysis