Command line: [exonerate -m p2g --showtargetgff --exhaustive yes -q selw_sili.aa.fa -t genome_selw_tri_subseq.fa]
Hostname: [luke]

C4 Alignment:
------------
         Query: W_E.siliculosus
        Target: gi|217405287|gb|CM000620.1|:subseq(321111,20000) Phaeodactylum tricornutum CCAP 1055/1 chromosome 18, whole genome shotgun sequence:[revcomp]
         Model: protein2genome:local
     Raw score: 155
   Query range: 47 -> 144
  Target range: 10390 -> 9999

    48 : SerThrThrThrAlaGlyArgGlyGluAspValAlaGluIleGluIleMetTyrCysGly :    67
         |||  !|||! !::!  !! !|||!  :!!  !  !...|||:!!:!:..!|||!! |||
         SerGlnThrLysSerThrMetGlyAlaAsnLysLeuArgIleGlnValGlnTyr***Gly
 10390 : TCTCAGACGAAGAGCACCATGGGCGCTAACAAGCTTCGCATTCAGGTGCAGTACTGAGGT : 10333

    68 : GlyUnk{G}  >>>> Target Intron 1 >>>>  {ly}TyrValArgPheTyrAsn :    76
         |||   {|}            94 bp            {||}|||! !|||!:!!:!|||
         Gly***{G}++                         ++{ly}TyrGlyArgTyrPheAsn
 10332 : GGATGA{G}gt.........................ag{GC}TACGGTCGCTATTTTAAC : 10212

    77 : ArgLeuThrProGlnLeu<-><->GluAlaValPheGlyAspArgIleSerPheIleGln :    94
          !!|||! !  !   |||      :!! !!  ! ! !.!:!!!:!|||..!.!.:!!   
         GlyLeuLysAspPheLeuSerSerGlnProSerValAlaAsnLysIleGluLeuValGly
 10211 : GGGTTGAAAGACTTTCTCTCCTCTCAACCCTCCGTAGCGAATAAGATTGAGCTAGTCGGT : 10152

    95 : ThrGluAspProGlyValThrGlyArgPheGluValValLeuArgAsnThrGlyGluLeu :   114
         ! !|||||||||  !|||||||||!.!|||||||||  !:!!  !   ||||||:!!|||
         LysGluAspProThrValThrGlyAsnPheGluValSerValAlaAlaThrGlyGlnLeu
 10151 : AAGGAAGATCCTACGGTGACGGGAAATTTTGAGGTGTCCGTCGCAGCGACTGGGCAACTT : 10092

   115 : LeuHisSerLysLysThrArgGlyGlnGlyArgCysGluAspProArgGluValGlnAla :   134
         ||||||||||||::!  !  !|||||||||!:!...  !..!||| !!|||  !  !|||
         LeuHisSerLysArgGlnAlaGlyGlnGlyLysAlaThrSerProGlyGluArgAlaAla
 10091 : TTGCATTCCAAGCGACAGGCCGGCCAAGGAAAGGCGACGTCCCCGGGCGAACGGGCGGCG : 10032

   135 : LeuValAspGlnIleGlnLysPheLeuAsp :   144
         :!!:!!!!:||||||:!!:!!.!.||||||
         IleLeuGluGlnIleLysGluLeuLeuAsp
 10031 : ATTTTGGAGCAAATAAAAGAGCTATTGGAT : 10000

vulgar: W_E.siliculosus 47 144 . gi|217405287|gb|CM000620.1|:subseq(321111,20000) 10390 9999 - 155 M 22 66 S 0 1 5 0 2 I 0 90 3 0 2 S 1 2 M 12 36 G 0 6 M 62 186
# --- START OF GFF DUMP ---
#
#
##gff-version 2
##source-version exonerate:protein2genome:local 2.2.0
##date 2011-03-11
##type DNA
#
#
# seqname source feature start end score strand frame attributes
#
gi|217405287|gb|CM000620.1|:subseq(321111,20000)	exonerate:protein2genome:local	gene	10000	10390	155	-	.	gene_id 1 ; sequence W_E.siliculosus ; gene_orientation +
gi|217405287|gb|CM000620.1|:subseq(321111,20000)	exonerate:protein2genome:local	cds	10324	10390	.	-	.	
gi|217405287|gb|CM000620.1|:subseq(321111,20000)	exonerate:protein2genome:local	exon	10324	10390	.	-	.	insertions 0 ; deletions 0
gi|217405287|gb|CM000620.1|:subseq(321111,20000)	exonerate:protein2genome:local	splice5	10322	10323	.	-	.	intron_id 1 ; splice_site "GT"
gi|217405287|gb|CM000620.1|:subseq(321111,20000)	exonerate:protein2genome:local	intron	10230	10323	.	-	.	intron_id 1
gi|217405287|gb|CM000620.1|:subseq(321111,20000)	exonerate:protein2genome:local	splice3	10230	10231	.	-	.	intron_id 0 ; splice_site "AG"
gi|217405287|gb|CM000620.1|:subseq(321111,20000)	exonerate:protein2genome:local	cds	10000	10229	.	-	.	
gi|217405287|gb|CM000620.1|:subseq(321111,20000)	exonerate:protein2genome:local	exon	10000	10229	.	-	.	insertions 6 ; deletions 0
gi|217405287|gb|CM000620.1|:subseq(321111,20000)	exonerate:protein2genome:local	similarity	10000	10390	155	-	.	alignment_id 1 ; Query W_E.siliculosus ; Align 10391 48 66 ; Align 10228 71 36 ; Align 10186 83 186
# --- END OF GFF DUMP ---
#
-- completed exonerate analysis