Command line: [exonerate -m p2g --showtargetgff -q ./proteins/SPP00000033_1.0.fa -t ./exonerate/SPP00000033_1.0_3_genomic.fa_exonerate.fa]
Hostname: [luke]

C4 Alignment:
------------
         Query: SPP00000033_1.0 # Protein # Selenophosphate synthetase 2 (SPS2) # Homo sapiens # Complete
        Target: gi|357434203|gb|JH378160.1|:subseq(5751640,301147) Saimiri boliviensis boliviensis unplaced genomic scaffold scaffold00056, whole genome shotgun sequence
         Model: protein2genome:local
     Raw score: 1161
   Query range: 171 -> 448
  Target range: 150342 -> 151177

    172 : AsnMetLeuMetLeuLeuSerValSerGlnSerMetSerGluGluGluArgGluLysV :    191
          ..!||||||||||||||||||||||||||||||:!!||||||||||||||||||||||
          GlyMetLeuMetLeuLeuSerValSerGlnSerValSerGluGluGluArgGluLysV
 150343 : GGCATGTTGATGTTACTCAGCGTCAGCCAGAGTGTGAGTGAGGAGGAACGCGAAAAGG : 150400

    192 : alThrProLeuMetValLysGlyPheArgAspAlaAlaGluGluGlyGlyThrAlaVa :    210
          ||||||||||||||||||||||||||||||||!.!||||||||||||||| !!|||||
          alThrProLeuMetValLysGlyPheArgAspValAlaGluGluGlyGlyProAlaVa
 150401 : TAACGCCGCTCATGGTCAAAGGCTTTCGGGATGTCGCTGAGGAAGGAGGGCCCGCAGT : 150457

    211 : lThrGlyGlyGlnThrValValAsnProTrpIleIleIleGlyGlyValAlaThrVal :    229
          |||||||||||||||||||! !:!!|||||||||! !|||||||||!.!|||! !|||
          lThrGlyGlyGlnThrValGlyHisProTrpIleArgIleGlyGlyAlaAlaIleVal
 150458 : GACCGGCGGGCAAACCGTGGGCCACCCTTGGATTAGAATCGGCGGAGCTGCCATTGTG : 150514

    230 : ValCysGlnProAsnGluPheIleMetProAspSerAlaValValGlyAspValLeuV :    249
          ||| !!|||||||||||||||||||||||| !!|||||||||||||||! !!.!|||!
          ValGlyGlnProAsnGluPheIleMetProHisSerAlaValValGlyGlyAlaLeuA
 150515 : GTGGGCCAGCCAAATGAATTCATAATGCCCCACAGCGCGGTGGTTGGGGGTGCGCTGG : 150574

    250 : alLeuThrLysProLeuGlyThrGlnValAlaValAsnAlaHisGlnTrpLeuAspAs :    268
          .!:!!|||||||||||||||||||||||||||||||||||||||||||||||||||||
          laValThrLysProLeuGlyThrGlnValAlaValAsnAlaHisGlnTrpLeuAspAs
 150575 : CGGTAACCAAACCGTTAGGAACCCAGGTTGCTGTCAATGCCCACCAATGGCTGGATAA : 150631

    269 : nProGluArgTrpAsnLysValLysMetValValSerArgGluGluValGluLeuAla :    287
          |||||||||| !!|||||||||||||||||||||||| !!|||:!!||||||||||||
          nProGluArgArgAsnLysValLysMetValValSerGlyGluGlnValGluLeuAla
 150632 : TCCTGAAAGACGGAATAAAGTAAAGATGGTGGTCTCCGGAGAACAGGTGGAGCTGGCA : 150688

    288 : TyrGlnGluAlaMetPheAsnMetAlaThrLeuAsnArgThrAlaAlaGlyLeuMetH :    307
          |||||||||||||||||||||||||||||||||||||||.!!||||||||||||||||
          TyrGlnGluAlaMetPheAsnMetAlaThrLeuAsnArgAlaAlaAlaGlyLeuMetH
 150689 : TACCAGGAAGCCATGTTCAATATGGCTACCCTCAACAGAGCTGCTGCAGGTTTAATGC : 150748

    308 : isThrPheAsnAlaHisAlaAlaThrAspIleThrGlyPheGlyIleLeuGlyHisSe :    326
          ||.!!|||||||||||||||||||||||||||||||||!!.|||||||||||||||||
          isAlaPheAsnAlaHisAlaAlaThrAspIleThrGlyLeuGlyIleLeuGlyHisSe
 150749 : ACGCATTTAATGCCCATGCGGCCACAGATATCACAGGCTTAGGCATTTTAGGACACTC : 150805

    327 : rGlnAsnLeuAlaLysGlnGlnArgAsnGluValSerPheValIleHisAsnLeuPro :    345
          | !!||||||||||||||||||||||||||||||||||||||||||||||||||||||
          r***AsnLeuAlaLysGlnGlnArgAsnGluValSerPheValIleHisAsnLeuPro
 150806 : CTAGAACCTTGCGAAACAACAAAGAAATGAAGTGTCTTTTGTTATTCATAATCTGCCA : 150862

    346 : IleIleAlaLysMetAlaAlaValSerLysAlaSerGlyArgPheGlyLeuLeuGlnG :    365
          ||||||||||||||||||!.!||||||||||||||||||!:!|||||||||! !||||
          IleIleAlaLysMetAlaValValSerLysAlaSerGlyGlnPheGlyLeuProGlnG
 150863 : ATAATTGCCAAGATGGCTGTCGTCAGCAAGGCCAGTGGACAGTTTGGGCTCCCTCAAG : 150922

    366 : lyThrSerAlaGluThrSerGlyGly----LeuLeuIleCysLeuProArgGluGlnA :    383
          ||||||||||||||||||||||||||####|||||||||||||||||| !!|||||||
          lyThrSerAlaGluThrSerGlyGly####LeuLeuIleCysLeuProGlyGluGlnA
 150923 : GAACCTCAGCTGAAACCTCTGGGGGATTGATTACTGATTTGTCTGCCTGGAGAACAGG : 150980

    384 : laAlaArgPheCysSerGluIleLysSerSerLysTyrGlyGluGlyHisGlnAlaTr :    402
          ||||||||||||||! !||||||||||||||||||||| !!|||.!!|||||||||||
          laAlaArgPheCysPheGluIleLysSerSerLysTyrArgGluSerHisGlnAlaTr
 150981 : CGGCTCGCTTTTGTTTCGAAATCAAATCCTCCAAGTACCGAGAGAGTCACCAAGCGTG : 151037

    403 : pIleValGlyIleValGluLysGlyAsnArgThrAlaArgIleIleAspLysProArg :    421
          |||||||||||||||||||||||||!:!!:!|||||||||||||||||||||||| !!
          pIleValGlyIleValGluLysGlySerGlnThrAlaArgIleIleAspLysPro***
 151038 : GATCGTTGGCATTGTGGAAAAGGGAAGCCAAACGGCCCGGATCATTGACAAGCCTTGA : 151094

    422 : ValIleGluValLeuProArgGlyAlaThrAlaAlaValLeuAlaProAspSerSerA :    441
          ||||||||||||||||||||||||:!!|||||||||!.!||| !!||||||!:!||||
          ValIleGluValLeuProArgGlySerThrAlaAlaAlaLeuProProAspAsnSerA
 151095 : GTTATTGAAGTCCTACCTCGTGGATCCACAGCTGCTGCTCTTCCTCCTGACAATTCAA : 151154

    442 : snAlaSerSerGluProSerSer :    448
          |||||||||||||||||!!!|||
          snAlaSerSerGluProCysSer
 151155 : ATGCCTCCTCTGAGCCTTGCTCG : 151177

vulgar: SPP00000033_1.0 171 448 . gi|357434203|gb|JH378160.1|:subseq(5751640,301147) 150342 151177 + 1161 M 202 606 F 0 4 M 75 225
# --- START OF GFF DUMP ---
#
#
##gff-version 2
##source-version exonerate:protein2genome:local 2.2.0
##date 2012-11-14
##type DNA
#
#
# seqname source feature start end score strand frame attributes
#
gi|357434203|gb|JH378160.1|:subseq(5751640,301147)	exonerate:protein2genome:local	gene	150343	151177	1161	+	.	gene_id 1 ; sequence SPP00000033_1.0 ; gene_orientation .
gi|357434203|gb|JH378160.1|:subseq(5751640,301147)	exonerate:protein2genome:local	cds	150343	151177	.	+	.	
gi|357434203|gb|JH378160.1|:subseq(5751640,301147)	exonerate:protein2genome:local	exon	150343	151177	.	+	.	insertions 0 ; deletions 0 ; frameshifts 4
gi|357434203|gb|JH378160.1|:subseq(5751640,301147)	exonerate:protein2genome:local	similarity	150343	151177	1161	+	.	alignment_id 1 ; Query SPP00000033_1.0 ; Align 150343 172 606 ; Align 150953 374 225
# --- END OF GFF DUMP ---
#