Command line: [exonerate -m p2g --showtargetgff -q selenou2_zf.aa.fa.canviat -t genomeselenou2_JXSJ01000066.1.fa] Hostname: [ws101073u] C4 Alignment: ------------ Query: SPP00000655_2.0 # Protein # Selenoprotein U2 (SELENOU2) # Zebrafish Target: JXSJ01000066.1:subseq(393031,101771) Miichthys miiuy scaffold66, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 683 Query range: 4 -> 222 Target range: 51771 -> 49999 5 : ProIleThrArgGlnIleSerSerSerValGlnAspAsnLysGlnAlaSerHisAsnIle : 24 |||||||||||||||:!!!!!!.!..! !!:!!!:! ! ! :!!||| !:!!||| ProIleThrArgGlnValThrLysGluAsnArgGluIleGlyThrSerSerValAspIle 51771 : CCGATAACTCGGCAAGTCACTAAAGAAAATCGAGAAATCGGAACCTCTTCTGTCGACATT : 51714 25 : CysLeuSerGluLeuLysAsnCysPheIlePheAspArgHisGlyAspLysLysSerPhe : 44 !|||!.!!!::!!:!!:!!|||.!.|||!:!||||||||||||! !...! !!||| ProLeuLysAspValGluAspCysIleIleTyrAspArgHisGlyValSerIleProPhe 51713 : CCTCTTAAAGACGTCGAGGATTGTATAATATACGACCGGCATGGAGTCTCTATTCCGTTC : 51654 45 : SerSerLeuPheGluHisAsnLysAlaIleValIlePheValArg >>>> Target I : 60 !.!!.!||| !:!! !!!..|||:!!:!::!!|||||||||||| 337 LysLysLeuHisGlnAspArgLysSerValIleIlePheValArg++ 51653 : AAGAAATTACACCAAGACAGGAAATCGGTGATAATTTTCGTGCGGgt............. : 51604 61 : ntron 1 >>>> HisPheLeuCysTyrThrCysLysGluTyrValGluAspLeuGlyL : 75 bp :!!||||||||||||!:!||||||||||||||||||||||||.!!| ++AsnPheLeuCysTyrSerCysLysGluTyrValGluAspLeuSerL 51603 : ............agAATTTCCTGTGCTATAGCTGTAAAGAATATGTGGAAGACTTGAGCA : 51224 76 : ysIleProGlnHisValLeuGln >>>> Target Intron 2 >>>> AspSerAs : 85 ||||||||:!!.!.!.!|||:!! 474 bp |||:!!:! ysIleProLysGluAlaLeuGlu++ ++AspAlaAs 51223 : AAATACCTAAAGAAGCCCTGGAGgt.........................agGATGCTGA : 50720 86 : nValArgLeuValValIleGlyGlnSerSerTyrSerHisIleGln >>>> Target : 101 !:!!||||||||||||||||||||||||:!!:!! !||||||:!: 93 pIleArgLeuValValIleGlyGlnSerAlaHisHisHisIleGlu++ 50719 : CATTAGACTGGTTGTGATTGGTCAGTCTGCTCACCATCATATAGAggt............ : 50670 102 : Intron 3 >>>> GlyPheCysSerLeuThrGlyPheProHisGluIleTyrValAsp : 115 bp !||||||||||||||||||!:!|||:!!|||!!:||||||||| ++ProPheCysSerLeuThrGlyTyrProTyrGluMetTyrValAsp 50669 : .............agCCATTCTGCTCACTGACGGGCTATCCCTATGAAATGTATGTGGAC : 50537 116 : ProGluArgGlnIleTyrLysArgLeuGlyLeuArgArgGlyGluThrTyrMetGluThr : 135 ||||||||| ||||||:!!!:!||||||::::::||| ||||||::: !!::!! ProGluArgCysIleTyrGlnLysLeuGlyMetLysArgGluGluThrPheThrAspSer 50536 : CCGGAGAGGTGCATTTATCAAAAGCTtgggatgaagagagaggaaacatttaCAGACTCA : 50477 136 : {P} >>>> Target Intron 4 >>>> {ro}SerValSerProHisValLysSer : 144 { } 95 bp {!!} ! !|||||||||||||||||| {A}++ ++{la}HisProSerProHisValLysSer 50476 : {G}gt.........................ag{CC}CACCCTAGTCCCCATGTGAAGTCT : 50355 145 : SerMetLeuSerGlySerLeuLysSerValTrpArgAlaMetThrSerProValPheAsp : 164 .!!|||!!!! !|||..!:!!||||||:!!|||||||||||||||||||||||||||||| GlyMetPheMetGlyGlnMetLysSerIleTrpArgAlaMetThrSerProValPheAsp 50354 : GGTATGTTTATGGGACAAATGAAGAGTATATGGAGGGCCATGACTAGCCCTGTATTTGAC : 50295 165 : PheGlnGlyAspProGlnGlnGlnGlyGlyAlaLeuIleValGlyPro{G} >>>> Ta : 181 ||||||||||||! !!!.|||||||||||||||:!!|||!.!||||||{|} PheGlnGlyAspLeuHisGlnGlnGlyGlyAlaIleIleAlaGlyPro{G}++ 50294 : TTCCAGGGTGACCTACATCAGCAAGGTGGAGCCATTATTGCAGGACCA{G}gt....... : 50241 182 : rget Intron 5 >>>> {ly}ProAspValHisPheAlaHisPheAspMetAsnArg : 193 116 bp {||} !!.!.|||||||||:!!|||||||||! |||||| ++{ly}SerGlnValHisPheSerHisPheAspThrAsnArg 50240 : ..................ag{GC}TCTCAAGTTCATTTTTCCCATTTTGACACAAACCGC : 50092 194 : LeuAspHisMetProIleAsnTrpLeuLeuGlnLeuAlaGlyLeu<->ProAlaLeuAsp : 212 ||||||||||||||||||||||||||||||||||||||||||:!! ! !.!!|||||| LeuAspHisMetProIleAsnTrpLeuLeuGlnLeuAlaGlyValGlnGlnThrLeuAsp 50091 : CTGGACCACATGCCCATTAACTGGCTCCTGCAGCTCGCCGGAGTTCAACAGACTCTGGAC : 50032 213 : PheSerAspSerProLysIleIleAsnIle : 222 |||||||||!.!||||||||||||:!!:!: PheSerAspLysProLysIleIleHisVal 50031 : TTCAGTGATAAGCCAAAGATCATCCATGTG : 50000 vulgar: SPP00000655_2.0 4 222 . JXSJ01000066.1:subseq(393031,101771) 51771 49999 - 683 M 55 165 5 0 2 I 0 333 3 0 2 M 23 69 5 0 2 I 0 470 3 0 2 M 18 54 5 0 2 I 0 89 3 0 2 M 35 105 S 0 1 5 0 2 I 0 91 3 0 2 S 1 2 M 44 132 S 0 1 5 0 2 I 0 112 3 0 2 S 1 2 M 27 81 G 0 3 M 14 42 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2016-11-21 ##type DNA # # # seqname source feature start end score strand frame attributes # JXSJ01000066.1:subseq(393031,101771) exonerate:protein2genome:local gene 50000 51771 683 - . gene_id 1 ; sequence SPP00000655_2.0 ; gene_orientation + JXSJ01000066.1:subseq(393031,101771) exonerate:protein2genome:local cds 51607 51771 . - . JXSJ01000066.1:subseq(393031,101771) exonerate:protein2genome:local exon 51607 51771 . - . insertions 0 ; deletions 0 JXSJ01000066.1:subseq(393031,101771) exonerate:protein2genome:local splice5 51605 51606 . - . intron_id 1 ; splice_site "GT" JXSJ01000066.1:subseq(393031,101771) exonerate:protein2genome:local intron 51270 51606 . - . intron_id 1 JXSJ01000066.1:subseq(393031,101771) exonerate:protein2genome:local splice3 51270 51271 . - . intron_id 0 ; splice_site "AG" JXSJ01000066.1:subseq(393031,101771) exonerate:protein2genome:local cds 51201 51269 . - . JXSJ01000066.1:subseq(393031,101771) exonerate:protein2genome:local exon 51201 51269 . - . insertions 0 ; deletions 0 JXSJ01000066.1:subseq(393031,101771) exonerate:protein2genome:local splice5 51199 51200 . - . intron_id 2 ; splice_site "GT" JXSJ01000066.1:subseq(393031,101771) exonerate:protein2genome:local intron 50727 51200 . - . intron_id 2 JXSJ01000066.1:subseq(393031,101771) exonerate:protein2genome:local splice3 50727 50728 . - . intron_id 1 ; splice_site "AG" JXSJ01000066.1:subseq(393031,101771) exonerate:protein2genome:local cds 50673 50726 . - . JXSJ01000066.1:subseq(393031,101771) exonerate:protein2genome:local exon 50673 50726 . - . insertions 0 ; deletions 0 JXSJ01000066.1:subseq(393031,101771) exonerate:protein2genome:local splice5 50671 50672 . - . intron_id 3 ; splice_site "gt" JXSJ01000066.1:subseq(393031,101771) exonerate:protein2genome:local intron 50580 50672 . - . intron_id 3 JXSJ01000066.1:subseq(393031,101771) exonerate:protein2genome:local splice3 50580 50581 . - . intron_id 2 ; splice_site "AG" JXSJ01000066.1:subseq(393031,101771) exonerate:protein2genome:local cds 50474 50579 . - . JXSJ01000066.1:subseq(393031,101771) exonerate:protein2genome:local exon 50474 50579 . - . insertions 0 ; deletions 0 JXSJ01000066.1:subseq(393031,101771) exonerate:protein2genome:local splice5 50472 50473 . - . intron_id 4 ; splice_site "GT" JXSJ01000066.1:subseq(393031,101771) exonerate:protein2genome:local intron 50379 50473 . - . intron_id 4 JXSJ01000066.1:subseq(393031,101771) exonerate:protein2genome:local splice3 50379 50380 . - . intron_id 3 ; splice_site "aG" JXSJ01000066.1:subseq(393031,101771) exonerate:protein2genome:local cds 50244 50378 . - . JXSJ01000066.1:subseq(393031,101771) exonerate:protein2genome:local exon 50244 50378 . - . insertions 0 ; deletions 0 JXSJ01000066.1:subseq(393031,101771) exonerate:protein2genome:local splice5 50242 50243 . - . intron_id 5 ; splice_site "GT" JXSJ01000066.1:subseq(393031,101771) exonerate:protein2genome:local intron 50128 50243 . - . intron_id 5 JXSJ01000066.1:subseq(393031,101771) exonerate:protein2genome:local splice3 50128 50129 . - . intron_id 4 ; splice_site "AG" JXSJ01000066.1:subseq(393031,101771) exonerate:protein2genome:local cds 50000 50127 . - . JXSJ01000066.1:subseq(393031,101771) exonerate:protein2genome:local exon 50000 50127 . - . insertions 3 ; deletions 0 JXSJ01000066.1:subseq(393031,101771) exonerate:protein2genome:local similarity 50000 51771 683 - . alignment_id 1 ; Query SPP00000655_2.0 ; Align 51772 5 165 ; Align 51270 60 69 ; Align 50727 83 54 ; Align 50580 101 105 ; Align 50377 137 132 ; Align 50126 182 81 ; Align 50042 209 42 # --- END OF GFF DUMP --- # -- completed exonerate analysis