Command line: [exonerate -m p2g --showtargetgff -q Sel.human/SPS1.fa -t gi|397458686|gb|JH736815.1|subseq.fa -E] Hostname: [ws101064] C4 Alignment: ------------ Query: SPP00000060_2.0 # Protein # Selenophosphate synthetase 1 (SPS1) # Human (Homo sapiens) Target: gi|397458686|gb|JH736815.1|:subseq(0,98669) Alligator mississippiensis unplaced genomic scaffold scaffold-9464, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 1635 Query range: 0 -> 321 Target range: 48849 -> 21200 1 : MetSerThrArgGluSerPheAsnProGluSerTyrGluLeuAspLysSerPheArgLeu : 20 ||||||..!||||||||||||||||||||||||||||||||||||||||||||||||||| MetSerValArgGluSerPheAsnProGluSerTyrGluLeuAspLysSerPheArgLeu 48849 : ATGTCTGTTCGGGAGTCCTTTAACCCAGAGAGCTATGAACTGGACAAGAGCTTCCGCCTG : 48792 21 : ThrArgPheThrGluLeuLysGlyThrGlyCysLysValProGlnAspValLeuGlnLys : 40 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| ThrArgPheThrGluLeuLysGlyThrGlyCysLysValProGlnAspValLeuGlnLys 48791 : ACCCGATTCACTGAACTGAAAGGCACAGGCTGCAAAGTGCCTCAGGATGTTTTACAGAAA : 48732 41 : LeuLeuGluSerLeuGlnGluAsnHisPheGlnGluAspGluGlnPheLeuGlyAlaVal : 60 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| LeuLeuGluSerLeuGlnGluAsnHisPheGlnGluAspGluGlnPheLeuGlyAlaVal 48731 : CTACTTGAATCCTTACAAGAAAACCATTTCCAAGAAGATGAACAGTTTCTTGGAGCTGTT : 48672 61 : MetProArgLeu{G} >>>> Target Intron 1 >>>> {ly}IleGlyMetAsp : 69 ||||||||||||{|} 2508 bp {||}|||||||||||| MetProArgLeu{G}++ ++{ly}IleGlyMetAsp 48671 : ATGCCCCGGTTG{G}gt.........................ag{GA}ATTGGGATGGAC : 46137 70 : ThrCysValIleProLeuArgHisGlyGlyLeuSerLeuValGlnThrThrAspTyrIle : 89 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| ThrCysValIleProLeuArgHisGlyGlyLeuSerLeuValGlnThrThrAspTyrIle 46136 : ACATGCGTTATTCCTTTAAGACATGGAGGTCTGTCTTTGGTCCAGACAACAGATTACATT : 46077 90 : TyrProIleValAspAspProTyrMetMet >>>> Target Intron 2 >>>> G : 100 |||||||||||||||||||||||||||||| 8456 bp | TyrProIleValAspAspProTyrMetMet++ ++G 46076 : TATCCAATTGTGGATGACCCATACATGATGgt.........................agG : 37588 101 : lyArgIleAlaCysAlaAsnValLeuSerAspLeuTyrAlaMetGlyValThrGluCysA : 120 ||||| !!|||||||||||||||||||||||||||||||||||||||||||||||||||| lyArgUnkAlaCysAlaAsnValLeuSerAspLeuTyrAlaMetGlyValThrGluCysA 37587 : GACGTNTAGCCTGTGCCAATGTCCTAAGTGACCTTTATGCCATGGGGGTCACAGAGTGCG : 37528 121 : spAsnMetLeuMetLeuLeuGlyValSerAsnLysMetThrAspArg >>>> Target : 136 |||||||||||||||||||||||:!!||||||||||||||||||||| 9 spAsnMetLeuMetLeuLeuGlyIleSerAsnLysMetThrAspArg++ 37527 : ACAACATGTTGATGCTCCTTGGAATCAGTAATAAAATGACGGACAGGgt........... : 37478 137 : Intron 3 >>>> GluArgAspLysValMetProLeuIleIleGlnGlyPheLysAs : 150 35 bp |||||||||||||||||||||||||||||||||||||||||||| ++GluArgAspLysValMetProLeuIleIleGlnGlyPheLysAs 37477 : ..............agGAAAGGGACAAAGTGATGCCTCTAATTATCCAGGGCTTCAAAGA : 36503 151 : pAlaAlaGluGluAlaGlyThrSerValThrGlyGlyGlnThrValLeuAsnProTrpIl : 170 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| pAlaAlaGluGluAlaGlyThrSerValThrGlyGlyGlnThrValLeuAsnProTrpIl 36502 : TGCAGCAGAAGAGGCAGGGACTTCTGTGACGGGTGGCCAAACGGTGTTAAATCCTTGGAT : 36443 171 : eValLeuGlyGlyValAlaThrThrValCysGlnProAsnGluPheIle{Me} >>>> : 187 |||||||||||||||||||||||||||||||||||||||||||||||||{||} eValLeuGlyGlyValAlaThrThrValCysGlnProAsnGluPheIle{Me}++ 36442 : TGTTCTTGGAGGAGTGGCCACAACAGTCTGTCAGCCTAATGAATTCATA{AT}gt..... : 36388 188 : Target Intron 4 >>>> {t}ProAspAsnAlaValProGlyAspValLeuValLe : 199 817 bp {|}||||||||||||||||||||||||||||||||||| ++{t}ProAspAsnAlaValProGlyAspValLeuValLe 36387 : ....................ag{G}CCAGACAATGCAGTGCCAGGAGATGTACTTGTATT : 35539 200 : uThrLysProLeuGlyThrGlnValAlaValAlaValHisGlnTrpLeuAspIle >>> : 218 ||||||||||||||||||||||||||||||||||||||||||||||||||||||| uThrLysProLeuGlyThrGlnValAlaValAlaValHisGlnTrpLeuAspIle++ 35538 : AACTAAACCCTTGGGAACACAAGTGGCTGTAGCTGTCCACCAGTGGCTAGATATTgt... : 35480 219 : > Target Intron 5 >>>> ProGluLysTrpAsnLysIleLysLeuValValThr : 229 3823 bp |||||||||||||||||||||||||||||||||||| ++ProGluLysTrpAsnLysIleLysLeuValValThr 35479 : ......................agCCAGAAAAATGGAATAAAATCAAACTGGTGGTAACG : 31626 230 : GlnGluAspValGluLeuAlaTyrGlnGluAlaMetMetAsnMetAlaArgLeuAsnArg : 249 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| GlnGluAspValGluLeuAlaTyrGlnGluAlaMetMetAsnMetAlaArgLeuAsnArg 31625 : CAAGAAGATGTAGAGCTAGCATACCAGGAGGCAATGATGAATATGGCAAGGCTCAACAGA : 31566 250 : Thr{G} >>>> Target Intron 6 >>>> {ly}GlyLeuLeuIleCysLeuPro : 258 |||{|} 10147 bp {||}||||||||||||||||||||| Thr{G}++ ++{ly}GlyLeuLeuIleCysLeuPro 31565 : ACA{G}gt.........................ag{GA}GGCCTCCTCATCTGTTTGCCA : 21392 259 : ArgGluGlnAlaAlaArgPheCysAlaGluIleLysSerProLysTyrGlyGluGlyHis : 278 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| ArgGluGlnAlaAlaArgPheCysAlaGluIleLysSerProLysTyrGlyGluGlyHis 21391 : CGTGAACAGGCTGCACGTTTCTGTGCTGAGATCAAATCTCCCAAGTACGGCGAGGGTCAC : 21332 279 : GlnAlaTrpIleIleGlyIleValGluLysGlyAsnArgThrAlaArgIleIleAspLys : 298 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| GlnAlaTrpIleIleGlyIleValGluLysGlyAsnArgThrAlaArgIleIleAspLys 21331 : CAAGCATGGATCATTGGGATCGTAGAGAAGGGGAACCGCACTGCCAGAATTATAGATAAA : 21272 299 : ProArgIleIleGluValAlaProGlnValAlaThrGlnAsnValAsnProThrProGly : 318 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| ProArgIleIleGluValAlaProGlnValAlaThrGlnAsnValAsnProThrProGly 21271 : CCACGAATCATTGAAGTTGCACCGCAAGTGGCCACTCAAAACGTGAACCCAACACCTGGT : 21212 319 : AlaThrSer : 321 ||||||||| AlaThrSer 21211 : GCCACTTCT : 21201 vulgar: SPP00000060_2.0 0 321 . gi|397458686|gb|JH736815.1|:subseq(0,98669) 48849 21200 - 1635 M 64 192 S 0 1 5 0 2 I 0 2504 3 0 2 S 1 2 M 34 102 5 0 2 I 0 8452 3 0 2 M 36 108 5 0 2 I 0 931 3 0 2 M 51 153 S 0 2 5 0 2 I 0 813 3 0 2 S 1 1 M 30 90 5 0 2 I 0 3819 3 0 2 M 33 99 S 0 1 5 0 2 I 0 10143 3 0 2 S 1 2 M 70 210 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2013-11-25 ##type DNA # # # seqname source feature start end score strand frame attributes # gi|397458686|gb|JH736815.1|:subseq(0,98669) exonerate:protein2genome:local gene 21201 48849 1635 - . gene_id 1 ; sequence SPP00000060_2.0 ; gene_orientation + gi|397458686|gb|JH736815.1|:subseq(0,98669) exonerate:protein2genome:local cds 48657 48849 . - . gi|397458686|gb|JH736815.1|:subseq(0,98669) exonerate:protein2genome:local exon 48657 48849 . - . insertions 0 ; deletions 0 gi|397458686|gb|JH736815.1|:subseq(0,98669) exonerate:protein2genome:local splice5 48655 48656 . - . intron_id 1 ; splice_site "GT" gi|397458686|gb|JH736815.1|:subseq(0,98669) exonerate:protein2genome:local intron 46149 48656 . - . intron_id 1 gi|397458686|gb|JH736815.1|:subseq(0,98669) exonerate:protein2genome:local splice3 46149 46150 . - . intron_id 0 ; splice_site "AG" gi|397458686|gb|JH736815.1|:subseq(0,98669) exonerate:protein2genome:local cds 46045 46148 . - . gi|397458686|gb|JH736815.1|:subseq(0,98669) exonerate:protein2genome:local exon 46045 46148 . - . insertions 0 ; deletions 0 gi|397458686|gb|JH736815.1|:subseq(0,98669) exonerate:protein2genome:local splice5 46043 46044 . - . intron_id 2 ; splice_site "GT" gi|397458686|gb|JH736815.1|:subseq(0,98669) exonerate:protein2genome:local intron 37589 46044 . - . intron_id 2 gi|397458686|gb|JH736815.1|:subseq(0,98669) exonerate:protein2genome:local splice3 37589 37590 . - . intron_id 1 ; splice_site "AG" gi|397458686|gb|JH736815.1|:subseq(0,98669) exonerate:protein2genome:local cds 37481 37588 . - . gi|397458686|gb|JH736815.1|:subseq(0,98669) exonerate:protein2genome:local exon 37481 37588 . - . insertions 0 ; deletions 0 gi|397458686|gb|JH736815.1|:subseq(0,98669) exonerate:protein2genome:local splice5 37479 37480 . - . intron_id 3 ; splice_site "GT" gi|397458686|gb|JH736815.1|:subseq(0,98669) exonerate:protein2genome:local intron 36546 37480 . - . intron_id 3 gi|397458686|gb|JH736815.1|:subseq(0,98669) exonerate:protein2genome:local splice3 36546 36547 . - . intron_id 2 ; splice_site "AG" gi|397458686|gb|JH736815.1|:subseq(0,98669) exonerate:protein2genome:local cds 36391 36545 . - . gi|397458686|gb|JH736815.1|:subseq(0,98669) exonerate:protein2genome:local exon 36391 36545 . - . insertions 0 ; deletions 0 gi|397458686|gb|JH736815.1|:subseq(0,98669) exonerate:protein2genome:local splice5 36389 36390 . - . intron_id 4 ; splice_site "GT" gi|397458686|gb|JH736815.1|:subseq(0,98669) exonerate:protein2genome:local intron 35574 36390 . - . intron_id 4 gi|397458686|gb|JH736815.1|:subseq(0,98669) exonerate:protein2genome:local splice3 35574 35575 . - . intron_id 3 ; splice_site "AG" gi|397458686|gb|JH736815.1|:subseq(0,98669) exonerate:protein2genome:local cds 35483 35573 . - . gi|397458686|gb|JH736815.1|:subseq(0,98669) exonerate:protein2genome:local exon 35483 35573 . - . insertions 0 ; deletions 0 gi|397458686|gb|JH736815.1|:subseq(0,98669) exonerate:protein2genome:local splice5 35481 35482 . - . intron_id 5 ; splice_site "GT" gi|397458686|gb|JH736815.1|:subseq(0,98669) exonerate:protein2genome:local intron 31660 35482 . - . intron_id 5 gi|397458686|gb|JH736815.1|:subseq(0,98669) exonerate:protein2genome:local splice3 31660 31661 . - . intron_id 4 ; splice_site "AG" gi|397458686|gb|JH736815.1|:subseq(0,98669) exonerate:protein2genome:local cds 31560 31659 . - . gi|397458686|gb|JH736815.1|:subseq(0,98669) exonerate:protein2genome:local exon 31560 31659 . - . insertions 0 ; deletions 0 gi|397458686|gb|JH736815.1|:subseq(0,98669) exonerate:protein2genome:local splice5 31558 31559 . - . intron_id 6 ; splice_site "GT" gi|397458686|gb|JH736815.1|:subseq(0,98669) exonerate:protein2genome:local intron 21413 31559 . - . intron_id 6 gi|397458686|gb|JH736815.1|:subseq(0,98669) exonerate:protein2genome:local splice3 21413 21414 . - . intron_id 5 ; splice_site "AG" gi|397458686|gb|JH736815.1|:subseq(0,98669) exonerate:protein2genome:local cds 21201 21412 . - . gi|397458686|gb|JH736815.1|:subseq(0,98669) exonerate:protein2genome:local exon 21201 21412 . - . insertions 0 ; deletions 0 gi|397458686|gb|JH736815.1|:subseq(0,98669) exonerate:protein2genome:local similarity 21201 48849 1635 - . alignment_id 1 ; Query SPP00000060_2.0 ; Align 48850 1 192 ; Align 46147 66 102 ; Align 37589 100 108 ; Align 36546 136 153 ; Align 35573 188 90 ; Align 31660 218 99 ; Align 21411 252 210 # --- END OF GFF DUMP --- # -- completed exonerate analysis