Command line: [exonerate -m p2g --showtargetgff -q GPx4/GPx4_U.fa -t GPx4/SUBSEQ/KB228918.1_1_subseq.txt --exhaustive yes] Hostname: [luke] C4 Alignment: ------------ Query: GPx4 Target: KB228918.1:subseq(4487468,25074) Odobenus rosmarus divergens isolate Ivan unplaced genomic scaffold Scaffold41, whole genome shotgun sequence:[revcomp] Model: protein2genome:local Raw score: 632 Query range: 0 -> 149 Target range: 25074 -> 24357 1 : ArgCysAlaGlnSerMetHisGluPheSerAlaLysAspIleAspGlyHisMetValAsn : 20 ||||||||||||||||||||||||||||||||||||||||||||||||!.!||||||||| ArgCysAlaGlnSerMetHisGluPheSerAlaLysAspIleAspGlyArgMetValAsn 25074 : CGCTGTGCTCAGTCTATGCACGAGTTCTCCGCCAAGGACATCGATGGACGCATGGTTAAC : 25017 21 : LeuAspLysTyr{Ar} >>>> Target Intron 1 >>>> {g}GlyPheValCys : 29 ||||||||||||{||} 83 bp {|}|||||||||||| LeuAspLysTyr{Ar}++ ++{g}GlyPheValCys 25016 : CTGGACAAGTAC{CG}gt.........................ag{G}GGCTTCGTGTGC : 24907 30 : IleValThrAsnValAlaSerGlnUnkGlyLysThrGluValAsnTyrThrGlnLeuVal : 49 |||||||||||||||||||||||| |||||||||!!:||||||||||||||||||||| IleValThrAsnValAlaSerGln***GlyLysThrAspValAsnTyrThrGlnLeuVal 24906 : ATCGTCACCAACGTGGCCTCGCAGTGAGGCAAAACAGACGTAAACTATACTCAGCTTGTC : 24847 50 : AspLeuHisAlaArgTyrAlaGluCysGlyLeuArgIleLeuAlaPheProCysAsnGln : 69 |||||||||||||||||||||||| !!||||||||||||||||||||||||||||||||| AspLeuHisAlaArgTyrAlaGluGlyGlyLeuArgIleLeuAlaPheProCysAsnGln 24846 : GACCTGCACGCCCGATACGCTGAGGGTGGTTTACGGATCCTGGCCTTCCCTTGCAACCAG : 24787 70 : PheGlyArgGln >>>> Target Intron 2 >>>> CysProCysLeuProGlnG : 80 |||||||||||| 84 bp |||||||||||||||| PheGlyArgGln++ --HisProCysLeuProGlnG 24786 : TTCGGGAGGCAGgt.........................cacacccctgcttgcCACAGG : 24670 81 : luProGlySerAsnAlaGluIleLysGluPheAlaAlaGlyTyrAsnValLysPheAspM : 100 |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||| luProGlySerAsnAlaGluIleLysGluPheAlaAlaGlyTyrAsnValLysPheAspM 24669 : AGCCAGGGAGCAATGCGGAGATCAAAGAGTTCGCTGCTGGCTATAATGTCAAATTCGATA : 24610 101 : etPheSerLysIleCysValAsnGlyAspAspAlaHisProLeuTrpLysTrpMetLysG : 120 ||!:!||||||||||||||||||||||||||||||||||||||||||||||||||||||! etTyrSerLysIleCysValAsnGlyAspAspAlaHisProLeuTrpLysTrpMetLysV 24609 : TGTACAGCAAGATCTGTGTGAATGGGGACGATGCTCACCCGCTGTGGAAGTGGATGAAAG : 24550 121 : lyProArgGly >>>> Target Intron 3 >>>> ValArgAsnArgHisArgPh : 130 !! !! ! ! 103 bp ||| ! !!|||||||| alGlnProLys+- ++AspArgLeuCysHisArgPh 24549 : TCCAGCCCAAGgg.........................aggatCGCCTTTGCCACCGATT : 24417 131 : eLeuMetIleSerHisLeuSerSerSerAlaIleLysTrpAsnPheThrLysValArg : 149 |||||||! !||||||! !! !!!!||||||||||||||||||||||||||||||||| eLeuMetThrSerHisCysPheCysSerAlaIleLysTrpAsnPheThrLysValArg 24416 : CCTCATGACTTCACATTGCTTTTGTAGTGCCATCAAATGGAATTTCACCAAGGTAAGA : 24358 vulgar: 0 149 . KB228918.1:subseq(4487468,25074) 25074 24357 - 632 M 24 72 S 0 2 5 0 2 I 0 79 3 0 2 S 1 1 M 48 144 5 0 2 I 0 80 3 0 2 M 50 150 5 0 2 I 0 99 3 0 2 M 26 78 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome:local 2.2.0 ##date 2015-11-26 ##type DNA # # # seqname source feature start end score strand frame attributes # KB228918.1:subseq(4487468,25074) exonerate:protein2genome:local gene 24358 25074 632 - . gene_id 1 ; sequence ; gene_orientation + KB228918.1:subseq(4487468,25074) exonerate:protein2genome:local cds 25001 25074 . - . KB228918.1:subseq(4487468,25074) exonerate:protein2genome:local exon 25001 25074 . - . insertions 0 ; deletions 0 KB228918.1:subseq(4487468,25074) exonerate:protein2genome:local splice5 24999 25000 . - . intron_id 1 ; splice_site "GT" KB228918.1:subseq(4487468,25074) exonerate:protein2genome:local intron 24918 25000 . - . intron_id 1 KB228918.1:subseq(4487468,25074) exonerate:protein2genome:local splice3 24918 24919 . - . intron_id 0 ; splice_site "AG" KB228918.1:subseq(4487468,25074) exonerate:protein2genome:local cds 24773 24917 . - . KB228918.1:subseq(4487468,25074) exonerate:protein2genome:local exon 24773 24917 . - . insertions 0 ; deletions 0 KB228918.1:subseq(4487468,25074) exonerate:protein2genome:local splice5 24771 24772 . - . intron_id 2 ; splice_site "GT" KB228918.1:subseq(4487468,25074) exonerate:protein2genome:local intron 24689 24772 . - . intron_id 2 KB228918.1:subseq(4487468,25074) exonerate:protein2genome:local splice3 24689 24690 . - . intron_id 1 ; splice_site "ca" KB228918.1:subseq(4487468,25074) exonerate:protein2genome:local cds 24539 24688 . - . KB228918.1:subseq(4487468,25074) exonerate:protein2genome:local exon 24539 24688 . - . insertions 0 ; deletions 0 KB228918.1:subseq(4487468,25074) exonerate:protein2genome:local splice5 24537 24538 . - . intron_id 3 ; splice_site "GG" KB228918.1:subseq(4487468,25074) exonerate:protein2genome:local intron 24436 24538 . - . intron_id 3 KB228918.1:subseq(4487468,25074) exonerate:protein2genome:local splice3 24436 24437 . - . intron_id 2 ; splice_site "ag" KB228918.1:subseq(4487468,25074) exonerate:protein2genome:local cds 24358 24435 . - . KB228918.1:subseq(4487468,25074) exonerate:protein2genome:local exon 24358 24435 . - . insertions 0 ; deletions 0 KB228918.1:subseq(4487468,25074) exonerate:protein2genome:local similarity 24358 25074 632 - . alignment_id 1 ; Query ; Align 25075 1 72 ; Align 24917 26 144 ; Align 24689 74 150 ; Align 24436 124 78 # --- END OF GFF DUMP --- # -- completed exonerate analysis