Command line: [exonerate -m p2g --showtargetgff -q ../SelU1-human.fa -t SelU1-NIVO01064043.1.subseq] Hostname: [sitdoc] C4 Alignment: ------------ Query: SPP00000026_1.0 # Protein # Selenoprotein U1 (SelU1) # Homo sapiens # Complete Target: NIVO01064043.1:subseq(0,374980) Ammotragus lervia isolate NAZ scaffold4425_len372517_cov72, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 1036 Query range: 0 -> 229 Target range: 30496 -> 22805 1 : MetSerPheLeuGlnAspProSerPhePheThrMetGlyMetTrpSerIleGlyAlaGly : 20 ||||||||||||||||||||||||||||||!:!||||||||||||||||||||||||||| MetSerPheLeuGlnAspProSerPhePheSerMetGlyMetTrpSerIleGlyAlaGly 30496 : ATGTCTTTCCTCCAGGATCCCAGTTTCTTCAGCATGGGAATGTGGTCCATCGGTGCGGGG : 30439 21 : AlaLeuGlyAlaAlaAlaLeuAlaLeuLeuLeuAlaAsnThrAspValPheLeuSerLys : 40 |||:!!|||!.!|||||||||||||||||||||||||||||||||:!!||||||:!!||| AlaValGlyValAlaAlaLeuAlaLeuLeuLeuAlaAsnThrAspMetPheLeuAlaLys 30438 : GCCGTAGGAGTCGCTGCCTTGGCACTGCTCCTGGCCAACACAGACATGTTTCTGGCCAAG : 30379 41 : ProGlnLysAlaAlaLeuGluTyrLeuGluAspIleAspLeuLysThrLeuGluLys{G} : 60 ||||||||||||||||||||||||||||||||||||||||||||||||:!!||||||{!} ProGlnLysAlaAlaLeuGluTyrLeuGluAspIleAspLeuLysThrMetGluLys{A} 30378 : CCCCAGAAAGCAGCATTGGAGTATCTGGAGGACATAGACCTGAAAACAATGGAGAAG{G} : 30319 61 : >>>> Target Intron 1 >>>> {lu}ProArgThrPheLysAlaLysGluLeu : 69 2018 bp {!:} !! !|||||||||||||||! !||| ++ ++{sp}AlaValThrPheLysAlaLysAlaLeu 30318 : gt.........................ag{AT}GCAGTAACTTTTAAAGCGAAGGCGCTC : 28274 70 : TrpGluLysAsnGlyAlaValIleMetAlaValArgArgProGlyCysPheLeuCysArg : 89 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| TrpGluLysAsnGlyAlaValIleMetAlaValArgArgProGlyCysPheLeuCysArg 28273 : TGGGAGAAGAATGGAGCTGTGATTATGGCTGTGCGGAGGCCAGGCTGCTTCCTCTGTCGG : 28214 90 : Glu >>>> Target Intron 2 >>>> GluAlaAlaAspLeuSerSerLeuLysS : 100 ||| 692 bp ||||||||||||||||||||||||||| Glu++ ++GluAlaAlaAspLeuSerSerLeuLysP 28213 : GAGgt.........................agGAGGCTGCAGATTTGTCCTCCCTGAAGC : 27489 101 : erMetLeuAspGlnLeuGlyValProLeuTyrAlaValValLysGluHisIleArgThrG : 120 !!! !||||||:!!|||||||||||||||||||||||||||||||||||||||!:!!.!| roArgLeuAspGluLeuGlyValProLeuTyrAlaValValLysGluHisIleLysAsnG 27488 : CCAGGTTGGACGAGCTGGGTGTCCCCCTCTATGCAGTGGTAAAGGAGCACATCAAGAATG : 27429 121 : luValLysAspPheGlnProTyrPheLysGlyGluIlePheLeuAspGluLys >>>> : 138 ||||||||||||||||||||||||||||||||||||||||||||||||||||| luValLysAspPheGlnProTyrPheLysGlyGluIlePheLeuAspGluLys++ 27428 : AAGTTAAGGACTTCCAGCCTTATTTCAAAGGAGAAATCTTCCTGGATGAAAAGgt..... : 27373 139 : Target Intron 3 >>>> LysLysPheTyrGlyProGlnArgArgLysMetMetPh : 150 564 bp |||||||||||||||||||||..!|||||||||:!!|| ++LysLysPheTyrGlyProGlnGluArgLysMetValPh 27372 : ....................agAAAAAGTTCTATGGCCCACAAGAACGGAAGATGGTGTT : 26775 151 : eMetGlyPheIleArgLeuGlyValTrpTyrAsnPhePheArgAlaTrpAsnGlyGlyPh : 170 ||||||||||:!!||||||||||||! ! ! ||||||||||||||||||!!.|||||||| eMetGlyPheValArgLeuGlyValLeuGlnAsnPhePheArgAlaTrpLysGlyGlyPh 26774 : CATGGGATTTGTCCGTCTGGGCGTCTTGCAAAACTTCTTCCGGGCCTGGAAAGGAGGCTT : 26715 171 : eSerGlyAsnLeuGluGlyGluGlyPheIleLeuGlyGlyValPheValValGlySerGl : 190 |||||||||||||!!:|||||||||||||||||||||||||||||||||:!!||| !!|| eSerGlyAsnLeuAspGlyGluGlyPheIleLeuGlyGlyValPheValMetGlyProGl 26714 : TTCTGGAAACCTGGATGGCGAGGGCTTCATCCTTGGGGGAGTCTTTGTGATGGGACCAGG : 26655 191 : yLysGln >>>> Target Intron 4 >>>> GlyIleLeuLeuGluHisArgGlu : 200 ||||||| 3730 bp |||||||||||||||||||||||| yLysGln++ ++GlyIleLeuLeuGluHisArgGlu 26654 : AAAGCAGgt.........................agGGCATCCTTCTTGAGCACCGAGAA : 22895 201 : LysGluPheGlyAspLysValAsnLeuLeuSerValLeuGluAlaAlaLysMetIleLys : 220 ||||||||||||||||||||||||||| !||||||||||||||||||!:!! !|||!:! LysGluPheGlyAspLysValAsnLeuThrSerValLeuGluAlaAlaArgLysIleArg 22894 : AAAGAATTTGGAGACAAAGTAAACCTGACTTCTGTTCTGGAAGCTGCTAGGAAGATCAGA : 22835 221 : ProGlnThrLeuAlaSerGluLysLys : 229 |||||||||! !||||||||||||:!! ProGlnThrSerAlaSerGluLysGln 22834 : CCACAGACTTCGGCCTCAGAGAAACAA : 22806 vulgar: SPP00000026_1.0 0 229 . NIVO01064043.1:subseq(0,374980) 30496 22805 - 1036 M 59 177 S 0 1 5 0 2 I 0 2014 3 0 2 S 1 2 M 30 90 5 0 2 I 0 688 3 0 2 M 47 141 5 0 2 I 0 560 3 0 2 M 55 165 5 0 2 I 0 3726 3 0 2 M 37 111 # --- 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 # NIVO01064043.1:subseq(0,374980) exonerate:protein2genome:local gene 22806 30496 1036 - . gene_id 1 ; sequence SPP00000026_1.0 ; gene_orientation + NIVO01064043.1:subseq(0,374980) exonerate:protein2genome:local cds 30319 30496 . - . NIVO01064043.1:subseq(0,374980) exonerate:protein2genome:local exon 30319 30496 . - . insertions 0 ; deletions 0 NIVO01064043.1:subseq(0,374980) exonerate:protein2genome:local splice5 30317 30318 . - . intron_id 1 ; splice_site "GT" NIVO01064043.1:subseq(0,374980) exonerate:protein2genome:local intron 28301 30318 . - . intron_id 1 NIVO01064043.1:subseq(0,374980) exonerate:protein2genome:local splice3 28301 28302 . - . intron_id 0 ; splice_site "ag" NIVO01064043.1:subseq(0,374980) exonerate:protein2genome:local cds 28209 28300 . - . NIVO01064043.1:subseq(0,374980) exonerate:protein2genome:local exon 28209 28300 . - . insertions 0 ; deletions 0 NIVO01064043.1:subseq(0,374980) exonerate:protein2genome:local splice5 28207 28208 . - . intron_id 2 ; splice_site "GT" NIVO01064043.1:subseq(0,374980) exonerate:protein2genome:local intron 27517 28208 . - . intron_id 2 NIVO01064043.1:subseq(0,374980) exonerate:protein2genome:local splice3 27517 27518 . - . intron_id 1 ; splice_site "AG" NIVO01064043.1:subseq(0,374980) exonerate:protein2genome:local cds 27376 27516 . - . NIVO01064043.1:subseq(0,374980) exonerate:protein2genome:local exon 27376 27516 . - . insertions 0 ; deletions 0 NIVO01064043.1:subseq(0,374980) exonerate:protein2genome:local splice5 27374 27375 . - . intron_id 3 ; splice_site "GT" NIVO01064043.1:subseq(0,374980) exonerate:protein2genome:local intron 26812 27375 . - . intron_id 3 NIVO01064043.1:subseq(0,374980) exonerate:protein2genome:local splice3 26812 26813 . - . intron_id 2 ; splice_site "AG" NIVO01064043.1:subseq(0,374980) exonerate:protein2genome:local cds 26647 26811 . - . NIVO01064043.1:subseq(0,374980) exonerate:protein2genome:local exon 26647 26811 . - . insertions 0 ; deletions 0 NIVO01064043.1:subseq(0,374980) exonerate:protein2genome:local splice5 26645 26646 . - . intron_id 4 ; splice_site "GT" NIVO01064043.1:subseq(0,374980) exonerate:protein2genome:local intron 22917 26646 . - . intron_id 4 NIVO01064043.1:subseq(0,374980) exonerate:protein2genome:local splice3 22917 22918 . - . intron_id 3 ; splice_site "AG" NIVO01064043.1:subseq(0,374980) exonerate:protein2genome:local cds 22806 22916 . - . NIVO01064043.1:subseq(0,374980) exonerate:protein2genome:local exon 22806 22916 . - . insertions 0 ; deletions 0 NIVO01064043.1:subseq(0,374980) exonerate:protein2genome:local similarity 22806 30496 1036 - . alignment_id 1 ; Query SPP00000026_1.0 ; Align 30497 1 177 ; Align 28299 61 90 ; Align 27517 91 141 ; Align 26812 138 165 ; Align 22917 193 111 # --- END OF GFF DUMP --- # -- completed exonerate analysis