Command line: [exonerate -m p2g --showtargetgff -q ../Sel.human/GPX4.fa -t gi|397462837|gb|JH732664.1|subseq.fa] Hostname: [ws101062] C4 Alignment: ------------ Query: SPP00000026_2.0 # Protein # Glutathione peroxidase 4 (GPx4) # Human (Homo sapiens) Target: gi|397462837|gb|JH732664.1|:subseq(220000,100000) Alligator mississippiensis unplaced genomic scaffold scaffold-2624, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 465 Query range: 15 -> 182 Target range: 51877 -> 51372 16 : CysGlyAlaLeuAlaAlaProGlyLeuAlaGlyThrMetCysAlaSerArgAspAspTrp : 35 |||!.!! !! !:!! !! !||| ! !!.! ! !||||||..! !||| !!||| CysAlaGluTrpSerTrpHisGlyAlaLeuAlaGlyAlaCysAlaGlnValAspTyrTrp 51877 : TGTGCGGAGTGGTCATGGCACGGTGCCTTGGCAGGAGCGTGTGCCCAGGTGGACTACTGG : 51820 36 : ArgCysAlaArgSerMetHisGluPheSerAlaLysAspIleAspGlyHisMetValAsn : 55 !!! !.!!!:!|||!!::!!|||! !!.!||||||||||||||| ! |||::: CysSerThrLysSerIleTyrGluSerHisValLysAspIleAspGlySerHisValSer 51819 : TGCTCTACCAAATCCATCTATGAATCACATGTCAAGGACATCGATGGCAGCCATGTTTCA : 51760 56 : LeuAspLysTyrArgGlyPheValCysIleValThrAsnValAlaSerGlnUnkGlyLys : 75 |||!!:|||||| !!|||!:!!.!||||||:!!||||||:!!||||||:!! |||||| LeuGluLysTyrTrpGlyTyrAlaCysIleIleThrAsnIleAlaSerLys***GlyLys 51759 : CTGGAAAAGTACTGGGGTTATGCCTGCATCATCACAAATATAGCTTCAAAATGAGGCAAA : 51700 76 : ThrGluValAsnTyrThrGlnLeuValAspLeuHisAlaArgTyrAlaGluCysGlyLeu : 95 |||!!::!!|||||||||||||||||||||:!!:!!|||! !!! .!!||| ! !!||| ThrAspIleAsnTyrThrGlnLeuValAspMetTyrAlaIle***ThrGluArgCysLeu 51699 : ACCGATATAAACTACACTCAGCTTGTCGATATGTATGCCATATAGACGGAGAGGTGTTTA : 51640 96 : ArgIleLeuAlaPheProCysAsnGlnPheGlyLysGlnGluProGlySerAsnGluGlu : 115 |||||||||||||||||||||||||||||| !|||||| !!||||||!:!|||! :!! ArgIleLeuAlaPheProCysAsnGlnPheLysLysGln***ProGlyAsnAsnAlaGln 51639 : CGCATCCTGGCCTTCCCCTGCAACCAGTTCAAGAAACAGTAACCAGGGAACAATGCCCAG : 51580 116 : IleLysGluPheAlaAlaGlyTyrAsnValLysPheAspMet----PheSerLysIleCy : 134 ||||||! |||!.!|||.!!|||...:!!|||!:!|||!!:####|||||||||! ! IleLysAlaPheValAlaSerTyrGlyMetLysTyrAspIle####PheSerLysThrAs 51579 : ATCAAAGCCTTTGTTGCAAGCTATGGAATGAAGTATGATATACACATTCAGCAAGACTGA : 51519 135 : sValAsnGlyAspAspAlaHisProLeuTrpLysTrpMetLysIleGlnProLysGlyLy : 154 !||| !!||||||:!!.!!! |||||||||||||||:!!:!! ||||||||||||!: pValTyrGlyAspAsnThrLeuProLeuTrpLysTrpLeuGlnGluGlnProLysGlyAr 51518 : TGTCTACGGGGACAACACCCTGCCTTTGTGGAAGTGGCTGCAGGAGCAGCCCAAAGGCAG : 51459 155 : sGlyIleLeuGlyAsnAlaIleLysTrpAsnPheThrLysPheLeuIleAspLysAsnGl : 174 !|||! !|||||||||||||||||||||||||||||||||||||||||||||!!..!.|| gGlyThrLeuGlyAsnAlaIleLysTrpAsnPheThrLysPheLeuIleAspAsnGluGl 51458 : AGGCACCCTAGGAAATGCAATAAAATGGAACTTCACAAAGTTCCTTATTGACAATGAAGG : 51399 175 : yCysValValLysArgTyrGlyPro : 182 |! !:!! ! ! !:! ! |||||| ySerLeuTrpAspGlnGluGlyPro 51398 : ATCCCTATGGGATCAAGAAGGACCT : 51373 vulgar: SPP00000026_2.0 15 182 . gi|397462837|gb|JH732664.1|:subseq(220000,100000) 51877 51372 - 465 M 114 342 F 0 4 M 53 159 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2013-11-15 ##type DNA # # # seqname source feature start end score strand frame attributes # gi|397462837|gb|JH732664.1|:subseq(220000,100000) exonerate:protein2genome:local gene 51373 51877 465 - . gene_id 1 ; sequence SPP00000026_2.0 ; gene_orientation . gi|397462837|gb|JH732664.1|:subseq(220000,100000) exonerate:protein2genome:local cds 51373 51877 . - . gi|397462837|gb|JH732664.1|:subseq(220000,100000) exonerate:protein2genome:local exon 51373 51877 . - . insertions 0 ; deletions 0 ; frameshifts 4 gi|397462837|gb|JH732664.1|:subseq(220000,100000) exonerate:protein2genome:local similarity 51373 51877 465 - . alignment_id 1 ; Query SPP00000026_2.0 ; Align 51878 16 342 ; Align 51532 130 159 # --- END OF GFF DUMP --- # -- completed exonerate analysis