Command line: [exonerate -m p2g --showtargetgff -q SPP00000015_1.0.fa -t SPP00000015_1.0.GL397313.1/GL397313.1.subseq] Hostname: [upf-OptiPlex-755] C4 Alignment: ------------ Query: SPP00000015_1.0 # Protein # Selenoprotein I (SelI) # Homo sapiens # Complete Target: GL397313.1:subseq(6485358,100140) dna:supercontig supercontig:Nleu1.0:GL397313.1:1:16672385:1 REF Model: protein2genome:local Raw score: 198 Query range: 6 -> 132 Target range: 38060 -> 50140 7 : ValSerProGluGlnLeuAlaGlyPheAspLysTyrLysTyrSerAlaValAspThrAsn : 26 :!!||| !.!.|||||| ! !!.!.!!::!!:!!!:!|||..!::!!.!! !::! LeuSerArgHisGlnLeuLysArgLeuGluGluHisArgTyrGlnSerAlaGlyArgSer 38061 : TTGTCAAGACACCAACTAAAGCGGCTAGAAGAACACAGATATCAAAGTGCTGGACGGTCC : 38118 27 : ProLeuSerLeuTyrValMetHisProPheTrpAsnThrIleValLysValPheProThr : 46 ! !|||..! !:!!|||!!. !!:!|||.!. !:!!|||!:! ! !!|||:!! LeuLeuGlu---ProLeuMetGlnGlyTyrTrpGluTrpLeuValArgArgValProSer 38119 : CTGCTTGAG---CCCTTAATGCAAGGGTATTGGGAATGGCTCGTTAGAAGAGTTCCCTCC : 38175 47 : TrpLeuAlaProAsnLeuIleThrPheSerGlyPheLeuLeuValValPheAsnPheLeu : 66 |||:!!||||||||||||||||||.!!! !|||.!.! !:!! !:!!! !!.. !:!! TrpIleAlaProAsnLeuIleThrIleIleGlyLeuSerIleAsnIleCysThrThrIle 38176 : TGGATTGCCCCAAATCTCATCACCATCATTGGACTGTCAATAAACATCTGTACAACTATT : 38235 67 : LeuMetAlaTyrPheAspProAspPheTyrAlaSer >>>> Target Intron 1 > : 79 |||:!:!.!!:!!:! !||| ! ! !..! 11705 bp LeuLeuValPheTyrCysProThrAlaThrGluGln++ 38236 : TTATTAGTCTTCTACTGCCCTACAGCTACAGAGCAGgt...................... : 38276 80 : >>> AlaProGlyHisLysHisValProAspTrpValTrpIleValValGlyIleLeuA : 97 :!! ! ! !..!!!.!.!||| |||!.!!::|||!.! !!.! ! ! --SerTyrPheIleSerGlnAlaProLeuTrpAlaTyrIleAlaCysAlaCysGlyL 38277 : ...ttTCTTATTTTATTTCACAGGCACCTCTGTGGGCATATATTGCTTGTGCCTGTGGCC : 50033 98 : snPheValAlaTyrThrLeuAspGlyValAspGlyLysGlnAlaArgArgThrAsnSerS : 117 !|||:!! ! ! :!!||||||!.!:!!||||||||||||||||||||||||||||||| euPheIleTyrGlnSerLeuAspAlaIleAspGlyLysGlnAlaArgArgThrAsnSerS 50034 : TTTTCATTTACCAGTCTTTGGATGCTATTGATGGGAAACAGGCAAGAAGAACCAATAGTA : 50093 118 : erThrProLeuGlyGluLeuPheAspHisGlyLeuAspSerTrpSer : 132 ||:!!|||||||||||||||||||||||||||! !|||||| ||| erSerProLeuGlyGluLeuPheAspHisGlyCysAspSerLeuSer 50094 : GTTCTCCTCTGGGAGAACTTTTTGATCATGGTTGTGATTCACTATCA : 50140 vulgar: SPP00000015_1.0 6 132 . GL397313.1:subseq(6485358,100140) 38060 50140 + 198 M 23 69 G 1 0 M 48 144 5 0 2 I 0 11701 3 0 2 M 54 162 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2012-11-27 ##type DNA # # # seqname source feature start end score strand frame attributes # GL397313.1:subseq(6485358,100140) exonerate:protein2genome:local gene 38061 50140 198 + . gene_id 1 ; sequence SPP00000015_1.0 ; gene_orientation + GL397313.1:subseq(6485358,100140) exonerate:protein2genome:local cds 38061 38273 . + . GL397313.1:subseq(6485358,100140) exonerate:protein2genome:local exon 38061 38273 . + . insertions 0 ; deletions 1 GL397313.1:subseq(6485358,100140) exonerate:protein2genome:local splice5 38274 38275 . + . intron_id 1 ; splice_site "GT" GL397313.1:subseq(6485358,100140) exonerate:protein2genome:local intron 38274 49978 . + . intron_id 1 GL397313.1:subseq(6485358,100140) exonerate:protein2genome:local splice3 49977 49978 . + . intron_id 0 ; splice_site "TT" GL397313.1:subseq(6485358,100140) exonerate:protein2genome:local cds 49979 50140 . + . GL397313.1:subseq(6485358,100140) exonerate:protein2genome:local exon 49979 50140 . + . insertions 0 ; deletions 0 GL397313.1:subseq(6485358,100140) exonerate:protein2genome:local similarity 38061 50140 198 + . alignment_id 1 ; Query SPP00000015_1.0 ; Align 38061 7 69 ; Align 38130 31 144 ; Align 49979 79 162 # --- END OF GFF DUMP --- # -- completed exonerate analysis