Command line: [exonerate -m p2g --showtargetgff -q secs_zf.aa.fa.canviat -t genomesecs_JXSJ01000055.1.fa] Hostname: [ws101071U] C4 Alignment: ------------ Query: SPP00000631_2.0 # Protein # Selenocysteine synthase (SecS) # Zebrafish Target: JXSJ01000055.1:subseq(1599664,6291) Miichthys miiuy scaffold55, whole genome shotgun sequence Model: protein2genome:local Raw score: 569 Query range: 0 -> 135 Target range: 2999 -> 3791 1 : MetAsnSerGluAsnPheArgLeuSerGluLysLeuValSerAlaSerTyrValArgGlnGl : 21 ||||||.!!|||||||||!!!||||||||||||:!!|||||| !!:!!|||:!!|||||||| MetAsnGlyGluAsnPheSerLeuSerGluLysIleValSerProAlaTyrIleArgGlnGl 3000 : aTGAACGGTGAAAACTTCAGCCTGAGTGAGAAGATCGTGTCACCCGCCTACATTCGTCAGGG : 3060 22 : yLeuAspAlaArgArgGlyHisGluGlnLeuIleArgArgLeuLeuAspGln >>>> Tar : 39 | !.!.|||||||||.!!||||||||||||||||||!.!||||||!!:||| ySerGlnAlaArgArgSerHisGluGlnLeuIleArgHisLeuLeuGluGln++ 3061 : GAGCCAGGCGCGTCGTAGCCATGAGCAGCTCATCAGACACTTACTGGAGCAGgt........ : 3116 40 : get Intron 1 >>>> GlyLysCysProGluGluGlyTrpSerGluSerThrIleGluL : 53 80 bp ||||||||||||||||||||||||||||||||||||||||||| ++GlyLysCysProGluGluGlyTrpSerGluSerThrIleGluL 3117 : .................agGGTAAGTGTCCGGAGGAGGGATGGAGCGAGAGCACCATAGAGC : 3236 54 : euPheLeuAsnGluLeuAlaValMetAspSerAsnAsnPheLeuGlyAsnCysGlyValGly : 73 ||||||||||||||||||||||||||||||||||||||||||||.!!||||||||||||||| euPheLeuAsnGluLeuAlaValMetAspSerAsnAsnPheLeuSerAsnCysGlyValGly 3237 : TCTTCTTGAATGAGCTGGCTGTGATGGACAGCAATAATTTCCTCAGCAACTGCGGcgtggga : 3296 74 : GluArgGluGlyArgValAlaSerSerLeuValAlaArgArgHisTyr{Ar} >>>> Tar : 90 ||||||||||||||||||||||||||||||||||||||||||||||||{||} GluArgGluGlyArgValAlaSerSerLeuValAlaArgArgHisTyr{Ar}++ 3297 : gagagggaaggcCGGGTGGCATCCAGCCTTGTGGCAAGACGACATTAC{AG}gt........ : 3351 91 : get Intron 2 >>>> {g}LeuIleHisGlyIleGlyArgSerGlyAspIleAlaAlaV : 104 307 bp {|}|||||||||||||||||||||||||||||||||||||||| ++{g}LeuIleHisGlyIleGlyArgSerGlyAspIleAlaAlaV 3352 : .................ag{G}TTGATCCATGGCATAGGTCGGTCAGGTGACATCGCTGCTG : 3696 105 : alGlnProLysAlaAlaGlySerSerLeuLeuAsnLysIleThrAsnSerValValLeuAsp : 124 ||||||||||||||||||||||||||||||||||||||:!:||||||||||||||||||||| alGlnProLysAlaAlaGlySerSerLeuLeuAsnLysLeuThrAsnSerValValLeuAsp 3697 : TTCAGCCCAAAGCAGCAGGATCCAGTCTGCTTAACAAGCTGACCAATTCAGTTGTGTTGGAC : 3756 125 : ValLeuLysLeuThrGlyValArgSerValSer : 135 :!!|||||||||:!!|||:!! !! !:!!||| IleLeuLysLeuSerGlyMetTyrMetLeuSer 3757 : ATCTTAAAGCTCTCAGGTATGTATATGTTATCA : 3791 vulgar: SPP00000631_2.0 0 135 . JXSJ01000055.1:subseq(1599664,6291) 2999 3791 + 569 M 38 114 5 0 2 I 0 76 3 0 2 M 51 153 S 0 2 5 0 2 I 0 303 3 0 2 S 1 1 M 45 135 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2016-11-28 ##type DNA # # # seqname source feature start end score strand frame attributes # JXSJ01000055.1:subseq(1599664,6291) exonerate:protein2genome:local gene 3000 3791 569 + . gene_id 1 ; sequence SPP00000631_2.0 ; gene_orientation + JXSJ01000055.1:subseq(1599664,6291) exonerate:protein2genome:local cds 3000 3113 . + . JXSJ01000055.1:subseq(1599664,6291) exonerate:protein2genome:local exon 3000 3113 . + . insertions 0 ; deletions 0 JXSJ01000055.1:subseq(1599664,6291) exonerate:protein2genome:local splice5 3114 3115 . + . intron_id 1 ; splice_site "GT" JXSJ01000055.1:subseq(1599664,6291) exonerate:protein2genome:local intron 3114 3193 . + . intron_id 1 JXSJ01000055.1:subseq(1599664,6291) exonerate:protein2genome:local splice3 3192 3193 . + . intron_id 0 ; splice_site "AG" JXSJ01000055.1:subseq(1599664,6291) exonerate:protein2genome:local cds 3194 3348 . + . JXSJ01000055.1:subseq(1599664,6291) exonerate:protein2genome:local exon 3194 3348 . + . insertions 0 ; deletions 0 JXSJ01000055.1:subseq(1599664,6291) exonerate:protein2genome:local splice5 3349 3350 . + . intron_id 2 ; splice_site "GT" JXSJ01000055.1:subseq(1599664,6291) exonerate:protein2genome:local intron 3349 3655 . + . intron_id 2 JXSJ01000055.1:subseq(1599664,6291) exonerate:protein2genome:local splice3 3654 3655 . + . intron_id 1 ; splice_site "AG" JXSJ01000055.1:subseq(1599664,6291) exonerate:protein2genome:local cds 3656 3791 . + . JXSJ01000055.1:subseq(1599664,6291) exonerate:protein2genome:local exon 3656 3791 . + . insertions 0 ; deletions 0 JXSJ01000055.1:subseq(1599664,6291) exonerate:protein2genome:local similarity 3000 3791 569 + . alignment_id 1 ; Query SPP00000631_2.0 ; Align 3000 1 114 ; Align 3194 39 153 ; Align 3657 91 135 # --- END OF GFF DUMP --- # -- completed exonerate analysis