Command line: [exonerate -m p2g --showtargetgff -q SEPHS2/SEPHS2_x.fa -t SEPHS2/SUBSEQ/CM004428.1_4_subseq.txt --exhaustive] Hostname: [ws101066u] C4 Alignment: ------------ Query: SPP00000630_2.0 Protein Selenophosphate synthetase 2 (SEPHS2) Zebrafish Target: CM004428.1:subseq(14101919,25248) Ictalurus punctatus breed USDA103 chromosome 15, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 1673 Query range: 16 -> 441 Target range: 25248 -> 19593 17 : TyrAlaValPheTyrProAlaGlyTyrGlnAlaLeuAsnProGluGluHisGlyLeuAsp : 36 :!!:!!!.!|||||||||||||||||||||||| |||||| !.!||||||||| HisSerAlaPheTyrProAlaGlyTyrGlnAla------ProGluThrArgGlyLeuAsp 25248 : CACTCGGCCTTTTACCCGGCCGGCTACCAGGCT------CCAGAAACTCGTGGCCTGGAC : 25197 37 : ArgGlyPheArgLeuThrAlaPheSerAspMetLysGlyUnkGlyCysLysValProGln : 56 |||..!||||||||||||||||||||||||||||||||| |||||||||||||||||| ArgAsnPheArgLeuThrAlaPheSerAspMetLysGly***GlyCysLysValProGln 25196 : CGCAATTTTCGGCTCACGGCCTTTTCCGATATGAAAGGATGAGGGTGCAAAGTCCCGCAG : 25137 57 : GluThrLeuLeuLysLeuLeuGlnGlyLeuGluProAspArgProProGlyGluAspGly : 76 ||||||||||||||||||||||||||||||||||||!!: !! !|||! !:!!||| GluThrLeuLeuLysLeuLeuGlnGlyLeuGluProGluValGlnProGluGlnAsp--- 25136 : GAGACACTGCTCAAACTCCTGCAGGGACTGGAGCCCGAGGTTCAGCCGGAACAGGAC--- : 25077 77 : GlyLeuGlyThrGlyValGlyAspGluThrAlaAspPheGly<->LeuValSerAlaAla : 95 |||:!!|||||||||!.!|||||||||.!!:!!!!:|||||| |||!.! !!!.!:!! GlyMetGlyThrGlyAlaGlyAspGluAlaSerGluPheGlyGlnLeuAlaProValSer 25076 : GGAATGGGGACAGGAGCCGGGGACGAGGCCTCGGAGTTCGGTCAGCTAGCTCCGGTCTCA : 25020 96 : GlnGlyProArgLeu{G} >>>> Target Intron 1 >>>> {ly}IleGlyMet : 104 !!.||||||||||||{|} 1328 bp {||}||||||||| HisGlyProArgLeu{G}++ ++{ly}IleGlyMet 25019 : CATGGTCCCCGGTTA{G}gt.........................ag{GC}ATTGGGATG : 23665 105 : AspSerCysValIleProLeuArgHisGlyGlyLeuSerLeuValGlnThrThrAspPhe : 124 ||||||||||||||||||||||||||||||||||||||||||:!!||||||||||||||| AspSerCysValIleProLeuArgHisGlyGlyLeuSerLeuIleGlnThrThrAspPhe 23664 : GATTCTTGTGTCATTCCCTTGAGGCATGGAGGTTTGTCCCTGATCCAGACTACAGATTTC : 23605 125 : PheTyrProLeuValGluAspProTyrMetMet >>>> Target Intron 2 >>>> : 136 ||||||||||||||||||||||||||||||||| 489 bp PheTyrProLeuValGluAspProTyrMetMet++ 23604 : TTCTATCCGTTAGTTGAAGATCCTTATATGATGgt......................... : 23567 137 : GlyArgIleAlaCysAlaAsnValLeuSerAspLeuTyrAlaMetGlyIleThrGluC : 155 ||||||||||||||||||||||||||||||||||||||||||||||||||||||!!:| ++GlyArgIleAlaCysAlaAsnValLeuSerAspLeuTyrAlaMetGlyIleThrAspC 23566 : agGGTCGTATAGCTTGTGCAAATGTGTTGAGCGATCTGTACGCCATGGGAATTACTGATT : 23023 156 : ysAspAsnMetLeuMetLeuLeuSerValSerGlnLysMetAsnGluLys >>>> Tar : 172 |||||||||||||||||||||||||||||||||||::!|||!:!|||||| ysAspAsnMetLeuMetLeuLeuSerValSerGlnArgMetSerGluLys++ 23022 : GTGACAACATGCTGATGCTATTGAGTGTCAGTCAGCGGATGAGCGAGAAGgt........ : 22970 173 : get Intron 3 >>>> GluArgAspLeuValLeuProLeuMetMetLysGlyPheAr : 185 1461 bp ||||||!!:! !|||:!!||||||:!!|||||||||||||| ++GluArgGluGlnValMetProLeuLeuMetLysGlyPheAr 22969 : .................agGAGAGAGAGCAGGTGATGCCCTTATTGATGAAGGGTTTCCG : 21472 186 : gAspAlaAlaGluGluGlyGlyThrSerValThrGlyGlyGlnThrValIleAsnProTr : 205 |||||||||||||||||||||||||||||||||||||||||||||||||:!:|||||||| gAspAlaAlaGluGluGlyGlyThrSerValThrGlyGlyGlnThrValValAsnProTr 21471 : TGATGCAGCTGAAGAGGGCGGGACTTCCGTGACAGGCGGACAGACAGTTGTGAACCCTTG : 21412 206 : pIleIleIleGlyGlyValAlaSerValValCysGlnProAsnAspPheIle{Le} >> : 223 ||||||||||||||||||||||!!!|||||||||||||||!.!|||||||||{:!} pIleIleIleGlyGlyValAlaThrValValCysGlnProThrAspPheIle{Me}++ 21411 : GATCATCATAGGGGGTGTGGCCACAGTAGTGTGCCAGCCTACTGACTTCATT{AT}gt.. : 21354 224 : >> Target Intron 4 >>>> {u}ProAspGlyAlaValProGlyAspValLeuVa : 234 401 bp {!}|||||||||||||||||||||||||||||||| ++{t}ProAspGlyAlaValProGlyAspValLeuVa 21353 : .......................ag{G}CCAGATGGTGCAGTACCAGGTGACGTTCTCGT : 20924 235 : lLeuThrLysProLeuGlyThrGlnValAlaValAsnAlaHisGlnTrpLeuAspIle : 254 ||||||||||||||||||||||||||||||||||||||||||||||||||||:!!||| lLeuThrLysProLeuGlyThrGlnValAlaValAsnAlaHisGlnTrpLeuAsnIle++ 20923 : TCTGACAAAACCTCTTGGAACCCAAGTAGCCGTCAACGCTCATCAGTGGCTGAATATAgt : 20864 255 : >>>> Target Intron 5 >>>> ProGluLysTrpAsnLysIleLysLeuValIle : 264 110 bp ||||||||||||||||||||||||||||||||| ++ProGluLysTrpAsnLysIleLysLeuValIle 20863 : .........................agCCTGAAAAATGGAACAAGATCAAGTTGGTGATC : 20724 265 : SerArgGluGluValGluGlnAlaTyrGlnGluAlaMetLeuAsnMetAlaThrLeuAsn : 284 |||!:!||||||||||||!!.||||||||||||||||||||||||||||||||||||||| SerLysGluGluValGluHisAlaTyrGlnGluAlaMetLeuAsnMetAlaThrLeuAsn 20723 : TCAAAGGAGGAGGTGGAGCATGCATATCAAGAGGCCATGCTCAACATGGCTACCCTTAAC : 20664 285 : ArgThr{A} >>>> Target Intron 6 >>>> {la}AlaAlaLeuMetHisLys : 293 ||||||{|} 522 bp {||}|||||||||||||||||| ArgThr{A}++ ++{la}AlaAlaLeuMetHisLys 20663 : CGCACT{G}gt.........................ag{CT}GCTGCTCTGATGCACAAG : 20115 294 : PheAsnAlaHisAlaAlaThrAspIleThrGlyPheGlyIleIleGlyHisAlaArgAsn : 313 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| PheAsnAlaHisAlaAlaThrAspIleThrGlyPheGlyIleIleGlyHisAlaArgAsn 20114 : TTCAATGCACATGCTGCCACGGACATCACCGGCTTCGGGATCATTGGGCATGCACGGAAC : 20055 314 : LeuAlaLysGlnGlnLysAsnAspValAlaPheValIleHisAsnLeuProIleIleSer : 333 ||||||||||||||| |||!!:||||||||||||||||||||||||||||||:!!:!! LeuAlaLysGlnGlnCysAsnGluValAlaPheValIleHisAsnLeuProIleLeuAla 20054 : CTGGCCAAGCAGCAGTGCAACGAAGTGGCCTTTGTTATCCACAACCTGCCAATCCTTGCT : 19995 334 : LysMetAlaAlaIleSerLysAlaGlyGlyAsnLeuPheGlyLeuLeuGlnGlyThrSer : 353 ||||||||||||:!!|||||||||||||||||||||||||||||||||!!.||||||||| LysMetAlaAlaValSerLysAlaGlyGlyAsnLeuPheGlyLeuLeuHisGlyThrSer 19994 : AAGATGGCTGCTGTCAGCAAAGCCGGGGGGAACCTCTTCGGCCTTCTGCACGGCACCTCC : 19935 354 : SerGluThrSer{G} >>>> Target Intron 7 >>>> {ly}GlyLeuLeuIle : 362 ||||||||||||{|} 96 bp {||}|||||||||||| SerGluThrSer{G}++ ++{ly}GlyLeuLeuIle 19934 : TCCGAGACTTCA{G}gt.........................ag{GT}GGTTTGCTGATC : 19812 363 : CysLeuProArgGluGlnAlaAlaArgPheCysAlaGluMetLysSerSerArgMetGly : 382 |||||||||||||||||||||||||||||||||||||||||||||:!!|||||| CysLeuProArgGluGlnAlaAlaArgPheCysAlaGluMetLysAlaSerArg------ 19811 : TGTCTTCCGCGGGAGCAGGCAGCACGTTTCTGTGCCGAGATGAAAGCGAGCCGG------ : 19755 383 : LeuLeuGlyAlaGlyGlnAspGlyGlyValGlyAspGlyGlnGlnAlaTrpIleIleGly : 402 .!!..!||||||||||||||||||||||||||||||||| ------------------ThrSerSerValGlyAspGlyGlnGlnAlaTrpIleIleGly 19754 : ------------------ACGAGCTCTGTGGGCGATGGACAGCAGGCCTGGATAATTGGA : 19716 403 : IleValGluLysGlyAsnArgCysAlaArgIleIleAspLysProArgIleIleGluVal : 422 |||||||||||||||!:!|||||||||||||||||||||||||||||||||||||||||| IleValGluLysGlySerArgCysAlaArgIleIleAspLysProArgIleIleGluVal 19715 : ATCGTGGAGAAGGGTAGCCGCTGTGCACGCATCATCGACAAGCCCCGCATCATTGAGGTG : 19656 423 : ProTyrArgGly<->SerValValSerValGlnGluGlySerAsnAsnAsnAlaSerPro : 441 |||||||||||| :!!!.!! !!:!!.!|||||||||:!!|||!.. !.!!!!!||| ProTyrArgGlyLeuAlaAlaGlyAsnAlaGlnGluGlyAlaAsnThrValThrThrPro 19655 : CCCTATCGTGGTCTCGCAGCAGGAAACGCACAGGAGGGTGCAAACACGGTCACCACTCCT : 19596 vulgar: SPP00000630_2.0 16 441 . CM004428.1:subseq(14101919,25248) 25248 19593 - 1673 M 11 33 G 2 0 M 46 138 G 1 0 M 14 42 G 0 3 M 10 30 S 0 1 5 0 2 I 0 1324 3 0 2 S 1 2 M 34 102 5 0 2 I 0 485 3 0 2 M 36 108 5 0 2 I 0 1457 3 0 2 M 51 153 S 0 2 5 0 2 I 0 397 3 0 2 S 1 1 M 30 90 5 0 2 I 0 106 3 0 2 M 33 99 S 0 1 5 0 2 I 0 518 3 0 2 S 1 2 M 70 210 S 0 1 5 0 2 I 0 92 3 0 2 S 1 2 M 22 66 G 8 0 M 38 114 G 0 3 M 15 45 # --- 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 # CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local gene 19594 25248 1673 - . gene_id 1 ; sequence SPP00000630_2.0 ; gene_orientation + CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local cds 25002 25248 . - . CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local exon 25002 25248 . - . insertions 3 ; deletions 3 CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local splice5 25000 25001 . - . intron_id 1 ; splice_site "GT" CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local intron 23674 25001 . - . intron_id 1 CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local splice3 23674 23675 . - . intron_id 0 ; splice_site "AG" CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local cds 23570 23673 . - . CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local exon 23570 23673 . - . insertions 0 ; deletions 0 CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local splice5 23568 23569 . - . intron_id 2 ; splice_site "GT" CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local intron 23081 23569 . - . intron_id 2 CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local splice3 23081 23082 . - . intron_id 1 ; splice_site "AG" CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local cds 22973 23080 . - . CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local exon 22973 23080 . - . insertions 0 ; deletions 0 CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local splice5 22971 22972 . - . intron_id 3 ; splice_site "GT" CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local intron 21512 22972 . - . intron_id 3 CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local splice3 21512 21513 . - . intron_id 2 ; splice_site "AG" CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local cds 21357 21511 . - . CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local exon 21357 21511 . - . insertions 0 ; deletions 0 CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local splice5 21355 21356 . - . intron_id 4 ; splice_site "GT" CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local intron 20956 21356 . - . intron_id 4 CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local splice3 20956 20957 . - . intron_id 3 ; splice_site "AG" CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local cds 20865 20955 . - . CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local exon 20865 20955 . - . insertions 0 ; deletions 0 CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local splice5 20863 20864 . - . intron_id 5 ; splice_site "GT" CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local intron 20755 20864 . - . intron_id 5 CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local splice3 20755 20756 . - . intron_id 4 ; splice_site "AG" CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local cds 20655 20754 . - . CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local exon 20655 20754 . - . insertions 0 ; deletions 0 CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local splice5 20653 20654 . - . intron_id 6 ; splice_site "GT" CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local intron 20133 20654 . - . intron_id 6 CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local splice3 20133 20134 . - . intron_id 5 ; splice_site "AG" CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local cds 19920 20132 . - . CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local exon 19920 20132 . - . insertions 0 ; deletions 0 CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local splice5 19918 19919 . - . intron_id 7 ; splice_site "GT" CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local intron 19824 19919 . - . intron_id 7 CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local splice3 19824 19825 . - . intron_id 6 ; splice_site "AG" CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local cds 19594 19823 . - . CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local exon 19594 19823 . - . insertions 3 ; deletions 8 CM004428.1:subseq(14101919,25248) exonerate:protein2genome:local similarity 19594 25248 1673 - . alignment_id 1 ; Query SPP00000630_2.0 ; Align 25249 17 33 ; Align 25216 30 138 ; Align 25078 77 42 ; Align 25033 91 30 ; Align 23672 102 102 ; Align 23081 136 108 ; Align 21512 172 153 ; Align 20955 224 90 ; Align 20755 254 99 ; Align 20131 288 210 ; Align 19822 359 66 ; Align 19756 389 114 ; Align 19639 427 45 # --- END OF GFF DUMP --- # -- completed exonerate analysis