Command line: [exonerate -m p2g --showtargetgff -q query_sacaromices.fa -t tricornatum1genome.fa] Hostname: [localhost.localdomain] C4 Alignment: ------------ Query: SPP00000077_1.0 # Protein # Selenoprotein R (SelR) # Saccharomyces cerevisiae # Complete Target: gi|217409007|gb|CM000609.1|:subseq(329661,2500) Phaeodactylum tricornutum CCAP 1055/1 chromosome 6, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 202 Query range: 34 -> 168 Target range: 1407 -> 999 35 : SerAspGluSerAsnAspValLysTrpAsnAspAlaLeuThrProLeuGlnLeuMetValLe : 55 !!!!!:.!...!!..!!: !! !|||..!! ! !||||||! ! !:!!!!! !||||| ThrGluHisAspThrGluSerMetTrpArgGlyLeuLeuThrArgGluGluPheArgValLe 1407 : ACCGAACACGACACGGAGTCTATGTGGCGCGGCCTATTGACACGGGAAGAATTCCGGGTCCT : 1347 56 : uArgAspLysAlaThrGluArgProAsnThrGlyAlaTyr<->LeuHisThrAsnGluSerG : 75 ||||... ! !.!||||||||||||..!:!! ! !||| ! ! !!!:!!!| uArgSerHisGlyThrGluArgProArgSerHisArgTyrAspThrTrpTyrProAspThrG 1346 : GCGATCGCACGGTACGGAACGACCCCGCTCCCACCGGTACGATACGTGGTACCCCGACACCG : 1284 76 : lyValTyrHisCysAlaAsnCysAspArgProLeuTyrSerSerLysAlaLysPheAspAla : 95 || !!:! !||| ! |||! !|||||||||:!!:!!::!||||||||||||:!! lyCysPheAlaCysArgAlaCysGlyLeuProLeuTyrAlaAlaArgAlaLysPheAspSer 1283 : GCTGCTTCGCCTGCCGTGCGTGCGGGTTGCCACTCTACGCGGCCCGGGCCAAATTTGATTCC : 1224 96 : ArgCysGlyTrpProAlaPheTyrGluGluValSerProGlyAlaIleThrTyrHisArgAs : 116 !!! |||||||||:!!||| !.!.|||..! !!.!!.! !||| ! ..!.! GlySerGlyTrpProSerPheGlyThrHisValGlnGlyAlaValAlaThrThrValGluGl 1223 : GGCTCGGGATGGCCCTCGTTTGGAACGCATGTACAAGGTGCCGTTGCCACGACCGTGGAACA : 1161 117 : pAsnSerLeuMetProAlaArgValGluIleCysCysAlaArgCysGlyGlyHisLeuGlyH : 137 .::! !!:!!||| ! !|||||||||||| !|||||||||||| !..!|||||||||| nSerProValMetGlyLysArgValGluIleHisCysAlaArgCysGlnSerHisLeuGlyH 1160 : GTCCCCGGTGATGGGGAAACGCGTCGAAATTCACTGCGCCCGCTGTCAGTCGCATTTGGGTC : 1098 138 : isValPheGluGlyGluGlyTrpLysGlnLeuLeuAsnLeuProLys<->AspThrArgHis : 156 ||||||||! ! ! ..! ! :!!! ! !...||| !! ... !|||||| isValPheAlaAspThrAsnThrAlaLysTrpAspArgLeuLysThrPheSerGluArgHis 1097 : ACGTGTTTGCCGACACCAATACGGCCAAATGGGATCGACTCAAAACCTTCTCGGAGCGTCAC : 1038 157 : UnkValAsnSerAlaSerLeuAsnLeuLysLysAsp : 168 ||||||.!! !||||||! ! !!..|||:!! CysValAsnGlyIleSerLeuMetTyrSerLysAsn 1037 : TGTGTGAACGGTATTAGTCTCATGTACAGTAAGAAC : 1000 vulgar: SPP00000077_1.0 34 168 . gi|217409007|gb|CM000609.1|:subseq(329661,2500) 1407 999 - 202 M 34 102 G 0 3 M 84 252 G 0 3 M 16 48 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2011-03-08 ##type DNA # # # seqname source feature start end score strand frame attributes # gi|217409007|gb|CM000609.1|:subseq(329661,2500) exonerate:protein2genome:local gene 1000 1407 202 - . gene_id 1 ; sequence SPP00000077_1.0 ; gene_orientation . gi|217409007|gb|CM000609.1|:subseq(329661,2500) exonerate:protein2genome:local cds 1000 1407 . - . gi|217409007|gb|CM000609.1|:subseq(329661,2500) exonerate:protein2genome:local exon 1000 1407 . - . insertions 6 ; deletions 0 gi|217409007|gb|CM000609.1|:subseq(329661,2500) exonerate:protein2genome:local similarity 1000 1407 202 - . alignment_id 1 ; Query SPP00000077_1.0 ; Align 1408 35 102 ; Align 1303 69 252 ; Align 1048 153 48 # --- END OF GFF DUMP --- # -- completed exonerate analysis