Command line: [exonerate -m p2g --showtargetgff -q Sel15.aa.fa.txt -t Sel15.subseq.txt] Hostname: [ws116044u] C4 Alignment: ------------ Query: SPP00000067_2.0 # Protein # 15 kDa selenoprotein (Sel15) # Human Target: JSZB01012017.1:subseq(0,120000) Manis javanica isolate MP_PG03-UM SCAFFOLD14246, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 719 Query range: 8 -> 162 Target range: 85564 -> 41462 9 : LeuValProAlaPheGlyLeuArgLeuLeuLeuAlaThrValLeuGlnAla >>>> Ta : 26 ||| !! !|||!!.||||||||||||||||||.!!||||||||||||||| LeuArgGlnAlaLeuGlyLeuArgLeuLeuLeuThrThrValLeuGlnAla++ 85564 : CTTAGGCAGGCGTTAGGGCTGCGGTTGCTATTGACCACTGTGCTTCAAGCGgt....... : 85511 27 : rget Intron 1 >>>> ValSerAlaPheGlyAlaGluPheSerSerGluAlaCysA : 39 7730 bp |||||||||!!.|||||||||||||||||||||||||||| ++ValSerAlaLeuGlyAlaGluPheSerSerGluAlaCysA 85510 : ..................aggtgTCTGCATTGGGGGCAGAGTTTTCATCAGAGGCTTGCA : 77744 40 : rgGluLeuGlyPheSerSerAsnLeuLeuCysSerSerCysAspLeuLeuGlyGlnPheA : 59 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| rgGluLeuGlyPheSerSerAsnLeuLeuCysSerSerCysAspLeuLeuGlyGlnPheA 77743 : GAGAATTAGGCTTCTCCAGCAACTTGCTGTGCAGCTCTTGTGATCTTCTTGGACAGTTCA : 77684 60 : snLeuLeuGlnLeuAspProAspCysArgGlyCysCysGlnGluGluAlaGlnPheGluT : 79 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| snLeuLeuGlnLeuAspProAspCysArgGlyCysCysGlnGluGluAlaGlnPheGluT 77683 : ATCTGCTTCAGCTGGATCCCGACTGTAGAGGATGCTGTCAGGAAGAAGCACAATTTGAAA : 77624 80 : hrLysLys >>>> Target Intron 2 >>>> LeuTyrAlaGlyAlaIleLeuGl : 89 |||||||| 21935 bp ||||||||||||||||||||||| hrLysLys++ ++LeuTyrAlaGlyAlaIleLeuGl 77623 : CCAAAAAGgt.........................agTTGTATGCAGGTGCTATTCTCGA : 55659 90 : uValCysGlyUnkLysLeuGlyArgPheProGlnValGln{A} >>>> Target Int : 103 |||||||||| |||||||||||||||||||||||||||{|} 11220 uValCysGly***LysLeuGlyArgPheProGlnValGln{A}++ 55658 : AGTCTGTGGATGAAAATTGGGAAGGTTCCCTCAAGTCCAA{G}gt............... : 55614 104 : ron 3 >>>> {la}PheValArgSerAspLysProLysLeuPheArgGlyLeuGlnIl : 118 bp {||}|||||||||||||||||||||||||||||||||||||||||||| ++{la}PheValArgSerAspLysProLysLeuPheArgGlyLeuGlnIl 55613 : ..........ag{cT}TTCGTCAGGAGTGATAAACCCAAACTCTTCAGAGGACTACAGAT : 44352 119 : eLys >>>> Target Intron 4 >>>> TyrValArgGlySerAspProValLeu : 128 |||| 2755 bp ||||||||||||||||||||||||||| eLys++ ++TyrValArgGlySerAspProValLeu 44351 : CAAGgt.........................agTATGTTCGTGGTTCAGACCCTGTACTC : 41567 129 : LysLeuLeuAspAspAsnGlyAsnIleAlaGluGluLeuSerIleLeuLysTrpAsnThr : 148 |||||||||||||||!:!|||||||||||||||||||||||||||||||||||||||||| LysLeuLeuAspAspSerGlyAsnIleAlaGluGluLeuSerIleLeuLysTrpAsnThr 41566 : AAGCTTCTGGATGACAGTGGGAACATTGCTGAGGAACTAAGCATCCTCAAATGGAACACA : 41507 149 : AspSerValGluGluPheLeuSerGluLysLeuGluArgIle : 162 |||||||||||||||||||||||||||||||||||||||||| AspSerValGluGluPheLeuSerGluLysLeuGluArgIle 41506 : GACAGCGTGGAAGAATTCCTCAGTGAAAAATTGGAACGCATA : 41463 vulgar: SPP00000067_2.0 8 162 . JSZB01012017.1:subseq(0,120000) 85564 41462 - 719 M 17 51 5 0 2 I 0 7726 3 0 2 M 56 168 5 0 2 I 0 21931 3 0 2 M 21 63 S 0 1 5 0 2 I 0 11216 3 0 2 S 1 2 M 16 48 5 0 2 I 0 2751 3 0 2 M 43 129 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2016-11-22 ##type DNA # # # seqname source feature start end score strand frame attributes # JSZB01012017.1:subseq(0,120000) exonerate:protein2genome:local gene 41463 85564 719 - . gene_id 1 ; sequence SPP00000067_2.0 ; gene_orientation + JSZB01012017.1:subseq(0,120000) exonerate:protein2genome:local cds 85514 85564 . - . JSZB01012017.1:subseq(0,120000) exonerate:protein2genome:local exon 85514 85564 . - . insertions 0 ; deletions 0 JSZB01012017.1:subseq(0,120000) exonerate:protein2genome:local splice5 85512 85513 . - . intron_id 1 ; splice_site "GT" JSZB01012017.1:subseq(0,120000) exonerate:protein2genome:local intron 77784 85513 . - . intron_id 1 JSZB01012017.1:subseq(0,120000) exonerate:protein2genome:local splice3 77784 77785 . - . intron_id 0 ; splice_site "ag" JSZB01012017.1:subseq(0,120000) exonerate:protein2genome:local cds 77616 77783 . - . JSZB01012017.1:subseq(0,120000) exonerate:protein2genome:local exon 77616 77783 . - . insertions 0 ; deletions 0 JSZB01012017.1:subseq(0,120000) exonerate:protein2genome:local splice5 77614 77615 . - . intron_id 2 ; splice_site "GT" JSZB01012017.1:subseq(0,120000) exonerate:protein2genome:local intron 55681 77615 . - . intron_id 2 JSZB01012017.1:subseq(0,120000) exonerate:protein2genome:local splice3 55681 55682 . - . intron_id 1 ; splice_site "AG" JSZB01012017.1:subseq(0,120000) exonerate:protein2genome:local cds 55617 55680 . - . JSZB01012017.1:subseq(0,120000) exonerate:protein2genome:local exon 55617 55680 . - . insertions 0 ; deletions 0 JSZB01012017.1:subseq(0,120000) exonerate:protein2genome:local splice5 55615 55616 . - . intron_id 3 ; splice_site "GT" JSZB01012017.1:subseq(0,120000) exonerate:protein2genome:local intron 44397 55616 . - . intron_id 3 JSZB01012017.1:subseq(0,120000) exonerate:protein2genome:local splice3 44397 44398 . - . intron_id 2 ; splice_site "ag" JSZB01012017.1:subseq(0,120000) exonerate:protein2genome:local cds 44347 44396 . - . JSZB01012017.1:subseq(0,120000) exonerate:protein2genome:local exon 44347 44396 . - . insertions 0 ; deletions 0 JSZB01012017.1:subseq(0,120000) exonerate:protein2genome:local splice5 44345 44346 . - . intron_id 4 ; splice_site "GT" JSZB01012017.1:subseq(0,120000) exonerate:protein2genome:local intron 41592 44346 . - . intron_id 4 JSZB01012017.1:subseq(0,120000) exonerate:protein2genome:local splice3 41592 41593 . - . intron_id 3 ; splice_site "AG" JSZB01012017.1:subseq(0,120000) exonerate:protein2genome:local cds 41463 41591 . - . JSZB01012017.1:subseq(0,120000) exonerate:protein2genome:local exon 41463 41591 . - . insertions 0 ; deletions 0 JSZB01012017.1:subseq(0,120000) exonerate:protein2genome:local similarity 41463 85564 719 - . alignment_id 1 ; Query SPP00000067_2.0 ; Align 85565 9 51 ; Align 77784 26 168 ; Align 55681 82 63 ; Align 44395 104 48 ; Align 41592 120 129 # --- END OF GFF DUMP --- # -- completed exonerate analysis