Command line: [exonerate -m p2g --showtargetgff -q ../Sel.human/GPX3.fa -t gi|397463514|gb|JH731987.1|subseq.fa] Hostname: [ws101062] C4 Alignment: ------------ Query: SPP00000021_2.0 # Protein # Glutathione peroxidase 3 (GPx3) # Human (Homo sapiens) Target: gi|397463514|gb|JH731987.1|:subseq(150000,100000) Alligator mississippiensis unplaced genomic scaffold scaffold-1586, whole genome shotgun sequence Model: protein2genome:local Raw score: 842 Query range: 0 -> 218 Target range: 57517 -> 60569 1 : MetAlaArgLeuLeuGlnAlaSerCysLeuLeuSerLeuLeuLeuAlaGlyPheValSer : 20 |||!.! !!! !|||!!.||||||!! :!!!!! !!|||:!!|||||||||.!.:!!..! MetGlyGlyGlnLeuHisAlaSerTrpIlePheProLeuValLeuAlaGlyLeuIleGln 57518 : ATGGGGGGCCAGCTGCATGCCTCCTGGATCTTCCCCTTGGTTTTGGCTGGCCTGATTCAG : 57575 21 : GlnSerArgGlyGlnGluLysSerLys >>>> Target Intron 1 >>>> MetA : 31 ! !!:! !!||||||||||||..!||| 1706 bp :!!: ProAsnTrpGlyGlnGluLysAspLys++ ++ValA 57576 : CCGAATTGGGGGCAGGAGAAAGACAAGgt.........................agGTGA : 59314 32 : spCysHisGlyGlyIleSerGlyThrIleTyrGluTyrGlyAlaLeuThrIleAspGlyG : 51 !!|||:!!.!!..!:!:.!!||||||:!!|||!!:||||||||||||||||||||||||| snCysTyrSerAsnValGlyGlyThrValTyrAspTyrGlyAlaLeuThrIleAspGlyG 59315 : ACTGCTACAGCAATGTGGGTGGCACAGTCTATGACTACGGGGCCCTCACCATCGACGGGG : 59374 52 : luGluTyrIleProPheLysGlnTyrAlaGlyLysTyrValLeuPheValAsnValAlaS : 71 ||||||||:!!|||||||||::!|||||||||||| :!!||||||||||||||||||! luGluTyrValProPheLysArgTyrAlaGlyLysMetIleLeuPheValAsnValAlaT 59375 : AAGAGTACGTCCCCTTCAAGAGGTACGCGGGGAAGATGATCCTCTTCGTCAACGTGGCCA : 59434 72 : erTyrUnkGlyLeuThrGlyGlnTyrIle{G} >>>> Target Intron 2 >>>> : 81 !!||| ||||||||| !||||||:!!{|} 471 bp hrTyr***GlyLeuThrMetGlnTyrVal{G}++ + 59435 : CGTACTGAGGGCTCACCATGCAGTACGTT{G}gt.........................a : 59934 82 : {lu}LeuAsnAlaLeuGlnGluGluLeuAlaProPheGlyLeuValIleLeuGlyPheP : 100 {||}|||||||||||||||...||||||!.!|||!:!|||||||||||||||||||||| +{lu}LeuAsnAlaLeuGlnSerGluLeuGlyProTyrGlyLeuValIleLeuGlyPheP 59935 : g{AA}CTGAATGCACTACAAAGCGAACTGGGGCCCTACGGACTCGTCATCCTGGGCTTCC : 59992 101 : roCysAsnGlnPheGlyLysGlnGluProGlyGluAsnSerGluIleLeuProThrLeu{ : 120 ||! !|||||||||||||||||||||||||||:!!|||:!!||||||||||||.!!|||{ roSerAsnGlnPheGlyLysGlnGluProGlyGlnAsnAlaGluIleLeuProAlaLeu{ 59993 : CCTCCAACCAGTTTGGGAAACAAGAACCCGGCCAGAACGCCGAGATCCTTCCTGCCCTG{ : 60052 121 : Ly} >>>> Target Intron 3 >>>> {s}TyrValArgProGlyGlyGlyPheV : 129 ||} 138 bp {|}||||||||||||||||||||||||! Ly}++ ++{s}TyrValArgProGlyGlyGlyPheA 60053 : AA}gt.........................ag{A}TACGTCCGCCCGGGGGGAGGCTTCG : 60217 130 : alProAsnPheGlnLeuPheGluLysGlyAspValAsnGlyGluLysGluGlnLysPheT : 149 .!||||||||||||||||||:!!||||||||||||||||||||||||||||||||| !!| laProAsnPheGlnLeuPheGlnLysGlyAspValAsnGlyGluLysGluGlnLysValT 60218 : CCCCCAATTTCCAGCTCTTCCAGAAGGGTGACGTGAACGGGGAGAAGGAACAGAAGGTCT : 60277 150 : yrThrPheLeuLys >>>> Target Intron 4 >>>> AsnSerCysProProTh : 159 |||||||||||||| 83 bp |||||||||||||||.. yrThrPheLeuLys++ ++AsnSerCysProProVa 60278 : ACACCTTCCTCAAGgt.........................agAACTCCTGCCCTCCCGT : 60390 160 : rSerGluLeuLeuGlyThrSerAspArgLeuPheTrpGluProMetLysValHisAspIl : 179 ! !||| !!!!||| ! !! ! !:!|||||||||||||||:!!|||:!!|||||||| lValGluAsnPheGlyAspProLysLysLeuPheTrpGluProLeuLysIleHisAspIl 60391 : GGTGGAAAACTTTGGAGACCCCAAGAAGCTTTTCTGGGAGCCCCTGAAGATCCACGACAT : 60450 180 : eArgTrpAsnPheGluLysPheLeuValGlyProAspGlyIleProIleMetArgTrpHi : 199 |!:!||||||||||||||||||||||||||||||||||||:!:|||:!:|||||||||:! eLysTrpAsnPheGluLysPheLeuValGlyProAspGlyValProValMetArgTrpTy 60451 : CAAGTGGAACTTCGAGAAGTTCCTGGTGGGGCCGGACGGGGTGCCGGTGATGCGGTGGTA : 60510 200 : sHisArgThrThrValSerAsnValLysMetAspIleLeuSerTyrMetArgArgGln : 218 !||||||.!!!.!:!!:!!!..||||||! ||||||:!!..!||||||!:! !||| rHisArgAlaAsnIleAlaThrValLysAsnAspIleMetGluTyrMetLysAlaGln 60511 : CCACCGCGCCAACATCGCCACGGTGAAGAACGACATCATGGAGTACATGAAGGCCCAG : 60569 vulgar: SPP00000021_2.0 0 218 . gi|397463514|gb|JH731987.1|:subseq(150000,100000) 57517 60569 + 842 M 29 87 5 0 2 I 0 1702 3 0 2 M 51 153 S 0 1 5 0 2 I 0 467 3 0 2 S 1 2 M 38 114 S 0 2 5 0 2 I 0 134 3 0 2 S 1 1 M 33 99 5 0 2 I 0 79 3 0 2 M 65 195 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2013-11-15 ##type DNA # # # seqname source feature start end score strand frame attributes # gi|397463514|gb|JH731987.1|:subseq(150000,100000) exonerate:protein2genome:local gene 57518 60569 842 + . gene_id 1 ; sequence SPP00000021_2.0 ; gene_orientation + gi|397463514|gb|JH731987.1|:subseq(150000,100000) exonerate:protein2genome:local cds 57518 57604 . + . gi|397463514|gb|JH731987.1|:subseq(150000,100000) exonerate:protein2genome:local exon 57518 57604 . + . insertions 0 ; deletions 0 gi|397463514|gb|JH731987.1|:subseq(150000,100000) exonerate:protein2genome:local splice5 57605 57606 . + . intron_id 1 ; splice_site "GT" gi|397463514|gb|JH731987.1|:subseq(150000,100000) exonerate:protein2genome:local intron 57605 59310 . + . intron_id 1 gi|397463514|gb|JH731987.1|:subseq(150000,100000) exonerate:protein2genome:local splice3 59309 59310 . + . intron_id 0 ; splice_site "AG" gi|397463514|gb|JH731987.1|:subseq(150000,100000) exonerate:protein2genome:local cds 59311 59464 . + . gi|397463514|gb|JH731987.1|:subseq(150000,100000) exonerate:protein2genome:local exon 59311 59464 . + . insertions 0 ; deletions 0 gi|397463514|gb|JH731987.1|:subseq(150000,100000) exonerate:protein2genome:local splice5 59465 59466 . + . intron_id 2 ; splice_site "GT" gi|397463514|gb|JH731987.1|:subseq(150000,100000) exonerate:protein2genome:local intron 59465 59935 . + . intron_id 2 gi|397463514|gb|JH731987.1|:subseq(150000,100000) exonerate:protein2genome:local splice3 59934 59935 . + . intron_id 1 ; splice_site "AG" gi|397463514|gb|JH731987.1|:subseq(150000,100000) exonerate:protein2genome:local cds 59936 60053 . + . gi|397463514|gb|JH731987.1|:subseq(150000,100000) exonerate:protein2genome:local exon 59936 60053 . + . insertions 0 ; deletions 0 gi|397463514|gb|JH731987.1|:subseq(150000,100000) exonerate:protein2genome:local splice5 60054 60055 . + . intron_id 3 ; splice_site "GT" gi|397463514|gb|JH731987.1|:subseq(150000,100000) exonerate:protein2genome:local intron 60054 60191 . + . intron_id 3 gi|397463514|gb|JH731987.1|:subseq(150000,100000) exonerate:protein2genome:local splice3 60190 60191 . + . intron_id 2 ; splice_site "AG" gi|397463514|gb|JH731987.1|:subseq(150000,100000) exonerate:protein2genome:local cds 60192 60291 . + . gi|397463514|gb|JH731987.1|:subseq(150000,100000) exonerate:protein2genome:local exon 60192 60291 . + . insertions 0 ; deletions 0 gi|397463514|gb|JH731987.1|:subseq(150000,100000) exonerate:protein2genome:local splice5 60292 60293 . + . intron_id 4 ; splice_site "GT" gi|397463514|gb|JH731987.1|:subseq(150000,100000) exonerate:protein2genome:local intron 60292 60374 . + . intron_id 4 gi|397463514|gb|JH731987.1|:subseq(150000,100000) exonerate:protein2genome:local splice3 60373 60374 . + . intron_id 3 ; splice_site "AG" gi|397463514|gb|JH731987.1|:subseq(150000,100000) exonerate:protein2genome:local cds 60375 60569 . + . gi|397463514|gb|JH731987.1|:subseq(150000,100000) exonerate:protein2genome:local exon 60375 60569 . + . insertions 0 ; deletions 0 gi|397463514|gb|JH731987.1|:subseq(150000,100000) exonerate:protein2genome:local similarity 57518 60569 842 + . alignment_id 1 ; Query SPP00000021_2.0 ; Align 57518 1 87 ; Align 59311 30 153 ; Align 59938 82 114 ; Align 60193 121 99 ; Align 60375 154 195 # --- END OF GFF DUMP --- # -- completed exonerate analysis