Command line: [exonerate -m p2g --showtargetgff -q SPP00000004_1.0.fa -t SPP00000004_1.0.GL397269.1/GL397269.1.subseq] Hostname: [upf-OptiPlex-755] C4 Alignment: ------------ Query: SPP00000004_1.0 # Protein # Glutathione peroxidase 1 (GPx1) # Homo sapiens # Complete Target: GL397269.1:subseq(37559082,100786) dna:supercontig supercontig:Nleu1.0:GL397269.1:1:47109586:1 REF:[revcomp] Model: protein2genome:local Raw score: 945 Query range: 0 -> 203 Target range: 50875 -> 49999 1 : MetCysAlaAlaArgLeuAlaAlaAlaAlaAlaAlaAlaGlnSerValTyrAlaPheSer : 20 |||||||||||||||||| |||||||||! !||||||||||||!! ||||||||| MetCysAlaAlaArgLeu------AlaAlaAlaAspAlaGlnSerVal***AlaPheSer 50875 : ATGTGTGCTGCTCGGCTA------GCGGCGGCGGATGCCCAGTCGGTGTAAGCCTTTTCC : 50824 21 : AlaArgProLeuAlaGlyGlyGluProValSerLeuGly--SerLeuArgGlyLysValL : 40 |||||||||||||||||||||||||||||||||||||||## ||||||||||||:!!| AlaArgProLeuAlaGlyGlyGluProValSerLeuGly##---LeuArgGlyLysIleL 50823 : GCGCGCCCGCTGGCCGGCGGGGAGCCTGTGAGCCTGGGCTC---CTGCGGGGCAAGATAC : 50765 41 : euLeuIleGluAsnValAlaSerLeuUnkGlyThrThrValArgAspTyrThrGlnMetA : 60 |||||||||||||||||||||||||| ||||||||||||||||||||||||||||||| euLeuIleGluAsnValAlaSerLeu***GlyThrThrValArgAspTyrThrGlnMetA 50764 : TACTTATCGAGAATGTGGCGTCCCTCTGAGGCACCACGGTCCGGGACTACACCCAGATGA : 50705 61 : snGluLeuGlnArgArgLeuGlyProArgGlyLeuValValLeuGlyPheProCysAsnG : 80 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| snGluLeuGlnArgArgLeuGlyProArgGlyLeuValValLeuGlyPheProCysAsnG 50704 : ACGAGCTGCAGCGGCGCCTCGGGCCCCGGGGCCTGGTGGTGCTCGGCTTCCCGTGCAACC : 50645 81 : lnPheGlyHisGln >>>> Target Intron 1 >>>> GluAsnAlaLysAsnGl : 90 |||||||||||||| 274 bp ||||||||||||||||| lnPheGlyHisGln++ ++GluAsnAlaLysAsnGl 50644 : AGTTTGGGCATCAGgt.........................agGAGAACGCCAAGAACGA : 50341 91 : uGluIleLeuAsnSerLeuLysTyrValArgProGlyGlyGlyPheGluProAsnPheMe : 110 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| uGluIleLeuAsnSerLeuLysTyrValArgProGlyGlyGlyPheGluProAsnPheMe 50340 : AGAGATTCTGAATTCCCTCAAGTACGTCCGACCTGGTGGTGGGTTCGAGCCCAACTTCAT : 50281 111 : tLeuPheGluLysCysGluValAsnGlyAlaGlyAlaHisProLeuPheAlaPheLeuAr : 130 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| tLeuPheGluLysCysGluValAsnGlyAlaGlyAlaHisProLeuPheAlaPheLeuAr 50280 : GCTCTTCGAGAAGTGCGAGGTGAACGGTGCGGGGGCGCACCCTCTCTTCGCCTTCCTGCG : 50221 131 : gGluAlaLeuProAlaProSerAspAspAlaThrAlaLeuMetThrAspProLysLeuIl : 150 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| gGluAlaLeuProAlaProSerAspAspAlaThrAlaLeuMetThrAspProLysLeuIl 50220 : GGAGGCCCTGCCAGCCCCCAGCGACGACGCCACTGCACTTATGACCGACCCCAAGCTCAT : 50161 151 : eThrTrpSerProValCysArgAsnAspValAlaTrpAsnPheGluLysPheLeuValGl : 170 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| eThrTrpSerProValCysArgAsnAspValAlaTrpAsnPheGluLysPheLeuValGl 50160 : CACCTGGTCTCCGGTGTGTCGCAACGATGTTGCCTGGAACTTTGAGAAGTTCCTGGTGGG : 50101 171 : yProAspGlyValProLeuArgArgTyrSerArgArgPheGlnThrIleAspIleGluPr : 190 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| yProAspGlyValProLeuArgArgTyrSerArgArgPheGlnThrIleAspIleGluPr 50100 : CCCTGACGGTGTGCCCCTACGCAGGTACAGCCGCCGCTTCCAGACCATTGACATCGAGCC : 50041 191 : oAspIleGluAlaLeuLeuSerGlnGlyProSerCysAla : 203 |||||||||||||||||||||||||||||||||||||||| oAspIleGluAlaLeuLeuSerGlnGlyProSerCysAla 50040 : TGACATCGAAGCCCTGCTCTCTCAAGGGCCCAGCTGTGCC : 50000 vulgar: SPP00000004_1.0 0 203 . GL397269.1:subseq(37559082,100786) 50875 49999 - 945 M 6 18 G 2 0 M 25 75 F 0 2 G 1 0 M 50 150 5 0 2 I 0 270 3 0 2 M 119 357 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2012-11-29 ##type DNA # # # seqname source feature start end score strand frame attributes # GL397269.1:subseq(37559082,100786) exonerate:protein2genome:local gene 50000 50875 945 - . gene_id 1 ; sequence SPP00000004_1.0 ; gene_orientation + GL397269.1:subseq(37559082,100786) exonerate:protein2genome:local cds 50631 50875 . - . GL397269.1:subseq(37559082,100786) exonerate:protein2genome:local exon 50631 50875 . - . insertions 0 ; deletions 3 ; frameshifts 2 GL397269.1:subseq(37559082,100786) exonerate:protein2genome:local splice5 50629 50630 . - . intron_id 1 ; splice_site "GT" GL397269.1:subseq(37559082,100786) exonerate:protein2genome:local intron 50357 50630 . - . intron_id 1 GL397269.1:subseq(37559082,100786) exonerate:protein2genome:local splice3 50357 50358 . - . intron_id 0 ; splice_site "AG" GL397269.1:subseq(37559082,100786) exonerate:protein2genome:local cds 50000 50356 . - . GL397269.1:subseq(37559082,100786) exonerate:protein2genome:local exon 50000 50356 . - . insertions 0 ; deletions 0 GL397269.1:subseq(37559082,100786) exonerate:protein2genome:local similarity 50000 50875 945 - . alignment_id 1 ; Query SPP00000004_1.0 ; Align 50876 1 18 ; Align 50858 9 75 ; Align 50781 35 150 ; Align 50357 85 357 # --- END OF GFF DUMP --- # -- completed exonerate analysis