Command line: [/Users/enekovillanuevaverdejo/Desktop/exonerate-1.4.0/src/program/exonerate -m p2g --showtargetgff --exhaustive yes -q /Users/enekovillanuevaverdejo/Desktop/msp_ddi_aa.fa -t /Users/enekovillanuevaverdejo/Desktop/scafold_17.fa] Hostname: [MBP-Eneko.local] C4 Alignment: ------------ Query: DDB0304589 |DDB_G0277015 |Protein| gene: msp on chromosome: 2 position 7225040 to 7225867 Target: scaffold_17 [revcomp] Model: protein2genome Raw score: 898 Query range: 4 -> 242 Target range: 39980 -> 39197 5 : AsnLeuProLysValAspLeuGlyGluAspCysGluGlyGlyValUnkAlaArgProSer : 24 |||||||||||||||!!::!!|||!!:|||||||||||||||||| |||||||||!!! AsnLeuProLysValGluValGlyAspAspCysGluGlyGlyVal***AlaArgProThr 39980 : AATTTACCAAAGGTAGAGGTTGGAGATGACTGTGAGGGCGGAGTTTGAGCAAGGCCAACT : 39923 25 : SerSerSerSerSerIleAsn<->AsnAlaSerAspGluSerThrProLeuIleSerLys : 43 ||||||:!!|||||||||||| !..!.!.!!||||||||||||||||||||||||||| SerSerAlaSerSerIleAsnLysThrValGlyAspGluSerThrProLeuIleSerLys 39922 : TCATCAGCATCATCAATTAATAAAACAGTTGGTGATGAGTCTACCCCATTAATTTCAAAG : 39863 44 : ThrAsnAspGluGluLys<->AlaAsnIleGlyIleSerSerThrSerAsnSerProGln : 62 |||::!|||||||||||| ! ! !..! !||| !!:!!|||:!!..!|||||| ThrSerAspGluGluLysAspIleAlaCysSerSerSerProSerSerAspGlnProGln 39862 : ACATCTGATGAAGAAAAAGATATTGCTTGCTCATCCTCACCATCCTCTGATCAACCACAA : 39803 63 : GluGluGlnThrLysLysProLeuPheIleSerIleLeuThrLeuLeuIleSerIlePro : 82 |||||||||!:!||||||||||||.!!|||||||||||||||:!!||||||||||||||| GluGluGlnSerLysLysProLeuIleIleSerIleLeuThrValLeuIleSerIlePro 39802 : GAAGAACAAAGTAAAAAGCCATTAATCATCTCAATTCTTACAGTTTTAATTTCAATTCCA : 39743 83 : AlaLeuValGlySerUnkCysTrpProValLeuIleAlaSerLeuSerGlyValAlaVal : 102 ||||||||||||||| |||||||||:!!|||:!!::!|||:!! !|||:!!|||||| AlaLeuValGlySer***CysTrpProIleLeuLeuSerSerIleLeuGlyIleAlaVal 39742 : GCTTTAGTTGGTTCTTGATGTTGGCCTATCCTTCTAAGCAGCATTCTCGGTATTGCTGTT : 39683 103 : SerAlaGlySerValGluLeuAlaHisSerLeuThrPheAlaIleThrLeuSerIleLeu : 122 |||!.!!.!|||||| !!!!||||||:!!:!!|||!:!||||||||||||:!!|||||| SerGlyAlaSerValThrPheAlaHisAlaIleThrTyrAlaIleThrLeuAlaIleLeu 39682 : TCAGGAGCTTCCGTCACATTTGCCCACGCCATTACTTATGCCATTACATTAGCTATTTTA : 39623 123 : SerAsnLeuAlaGlnTyrHisPheHisLysCysLysLysArgProSerAspArgGlyHis : 142 |||||||||||||||||||||!:!|||||||||:!!...||||||!!!!!:!:!|||||| SerAsnLeuAlaGlnTyrHisTyrHisLysCysGlnSerArgProArgGluLysGlyHis 39622 : TCCAATTTAGCTCAATACCATTACCACAAATGTCAATCTAGACCAAGAGAAAAAGGCCAT : 39563 143 : TrpIleLysPheGlyProPheTyrLeuThrAlaIleAlaValProLeuAlaThrPheAsp : 162 ||| !|||!:!|||||||||||||||||| !:!:|||||||||||||||||||||||| TrpHisLysTyrGlyProPheTyrLeuThrLeuLeuAlaValProLeuAlaThrPheAsp 39562 : TGGCACAAATATGGACCATTTTATCTTACTTTATTGGCTGTTCCTTTAGCAACTTTTGAT : 39503 163 : IleLeuArgHisIleLeuValAspAsnSerIleTrp >>>> Target Intron 1 > : 175 |||||||||||||||||||||||||||..!:!!||| 63 bp IleLeuArgHisIleLeuValAspAsnGlnLeuTrp++ 39502 : ATTCTAAGACATATTTTAGTAGATAACCAACTTTGGgt...................... : 39462 176 : >>> ThrIleHisSerPheIleSerProAlaAlaTyrArgProGlyCysGluAsnGluA : 193 |||:!!|||||||||:!!||||||||||||||||||||||||||| !:!!|||| ++ThrValHisSerPheValSerProAlaAlaTyrArgProGlyCysThrHisGluA 39461 : ...agACTGTACATTCATTTGTTTCACCAGCAGCATATAGACCAGGTTGCACACACGAAA : 39347 194 : snIleThrCysLeuSerValMetGlyTrpPheSerAlaIleValPheThrTyrThrGlyT : 213 |||||! !|||||||||||||||||||||||||||||||||:!!|||||||||||||||| snIleArgCysLeuSerValMetGlyTrpPheSerAlaIleIlePheThrTyrThrGlyT 39346 : ATATTAGATGCCTTTCAGTCATGGGTTGGTTTTCAGCAATTATTTTCACATATACAGGTT : 39287 214 : yrAlaCysLeuLeuValGlyThrIleTrpAlaAlaAspLeuIleProLysIleLysLysV : 233 ||||| !||||||..!||||||||||||..!|||!!:||| ! !||||||||||||: yrAlaLeuLeuLeuThrGlyThrIleTrpCysAlaGluLeuHisIleLysIleLysLysL 39286 : ATGCCCTTTTATTAACTGGTACAATTTGGTGTGCTGAACTCCACATAAAGATCAAAAAGT : 39227 234 : alTrpThrGlnLeuArgProSerLysLys : 242 !!||| !||||||||||||!.!|||||| euTrpGlnGlnLeuArgProLysLysLys 39226 : TATGGCAACAACTCCGTCCAAAAAAGAAA : 39198 vulgar: DDB0304589 4 242 . scaffold_17 39980 39197 - 898 M 27 81 G 0 3 M 18 54 G 0 3 M 125 375 5 0 2 I 0 59 3 0 2 M 68 204 # --- START OF GFF DUMP --- # # ##gff-version 2 ##source-version exonerate:protein2genome 1.4.0 ##date 2009-03-06 ##type DNA # # # seqname source feature start end score strand frame attributes # scaffold_17 exonerate:protein2genome gene 39198 39980 898 - . gene_id 1 ; sequence DDB0304589 ; gene_orientation + scaffold_17 exonerate:protein2genome cds 39465 39980 . - . scaffold_17 exonerate:protein2genome exon 39465 39980 . - . insertions 6 ; deletions 0 scaffold_17 exonerate:protein2genome splice5 39463 39464 . - . intron_id 1 ; splice_site "GT" scaffold_17 exonerate:protein2genome intron 39402 39464 . - . intron_id 1 scaffold_17 exonerate:protein2genome splice3 39402 39403 . - . intron_id 0 ; splice_site "AG" scaffold_17 exonerate:protein2genome cds 39198 39401 . - . scaffold_17 exonerate:protein2genome exon 39198 39401 . - . insertions 0 ; deletions 0 scaffold_17 exonerate:protein2genome similarity 39198 39980 898 - . alignment_id 1 ; Query DDB0304589 ; Align 39981 5 81 ; Align 39897 32 54 ; Align 39840 50 375 ; Align 39402 175 204 # --- END OF GFF DUMP --- # -- completed exonerate analysis