Command line: [exonerate -m p2g --showtargetgff -q gpx1a_zf.aa.fa.canviat -t genomegpx1a_JXSJ01000122.1.fa] Hostname: [ws101097u] C4 Alignment: ------------ Query: SPP00000621_2.0 # Protein # Glutathione peroxidase 1a (GPx1a) # Zebrafish Target: JXSJ01000122.1:subseq(419000,101574) Miichthys miiuy scaffold122, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 794 Query range: 5 -> 187 Target range: 51574 -> 49999 6 : SerPheTyrAspIleThrAlaLysThrLeuThrGlyGluGluPheLysPheSerSerLeu : 25 !!!|||||||||:!!! !|||!!. !||||||||||||! !|||!!.|||||||||||| ThrPheTyrAspLeuLysAlaAsnLeuLeuThrGlyGluAlaPheAsnPheSerSerLeu 51574 : ACATTTTACGACCTCAAGGCCAACCTCTTAACAGGAGAAGCGTTTAATTTCTCGTCCCTG : 51517 26 : GlnGlyLysValValLeuIleGluAsnValAlaSerLeuUnkGlyThrThrSerArgAsp : 45 ||||||||||||||||||||||||||||||||||||||| |||||||||!!!|||||| GlnGlyLysValValLeuIleGluAsnValAlaSerLeu***GlyThrThrThrArgAsp 51516 : CAGGGAAAAGTTGTCCTCATTGAGAATGTCGCGTCTCTCTGAGGTACGACCACCAGGGAT : 51457 46 : TyrThrGlnMetAsnGluLeuHisGluArgPheAlaGluLysGlyLeuValValLeuGly : 65 ||||||||||||||||||||||||!!:|||!:!|||...||||||||||||:!!|||||| TyrThrGlnMetAsnGluLeuHisAspArgTyrAlaSerLysGlyLeuValIleLeuGly 51456 : TACACCCAGATGAACGAGCTCCATGACCGCTATGCCAGTAAGGGGCTCGTTATCCTGGGA : 51397 66 : ValProCysAsnGlnPheGlyTyrGln >>>> Target Intron 1 >>>> GluA : 76 |||||||||||||||||||||:!!||| 1029 bp |||| ValProCysAsnGlnPheGlyHisGln++ ++GluA 51396 : GTGCCCTGCAACCAGTTCGGCCATCAGgt.........................agGAGA : 50335 77 : snCysThrAsnGluGluIleLeuLeuSerLeuLysTyrValArgProGlyAsnGlyTyrG : 96 |||||! !||||||||||||||||||||||||||||||||||||||||||||||||!:!| snCysLysAsnGluGluIleLeuLeuSerLeuLysTyrValArgProGlyAsnGlyPheG 50334 : ACTGTAAGAATGAAGAAATCCTCTTGTCTCTGAAGTACGTCCGTCCTGGAAATGGCTTTG : 50275 97 : luProLysPheGlnLeuLeuGluLysValAspValAsnGlyLysAsnAlaHisProLeuP : 116 ||||||||||||||||||||||||||||||||||||||||||||:!!|||||| !!|||| luProLysPheGlnLeuLeuGluLysValAspValAsnGlyLysAspAlaHisSerLeuP 50274 : AGCCAAAGTTTCAGCTCCTGGAGAAGGTGGATGTGAATGGAAAGGATGCCCATTCCCTGT : 50215 117 : heThrPheLeuLysGluLysLeuProPheProSerAspGluProMetProPheMetSerA : 136 ||..!||||||!:!|||!..||||||! !||||||||||||||| !!|||!!:!!!| heValPheLeuArgGluSerLeuProSerProSerAspGluProCysAlaPheIleThrA 50214 : TTGTGTTCCTGAGAGAAAGTCTCCCATCCCCCAGCGACGAGCCGTGCGCCTTCATTACTG : 50155 137 : spProLysPheIleIleTrpSerProValCysArgAsnAspIleAlaTrpAsnPheGluL : 156 ||||||||!!.:!!|||||||||||||||||||||||||||:!:|||||||||||||||| spProLysLeuValIleTrpSerProValCysArgAsnAspValAlaTrpAsnPheGluL 50154 : ACCCCAAGTTAGTCATCTGGAGCCCGGTCTGCAGGAACGATGTGGCGTGGAACTTTGAGA : 50095 157 : ysPheLeuIleGlySerAspGlyValProPheLysArgTyrSerArgArgPheLeuThrS : 176 |||||||||||||| !!||||||||||||||||||||||||||||||!:!|||||||||! ysPheLeuIleGlyProAspGlyValProPheLysArgTyrSerArgLysPheLeuThrA 50094 : AGTTCCTCATCGGGCCCGATGGGGTTCCCTTCAAGCGTTACAGCAGGAAGTTCCTCACCA : 50035 177 : erGlyIleAspGlyAspIleLysLysLeuLeuSer : 187 :!! !|||!!:|||||||||||||||||||||!:! snAspIleGluGlyAspIleLysLysLeuLeuAsn 50034 : ATGACATCGAGGGAGACATCAAGAAGCTCCTCAAC : 50000 vulgar: SPP00000621_2.0 5 187 . JXSJ01000122.1:subseq(419000,101574) 51574 49999 - 794 M 69 207 5 0 2 I 0 1025 3 0 2 M 113 339 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2016-11-16 ##type DNA # # # seqname source feature start end score strand frame attributes # JXSJ01000122.1:subseq(419000,101574) exonerate:protein2genome:local gene 50000 51574 794 - . gene_id 1 ; sequence SPP00000621_2.0 ; gene_orientation + JXSJ01000122.1:subseq(419000,101574) exonerate:protein2genome:local cds 51368 51574 . - . JXSJ01000122.1:subseq(419000,101574) exonerate:protein2genome:local exon 51368 51574 . - . insertions 0 ; deletions 0 JXSJ01000122.1:subseq(419000,101574) exonerate:protein2genome:local splice5 51366 51367 . - . intron_id 1 ; splice_site "GT" JXSJ01000122.1:subseq(419000,101574) exonerate:protein2genome:local intron 50339 51367 . - . intron_id 1 JXSJ01000122.1:subseq(419000,101574) exonerate:protein2genome:local splice3 50339 50340 . - . intron_id 0 ; splice_site "AG" JXSJ01000122.1:subseq(419000,101574) exonerate:protein2genome:local cds 50000 50338 . - . JXSJ01000122.1:subseq(419000,101574) exonerate:protein2genome:local exon 50000 50338 . - . insertions 0 ; deletions 0 JXSJ01000122.1:subseq(419000,101574) exonerate:protein2genome:local similarity 50000 51574 794 - . alignment_id 1 ; Query SPP00000621_2.0 ; Align 51575 6 207 ; Align 50339 75 339 # --- END OF GFF DUMP --- # -- completed exonerate analysis