Command line: [exonerate -m p2g --showtargetgff -q SELENOU.fa -t SELENOU.subseq] Hostname: [ws101098U] C4 Alignment: ------------ Query: SPP00000109_2.0 # Protein # Selenoprotein U1 (SELENOU1) # Human Target: JSZB01001214.1:subseq(190000,100000) Manis javanica isolate MP_PG03-UM SCAFFOLD1338, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 941 Query range: 0 -> 229 Target range: 49995 -> 49254 1 : MetSerPheLeuGlnAspProSerPhePheThrMetGlyMetTrpSerIleGlyAlaGly : 20 |||||||||||||||||| !!||||||||||||!!:||||||||||||||||||!.!||| MetSerPheLeuGlnAspThrSerPhePheThrIleGlyMetTrpSerIleGlyValGly 49995 : ATGTCTTTCCTCCAGGATACAAGTTTCTTCACCATAGGAATGTGGTCCATTGGTGTGGGG : 49938 21 : AlaLeuGlyAlaAlaAlaLeuAlaLeuLeuLeuAlaAsnThrAspValPheLeuSerLys : 40 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| AlaLeuGlyAlaAlaAlaLeuAlaLeuLeuLeuAlaAsnThrAspValPheLeuSerLys 49937 : GCTCTTGGGGCTGCTGCCTTAGCACTGCTCCTGGCCAACACAGACGTGTTTCTGTCCAAG : 49878 41 : ProGlnLysAlaAlaLeuGluTyrLeuGluAspIleAspLeuLysThrLeuGluLysGlu : 60 !!|||||||||.!!||||||||||||||||||||||||||||||||||||||||||||| ThrGlnLysAlaThrLeuGluTyrLeuGluAspIleAspLeuLysThrLeuGluLysGlu 49877 : ACCCAGAAAGCAACGCTGGAGTATCTGGAGGATATAGACCTGAAAACTCTGGAGAAGGAA : 49818 61 : ProArgThrPheLysAlaLysGluLeuTrpGluLysAsnGlyAlaValIleMetAlaVal : 80 ||||||||||||||||||||||||||||||||||||!:!||||||||||||||||||||| ProArgThrPheLysAlaLysGluLeuTrpGluLysSerGlyAlaValIleMetAlaVal 49817 : CCAAGGACTTTCAAAGCAAAGGAGCTCTGGGAAAAGAGTGGGGCTGTGATTATGGCTGTG : 49758 81 : ArgArgProGlyCysPheLeuCysArgGluGluAlaAlaAspLeuSerSerLeuLysSer : 100 !!|||||||||||||||||||||||||||||||||! !|||||||||||||||||| !! TrpArgProGlyCysPheLeuCysArgGluGluAlaGluAspLeuSerSerLeuLysPro 49757 : TGGAGGCCAGGCTGTTTCCTATGTCGAGAGGAGGCTGAGGACCTGTCCTCCCTGAAGCCC : 49698 101 : MetLeuAspGlnLeuGlyValProLeuTyrAla{V} >>>> Target Intron 1 > : 112 ! !||||||:!!|||||||||||||||||||||{|} 55 bp LysLeuAspGluLeuGlyValProLeuTyrAla{V}-- 49697 : AAGTTGGATGAGCTGGGTGTCCCCCTCTATGCA{G}tg...................... : 49659 113 : >>> {al}ValLysGluHisIleArgThrGluValLysAspPheGlnProTyrPheLys : 129 {||}|||||||||!!.|||||||||||||||:!!||||||||||||||||||||| ++{al}ValLysGluGlnIleArgThrGluValGluAspPheGlnProTyrPheLys 49658 : ...ag{TG}GTGAAGGAGCAGATTAGGACTGAAGTGGAGGACTTCCAGCcttatttcaaa : 49556 130 : GlyGluIlePheLeuAspGluLysLysLysPheTyrGlyProGlnArgArgLysMetMet : 149 ||||||||||||||||||||||||||||||||||||! !|||||||||!:!|||||| GluGluIlePheLeuAspGluLysLysLysPheTyrGlyLeuGlnArgArgArgMetMet 49555 : gaagaaatcttcctggatgaaaagaaaaagttctaTGGTCTCCAAAGGCGGAGGATGATG : 49496 150 : PheMetGlyPheIleArgLeuGlyValTrpTyrAsnPhePheArgAlaTrpAsnGlyGly : 169 |||:!!||||||:!!|||||||||||||||||||||||||||!:!||||||||||||||| PheLeuGlyPheValArgLeuGlyValTrpTyrAsnPhePheGlnAlaTrpAsnGlyGly 49495 : TTTCTGGGATTTGTCCGTCTGGGTGTCTGGTACAACTTCTTCCAGGCCTGGAATGGAGGC : 49436 170 : PheSer--GlyAsnLeuGluGlyGluGlyPheIleLeuGlyGlyValPheValValGlyS : 189 .!!|||## |||||||||||||||||||||||||||||||||||||||:!!||||||| IleSer##---AsnLeuGluGlyGluGlyPheIleLeuGlyGlyValPheMetValGlyS 49435 : ATCTCTGG---AACCTGGAAGGCGAAGGCTTCATCCTTGGGGGAGTTTTTATGGTGGGAT : 49377 190 : erGlyLysGlnGlyIleLeuLeuGluHisArgGluLysGluPheGlyAspLysValAsnL : 209 ||||||||:!!||||||||||||||||||:::|||||||||||||||||||||!.!!!.! erGlyLysGluGlyIleLeuLeuGluHisGlnGluLysGluPheGlyAspLysAlaLysP 49376 : CAGGAAAGGAGGGCATTCTTCTTGAGcatcaagaaaaagaatttgGAGACAAAGCAAAGC : 49317 210 : euLeuSerValLeuGluAlaAlaLysMetIleLysProGlnThrLeuAlaSerGluLysL : 229 !:!!||||||||||||!.!|||!:!! !|||||| !||||||! !||||||||||||| roIleSerValLeuGluValAlaArgLysIleLysValGlnThrSerAlaSerGluLysL 49316 : CAATTTCTGTTCTGGAAGTTGCTAGGAAGATCAAAGTACAGACTTCAGCCTCAGAGAAAA : 49257 230 : ys : 229 || ys 49256 : AA : 49255 vulgar: SPP00000109_2.0 0 229 . JSZB01001214.1:subseq(190000,100000) 49995 49254 - 941 M 111 333 S 0 1 5 0 2 I 0 51 3 0 2 S 1 2 M 59 177 F 0 2 G 1 0 M 57 171 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2016-11-09 ##type DNA # # # seqname source feature start end score strand frame attributes # JSZB01001214.1:subseq(190000,100000) exonerate:protein2genome:local gene 49255 49995 941 - . gene_id 1 ; sequence SPP00000109_2.0 ; gene_orientation + JSZB01001214.1:subseq(190000,100000) exonerate:protein2genome:local cds 49662 49995 . - . JSZB01001214.1:subseq(190000,100000) exonerate:protein2genome:local exon 49662 49995 . - . insertions 0 ; deletions 0 JSZB01001214.1:subseq(190000,100000) exonerate:protein2genome:local splice5 49660 49661 . - . intron_id 1 ; splice_site "TG" JSZB01001214.1:subseq(190000,100000) exonerate:protein2genome:local intron 49607 49661 . - . intron_id 1 JSZB01001214.1:subseq(190000,100000) exonerate:protein2genome:local splice3 49607 49608 . - . intron_id 0 ; splice_site "AG" JSZB01001214.1:subseq(190000,100000) exonerate:protein2genome:local cds 49255 49606 . - . JSZB01001214.1:subseq(190000,100000) exonerate:protein2genome:local exon 49255 49606 . - . insertions 0 ; deletions 1 ; frameshifts 2 JSZB01001214.1:subseq(190000,100000) exonerate:protein2genome:local similarity 49255 49995 941 - . alignment_id 1 ; Query SPP00000109_2.0 ; Align 49996 1 333 ; Align 49605 113 177 ; Align 49426 173 171 # --- END OF GFF DUMP --- # -- completed exonerate analysis