Command line: [exonerate -m p2g --showtargetgff -q gpx3b.fa -t genomic.fa] Hostname: [ws101183u] C4 Alignment: ------------ Query: SPP00000623_2.0 # Protein # Glutathione peroxidase 3b (GPx3b) # Zebrafish Target: CM004426.1:subseq(6030000,150000) Ictalurus punctatus breed USDA103 chromosome 13, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 607 Query range: 0 -> 162 Target range: 5515 -> 1568 1 : MetAlaAlaAsnSerValTyrAspPheThrValGluThrLeuGluGluLysSerValSerLe : 21 ||||||::!!:!|||||||||:!!||||||:!!!!::!!||||||..!|||||||||||||| MetAlaSerSerSerValTyrAsnPheThrMetAspSerLeuGluArgLysSerValSerLe 5515 : ATGGCGAGCAGTTCAGTTTACAACTTCACCATGGACTCACTGGAAAGGAAATCTGTTTCTCT : 5455 22 : uAspIlePheArgGlyLysValLeuLeuIleMetAsnValAlaThrPheUnkGlySerThr : 42 |! :!!!:!|||||||||||||||||||||:!:||||||||||||||| ||||||||| uGlyValTyrArgGlyLysValLeuLeuIleValAsnValAlaThrPhe***GlySerThr+ 5454 : GGGGGTTTACAGGGGGAAAGTGCTTCTCATCGTCAACGTGGCCACCTTCTGAGGCTCCACTg : 5392 43 : >>>> Target Intron 1 >>>> GlnTyrHisArgMetAsnAlaLeuMetAspMetT : 53 358 bp ||||||||||||:!!|||||||||||||||! ! + --GlnTyrHisArgLeuAsnAlaLeuMetAspArgP 5391 : t.........................gtCAGTACCACAGGCTGAATGCACTCATGGACAGAT : 5001 54 : yrGlyGlyGlnArgPheThrValLeuGlyPheProCysAsnGlnPheGlyLeuGlnSerPro : 73 :!.!!|||:!!!.!||||||:!!||||||||| !!|||||||||||||||||||||:!!||| heSerGlyGluAsnPheThrIleLeuGlyPheSerCysAsnGlnPheGlyLeuGlnAlaPro 5000 : TCAGTGGTGAGAACTTCACTATTCTGGGATTCTCATGTAATCAGTTTGGACTACAGGCACCA : 4941 74 : {G} >>>> Target Intron 2 >>>> {ly}GluAsnHisGluThrLeuAsnValLe : 83 {!} 2647 bp { !}||||||||||||||||||:!!:!!|| {G}++ ++{lu}GluAsnHisGluThrLeuAspIleLe 4940 : {G}gt.........................ag{aG}GAGAATCACGAGACTTTGGACATTCT : 2264 84 : uGlnHisValArgProGlySerGlyPheLeuProLysPheProIlePheSerArgIleGluV : 104 |!:! !||||||||||||.!!.!!|||:!!||||||||||||!!:||||||!:!||||||| uArgCysValArgProGlyGlySerPheIleProLysPheProMetPheSerLysIleGluV 2263 : GCGGTGTGTCAGACCTGGAGGCAGTTTCATCCCAAAGTTCCCCATGTTCAGTAAAATCGAGG : 2201 105 : alAsnGlySerAspGluAspProLeuTyrAlaTyrLeuLys >>>> Target Intron : 118 |||||||| !:!!|||!!:|||||||||:!!!:!|||||| 456 bp alAsnGlyLeuAsnGluGluProLeuTyrSerPheLeuLys+- 2200 : TTAATGGGCTGAATGAGGAGCCTCTGTACTCCTTCCTGAAGgc................... : 2157 119 : 3 >>>> GluSerLeuProPheValAsnProValIleGlyAspIleArgLysLeuTyrTrp : 135 ! !|||||||||! !||||||||||||||||||||||||!:!|||||||||||| ++AlaSerLeuProSerValAsnProValIleGlyAspIleLysLysLeuTyrTrp 2156 : ......agGCGTCTCTGCCTTCTGTGAACCCTGTCATTGGAGACATCAAGAAACTCTATTGG : 1652 136 : SerProIleLysAlaAsnAspValArgTrpAsnPheGluLysPheLeuIleThrAlaAspGl : 156 ||||||||||||!.!||||||:!!|||||||||||||||||||||||||||||||||||||| SerProIleLysValAsnAspIleArgTrpAsnPheGluLysPheLeuIleThrAlaAspGl 1651 : TCTCCGATTAAAGTCAACGACATCCGCTGGAACTTTGAGAAGTTTCTTATCACTGCAGACGG : 1589 157 : yArgProTyrLysArgTyr : 162 |! !||||||!:!|||||| yLeuProTyrArgArgTyr 1588 : GCTGCCATACAGACGGTAT : 1569 vulgar: SPP00000623_2.0 0 162 . CM004426.1:subseq(6030000,150000) 5515 1568 - 607 M 41 123 5 0 2 I 0 354 3 0 2 M 32 96 S 0 1 5 0 2 I 0 2643 3 0 2 S 1 2 M 43 129 5 0 2 I 0 452 3 0 2 M 45 135 # --- 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 # CM004426.1:subseq(6030000,150000) exonerate:protein2genome:local gene 1569 5515 607 - . gene_id 1 ; sequence SPP00000623_2.0 ; gene_orientation + CM004426.1:subseq(6030000,150000) exonerate:protein2genome:local cds 5393 5515 . - . CM004426.1:subseq(6030000,150000) exonerate:protein2genome:local exon 5393 5515 . - . insertions 0 ; deletions 0 CM004426.1:subseq(6030000,150000) exonerate:protein2genome:local splice5 5391 5392 . - . intron_id 1 ; splice_site "GT" CM004426.1:subseq(6030000,150000) exonerate:protein2genome:local intron 5035 5392 . - . intron_id 1 CM004426.1:subseq(6030000,150000) exonerate:protein2genome:local splice3 5035 5036 . - . intron_id 0 ; splice_site "GT" CM004426.1:subseq(6030000,150000) exonerate:protein2genome:local cds 4938 5034 . - . CM004426.1:subseq(6030000,150000) exonerate:protein2genome:local exon 4938 5034 . - . insertions 0 ; deletions 0 CM004426.1:subseq(6030000,150000) exonerate:protein2genome:local splice5 4936 4937 . - . intron_id 2 ; splice_site "GT" CM004426.1:subseq(6030000,150000) exonerate:protein2genome:local intron 2291 4937 . - . intron_id 2 CM004426.1:subseq(6030000,150000) exonerate:protein2genome:local splice3 2291 2292 . - . intron_id 1 ; splice_site "ag" CM004426.1:subseq(6030000,150000) exonerate:protein2genome:local cds 2160 2290 . - . CM004426.1:subseq(6030000,150000) exonerate:protein2genome:local exon 2160 2290 . - . insertions 0 ; deletions 0 CM004426.1:subseq(6030000,150000) exonerate:protein2genome:local splice5 2158 2159 . - . intron_id 3 ; splice_site "GC" CM004426.1:subseq(6030000,150000) exonerate:protein2genome:local intron 1704 2159 . - . intron_id 3 CM004426.1:subseq(6030000,150000) exonerate:protein2genome:local splice3 1704 1705 . - . intron_id 2 ; splice_site "AG" CM004426.1:subseq(6030000,150000) exonerate:protein2genome:local cds 1569 1703 . - . CM004426.1:subseq(6030000,150000) exonerate:protein2genome:local exon 1569 1703 . - . insertions 0 ; deletions 0 CM004426.1:subseq(6030000,150000) exonerate:protein2genome:local similarity 1569 5515 607 - . alignment_id 1 ; Query SPP00000623_2.0 ; Align 5516 1 123 ; Align 5035 42 96 ; Align 2289 75 129 ; Align 1704 118 135 # --- END OF GFF DUMP --- # -- completed exonerate analysis