Command line: [/mnt/NFS_UPF/soft/exonerate/exonerate-2.2.0-x86_64/bin/exonerate -m p2g --showtargetgff --exhaustive yes -q query15 -t q_subseq] Hostname: [am30751-101042u] C4 Alignment: ------------ Query: SPP00000109_1.0 # Protein # Selenoprotein 15 (Sel15) # Mus musculus # Complete Target: PVKD010011498.1:subseq(1,48457) Cricetomys gambianus isolate US082 CriGam_scaffold_22982, whole genome shotgun sequence Model: protein2genome:local Raw score: 770 Query range: 0 -> 162 Target range: 10110 -> 41086 1 : MetAlaAlaGlyGlnGlyGlyTrpLeuArgProAlaLeuGlyLeuArgLeuLeuLeuAla : 20 |||||||||!.!|||! !|||||||||||||||||||||||||||||||||||||||||| MetAlaAlaAlaGlnAspGlyTrpLeuArgProAlaLeuGlyLeuArgLeuLeuLeuAla 10111 : ATGGCGGCCGCGCAGGACGGGTGGCTGCGGCCGGCTCTGGGGCTGCGCTTGCTGCTGGCG : 10168 21 : ThrAlaPheGlnAla >>>> Target Intron 1 >>>> AlaSerAlaLeuGlyA : 31 ||||||||||||||| 10237 bp !.!||||||||||||| ThrAlaPheGlnAla++ ++ValSerAlaLeuGlyA 10169 : ACTGCGTTTCAAGCGgt.........................agGTGTCTGCTCTTGGTG : 20438 32 : laGluPheAlaSerGluAlaCysArgGluLeuGlyPheSerSerAsnLeuLeuCysSerS : 51 ||||||||:!!||||||||||||||||||||||||||||||||||||||||||||||||| laGluPheSerSerGluAlaCysArgGluLeuGlyPheSerSerAsnLeuLeuCysSerS 20439 : CAGAGTTTTCATCAGAGGCATGTAGAGAGTTGGGCTTCTCCAGCAACTTGCTCTGCAGTT : 20498 52 : erCysAspLeuLeuGlyGlnPheAsnLeuLeuProLeuAspProValCysArgGlyCysC : 71 ||||||||||||||||||||||||||||||||! !||||||||||||||||||||||||| erCysAspLeuLeuGlyGlnPheAsnLeuLeuGlnLeuAspProValCysArgGlyCysC 20499 : CTTGTGATCTCCTAGGACAGTTTAACCTGCTGCAGCTGGATCCTGTTTGCAGAGGATGCT : 20558 72 : ysGlnGluGluAlaGlnPheGluThrLysLys >>>> Target Intron 2 >>>> : 82 |||||||||||||||||||||||||||||||| 15533 bp ysGlnGluGluAlaGlnPheGluThrLysLys++ + 20559 : GTCAGGAAGAAGCACAGTTTGAAACCAAGAAGgt.........................a : 36122 83 : LeuTyrAlaGlyAlaIleLeuGluValCysGlyUnkLysLeuGlyArgPheProGlnVa : 101 ||||||||||||||||||||||||||||||||| ||||||||||||||||||||||| +LeuTyrAlaGlyAlaIleLeuGluValCysGly***LysLeuGlyArgPheProGlnVa 36123 : gCTATATGCAGGAGCCATCCTTGAAGTCTGTGGATGAAAATTGGGGAGGTTCCCTCAAGT : 36181 102 : lGln{A} >>>> Target Intron 3 >>>> {la}PheValArgSerAspLysPr : 110 ||||{|} 2997 bp {||}|||||||||||||||||||| lGln{A}++ ++{la}PheValArgSerAspLysPr 36182 : CCAA{G}gt.........................ag{ct}tttgTCAGAAGTGATAAACC : 39205 111 : oLysLeuPheArgGlyLeuGlnIleLys >>>> Target Intron 4 >>>> Tyr : 120 |||||||||||||||||||||||||||| 1723 bp ||| oLysLeuPheArgGlyLeuGlnIleLys++ ++Tyr 39206 : CAAGCTGTTCAGAGGTCTGCAGATTAAGgt.........................agTAC : 40958 121 : ValArgGlySerAspProValLeuLysLeuLeuAspAspAsnGlyAsnIleAlaGluGlu : 140 ||||||||||||||||||:!!||||||||||||||||||||||||||||||||||||||| ValArgGlySerAspProIleLeuLysLeuLeuAspAspAsnGlyAsnIleAlaGluGlu 40959 : GTTCGTGGTTCAGATCCTATACTAAAGCTTTTGGACGACAACGGGAACATTGCTGAAGAA : 41018 141 : LeuSerIleLeuLysTrpAsnThrAspSerValGluGluPheLeuSerGluLysLeuGlu : 160 |||!:!|||||||||||||||||||||||||||||||||||||||||||||||||||||| LeuAsnIleLeuLysTrpAsnThrAspSerValGluGluPheLeuSerGluLysLeuGlu 41019 : CTGAACATCCTCAAATGGAACACAGATAGTGTGGAAGAGTTTCTGAGCGAGAAGTTGGAA : 41078 161 : ArgIle : 162 |||||| ArgIle 41079 : CGCATA : 41086 vulgar: SPP00000109_1.0 0 162 . PVKD010011498.1:subseq(1,48457) 10110 41086 + 770 M 25 75 5 0 2 I 0 10233 3 0 2 M 56 168 5 0 2 I 0 15529 3 0 2 M 21 63 S 0 1 5 0 2 I 0 2993 3 0 2 S 1 2 M 16 48 5 0 2 I 0 1719 3 0 2 M 43 129 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2019-11-20 ##type DNA # # # seqname source feature start end score strand frame attributes # PVKD010011498.1:subseq(1,48457) exonerate:protein2genome:local gene 10111 41086 770 + . gene_id 1 ; sequence SPP00000109_1.0 ; gene_orientation + PVKD010011498.1:subseq(1,48457) exonerate:protein2genome:local cds 10111 10185 . + . PVKD010011498.1:subseq(1,48457) exonerate:protein2genome:local exon 10111 10185 . + . insertions 0 ; deletions 0 PVKD010011498.1:subseq(1,48457) exonerate:protein2genome:local splice5 10186 10187 . + . intron_id 1 ; splice_site "GT" PVKD010011498.1:subseq(1,48457) exonerate:protein2genome:local intron 10186 20422 . + . intron_id 1 PVKD010011498.1:subseq(1,48457) exonerate:protein2genome:local splice3 20421 20422 . + . intron_id 0 ; splice_site "AG" PVKD010011498.1:subseq(1,48457) exonerate:protein2genome:local cds 20423 20590 . + . PVKD010011498.1:subseq(1,48457) exonerate:protein2genome:local exon 20423 20590 . + . insertions 0 ; deletions 0 PVKD010011498.1:subseq(1,48457) exonerate:protein2genome:local splice5 20591 20592 . + . intron_id 2 ; splice_site "GT" PVKD010011498.1:subseq(1,48457) exonerate:protein2genome:local intron 20591 36123 . + . intron_id 2 PVKD010011498.1:subseq(1,48457) exonerate:protein2genome:local splice3 36122 36123 . + . intron_id 1 ; splice_site "AG" PVKD010011498.1:subseq(1,48457) exonerate:protein2genome:local cds 36124 36187 . + . PVKD010011498.1:subseq(1,48457) exonerate:protein2genome:local exon 36124 36187 . + . insertions 0 ; deletions 0 PVKD010011498.1:subseq(1,48457) exonerate:protein2genome:local splice5 36188 36189 . + . intron_id 3 ; splice_site "GT" PVKD010011498.1:subseq(1,48457) exonerate:protein2genome:local intron 36188 39184 . + . intron_id 3 PVKD010011498.1:subseq(1,48457) exonerate:protein2genome:local splice3 39183 39184 . + . intron_id 2 ; splice_site "ag" PVKD010011498.1:subseq(1,48457) exonerate:protein2genome:local cds 39185 39234 . + . PVKD010011498.1:subseq(1,48457) exonerate:protein2genome:local exon 39185 39234 . + . insertions 0 ; deletions 0 PVKD010011498.1:subseq(1,48457) exonerate:protein2genome:local splice5 39235 39236 . + . intron_id 4 ; splice_site "GT" PVKD010011498.1:subseq(1,48457) exonerate:protein2genome:local intron 39235 40957 . + . intron_id 4 PVKD010011498.1:subseq(1,48457) exonerate:protein2genome:local splice3 40956 40957 . + . intron_id 3 ; splice_site "AG" PVKD010011498.1:subseq(1,48457) exonerate:protein2genome:local cds 40958 41086 . + . PVKD010011498.1:subseq(1,48457) exonerate:protein2genome:local exon 40958 41086 . + . insertions 0 ; deletions 0 PVKD010011498.1:subseq(1,48457) exonerate:protein2genome:local similarity 10111 41086 770 + . alignment_id 1 ; Query SPP00000109_1.0 ; Align 10111 1 75 ; Align 20423 26 168 ; Align 36124 82 63 ; Align 39187 104 48 ; Align 40958 120 129 # --- END OF GFF DUMP --- # -- completed exonerate analysis