Command line: [exonerate -m p2g --showtargetgff -q gpx4a.fa -t genomic.fa] Hostname: [ws101183u] C4 Alignment: ------------ Query: SPP00000618_2.0 # Protein # Glutathione peroxidase 4a (GPx4a) # Zebrafish Target: CM004423.1:subseq(16400000,150000) Ictalurus punctatus breed USDA103 chromosome 10, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 460 Query range: 0 -> 113 Target range: 72556 -> 71233 1 : PheThrTyrPheArgGlyLysValValIleIleThrAsnValAlaSerLysUnkGlyLys : 20 ::: ::: |||||||||||||||||||||||||||||||||!!!:!! !.!||| TyrPhePheCysArgGlyLysValValIleIleThrAsnValAlaThrGln***AlaLys 72556 : tattttttctgtagaGGCAAGGTGGTTATCATCACCAACGTTGCCACCCAATGAGCCAAA : 72499 21 : ThrProValAsnTyrSerGlnPheAlaGluMetHisAlaLysTyrSerGluArgGlyLeu : 40 .!!||||||||||||||||||||||||:!!|||||||||||||||:!!|||||||||||| AlaProValAsnTyrSerGlnPheAlaGlnMetHisAlaLysTyrAlaGluArgGlyLeu 72498 : GCCCCAGTGAACTACTCTCAGTTTGCACAGATGCACGCCAAGTACGCTGAGAGAGGTTTA : 72439 41 : ArgIleLeuAlaPheProSerAsnGlnPheGlyArgGln >>>> Target Intron : 54 !|||||||||||||||!!!||||||||||||!.!||| 697 bp SerIleLeuAlaPheProCysAsnGlnPheGlyAsnGln++ 72438 : TCCATCCTGGCTTTCCCATGCAATCAGTTTGGCAACCAGgt................... : 72395 55 : 1 >>>> GluProGlyThrAsnSerGlnIleLysGluPheAlaLysSerTyrAsnAlaG : 71 ||||||||||||:!!..!:!!|||||||||!:!||||||||||||||||||: ++GluProGlyThrHisGluGluIleLysGluTyrAlaLysSerTyrAsnAlaG 72394 : ......agGAACCTGGCACCCATGAAGAGATCAAAGAATATGCCAAATCCTACAATGCGC : 71649 72 : luPheAspMetPheSerLysIleAspValAsnGlyAspGlyAlaHisProLeuTrpLysT : 91 !!||||||:!:|||||||||||| ! ||||||||| ! !|||||| !|||||||||| lnPheAspLeuPheSerLysIleLysValAsnGlyProGluAlaHisAspLeuTrpLysT 71648 : AATTTGACTTATTCAGCAAGATTAAAGTAAACGGACCAGAAGCTCATGACCTGTGGAAGT : 71589 92 : rpLeuLysAspGlnProAsnGlyLysGlyPheLeuGly{As} >>>> Target Intr : 104 ||||||||! |||||||||||||||||||||||||||{||} 287 bp rpLeuLysAlaGlnProAsnGlyLysGlyPheLeuGly{As}++ 71588 : GGTTGAAAGCGCAACCCAATGGCAAAGGATTCCTAGGA{AA}gt................ : 71546 105 : on 2 >>>> {n}GlyIleLysTrpAsnPheThrLysVal : 113 {|}||||||||||||||||||||||||||| ++{n}GlyIleLysTrpAsnPheThrLysVal 71545 : .........ag{T}GGAATCAAATGGAATTTCACCAAGGTA : 71234 vulgar: SPP00000618_2.0 0 113 . CM004423.1:subseq(16400000,150000) 72556 71233 - 460 M 53 159 5 0 2 I 0 693 3 0 2 M 50 150 S 0 2 5 0 2 I 0 283 3 0 2 S 1 1 M 9 27 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2016-11-14 ##type DNA # # # seqname source feature start end score strand frame attributes # CM004423.1:subseq(16400000,150000) exonerate:protein2genome:local gene 71234 72556 460 - . gene_id 1 ; sequence SPP00000618_2.0 ; gene_orientation + CM004423.1:subseq(16400000,150000) exonerate:protein2genome:local cds 72398 72556 . - . CM004423.1:subseq(16400000,150000) exonerate:protein2genome:local exon 72398 72556 . - . insertions 0 ; deletions 0 CM004423.1:subseq(16400000,150000) exonerate:protein2genome:local splice5 72396 72397 . - . intron_id 1 ; splice_site "GT" CM004423.1:subseq(16400000,150000) exonerate:protein2genome:local intron 71701 72397 . - . intron_id 1 CM004423.1:subseq(16400000,150000) exonerate:protein2genome:local splice3 71701 71702 . - . intron_id 0 ; splice_site "AG" CM004423.1:subseq(16400000,150000) exonerate:protein2genome:local cds 71549 71700 . - . CM004423.1:subseq(16400000,150000) exonerate:protein2genome:local exon 71549 71700 . - . insertions 0 ; deletions 0 CM004423.1:subseq(16400000,150000) exonerate:protein2genome:local splice5 71547 71548 . - . intron_id 2 ; splice_site "GT" CM004423.1:subseq(16400000,150000) exonerate:protein2genome:local intron 71262 71548 . - . intron_id 2 CM004423.1:subseq(16400000,150000) exonerate:protein2genome:local splice3 71262 71263 . - . intron_id 1 ; splice_site "AG" CM004423.1:subseq(16400000,150000) exonerate:protein2genome:local cds 71234 71261 . - . CM004423.1:subseq(16400000,150000) exonerate:protein2genome:local exon 71234 71261 . - . insertions 0 ; deletions 0 CM004423.1:subseq(16400000,150000) exonerate:protein2genome:local similarity 71234 72556 460 - . alignment_id 1 ; Query SPP00000618_2.0 ; Align 72557 1 159 ; Align 71701 54 150 ; Align 71261 105 27 # --- END OF GFF DUMP --- # C4 Alignment: ------------ Query: SPP00000618_2.0 # Protein # Glutathione peroxidase 4a (GPx4a) # Zebrafish Target: CM004423.1:subseq(16400000,150000) Ictalurus punctatus breed USDA103 chromosome 10, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 449 Query range: 0 -> 113 Target range: 91784 -> 88180 1 : PheThrTyrPheArgGlyLysValValIleIleThrAsnValAlaSerLysUnkGlyLys : 20 ::: ::: |||||||||||||||||||||||||||||||||!!!:!! !.!||| TyrPhePheCysArgGlyLysValValIleIleThrAsnValAlaThrGln***AlaLys 91784 : tattttttctgtagAGGCAAGGTGGTTATCATCACCAACGTTGCCACCCAATGAGCCAAA : 91727 21 : ThrProValAsnTyrSerGlnPheAlaGluMetHisAlaLysTyrSerGluArgGlyLeu : 40 .!!||||||||||||||||||||||||:!!|||||||||||||||:!!|||||||||||| AlaProValAsnTyrSerGlnPheAlaGlnMetHisAlaLysTyrAlaGluArgGlyLeu 91726 : GCCCCAGTGAACTACTCTCAGTTTGCACAGATGCACGCCAAGTACGCTGAGAGAGGTTTA : 91667 41 : ArgIleLeuAlaPheProSerAsnGlnPheGlyArgGln >>>> Target Intron : 54 !|||||||||||||||!!!||||||||||||!.!||| 2543 bp SerIleLeuAlaPheProCysAsnGlnPheGlyAsnGln++ 91666 : TCCATCCTGGCTTTCCCATGCAATCAGTTTGGCAACCAGgt................... : 91623 55 : 1 >>>> GluProGlyThrAsnSerGlnIleLysGluPheAlaLysSerTyrAsnAlaG : 71 ||||||||||||:!!:!!|||||||||.!. !!.!||||||||||||!.!: ++GluProGlyThrHisAlaGlnIleLysHisHisValLysSerTyrAsnValG 91622 : ......agGAACCTGGCACCCATGCACAGATCAAACACCATGTCAAATCCTACAATGTGC : 89031 72 : luPheAspMetPheSerLysIleAspValAsnGlyAspGlyAlaHisProLeuTrpLysT : 91 !!||||||:!:|||||||||||| ! |||||||||.!.! !|||||| !|||||||||| lnPheAspLeuPheSerLysIleLysValAsnGlyGlnGluAlaHisAspLeuTrpLysT 89030 : AATTTGATTTATTCAGTAAGATCAAAGTAAACGGACAAGAAGCTCATGACCTGTGGAAGT : 88971 92 : rpLeuLysAspGlnProAsnGlyLysGlyPheLeuGly{As} >>>> Target Intr : 104 ||:!!|||! |||||||||||||||||||||||||||{||} 722 bp rpMetLysAlaGlnProAsnGlyLysGlyPheLeuGly{As}++ 88970 : GGATGAAGGCGCAACCCAATGGCAAAGGATTCCTAGGA{AA}gt................ : 88928 105 : on 2 >>>> {n}GlyIleLysTrpAsnPheThrLysVal : 113 {|}||||||||||||||||||||||||||| ++{n}GlyIleLysTrpAsnPheThrLysVal 88927 : .........ag{T}GGAATCAAATGGAATTTCACCAAGGTA : 88181 vulgar: SPP00000618_2.0 0 113 . CM004423.1:subseq(16400000,150000) 91784 88180 - 449 M 53 159 5 0 2 I 0 2539 3 0 2 M 50 150 S 0 2 5 0 2 I 0 718 3 0 2 S 1 1 M 9 27 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2016-11-14 ##type DNA # # # seqname source feature start end score strand frame attributes # CM004423.1:subseq(16400000,150000) exonerate:protein2genome:local gene 88181 91784 449 - . gene_id 2 ; sequence SPP00000618_2.0 ; gene_orientation + CM004423.1:subseq(16400000,150000) exonerate:protein2genome:local cds 91626 91784 . - . CM004423.1:subseq(16400000,150000) exonerate:protein2genome:local exon 91626 91784 . - . insertions 0 ; deletions 0 CM004423.1:subseq(16400000,150000) exonerate:protein2genome:local splice5 91624 91625 . - . intron_id 1 ; splice_site "GT" CM004423.1:subseq(16400000,150000) exonerate:protein2genome:local intron 89083 91625 . - . intron_id 1 CM004423.1:subseq(16400000,150000) exonerate:protein2genome:local splice3 89083 89084 . - . intron_id 0 ; splice_site "AG" CM004423.1:subseq(16400000,150000) exonerate:protein2genome:local cds 88931 89082 . - . CM004423.1:subseq(16400000,150000) exonerate:protein2genome:local exon 88931 89082 . - . insertions 0 ; deletions 0 CM004423.1:subseq(16400000,150000) exonerate:protein2genome:local splice5 88929 88930 . - . intron_id 2 ; splice_site "GT" CM004423.1:subseq(16400000,150000) exonerate:protein2genome:local intron 88209 88930 . - . intron_id 2 CM004423.1:subseq(16400000,150000) exonerate:protein2genome:local splice3 88209 88210 . - . intron_id 1 ; splice_site "AG" CM004423.1:subseq(16400000,150000) exonerate:protein2genome:local cds 88181 88208 . - . CM004423.1:subseq(16400000,150000) exonerate:protein2genome:local exon 88181 88208 . - . insertions 0 ; deletions 0 CM004423.1:subseq(16400000,150000) exonerate:protein2genome:local similarity 88181 91784 449 - . alignment_id 2 ; Query SPP00000618_2.0 ; Align 91785 1 159 ; Align 89083 54 150 ; Align 88208 105 27 # --- END OF GFF DUMP --- # -- completed exonerate analysis