Command line: [exonerate -m p2g --showtargetgff -q SELENOU1a.x.fa -t SELENOU1a.extraction.KV884748.1.fa] Hostname: [sitdoc] C4 Alignment: ------------ Query: SELENOU1a Protein # Selenoprotein U1a (SELENOU1a) # Zebrafish Target: KV884748.1:subseq(1884943,100000) Monopterus albus unplaced genomic scaffold scaffold60.1, whole genome shotgun sequence Model: protein2genome:local Raw score: 781 Query range: 7 -> 222 Target range: 49644 -> 51035 8 : PheGluGlyGluThrLeuAlaMetGlyMetTrpSerLeuGlyLeuGlyAlaValGlyAla : 27 .!!||||||||| !|||! !|||! !:!!|||||||||||||||||||||||||||||| LeuGluGlyGluLeuLeuGluMetValValTrpSerLeuGlyLeuGlyAlaValGlyAla 49645 : CTCGAGGGGGAGCTGCTTGAGATGGTGGTGTGGTCTCTGGGACTGGGTGCAGTAGGTGCT : 49702 28 : AlaIleAlaGlyLeuIleLeuAlaAsnThrAspPheLeuLeuThrLysSerAlaProAla : 47 |||:!:||||||:!!.!!||||||||||||:!!.!.! !||| !!||| !||| !||| AlaLeuAlaGlyIlePheLeuAlaAsnThrAsnLeuCysLeuProLysValAlaLysAla 49703 : GCGCTGGCTGGGATCTTTTTGGCCAACACTAACCTGTGTCTGCCCAAAGTTGCCAAGGCA : 49762 48 : ThrValAspTyrLeuAlaAsnAlaAspLeu<->LysThrIleAspGly{A} >>>> Ta : 63 :!!:!!.!.||||||! !:!!||||||||| ...|||:!!|||.!!{|} SerLeuGlnTyrLeuGluAspAlaAspLeuCysSerThrValAspSer{A}++ 49763 : TCGCTGCAGTACCTCGAAGACGCTGACCTGTGCTCCACTGTAGATAGT{G}gt....... : 49816 64 : rget Intron 1 >>>> {sp}GluArgSerLeuLysAlaLysAlaLeuTrpGluLys : 75 199 bp {||}!!:!:! !:!!|||||||||::!|||||||||||| ++{sp}AspLysValIleLysAlaLysSerLeuTrpGluLys 49817 : ..................ag{AT}GATAAGGTCATCAAAGCAAAGAGCCTGTGGGAGAAG : 50048 76 : SerGlyAlaValIleMetAlaValArgArgProGlyUnkPheLeuCysArgGlu >>>> : 94 !!!|||||||||:!!||||||||||||||||||||| ||||||||||||||| ThrGlyAlaValValMetAlaValArgArgProGly***PheLeuCysArgGlu++ 50049 : ACTGGGGCCGTGGTCATGGCCGTTCGACGACCTGGATGATTTTTGTGCAGAGAGgt.... : 50107 95 : Target Intron 2 >>>> GluAlaSerGluLeuSerSerLeuLysProGlnLeuA : 106 81 bp ||||||||||||||||||||||||||||||||||||! ++GluAlaSerGluLeuSerSerLeuLysProGlnLeuG 50108 : .....................agGAGGCCTCTGAGCTGTCCTCTCTAAAGCCCCAGCTGG : 50222 107 : spGluLeuGlyValProLeuTyrAlaValValLysGluAsnValGlyThrGluIleGlnA : 126 !:|||||||||||||||||| !|||||||||||||||:!!||||||.!!||||||::!| luGluLeuGlyValProLeuValAlaValValLysGluAspValGlyAlaGluIleArgA 50223 : AAGAGCTCGGGGTCCCTCTGGTTGCTGTGGTGAAGGAGGACGTTGGCGCCGAGATCAGGG : 50282 127 : spPheArgProHisPheAlaGlyGluIlePheLeuAspGluLys >>>> Target In : 141 |||||||||||||||||||||||!!:|||!:!:!!|||:!!||| 368 spPheArgProHisPheAlaGlyAspIleTyrValAspLysLys++ 50283 : ACTTCAGACCGCACTTTGCTGGGGACATCTATGTAGACAAAAAGgt.............. : 50329 142 : tron 3 >>>> GlnAlaPheTyrGlyProGlnGlnArgLysMetGlyGlyLeuGlyPh : 156 bp :!! !||||||||||||! !||||||!:!||||||!.!||||||.! ++LysHisPheTyrGlyProLeuGlnArgArgMetGlyAlaLeuGlyLe 50330 : ...........agAAACATTTCTACGGCCCGCTGCAGAGGAGGATGGGGGCTCTGGGGCT : 50740 157 : eIleArgLeuGlyValTrpGlnAsnPheValArgAlaTrpArgAlaGlyTyrGlnGlyAs : 176 !|||||||||||||||||||||||||||:!!||||||||||||:!!|||||||||||||| uIleArgLeuGlyValTrpGlnAsnPheIleArgAlaTrpArgSerGlyTyrGlnGlyAs 50741 : CATCCGCCTCGGGGTGTGGCAGAACTTCATACGGGCCTGGAGGTCTGGTTACCAGGGCAA : 50800 177 : nMetAsnGlyGluGlyPheIleLeuGlyGlyValPheValMetGlySerGlyGlyGln : 196 |:!!||||||||||||||||||||||||||||||||||||!!:||| !!|||! !||| nLeuAsnGlyGluGlyPheIleLeuGlyGlyValPheValIleGlyProGlyAspGln++ 50801 : CTTGAACGGCGAGGGCTTCATCCTCGGGGGAGTGTTTGTCATTGGACCAGGAGACCAGgt : 50860 197 : >>>> Target Intron 4 >>>> GlyValLeuLeuGluHisArgGluLysGluPhe : 206 95 bp |||:!!||||||||||||||||||||||||||| ++GlyIleLeuLeuGluHisArgGluLysGluPhe 50861 : .........................agGGGATTCTCCTGGAACACCGTGAGAAGGAGTTT : 50985 207 : GlyAspLysValSerLeuGluSerValLeuGluAlaAlaLysLysVal : 222 ||||||||||||!:! !! !..!||||||||||||!.!||||||:!! GlyAspLysValAsnThrAlaAspValLeuGluAlaValLysLysIle 50986 : GGTGATAAGGTGAACACTGCAGATGTTTTAGAGGCTGTCAAGAAAATA : 51035 vulgar: SELENOU1a 7 222 . KV884748.1:subseq(1884943,100000) 49644 51035 + 781 M 50 150 G 0 3 M 5 15 S 0 1 5 0 2 I 0 195 3 0 2 S 1 2 M 30 90 5 0 2 I 0 77 3 0 2 M 47 141 5 0 2 I 0 364 3 0 2 M 55 165 5 0 2 I 0 91 3 0 2 M 27 81 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2017-11-24 ##type DNA # # # seqname source feature start end score strand frame attributes # KV884748.1:subseq(1884943,100000) exonerate:protein2genome:local gene 49645 51035 781 + . gene_id 1 ; sequence SELENOU1a ; gene_orientation + KV884748.1:subseq(1884943,100000) exonerate:protein2genome:local cds 49645 49813 . + . KV884748.1:subseq(1884943,100000) exonerate:protein2genome:local exon 49645 49813 . + . insertions 3 ; deletions 0 KV884748.1:subseq(1884943,100000) exonerate:protein2genome:local splice5 49814 49815 . + . intron_id 1 ; splice_site "GT" KV884748.1:subseq(1884943,100000) exonerate:protein2genome:local intron 49814 50012 . + . intron_id 1 KV884748.1:subseq(1884943,100000) exonerate:protein2genome:local splice3 50011 50012 . + . intron_id 0 ; splice_site "AG" KV884748.1:subseq(1884943,100000) exonerate:protein2genome:local cds 50013 50104 . + . KV884748.1:subseq(1884943,100000) exonerate:protein2genome:local exon 50013 50104 . + . insertions 0 ; deletions 0 KV884748.1:subseq(1884943,100000) exonerate:protein2genome:local splice5 50105 50106 . + . intron_id 2 ; splice_site "GT" KV884748.1:subseq(1884943,100000) exonerate:protein2genome:local intron 50105 50185 . + . intron_id 2 KV884748.1:subseq(1884943,100000) exonerate:protein2genome:local splice3 50184 50185 . + . intron_id 1 ; splice_site "AG" KV884748.1:subseq(1884943,100000) exonerate:protein2genome:local cds 50186 50326 . + . KV884748.1:subseq(1884943,100000) exonerate:protein2genome:local exon 50186 50326 . + . insertions 0 ; deletions 0 KV884748.1:subseq(1884943,100000) exonerate:protein2genome:local splice5 50327 50328 . + . intron_id 3 ; splice_site "GT" KV884748.1:subseq(1884943,100000) exonerate:protein2genome:local intron 50327 50694 . + . intron_id 3 KV884748.1:subseq(1884943,100000) exonerate:protein2genome:local splice3 50693 50694 . + . intron_id 2 ; splice_site "AG" KV884748.1:subseq(1884943,100000) exonerate:protein2genome:local cds 50695 50859 . + . KV884748.1:subseq(1884943,100000) exonerate:protein2genome:local exon 50695 50859 . + . insertions 0 ; deletions 0 KV884748.1:subseq(1884943,100000) exonerate:protein2genome:local splice5 50860 50861 . + . intron_id 4 ; splice_site "GT" KV884748.1:subseq(1884943,100000) exonerate:protein2genome:local intron 50860 50954 . + . intron_id 4 KV884748.1:subseq(1884943,100000) exonerate:protein2genome:local splice3 50953 50954 . + . intron_id 3 ; splice_site "AG" KV884748.1:subseq(1884943,100000) exonerate:protein2genome:local cds 50955 51035 . + . KV884748.1:subseq(1884943,100000) exonerate:protein2genome:local exon 50955 51035 . + . insertions 0 ; deletions 0 KV884748.1:subseq(1884943,100000) exonerate:protein2genome:local similarity 49645 51035 781 + . alignment_id 1 ; Query SELENOU1a ; Align 49645 8 150 ; Align 49798 58 15 ; Align 50015 64 90 ; Align 50186 94 141 ; Align 50695 141 165 ; Align 50955 196 81 # --- END OF GFF DUMP --- # -- completed exonerate analysis