Command line: [exonerate -m p2g --showtargetgff -q gpx2.fa -t genomic.fa] Hostname: [ws101183u] C4 Alignment: ------------ Query: SPP00000616_2.0 # Protein # Glutathione peroxidase 2 (GPx2) # Zebrafish Target: CM004422.1:subseq(9320000,150000) Ictalurus punctatus breed USDA103 chromosome 9, whole genome shotgun sequence Model: protein2genome:local Raw score: 873 Query range: 0 -> 189 Target range: 1381 -> 2729 1 : MetThrPheIleAlaLysThrPheTyrAspLeuHisAlaThrThrLeuGluGlyAspThrIl : 21 ||||||||||||||||||! !||||||||||||!.!|||||||||||||||||||||:!!|| MetThrPheIleAlaLysMetPheTyrAspLeuArgAlaThrThrLeuGluGlyAspSerIl 1382 : ATGACCTTCATAGCCAAGATGTTTTATGACCTGCGTGCTACCACGCTGGAAGGGGACTCTAT : 1442 22 : eAspPheAsnIleTyrArgGlyArgValValLeuIleGluAsnValAlaSerLeuUnkGlyT : 42 ||||||||||:!!|||||||||||||||||||||||||||||||||||||||||| |||| eAspPheAsnValTyrArgGlyArgValValLeuIleGluAsnValAlaSerLeu***GlyT 1443 : CGATTTTAATGTCTACCGAGGACGAGTGGTACTTATCGAGAATGTTGCATCGCTATGAGGCA : 1505 43 : hrThrThrGlnAspTyrThrGlnLeuAsnGluLeuGlnSerArgTyrProHisArgLeuVal : 62 ||||||||!:!||||||!:!|||||||||:!!||||||||||||||||||||||||!!!||| hrThrThrArgAspTyrSerGlnLeuAsnGlnLeuGlnSerArgTyrProHisArgPheVal 1506 : CCACCACCCGGGACTACAGCCAGCTCAACCAGCTCCAGAGCAGGTACCCGCACCGGTTTGTG : 1565 63 : ValLeuGlyPheProCysAsnGlnPheGlyTyrGln >>>> Target Intron 1 >>> : 75 ||||||!.!||||||||||||||||||||||||||| 781 bp ValLeuAlaPheProCysAsnGlnPheGlyTyrGln++ 1566 : GTCCTGGCCTTCCCCTGTAATCAGTTTGGCTACCAGgt........................ : 1606 76 : > GluAsnCysSerAspGlyGluIleLeuAsnSerLeuLysTyrValArgProGlyGluGl : 94 |||||||||!!!:!!||||||||||||||||||||||||:!!||||||||||||! !|| ++GluAsnCysArgAsnGlyGluIleLeuAsnSerLeuLysHisValArgProGlyGlyGl 1607 : .agGAAAACTGCAGAAATGGAGAGATCCTGAACTCCTTGAAGCACGTACGACCAGGTGGAGG : 2442 95 : yTyrLysProSerPheThrIlePheGlnLysCysValValAsnGlySerAspAlaHisProV : 115 |!:!:!!|||! !||||||||||||:!!||||||! !|||||||||:!!..!.!!||||||| yPheGluProLeuPheThrIlePheGluLysCysAspValAsnGlyAlaSerThrHisProV 2443 : GTTTGAGCCTTTGTTCACCATTTTTGAGAAGTGCGATGTGAATGGGGCTAGCACACATCCTG : 2505 116 : alPheSerTyrLeuLysAspLysLeuProTyrProAspAspAspProValThrLeuIleGln : 135 |||||:!!||||||||||||||||||||||||||||||||||||||||||:!!||||||||| alPheAlaTyrLeuLysAspLysLeuProTyrProAspAspAspProValSerLeuIleGln 2506 : TCTTCGCATACCTGAAAGACAAACTTCCGTACCCCGATGATGATCCAGTCTCTCTTATTCAG : 2565 136 : AspProLysTyrLeuValTrpAsnProValSerArgAsnAspIleSerTrpAsnPheGluLy : 156 |||||||||!:!|||||||||!:!|||:!!||||||!.!||||||||||||||||||||||| AspProLysPheLeuValTrpSerProIleSerArgThrAspIleSerTrpAsnPheGluLy 2566 : GACCCCAAGTTCTTAGTGTGGAGTCCCATCAGCAGAACTGATATCTCATGGAATTTTGAGAA : 2628 157 : sPheLeuIleGlyProGluGlyGluProPheLysArgTyrSerLysLysPheGlnThrIleA : 177 ||||||||||||||||||||||||||||||||||||||||||||||!!.|||:!!||||||| sPheLeuIleGlyProGluGlyGluProPheLysArgTyrSerLysAsnPheGluThrIleA 2629 : GTTCCTGATTGGCCCAGAGGGCGAGCCCTTCAAGAGATACAGCAAAAACTTCGAGACCATTA : 2691 178 : snIleGluProAspIleGlnArgLeuLeuLysLeuThr : 189 |||||||||||||||||||||||||||||!:!|||||| snIleGluProAspIleGlnArgLeuLeuArgLeuThr 2692 : ACATCGAACCTGACATTCAAAGGCTATTGAGGCTGACC : 2729 vulgar: SPP00000616_2.0 0 189 . CM004422.1:subseq(9320000,150000) 1381 2729 + 873 M 74 222 5 0 2 I 0 777 3 0 2 M 115 345 # --- 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 # CM004422.1:subseq(9320000,150000) exonerate:protein2genome:local gene 1382 2729 873 + . gene_id 1 ; sequence SPP00000616_2.0 ; gene_orientation + CM004422.1:subseq(9320000,150000) exonerate:protein2genome:local cds 1382 1603 . + . CM004422.1:subseq(9320000,150000) exonerate:protein2genome:local exon 1382 1603 . + . insertions 0 ; deletions 0 CM004422.1:subseq(9320000,150000) exonerate:protein2genome:local splice5 1604 1605 . + . intron_id 1 ; splice_site "GT" CM004422.1:subseq(9320000,150000) exonerate:protein2genome:local intron 1604 2384 . + . intron_id 1 CM004422.1:subseq(9320000,150000) exonerate:protein2genome:local splice3 2383 2384 . + . intron_id 0 ; splice_site "AG" CM004422.1:subseq(9320000,150000) exonerate:protein2genome:local cds 2385 2729 . + . CM004422.1:subseq(9320000,150000) exonerate:protein2genome:local exon 2385 2729 . + . insertions 0 ; deletions 0 CM004422.1:subseq(9320000,150000) exonerate:protein2genome:local similarity 1382 2729 873 + . alignment_id 1 ; Query SPP00000616_2.0 ; Align 1382 1 222 ; Align 2385 75 345 # --- END OF GFF DUMP --- # -- completed exonerate analysis