Command line: [exonerate -m p2g --showtargetgff -q SelU1a/SelU1a_x.fa -t SelU1a/SUBSEQ/CM004426.1_5_subseq.txt --exhaustive] Hostname: [ws101122u] C4 Alignment: ------------ Query: SPP00000652_2.0 Protein Selenoprotein X1a (SELENOX1a) Zebrafish Target: CM004426.1:subseq(2597757,25152) Ictalurus punctatus breed USDA103 chromosome 13, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 559 Query range: 56 -> 222 Target range: 25110 -> 22547 57 : LeuLysThrIleAspGlyAspGluArgSerLeuLysAlaLysAlaLeuTrpGluLysSer : 76 :::!!.||| !|||!.!!!:|||!:!|||!!!|||||| .!!||||||||||||!!! ValAsnThrGlyAspAlaGluGluLysSerPheLysAlaPheThrLeuTrpGluLysThr 25110 : gtgAACACTGGAGACGCAGAGGAGAAGAGCTTTAAAGCGTTCACACTCTGGGAGAAAACC : 25053 77 : GlyAlaValIleMetAlaValArgArgProGlyUnkPheLeuCysArgGlu >>>> Ta : 94 ||||||||||||||||||||||||||||||||| ! !|||||||||||| GlyAlaValIleMetAlaValArgArgProGly***SerLeuCysArgGlu++ 25052 : GGAGCCGTGATCATGGCCGTCAGAAGACCTGGATGATCTTTGTGCAGagaggt....... : 24997 95 : rget Intron 1 >>>> GluAlaSerGluLeuSerSerLeuLysProGlnLeuAspG : 107 1292 bp |||||||||||||||||||||||||||||||||||||||| ++GluAlaSerGluLeuSerSerLeuLysProGlnLeuAspG 24996 : ..................agGAAGCCTCTGAGCTGTCCTCTCTGAAGCCTCAGCTGGATG : 23668 108 : luLeuGlyValProLeuTyrAlaValValLysGluAsnValGlyThrGluIleGlnAspP : 127 ||:!!|||||||||||||||||||||||||||||||||:!!|||||||||||||||:!!| luValGlyValProLeuTyrAlaValValLysGluAsnIleGlyThrGluIleGlnAsnP 23667 : AGGTCGGTGTTCCTCTGTACGCGGTGGTGAAAGAGAACATCGGTACAGAGATACAAAACT : 23608 128 : heArgProHisPheAlaGlyGluIlePheLeuAspGluLys >>>> Target Intro : 141 ||||||||:!!|||::!|||||||||||||||||||||||| 687 bp heArgProTyrPheSerGlyGluIlePheLeuAspGluLys++ 23607 : TCAGACCCTACTTCAGCGGAGAAATCTTCCTGGATGAGAAGgt................. : 23564 142 : n 2 >>>> GlnAlaPheTyrGlyProGlnGlnArgLysMetGlyGlyLeuGlyPheIl : 157 ||| !|||||| !|||::!:!!|||! !|||! ! !|||||||||! ++GlnLeuPheTyrThrProArgLysArgMetMetValLysLeuGlyPheTh 23563 : ........agCAATTATTTTACACTCCCAGGAAGAGGATGATGGTGAAACTGGGCTTCAC : 22831 158 : eArgLeuGlyValTrpGlnAsnPheValArgAlaTrpArgAlaGlyTyrGlnGlyAsnMe : 177 !|||:!!.!!||||||::!|||.!!:!!|||||| !! !||||||||||||||||| rArgIleSerValTrpArgAsnIleLeuArgAlaAsnTrpLysGlyTyrGlnGlyAsnMe 22830 : TCGTATCAGCGTGTGGAGGAACATCCTGCGCGCAAATTGGAAAGGTTACCAGGGCAACAT : 22771 178 : tAsnGlyGluGlyPheIleLeuGlyGlyValPheValMetGlySerGlyGlyGln >>> : 196 |!!.||||||||||||||||||||||||:!!!:!|||!!:||| !!! !..!||| tLysGlyGluGlyPheIleLeuGlyGlyLeuTyrValIleGlyProGluAsnGln++ 22770 : GAAAGGGGAGGGGTTTATACTTGGAGGCCTTTATGTGATTGGTCCAGAGAATCAGgt... : 22712 197 : > Target Intron 3 >>>> GlyValLeuLeuGluHisArgGluLysGluPheGly : 207 86 bp |||:!!|||:!!||||||!:!||||||||||||||| ++GlyIleLeuMetGluHisGlnGluLysGluPheGly 22711 : ......................agGGTATTTTAATGGAGCACCAGGAGAAGGAGTTTGGC : 22595 208 : AspLysValSerLeuGluSerValLeuGluAlaAlaLysLysVal : 222 ||||||!.!..!|||.!.|||||||||.!.|||!.!! !|||||| AspLysAlaAspLeuAsnSerValLeuHisAlaValMetLysVal 22594 : GATAAAGCTGATCTCAATTCAGTGTTACACGCAGTAATGAAAGTG : 22548 vulgar: SPP00000652_2.0 56 222 . CM004426.1:subseq(2597757,25152) 25110 22547 - 559 M 37 111 5 0 2 I 0 1288 3 0 2 M 47 141 5 0 2 I 0 683 3 0 2 M 55 165 5 0 2 I 0 82 3 0 2 M 27 81 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2016-11-15 ##type DNA # # # seqname source feature start end score strand frame attributes # CM004426.1:subseq(2597757,25152) exonerate:protein2genome:local gene 22548 25110 559 - . gene_id 1 ; sequence SPP00000652_2.0 ; gene_orientation + CM004426.1:subseq(2597757,25152) exonerate:protein2genome:local cds 25000 25110 . - . CM004426.1:subseq(2597757,25152) exonerate:protein2genome:local exon 25000 25110 . - . insertions 0 ; deletions 0 CM004426.1:subseq(2597757,25152) exonerate:protein2genome:local splice5 24998 24999 . - . intron_id 1 ; splice_site "gt" CM004426.1:subseq(2597757,25152) exonerate:protein2genome:local intron 23708 24999 . - . intron_id 1 CM004426.1:subseq(2597757,25152) exonerate:protein2genome:local splice3 23708 23709 . - . intron_id 0 ; splice_site "ag" CM004426.1:subseq(2597757,25152) exonerate:protein2genome:local cds 23567 23707 . - . CM004426.1:subseq(2597757,25152) exonerate:protein2genome:local exon 23567 23707 . - . insertions 0 ; deletions 0 CM004426.1:subseq(2597757,25152) exonerate:protein2genome:local splice5 23565 23566 . - . intron_id 2 ; splice_site "GT" CM004426.1:subseq(2597757,25152) exonerate:protein2genome:local intron 22880 23566 . - . intron_id 2 CM004426.1:subseq(2597757,25152) exonerate:protein2genome:local splice3 22880 22881 . - . intron_id 1 ; splice_site "AG" CM004426.1:subseq(2597757,25152) exonerate:protein2genome:local cds 22715 22879 . - . CM004426.1:subseq(2597757,25152) exonerate:protein2genome:local exon 22715 22879 . - . insertions 0 ; deletions 0 CM004426.1:subseq(2597757,25152) exonerate:protein2genome:local splice5 22713 22714 . - . intron_id 3 ; splice_site "GT" CM004426.1:subseq(2597757,25152) exonerate:protein2genome:local intron 22629 22714 . - . intron_id 3 CM004426.1:subseq(2597757,25152) exonerate:protein2genome:local splice3 22629 22630 . - . intron_id 2 ; splice_site "AG" CM004426.1:subseq(2597757,25152) exonerate:protein2genome:local cds 22548 22628 . - . CM004426.1:subseq(2597757,25152) exonerate:protein2genome:local exon 22548 22628 . - . insertions 0 ; deletions 0 CM004426.1:subseq(2597757,25152) exonerate:protein2genome:local similarity 22548 25110 559 - . alignment_id 1 ; Query SPP00000652_2.0 ; Align 25111 57 111 ; Align 23708 94 141 ; Align 22880 141 165 ; Align 22629 196 81 # --- END OF GFF DUMP --- # -- completed exonerate analysis