Command line: [exonerate -m p2g --showtargetgff -q DIO3a.x.fa -t DIO3a.extraction.KV884750.1.fa]
Hostname: [sitdoc]

C4 Alignment:
------------
         Query: DIO3a # Protein # Iodothyronine deiodinase 3a (DIO3a) # Zebrafish
        Target: KV884750.1:subseq(768675,100000) Monopterus albus unplaced genomic scaffold scaffold62.1, whole genome shotgun sequence
         Model: protein2genome:local
     Raw score: 919
   Query range: 5 -> 265
  Target range: 50020 -> 50791

     6 : GlnGlySerAlaGlyValGlnSerAlaLeuLysAsnAlaAlaValCysValLeuLeuLeu :    25
         ...|||.!!!.!  !:!!  !!!!|||||||||:!!||||||:!!|||:!!:!!||||||
         SerGlyGlyValGlnMetAlaArgAlaLeuLysHisAlaAlaLeuCysLeuMetLeuLeu
 50021 : TCTGGCGGTGTCCAAATGGCGAGGGCGCTTAAACATGCAGCCCTGTGCCTGATGCTGCTG : 50078

    26 : ProArgPheLeuLeuAlaAlaLeuMetLeuCysLeuLeuAspPheLeuCysIleArgArg :    45
         |||||||||||||||||||||:!!||||||!! !!!|||||||||||||||||||||!:!
         ProArgPheLeuLeuAlaAlaValMetLeuTrpPheLeuAspPheLeuCysIleArgLys
 50079 : CCGCGCTTCCTCTTGGCTGCAGTCATGCTATGGTTCTTGGATTTCTTGTGCATCAGGAAA : 50138

    46 : LysLeuLeuLeuLysMetGlnGluGlyAlaPheSerSerProAspAspProProLeuArg :    65
         |||:!!||||||||||||  !|||     !.!...!||||||||||||||||||:!! !!
         LysMetLeuLeuLysMetGlyGlu---LysLeuAspSerProAspAspProProValCys
 50139 : AAAATGCTGCTGAAAATGGGAGAG---AAGCTGGACAGCCCGGACGACCCGCCGGTGTGT : 50195

    66 : ValSerAspSerAsnLysMetPheThrLeuGluSerLeuArgAlaValTrpTyrGlyGln :    85
         ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
         ValSerAspSerAsnLysMetPheThrLeuGluSerLeuArgAlaValTrpTyrGlyGln
 50196 : GTCTCCGACTCTAATAAGATGTTCACCTTGGAGTCCCTCAGGGCCGTGTGGTACGGCCAG : 50255

    86 : LysLeuAspPhePheLysSerAlaArgLeuGlyGlyAlaAlaProAsnThrGluValPhe :   105
         ||||||||||||.!!|||||||||!.!||||||  !.!!|||||||||||||||||| ! 
         LysLeuAspPheLeuLysSerAlaHisLeuGlyHisThrAlaProAsnThrGluValVal
 50256 : AAATTAGACTTTCTCAAATCTGCGCACCTCGGGCACACTGCGCCCAACACCGAGGTGGTG : 50315

   106 : ProLeuAspGlyAspAlaArgAlaAlaGluArgIleLeuAspTyrAlaArgGlyArgArg :   125
         ! !:!!      .!.! !|||  !  !!  ||||||||||||! !  !!:!! !!:!|||
         LeuVal------GlnGluArgArgGlnValArgIleLeuAspCysMetLysGluLysArg
 50316 : CTGGTC------CAGGAGCGGAGGCAGGTCCGGATCCTGGACTGTATGAAAGAGAAGAGA : 50369

   126 : ProLeuIleLeuAsnPheGlySerCysSerUnkProProPheMetThrArgLeuSerAla :   145
         ||||||||||||||||||||||||||||||   |||||||||||||||||||||:!!|||
         ProLeuIleLeuAsnPheGlySerCysSer***ProProPheMetThrArgLeuAlaAla
 50370 : CCGCTCATTCTTAACTTTGGCAGCTGCTCCTGACCGCCATTCATGACGCGCCTGGCGGCG : 50429

   146 : PheGlnArgValAlaArgGlnTyrAlaAspIleAlaAspSerLeuLeuValTyrIleGlu :   165
         |||||||||:!!!.!!!!|||||||||||||||||||||! !|||:!!||||||||||||
         PheGlnArgLeuValSerGlnTyrAlaAspIleAlaAspPheLeuValValTyrIleGlu
 50430 : TTTCAGCGCCTCGTAAGCCAGTACGCAGATATTGCGGACTTTTTAGTAGTATATATCGAG : 50489

   166 : GluAlaHisProSerAspGlyTrpValSerSerAspAlaProValGlnIleProArgHis :   185
         ||||||||||||||||||||||||||||||||||||||||||  !||||||! !!:!|||
         GluAlaHisProSerAspGlyTrpValSerSerAspAlaProTyrGlnIleArgLysHis
 50490 : GAGGCGCATCCGTCGGACGGCTGGGTGAGCTCGGACGCACCGTATCAGATCCGCAAGCAC : 50549

   186 : ArgCysLeuGluAspArgLeuArgAlaAlaGlnMetLeuHisArgGluAlaAlaGlySer :   205
         |||||||||||||||||||||||||||||||||:!!:!!! !  !|||!.! !!||||||
         ArgCysLeuGluAspArgLeuArgAlaAlaGlnLeuMetLeuSerGluValProGlySer
 50550 : CGCTGCTTGGAGGACAGACTGAGAGCCGCTCAGCTGATGCTCTCCGAGGTGCCGGGGAGC : 50609

   206 : AlaValValValAspSerMetGlnAsnSerCysAsnAlaAlaTyrGlyAlaTyrPheGlu :   225
           !||||||||||||!:!|||.!.||||||! !|||||||||||||||||||||||||||
         AsnValValValAspAsnMetAspAsnSerSerAsnAlaAlaTyrGlyAlaTyrPheGlu
 50610 : AACGTAGTGGTGGATAATATGGACAACTCGTCCAACGCCGCGTACGGAGCCTACTTTGAG : 50669

   226 : ArgLeuTyrIleValLysAspAlaThrValValTyrGlnGlyGlyArgGlyProGluGly :   245
         |||||||||||||||!:!|||! !! !|||||||||||||||||||||||||||||||||
         ArgLeuTyrIleValArgAspGluArgValValTyrGlnGlyGlyArgGlyProGluGly
 50670 : AGACTTTACATCGTGAGGGATGAAAGAGTGGTGTATCAGGGAGGCAGAGGTCCAGAGGGA : 50729

   246 : TyrArgIleAlaGluLeuArgAspTrpLeuGluArgTyrArgSerGlyLeuGlnGluSer :   265
         |||||||||:!!|||||||||:!!|||||||||!:!||||||!.!! !|||  !.!.|||
         TyrArgIleSerGluLeuArgAsnTrpLeuGluGlnTyrArgLysAspLeuValAsnSer
 50730 : TACCGGATTTCCGAGCTCAGAAACTGGCTGGAGCAATACAGGAAGGATTTGGTGAATTCC : 50789

vulgar: DIO3a 5 265 . KV884750.1:subseq(768675,100000) 50020 50791 + 919 M 48 144 G 1 0 M 53 159 G 2 0 M 156 468
# --- START OF GFF DUMP ---
#
#
##gff-version 2
##source-version exonerate:protein2genome:local 2.2.0
##date 2017-11-23
##type DNA
#
#
# seqname source feature start end score strand frame attributes
#
KV884750.1:subseq(768675,100000)	exonerate:protein2genome:local	gene	50021	50791	919	+	.	gene_id 1 ; sequence DIO3a ; gene_orientation .
KV884750.1:subseq(768675,100000)	exonerate:protein2genome:local	cds	50021	50791	.	+	.	
KV884750.1:subseq(768675,100000)	exonerate:protein2genome:local	exon	50021	50791	.	+	.	insertions 0 ; deletions 3
KV884750.1:subseq(768675,100000)	exonerate:protein2genome:local	similarity	50021	50791	919	+	.	alignment_id 1 ; Query DIO3a ; Align 50021 6 144 ; Align 50165 55 159 ; Align 50324 110 468
# --- END OF GFF DUMP ---
#
-- completed exonerate analysis