Command line: [exonerate -m p2g --showtargetgff -q ../SPS2-human.fa -t SPS2-NIVO01016939.1.subseq] Hostname: [sitdoc] C4 Alignment: ------------ Query: SPP00000033_1.0 # Protein # Selenophosphate synthetase 2 (SPS2) # Homo sapiens # Complete Target: NIVO01016939.1:subseq(853875,54367) Ammotragus lervia isolate NAZ scaffold1305_len909451_cov72, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 2079 Query range: 0 -> 447 Target range: 51319 -> 49972 1 : MetAlaGluAlaSerAlaThrGlyAlaCysGlyGluAlaMet<-><-><->AlaAlaAla : 17 |||||||||!.!:!!|||||||||||| !!|||||||||||| ||||||!.! MetAlaGluValAlaAlaThrGlyAlaGlyGlyGluAlaMetAlaAlaValAlaAlaGly 51319 : ATGGCGGAAGTCGCGGCGACTGGCGCCGGCGGAGAGGCTATGGCGGCAGTGGCGGCGGGG : 51262 18 : GluGlySerSerGlyProAlaGlyLeuThrLeuGlyArgSerPheSerAsnTyrArgPro : 37 ||||||||||||||||||!.!|||||| !!|||||||||.!!||||||!:!||||||||| GluGlySerSerGlyProValGlyLeuProLeuGlyArgGlyPheSerSerTyrArgPro 51261 : GAAGGCTCCTCGGGCCCGGTGGGCTTGCCTCTAGGCCGGGGCTTCTCGAGCTACCGGCCC : 51202 38 : PheGluProGlnAlaLeuGlyLeuSerProSerTrpArgLeuThrGlyPheSerGlyMet : 57 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| PheGluProGlnAlaLeuGlyLeuSerProSerTrpArgLeuThrGlyPheSerGlyMet 51201 : TTCGAGCCCCAGGCGTTGGGTCTCAGCCCGAGCTGGCGGCTCACCGGCTTCTCTGGCATG : 51142 58 : LysGlyUnkGlyCysLysValProGlnGluAlaLeuLeuLysLeuLeuAlaGlyLeuThr : 77 |||||| |||||||||||||||||||||.!!||||||||||||||||||||||||||| LysGly***GlyCysLysValProGlnGluThrLeuLeuLysLeuLeuAlaGlyLeuThr 51141 : AAGGGCTGAGGCTGCAAGGTCCCCCAGGAGACGCTGCTCAAACTCCTGGCGGGACTGACG : 51082 78 : ArgProAspValArgProProLeuGlyArgGlyLeuValGlyGlyGlnGluGluAlaSer : 97 ||||||!!:|||! !||| !!! !|||||||||||||||||||||! ||||||.!!:!! ArgProGluValLeuProAlaProGlyArgGlyLeuValGlyGlyLeuGluGluThrAla 51081 : CGACCGGAGGTGCTGCCCGCGCCGGGCCGGGGCCTGGTCGGAGGCCTTGAAGAGACGGCC : 51022 98 : GlnGluAlaGlyLeuProAlaGlyAlaGlyProSerProThrPheProAlaLeuGlyIle : 117 ||||||||||||||||||!.! !!|||! !||||||||| ||||||:::||||||||| GlnGluAlaGlyLeuProValArgAlaGluProSerProProPheProSerLeuGlyIle 51021 : CAGGAAGCCGGTCTGCCGGTGAGGGCGGAACCCAGcccccccttccccagcctgggCATT : 50962 118 : GlyMetAspSerCysValIleProLeuArgHisGlyGlyLeuSerLeuValGlnThrThr : 137 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| GlyMetAspSerCysValIleProLeuArgHisGlyGlyLeuSerLeuValGlnThrThr 50961 : GGGATGGACTCCTGCGTCATACCCCTGAGACACGGGGGCCTGTCCTTGGTGCAGACCACC : 50902 138 : AspPhePheTyrProLeuValGluAspProTyrMetMetGlyArgIleAlaCysAlaAsn : 157 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| AspPhePheTyrProLeuValGluAspProTyrMetMetGlyArgIleAlaCysAlaAsn 50901 : GACTTCTTTTACCCCCTGGTGGAAGATCCCTATATGATGGGGCGCATCGCATGTGCCAAC : 50842 158 : ValLeuSerAspLeuTyrAlaMetGlyIleThrGluCysAspAsnMetLeuMetLeuLeu : 177 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| ValLeuSerAspLeuTyrAlaMetGlyIleThrGluCysAspAsnMetLeuMetLeuLeu 50841 : GTGCTGAGTGACCTCTACGCGATGGGCATTACTGAGTGTGATAACATGTTGATGCTACTC : 50782 178 : SerValSerGlnSerMetSerGluGluGluArgGluLysValThrProLeuMetValLys : 197 |||||||||||||||||| !!||||||||||||||||||:!!||||||||||||:!!||| SerValSerGlnSerMetProGluGluGluArgGluLysIleThrProLeuMetIleLys 50781 : AGCGTCAGCCAGAGTATGCCTGAGGAGGAGCGAGAAAAGATAACACCACTCATGATCAAA : 50722 198 : GlyPheArgAspAlaAlaGluGluGlyGlyThrAlaValThrGlyGlyGlnThrValVal : 217 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| GlyPheArgAspAlaAlaGluGluGlyGlyThrAlaValThrGlyGlyGlnThrValVal 50721 : GGCTTTCGGGATGCTGCCGAGGAAGGAGGGACTGCAGTGACTGGTGGGCAAACGGTGGTA : 50662 218 : AsnProTrpIleIleIleGlyGlyValAlaThrValValCysGlnProAsnGluPheIle : 237 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| AsnProTrpIleIleIleGlyGlyValAlaThrValValCysGlnProAsnGluPheIle 50661 : AACCCTTGGATTATCATCGGTGGAGTTGCCACTGTAGTATGTCAACCAAATGAATTCATC : 50602 238 : MetProAspSerAlaValValGlyAspValLeuValLeuThrLysProLeuGlyThrGln : 257 :!!||||||||||||||||||||||||||||||||||||||||||||||||||||||||| LeuProAspSerAlaValValGlyAspValLeuValLeuThrLysProLeuGlyThrGln 50601 : TTGCCTGACAGCGCCGTTGTTGGGGATGTGCTCGTGTTAACTAAACCCTTAGGAACCCAG : 50542 258 : ValAlaValAsnAlaHisGlnTrpLeuAspAsnProGluArgTrpAsnLysValLysMet : 277 |||||||||||||||||||||||||||||||||||||||||||||||||||:!!|||||| ValAlaValAsnAlaHisGlnTrpLeuAspAsnProGluArgTrpAsnLysIleLysMet 50541 : GTAGCTGTCAATGCCCACCAATGGCTGGATAATCCTGAACGATGGAATAAGATAAAGATG : 50482 278 : ValValSerArgGluGluValGluLeuAlaTyrGlnGluAlaMetPheAsnMetAlaThr : 297 |||||||||||||||||||||||||||!.!|||||||||||||||||||||||||||||| ValValSerArgGluGluValGluLeuGlyTyrGlnGluAlaMetPheAsnMetAlaThr 50481 : GTGGTCTCCAGAGAAGAGGTAGAGCTGGGTTATCAGGAAGCCATGTTCAACATGGCTACC : 50422 298 : LeuAsnArgThrAlaAlaGlyLeuMetHisThrPheAsnAlaHisAlaAlaThrAspIle : 317 |||||||||! !|||||||||||||||||||||||||||||||||||||||||||||||| LeuAsnArgIleAlaAlaGlyLeuMetHisThrPheAsnAlaHisAlaAlaThrAspIle 50421 : CTCAATAGAATTGCTGCCGGGCTGATGCACACGTTTAATGCCCATGCAGCCACAGATATC : 50362 318 : ThrGlyPheGlyIleLeuGlyHisSerGlnAsnLeuAlaLysGlnGlnArgAsnGluVal : 337 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| ThrGlyPheGlyIleLeuGlyHisSerGlnAsnLeuAlaLysGlnGlnArgAsnGluVal 50361 : ACCGGCTTTGGCATTCTGGGTCACTCTCAGAACCTTgcaaaacaacaaagaaatgagGTC : 50302 338 : SerPheValIleHisAsnLeuProIleIleAlaLysMetAlaAlaValSerLysAlaSer : 357 |||||||||||||||||||||||||||||||||||||||||||||:!!|||||||||||| SerPheValIleHisAsnLeuProIleIleAlaLysMetAlaAlaIleSerLysAlaSer 50301 : TCGTTTGTCATTCATAATCTGCCCATCATTGCCAAGATGGCTGCCATCAGCAAGGCCAGC : 50242 358 : GlyArgPheGlyLeuLeuGlnGlyThrSerAlaGluThrSerGlyGlyLeuLeuIleCys : 377 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| GlyArgPheGlyLeuLeuGlnGlyThrSerAlaGluThrSerGlyGlyLeuLeuIleCys 50241 : GGGCGCTTTGGTCTCCTGCAAGGAACCTCAGCAGAAACTTCTGGGGGGCTACTGATTTGT : 50182 378 : LeuProArgGluGlnAlaAlaArgPheCysSerGluIleLysSerSerLysTyrGlyGlu : 397 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| LeuProArgGluGlnAlaAlaArgPheCysSerGluIleLysSerSerLysTyrGlyGlu 50181 : CTGCCAAGAGAACAAGCAGCTCGCTTTTGTTCAGAAATCAAATCTTCAAAGTATGGAGAG : 50122 398 : GlyHisGlnAlaTrpIleValGlyIleValGluLysGlyAsnArgThrAlaArgIleIle : 417 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| GlyHisGlnAlaTrpIleValGlyIleValGluLysGlyAsnArgThrAlaArgIleIle 50121 : GGTCATCAAGCCTGGATCGTCGGTATTGTGGAAAAAGGAAACCGGACAGCCCGGATCATT : 50062 418 : AspLysProArgValIleGluValLeuProArgGlyAlaThrAlaAlaValLeuAlaPro : 437 |||||||||||||||||||||||||||||||||||||||||| |||..!||||||||| AspLysProArgValIleGluValLeuProArgGlyAlaThr---AlaThrLeuAlaPro 50061 : GACAAGCCTCGCGTTATTGAGGTCCTACCTCGTGGGGCCACT---GCTACTCTTGCTCCT : 50005 438 : AspSerSerAsnAlaSerSerGluProSer : 447 !!:!:!|||!:!|||||||||||||||||| GluAsnSerSerAlaSerSerGluProSer 50004 : GAGAACTCCAGTGCCTCCTCTGAGCCTAGC : 49973 vulgar: SPP00000033_1.0 0 447 . NIVO01016939.1:subseq(853875,54367) 51319 49972 - 2079 M 14 42 G 0 9 M 417 1251 G 1 0 M 15 45 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2017-11-21 ##type DNA # # # seqname source feature start end score strand frame attributes # NIVO01016939.1:subseq(853875,54367) exonerate:protein2genome:local gene 49973 51319 2079 - . gene_id 1 ; sequence SPP00000033_1.0 ; gene_orientation . NIVO01016939.1:subseq(853875,54367) exonerate:protein2genome:local cds 49973 51319 . - . NIVO01016939.1:subseq(853875,54367) exonerate:protein2genome:local exon 49973 51319 . - . insertions 9 ; deletions 1 NIVO01016939.1:subseq(853875,54367) exonerate:protein2genome:local similarity 49973 51319 2079 - . alignment_id 1 ; Query SPP00000033_1.0 ; Align 51320 1 42 ; Align 51269 15 1251 ; Align 50018 433 45 # --- END OF GFF DUMP --- # -- completed exonerate analysis