Command line: [exonerate -m p2g --showtargetgff -q gpx1a.fa -t genomic.fa] Hostname: [ws101183u] C4 Alignment: ------------ Query: SPP00000621_2.0 # Protein # Glutathione peroxidase 1a (GPx1a) # Zebrafish Target: CM004424.1:subseq(15520000,150000) Ictalurus punctatus breed USDA103 chromosome 11, whole genome shotgun sequence Model: protein2genome:local Raw score: 834 Query range: 4 -> 190 Target range: 3229 -> 5291 5 : LysSerPheTyrAspIleThrAlaLysThrLeuThrGlyGluGluPheLysPheSerSerLe : 25 |||||||||||||||:!!:!!!.!|||||||||||||||||||||!!.|||!!.|||!!!|| LysSerPheTyrAspLeuSerValLysThrLeuThrGlyGluGluLeuLysLeuSerCysLe 3230 : AAGTCGTTTTACGACCTCTCGGTGAAGACGCTGACAGGAGAGGAGTTGAAGTTGTCGTGTTT : 3290 26 : uGlnGlyLysValValLeuIleGluAsnValAlaSerLeuUnkGlyThrThrSerArgAspT : 46 |||||||||||||||||||||||||||||||||||||||| |||||||||!!!||||||| uGlnGlyLysValValLeuIleGluAsnValAlaSerLeu***GlyThrThrThrArgAspT 3291 : ACAGGGTAAGGTAGTCCTGATCGAGAACGTTGCGTCTCTCTGAGGCACGACCACCAGGGATT : 3353 47 : yrThrGlnMetAsnGluLeuHisGluArgPheAlaGluLysGlyLeuValValLeuGlyVal : 66 |||||||||||||||||||||||||||||||||||! !:!!|||||||||:!!||||||||| yrThrGlnMetAsnGluLeuHisGluArgPheAlaAlaGluGlyLeuValIleLeuGlyVal 3354 : ACACCCAGATGAACGAGCTCCATGAGCGCTTCGCCGCGGAGGGACTCGTCATTCTGGGGGTG : 3413 67 : ProCysAsnGlnPheGlyTyrGln >>>> Target Intron 1 >>>> GluAsnCys : 77 ||||||||||||||||||:!!||| 1504 bp ||||||||| ProCysAsnGlnPheGlyHisGln++ ++GluAsnCys 3414 : CCCTGCAACCAGTTTGGCCATCAGgt.........................aggaaaactgT : 4950 78 : ThrAsnGluGluIleLeuLeuSerLeuLysTyrValArgProGlyAsnGlyTyrGluProLy : 98 ! !|||||||||||||||:!!||||||||||||:!!|||||||||||||||!:!|||||||| LysAsnGluGluIleLeuValSerLeuLysTyrIleArgProGlyAsnGlyPheGluProLy 4951 : AAAAACGAGGAGATCCTGGTGTCGCTGAAATACATCCGTCCTGGAAATGGCTTCGAACCAAA : 5013 99 : sPheGlnLeuLeuGluLysValAspValAsnGlyLysAsnAlaHisProLeuPheThrPheL : 119 |||||||||||||||||||||||||||||||||||||:!!|||||||||||||||..!|||| sPheGlnLeuLeuGluLysValAspValAsnGlyLysAspAlaHisProLeuPheValPheL 5014 : ATTCCAGCTCTTGGAGAAGGTCGATGTGAACGGGAAAGACGCTCACCCGCTGTTTGTCTTTC : 5076 120 : euLysGluLysLeuProPheProSerAspGluProMetProPheMetSerAspProLysPhe : 139 |||||||||||||||||||||||||||||||||||! ! !!!!.|||.!!|||||||||! ! euLysGluLysLeuProPheProSerAspGluProThrAlaLeuMetGlyAspProLysCys 5077 : TCAAGGAGAAGCTTCCTTTCCCAAGCGATGAGCCGACGGCATTAATGGGCGACCCCAAATGC : 5136 140 : IleIleTrpSerProValCysArgAsnAspIleAlaTrpAsnPheGluLysPheLeuIleGl : 160 |||||||||||||||||||||||||||||||||:!!|||||||||||||||||||||||||| IleIleTrpSerProValCysArgAsnAspIleSerTrpAsnPheGluLysPheLeuIleGl 5137 : ATCATCTGGAGTCCTGTGTGCAGAAACGACATCTCCTGGAACTTTGAGAAGTTCCTTATCGG : 5199 161 : ySerAspGlyValProPheLysArgTyrSerArgArgPheLeuThrSerGlyIleAspGlyA : 181 | !!|||||||||||||||||||||||||||||||||!:!||||||:!!! !|||!!:|||| yProAspGlyValProPheLysArgTyrSerArgArgTyrLeuThrAlaAspIleGluGlyA 5200 : CCCAGATGGCGTCCCATTCAAGCGTTACAGCAGAAGATACCTGACCGCTGACATAGAGGGAG : 5262 182 : spIleLysLysLeuLeuSerIleProLys : 190 |||||||||||||||||.!!||| !!||| spIleLysLysLeuLeuGlyIleAlaLys 5263 : ACATAAAGAAGCTTCTGGGCATTGCCAAA : 5291 vulgar: SPP00000621_2.0 4 190 . CM004424.1:subseq(15520000,150000) 3229 5291 + 834 M 70 210 5 0 2 I 0 1500 3 0 2 M 116 348 # --- 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 # CM004424.1:subseq(15520000,150000) exonerate:protein2genome:local gene 3230 5291 834 + . gene_id 1 ; sequence SPP00000621_2.0 ; gene_orientation + CM004424.1:subseq(15520000,150000) exonerate:protein2genome:local cds 3230 3439 . + . CM004424.1:subseq(15520000,150000) exonerate:protein2genome:local exon 3230 3439 . + . insertions 0 ; deletions 0 CM004424.1:subseq(15520000,150000) exonerate:protein2genome:local splice5 3440 3441 . + . intron_id 1 ; splice_site "GT" CM004424.1:subseq(15520000,150000) exonerate:protein2genome:local intron 3440 4943 . + . intron_id 1 CM004424.1:subseq(15520000,150000) exonerate:protein2genome:local splice3 4942 4943 . + . intron_id 0 ; splice_site "ag" CM004424.1:subseq(15520000,150000) exonerate:protein2genome:local cds 4944 5291 . + . CM004424.1:subseq(15520000,150000) exonerate:protein2genome:local exon 4944 5291 . + . insertions 0 ; deletions 0 CM004424.1:subseq(15520000,150000) exonerate:protein2genome:local similarity 3230 5291 834 + . alignment_id 1 ; Query SPP00000621_2.0 ; Align 3230 5 210 ; Align 4944 75 348 # --- END OF GFF DUMP --- # -- completed exonerate analysis