Command line: [exonerate -m p2g --showtargetgff -q SPP00000008_1.0.fa -t SPP00000008_1.0.GL397305.1/GL397305.1.subseq] Hostname: [upf-OptiPlex-755] C4 Alignment: ------------ Query: SPP00000008_1.0 # Protein # Glutathione peroxidase 5 (GPx5) # Homo sapiens # Complete Target: GL397305.1:subseq(18711209,102757) dna:supercontig supercontig:Nleu1.0:GL397305.1:1:19291937:1 REF Model: protein2genome:local Raw score: 744 Query range: 27 -> 217 Target range: 49996 -> 52757 28 : MetLysMetAspCysHisLysAspGluLysGlyThrIleTyrAspTyrGluAlaIleAla : 47 .!.:!!|||||||||||| ! ! !!...!!|||||||||!!:|||! !|||:!!.!! PheGlnMetAspCysHisGlyGlyMetSerSerThrIleTyrGluTyrGlyAlaLeuThr 49997 : TTCCAGATGGACTGCCATGGTGGCATGAGCAGCACCATTTATGAGTACGGAGCCCTCACC : 50054 48 : LeuAsnLysAsnGluTyrValSerPheLysGlnTyrValGlyLysHisIleLeuPheVal : 67 :!!:!! !.!.||||||:!! !!||||||||||||!.!||||||:!!:!!||||||||| IleAspGlyGluGluTyrIleProPheLysGlnTyrAlaGlyLysTyrValLeuPheVal 50055 : ATTGATGGGGAGGAGTACATCCCCTTCAAGCAGTATGCTGGCAAATACGTCCTCTTTGTC : 50114 68 : AsnValAlaThrTyrCysGlyLeuThrAlaGlnTyrPro{G} >>>> Target Intr : 81 |||||||||!:!|||!! |||||||||!.!|||||| !{|} 1390 bp AsnValAlaSerTyr***GlyLeuThrGlyGlnTyrIle{G}++ 50115 : AACGTGGCCAGCTACTGAGGCCTGACAGGCCAGTACATT{G}gt................ : 50159 82 : on 1 >>>> {lu}LeuAsnAlaLeuGlnGluGluLeuLysProTyrGlyLeuValVal : 96 {||}|||||||||||||||||||||||| !|||!:!|||||||||:!! ++{lu}LeuAsnAlaLeuGlnGluGluLeuAlaProPheGlyLeuValIle 50160 : .........ag{AA}CTGAATGCACTACAGGAAGAGCTTGCACCATTTGGTCTGGTCATT : 51591 97 : LeuGlyPheProCysAsnGlnPheGlyLysGlnGluProGlyAspAsnLysGluIleLeu : 116 ||||||||||||||||||||||||||||||||||||||||||!!:|||..!||||||||| LeuGlyPheProCysAsnGlnPheGlyLysGlnGluProGlyGluAsnSerGluIleLeu 51592 : CTGGGCTTTCCCTGCAACCAGTTTGGAAAACAGGAACCAGGAGAGAACTCAGAGATCCTT : 51651 117 : ProGlyLeu{Ly} >>>> Target Intron 2 >>>> {s}TyrValArgProGly : 125 ||| !|||{||} 415 bp {|}||||||||||||||| ProThrLeu{Ly}++ ++{s}TyrValArgProGly 51652 : CCCACCCTC{AA}gt.........................ag{G}TATGTCCGACCAGGT : 52093 126 : GlyGlyPheValProSerPheGlnLeuPheGluLysGlyAspValAsnGlyGluLysGlu : 145 |||||||||||||||!:!|||||||||||||||||||||||||||||||||||||||||| GlyGlyPheValProAsnPheGlnLeuPheGluLysGlyAspValAsnGlyGluLysGlu 52094 : GGAGGCTTTGTCCCTAATTTCCAGCTCTTTGAGAAAGGGGATGTCAATGGAGAGAAAGAG : 52153 146 : GlnLysValPheSerPheLeuLys >>>> Target Intron 3 >>>> HisSerC : 156 |||||| !!!:!!!!||||||||| 386 bp :!!|||| GlnLysPheTyrThrPheLeuLys++ ++AsnSerC 52154 : CAGAAATTCTACACTTTCCTAAAGgt.........................agAACTCCT : 52572 157 : ysProHisProSerGluIleLeuGlyThrPheLysSerIleSerTrpAspProValLysV : 176 |||||! ! !!||||||:!!|||||||||! ! ! !!:!!! !|||!!:|||:!!|||| ysProProThrSerGluLeuLeuGlyThrSerAspArgLeuPheTrpGluProMetLysV 52573 : GTCCTCCCACCTCGGAGCTCCTGGGTACATCTGACCGCCTCTTCTGGGAACCCATGAAGG : 52632 177 : alHisAspIleArgTrpAsnPheGluLysPheLeuValGlyProAspGlyIleProValM : 196 ||||||||||||||||||||||||||||||||||||||||||||||||||! !|||:!!| alHisAspIleArgTrpAsnPheGluLysPheLeuValGlyProAspGlyThrProIleM 52633 : TTCACGACATCCGCTGGAACTTTGAGAAGTTCCTGGTGGGGCCAGATGGTACACCCATCA : 52692 197 : etArgTrpSerHisArgAlaThrValSerSerValLysThrAspIleLeuAlaTyrLeuL : 216 |||||||| !||||||.!!|||||||||!:!||||||! !|||||||||:!!|||:!!! etArgTrpHisHisArgThrThrValSerAsnValLysMetAspIleLeuSerTyrMetA 52693 : TGCGCTGGCACCACCGGACCACGGTCAGCAACGTCAAGATGGACATCCTGTCCTACATGA : 52752 217 : ysGln : 217 :!!:! rgArg 52753 : GGCGG : 52757 vulgar: SPP00000008_1.0 27 217 . GL397305.1:subseq(18711209,102757) 49996 52757 + 744 M 53 159 S 0 1 5 0 2 I 0 1386 3 0 2 S 1 2 M 38 114 S 0 2 5 0 2 I 0 411 3 0 2 S 1 1 M 33 99 5 0 2 I 0 382 3 0 2 M 64 192 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2012-11-30 ##type DNA # # # seqname source feature start end score strand frame attributes # GL397305.1:subseq(18711209,102757) exonerate:protein2genome:local gene 49997 52757 744 + . gene_id 1 ; sequence SPP00000008_1.0 ; gene_orientation + GL397305.1:subseq(18711209,102757) exonerate:protein2genome:local cds 49997 50156 . + . GL397305.1:subseq(18711209,102757) exonerate:protein2genome:local exon 49997 50156 . + . insertions 0 ; deletions 0 GL397305.1:subseq(18711209,102757) exonerate:protein2genome:local splice5 50157 50158 . + . intron_id 1 ; splice_site "GT" GL397305.1:subseq(18711209,102757) exonerate:protein2genome:local intron 50157 51546 . + . intron_id 1 GL397305.1:subseq(18711209,102757) exonerate:protein2genome:local splice3 51545 51546 . + . intron_id 0 ; splice_site "AG" GL397305.1:subseq(18711209,102757) exonerate:protein2genome:local cds 51547 51664 . + . GL397305.1:subseq(18711209,102757) exonerate:protein2genome:local exon 51547 51664 . + . insertions 0 ; deletions 0 GL397305.1:subseq(18711209,102757) exonerate:protein2genome:local splice5 51665 51666 . + . intron_id 2 ; splice_site "GT" GL397305.1:subseq(18711209,102757) exonerate:protein2genome:local intron 51665 52079 . + . intron_id 2 GL397305.1:subseq(18711209,102757) exonerate:protein2genome:local splice3 52078 52079 . + . intron_id 1 ; splice_site "AG" GL397305.1:subseq(18711209,102757) exonerate:protein2genome:local cds 52080 52179 . + . GL397305.1:subseq(18711209,102757) exonerate:protein2genome:local exon 52080 52179 . + . insertions 0 ; deletions 0 GL397305.1:subseq(18711209,102757) exonerate:protein2genome:local splice5 52180 52181 . + . intron_id 3 ; splice_site "GT" GL397305.1:subseq(18711209,102757) exonerate:protein2genome:local intron 52180 52565 . + . intron_id 3 GL397305.1:subseq(18711209,102757) exonerate:protein2genome:local splice3 52564 52565 . + . intron_id 2 ; splice_site "AG" GL397305.1:subseq(18711209,102757) exonerate:protein2genome:local cds 52566 52757 . + . GL397305.1:subseq(18711209,102757) exonerate:protein2genome:local exon 52566 52757 . + . insertions 0 ; deletions 0 GL397305.1:subseq(18711209,102757) exonerate:protein2genome:local similarity 49997 52757 744 + . alignment_id 1 ; Query SPP00000008_1.0 ; Align 49997 28 159 ; Align 51549 82 114 ; Align 52081 121 99 ; Align 52566 154 192 # --- END OF GFF DUMP --- # -- completed exonerate analysis