Command line: [exonerate --exhaustive yes -m p2g --showtargetgff -q ../GPx5-human.fa -t GPx5-NIVO01038789.1.subseq] Hostname: [sitdoc] C4 Alignment: ------------ Query: SPP00000008_1.0 # Protein # Glutathione peroxidase 5 (GPx5) # Homo sapiens # Complete Target: NIVO01038789.1:subseq(406726,133103) Ammotragus lervia isolate NAZ scaffold232_len768847_cov70, whole genome shotgun sequence Model: protein2genome:local Raw score: 906 Query range: 0 -> 221 Target range: 75124 -> 83103 1 : MetThrThrGlnLeuArgValValHisLeuLeuProLeuLeuLeuAlaCysPheValGln : 20 ||||||..!|||||||||!.! |||... |||...|||||| !!||||||||| MetThrValGlnLeuArgAlaSerCysLeuPheLeuLeuPheLeuAlaGlyPheValGln 75125 : ATGACTGTACAGTTAAGGGCCtcttgtctttttctccttttcctagcTGGCTTTGTGCAA : 75182 21 : ThrSerProLysGlnGluLys >>>> Target Intron 1 >>>> MetLys >> : 30 |||!:! !!..! !|||||| 230 bp :!::!! ThrAsnTyrSerLeuGluLys++ ++ValGlu++ 75183 : ACAAATTACAGTCTAGAGAAGgt.........................agGTAGAGgt.. : 75444 31 : >> Target Intron 2 >>>> MetAspCysHisLysAspGluLysGlyThrIleTy : 41 2570 bp |||||||||:!!||||||! !|||||||||||||| ++MetAspCysTyrLysAspValLysGlyThrIleTy 75445 : .......................agATGGATTGCTACAAAGATGTGAAAGGCACCATCTA : 78045 42 : rAspTyrGluAlaIleAlaLeuAsnLysAsnGluTyrValSerPheLysGlnTyrValGl : 61 |||||||!!:|||.!!.!!|||||| !!!.|||:!!:!!..!||||||||||||!.!|| rAspTyrAspAlaPheThrLeuAsnGlyLysGluHisIleGlnPheLysGlnTyrAlaGl 78046 : TGATTATGACGCTTTTACTCTTAATGGAAAGGAACACATTCAGTTCAAGCAGTATGCGGG : 78105 62 : yLysHisIleLeuPheValAsnValAlaThrTyrCysGlyLeuThrAlaGlnTyrPro{G : 81 |||||||:!!||||||||||||||||||||||||||||||||||||||||||||||||{| yLysHisValLeuPheValAsnValAlaThrTyrCysGlyLeuThrAlaGlnTyrPro{G 78106 : CAAACATGTCCTTTTTGTCAATGTGGCCACCTATTGTGGTCTGACAGCTCAATATCCT{G : 78165 82 : } >>>> Target Intron 3 >>>> {lu}LeuAsnAlaLeuGlnGluGluLeuLy : 90 } 2920 bp {||}|||||||||||||||||||||||||| }++ ++{lu}LeuAsnAlaLeuGlnGluGluLeuLy 78166 : }gt.........................ag{AA}CTGAATGCACTACAGGAGGAGCTGAA : 81112 91 : sProTyrGlyLeuValValLeuGlyPheProCysAsnGlnPheGlyLysGlnGluProGl : 110 ||||!:!||||||||||||||||||||||||||||||||||||||||||||||||||||| sProPheGlyLeuValValLeuGlyPheProCysAsnGlnPheGlyLysGlnGluProGl 81113 : GCCCTTTGGCTTAGTTGTGTTGGGATTTCCCTGTAACCAATTCGGAAAGCAAGAACCAGG : 81172 111 : yAspAsnLysGluIleLeuProGlyLeu{Ly} >>>> Target Intron 4 >>>> : 120 |!!:|||..!||||||||||||||||||{||} 481 bp yGluAsnSerGluIleLeuProGlyLeu{Ly}++ + 81173 : AGAAAACTCAGAAATCCTTCCAGGACTG{AA}gt.........................a : 81683 121 : {s}TyrValArgProGlyGlyGlyPheValProSerPheGlnLeuPheGluLysGlyAs : 139 {|}|||||||||||||||||||||!:!||||||!:!||||||||||||||||||||||| +{s}TyrValArgProGlyGlyGlyTyrValProAsnPheGlnLeuPheGluLysGlyAs 81684 : g{G}TATGTCCGTCCAGGAGGAGGATACGTACCTAATTTCCAGCTGTTTGAGAAAGGGGA : 81740 140 : pValAsnGlyGluLysGluGlnLysValPheSerPheLeuLys >>>> Target Int : 154 |||||||||||||! !|||||||||||||||!!!||||||||| 1115 b pValAsnGlyGluThrGluGlnLysValPheThrPheLeuLys++ 81741 : TGTGAATGGTGAAACAGAGCAGAAAGTCTTCACCTTCTTAAAGgt............... : 81787 155 : ron 5 >>>> HisSerCysProHisProSerGluIleLeuGlyThrPheLysSerIle : 169 p !!.|||||||||||||||||||||.!!:!!|||:!!.!!||| !||| ++GlnSerCysProHisProSerGluPheMetGlySerIleLysHisIle 81788 : ..........agCAATCCTGTCCTCACCCTTCTGAGTTTATGGGCTCAATCAAACATATA : 82945 170 : SerTrpAspProValLysValHisAspIleArgTrpAsnPheGluLysPheLeuValGly : 189 ||||||!!:|||:!!! !|||!.!|||||||||||||||||||||||||||||||||||| SerTrpGluProIleMetValArgAspIleArgTrpAsnPheGluLysPheLeuValGly 82946 : TCCTGGGAACCCATCATGGTCCGTGACATCCGCTGGAATTTTGAAAAGTTCCTAGTGGGA : 83005 190 : ProAspGlyIleProValMetArgTrpSerHisArgAlaThrValSerSerValLysThr : 209 |||||||||:!!|||||||||||||||! !||||||.!! !!||||||!!!||||||:!! ProAspGlyValProValMetArgTrpPheHisArgThrProValSerThrValLysSer 83006 : CCTGATGGCGTCCCTGTCATGCGCTGGTTTCATCGGACTCCAGTCAGTACAGTCAAGTCA : 83065 210 : AspIleLeuAlaTyrLeuLysGlnPheLysThrLys : 221 |||||||||||||||:!!|||||||||:!!|||||| AspIleLeuAlaTyrMetLysGlnPheGlnThrLys 83066 : GATATTCTGGCATACATGAAACAGTTCCAAACCAAA : 83103 vulgar: SPP00000008_1.0 0 221 . NIVO01038789.1:subseq(406726,133103) 75124 83103 + 906 M 27 81 5 0 2 I 0 226 3 0 2 M 2 6 5 0 2 I 0 2566 3 0 2 M 51 153 S 0 1 5 0 2 I 0 2916 3 0 2 S 1 2 M 38 114 S 0 2 5 0 2 I 0 477 3 0 2 S 1 1 M 33 99 5 0 2 I 0 1111 3 0 2 M 68 204 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2017-11-17 ##type DNA # # # seqname source feature start end score strand frame attributes # NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local gene 75125 83103 906 + . gene_id 1 ; sequence SPP00000008_1.0 ; gene_orientation + NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local cds 75125 75205 . + . NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local exon 75125 75205 . + . insertions 0 ; deletions 0 NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local splice5 75206 75207 . + . intron_id 1 ; splice_site "GT" NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local intron 75206 75435 . + . intron_id 1 NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local splice3 75434 75435 . + . intron_id 0 ; splice_site "aG" NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local cds 75436 75441 . + . NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local exon 75436 75441 . + . insertions 0 ; deletions 0 NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local splice5 75442 75443 . + . intron_id 2 ; splice_site "GT" NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local intron 75442 78011 . + . intron_id 2 NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local splice3 78010 78011 . + . intron_id 1 ; splice_site "ag" NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local cds 78012 78165 . + . NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local exon 78012 78165 . + . insertions 0 ; deletions 0 NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local splice5 78166 78167 . + . intron_id 3 ; splice_site "GT" NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local intron 78166 81085 . + . intron_id 3 NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local splice3 81084 81085 . + . intron_id 2 ; splice_site "AG" NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local cds 81086 81203 . + . NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local exon 81086 81203 . + . insertions 0 ; deletions 0 NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local splice5 81204 81205 . + . intron_id 4 ; splice_site "GT" NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local intron 81204 81684 . + . intron_id 4 NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local splice3 81683 81684 . + . intron_id 3 ; splice_site "AG" NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local cds 81685 81784 . + . NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local exon 81685 81784 . + . insertions 0 ; deletions 0 NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local splice5 81785 81786 . + . intron_id 5 ; splice_site "GT" NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local intron 81785 82899 . + . intron_id 5 NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local splice3 82898 82899 . + . intron_id 4 ; splice_site "AG" NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local cds 82900 83103 . + . NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local exon 82900 83103 . + . insertions 0 ; deletions 0 NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local similarity 75125 83103 906 + . alignment_id 1 ; Query SPP00000008_1.0 ; Align 75125 1 81 ; Align 75436 28 6 ; Align 78012 30 153 ; Align 81088 82 114 ; Align 81686 121 99 ; Align 82900 154 204 # --- END OF GFF DUMP --- # C4 Alignment: ------------ Query: SPP00000008_1.0 # Protein # Glutathione peroxidase 5 (GPx5) # Homo sapiens # Complete Target: NIVO01038789.1:subseq(406726,133103) Ammotragus lervia isolate NAZ scaffold232_len768847_cov70, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 869 Query range: 0 -> 221 Target range: 56729 -> 49999 1 : MetThrThrGlnLeuArgValValHisLeuLeuProLeuLeuLeuAlaCysPheValGln : 20 |||! ! !|||||| !!!.! ! !|||!!!||||||!!!|||!.! !!!!.!.!||| MetIleLeuGlnLeuTrpAlaSerCysLeuPheProLeuPheLeuValSerLeuAlaGln 56729 : ATGATCCTGCAGCTCTGGGCCTCCTGTCTTTTCCCACTGTTCCTGGTCAGCTTAGCTCAG : 56672 21 : ThrSerProLysGlnGluLysMetLys >>>> Target Intron 1 >>>> MetA : 31 !!!!|||||||||:!!:!!|||||| 2330 bp :!!| LeuThrProLysGlnGlnGlnMetLys++ ++ValA 56671 : CTGACTCCGAAGCAGCAACAGATGAAGgt.........................agGTGG : 54309 32 : spCysHisLysAspGluLysGlyThrIleTyrAspTyrGluAlaIleAlaLeuAsnLysA : 51 |||||:!!|||! ! !! !||||||||||||!!:|||! !|||:!!.!!|||||| . spCysTyrLysGlyValThrGlyThrIleTyrGluTyrGlyAlaLeuThrLeuAsnGlyG 54308 : ATTGCTATAAGGGGGTGACAGGCACCATCTATGAGTATGGAGCCCTCACCCTCAATGGTG : 54249 52 : snGluTyrValSerPheLysGlnTyrValGlyLysHisIleLeuPheValAsnValAlaT : 71 !.||||||:!!..!||||||||||||!.!|||||||||:!!||||||||||||||||||| luGluTyrIleGlnPheLysGlnTyrAlaGlyLysHisValLeuPheValAsnValAlaT 54248 : AGGAGTATATCCAGTTCAAGCAGTATGCAGGCAAGCATGTCCTGTTTGTCAATGTGGCCA : 54189 72 : hrTyrCysGlyLeuThrAlaGlnTyrPro{G} >>>> Target Intron 2 >>>> : 81 |||||!! |||||||||||||||||||||{|} 1830 bp hrTyr***GlyLeuThrAlaGlnTyrPro{G}++ + 54188 : CCTATTGAGGCTTGACAGCTCAGTATCCA{G}gt.........................a : 52330 82 : {lu}LeuAsnAlaLeuGlnGluGluLeuLysProTyrGlyLeuValValLeuGlyPheP : 100 {||}||||||||||||||||||||||||||||||!:!|||:!!|||||||||||||||| +{lu}LeuAsnAlaLeuGlnGluGluLeuLysProPheGlyValValValLeuGlyPheP 52329 : g{AA}CTGAATGCACTACAGGAGGAGCTGAAGCCTTTCGGTGTAGTTGTGTTGGGCTTTC : 52272 101 : roCysAsnGlnPheGlyLysGlnGluProGlyAspAsnLysGluIleLeuProGlyLeu{ : 120 |||||||||||||||||||||||||||||!.! ! |||..!||||||||| !||||||{ roCysAsnGlnPheGlyLysGlnGluProAlaLysAsnSerGluIleLeuMetGlyLeu{ 52271 : CCTGCAACCAATTTGGAAAACAAGAGCCAGCAAAGAACTCAGAGATCCTTATGGGGCTC{ : 52212 121 : Ly} >>>> Target Intron 3 >>>> {s}TyrValArgProGlyGlyGlyPheV : 129 ||} 1156 bp {|}||||||||||||||||||||||||| Ly}++ ++{s}TyrValArgProGlyGlyGlyPheV 52211 : AA}gt.........................ag{G}TATGTACGTCCAGGTGGTGGCTTTG : 51029 130 : alProSerPheGlnLeuPheGluLysGlyAspValAsnGlyGluLysGluGlnLysValP : 149 |||||!:!|||||||||||||||||||||||||||||||||||||||||||||||||||| alProAsnPheGlnLeuPheGluLysGlyAspValAsnGlyGluLysGluGlnLysValP 51028 : TTCCCAATTTTCAGCTCTTTGAGAAAGGGGATGTGaatggagaaaaagaacagaaggtCT : 50969 150 : heSerPheLeuLys >>>> Target Intron 4 >>>> HisSerCysProHisPr : 159 ||!!!||||||||| 751 bp :!!|||||||||! ! heThrPheLeuLys++ ++AsnSerCysProProTh 50968 : TTACCTTCCTGAAGgt.........................agAACTCCTGCCCTCCAAC : 50188 160 : oSerGluIleLeuGlyThrPheLysSerIleSerTrpAspProValLysValHisAspIl : 179 !|||!!::!!||||||:!!! !....!:!!! !|||!!:|||:!!|||||||||||||| rSerAspLeuLeuGlySerSerSerGlnLeuPheTrpGluProMetLysValHisAspIl 50187 : CTCTGATCTTTTGGGCTCATCAAGTCAGCTCTTCTGGGAGCCCATGAAGGTCCATGACAT : 50128 180 : eArgTrpAsnPheGluLysPheLeuValGlyProAspGlyIleProValMetArgTrpSe : 199 ||||||||||||||||||||||||||||||||||||||||:!!|||||||||||||||! eArgTrpAsnPheGluLysPheLeuValGlyProAspGlyValProValMetArgTrpTy 50127 : TCGCTGGAACTTTGAGAAGTTTCTAGTGGGGCCTGATGGTGTCCCTGTCATGCGTTGGTA : 50068 200 : rHisArgAlaThrValSerSerValLysThrAspIleLeuAlaTyrLeuLysGlnPheLy : 219 !|||||||||:!!||||||!!!||||||:!!|||!!:|||! !||||||||||||||||| rHisArgAlaSerValSerThrValLysSerAspMetLeuGluTyrLeuLysGlnPheLy 50067 : CCACCGGGCTTCAGTCAGCACAGTTAAATCAGACATGCTGGAGTACCTGAAGCAGTTCAA : 50008 220 : sThrLys : 221 |!:!:!! sSerGlu 50007 : AAGCGAA : 50000 vulgar: SPP00000008_1.0 0 221 . NIVO01038789.1:subseq(406726,133103) 56729 49999 - 869 M 29 87 5 0 2 I 0 2326 3 0 2 M 51 153 S 0 1 5 0 2 I 0 1826 3 0 2 S 1 2 M 38 114 S 0 2 5 0 2 I 0 1152 3 0 2 S 1 1 M 33 99 5 0 2 I 0 747 3 0 2 M 68 204 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2017-11-17 ##type DNA # # # seqname source feature start end score strand frame attributes # NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local gene 50000 56729 869 - . gene_id 1 ; sequence SPP00000008_1.0 ; gene_orientation + NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local cds 56643 56729 . - . NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local exon 56643 56729 . - . insertions 0 ; deletions 0 NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local splice5 56641 56642 . - . intron_id 1 ; splice_site "GT" NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local intron 54313 56642 . - . intron_id 1 NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local splice3 54313 54314 . - . intron_id 0 ; splice_site "ag" NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local cds 54159 54312 . - . NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local exon 54159 54312 . - . insertions 0 ; deletions 0 NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local splice5 54157 54158 . - . intron_id 2 ; splice_site "GT" NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local intron 52329 54158 . - . intron_id 2 NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local splice3 52329 52330 . - . intron_id 1 ; splice_site "AG" NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local cds 52211 52328 . - . NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local exon 52211 52328 . - . insertions 0 ; deletions 0 NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local splice5 52209 52210 . - . intron_id 3 ; splice_site "GT" NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local intron 51055 52210 . - . intron_id 3 NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local splice3 51055 51056 . - . intron_id 2 ; splice_site "AG" NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local cds 50955 51054 . - . NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local exon 50955 51054 . - . insertions 0 ; deletions 0 NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local splice5 50953 50954 . - . intron_id 4 ; splice_site "GT" NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local intron 50204 50954 . - . intron_id 4 NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local splice3 50204 50205 . - . intron_id 3 ; splice_site "AG" NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local cds 50000 50203 . - . NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local exon 50000 50203 . - . insertions 0 ; deletions 0 NIVO01038789.1:subseq(406726,133103) exonerate:protein2genome:local similarity 50000 56729 869 - . alignment_id 1 ; Query SPP00000008_1.0 ; Align 56730 1 87 ; Align 54313 30 153 ; Align 52327 82 114 ; Align 51054 121 99 ; Align 50204 154 204 # --- END OF GFF DUMP --- # -- completed exonerate analysis