Resultados del exonerate
Command line: [exonerate -m p2g -q eEFSec.fa -t rblast3160_15500-18500.fa]
Hostname: [am75211]
C4 Alignment:
------------
Query: SPP00000038_1.0 # Protein # Eukaryotic elongation factor (eEFSec) # Homo sapiens # Complete
Target: Contig3160:subseq(15500,3000) [revcomp]
Model: protein2genome
Raw score: 111
Query range: 90 -> 149
Target range: 985 -> 1165
91 : ValAspCysProGlyHisAlaSerLeuIleArgThrIleIleGlyGlyAlaGlnIleIleAs : 111
:!!||||||||||||||||||..!! !|||!:!!.!!!:||| !|||||| !!!:||
IleAspCysProGlyHisAlaAspTyrIleLysAsnMetIleThrGlyAlaValGlnMetAs
986 : ATAGACTGTCCAGGACACGCAGATTATATAAAAAATATGATTACTGGTGCTGTTCAAATGGA : 1046
112 : pLeuMetMetLeuValIleAspValThrLysGlyMetGlnThrGlnSerAlaGlu<->CysL : 131
| ! !!:|||||||||...:!!||| ! ||| ..! !!|||!!! !||| ! |
pGlyAlaIleLeuValIleSerLeuThrAspGlyProMetProGlnThrIleGluHisLeuL
1047 : CGGGGCTATATTAGTTATTTCATTGACAGACGGTCCAATGCCGCAAACAATAGAACATTTAT : 1109
132 : euValIleGlyGlnIleAlaCysGlnLysLeuValValValLeuAsnLysIleAsp : 149
||:!!||| !||||||!.! :!!!!.:!!:!!:!! !!||||||||| !|||
euLeuIleLysGlnIleGlyIleGluAsnValIleIlePheLeuAsnLysGluAsp
1110 : TACTTATAAAACAAATAGGTATAGAAAATGTTATAATATTTTTAAATAAAGAAGAT : 1165
vulgar: SPP00000038_1.0 90 149 . Contig3160:subseq(15500,3000) 985 1165 + 111 M 39 117 G 0 3 M 20 60
-- completed exonerate analysis
Command line: [exonerate -m p2g -q GPx2.fa -t blast2766_4000-7000.fa]
Hostname: [am75250]
-- completed exonerate analysis
Command line: [exonerate -m p2g -q DI2.fa -t rblast2984_3000-4500.fa]
Hostname: [am75211]
-- completed exonerate analysis
Command line: [exonerate -m p2g -q SPS1.fa -t rblast3973_25000-5500.fa]
Hostname: [am75211]
Command line: [exonerate -m p2g -q Sel15.fa -t rblast4131.0_9500-12500.fa]
Hostname: [am75211]
-- completed exonerate analysis
Command line: [exonerate -m p2g -q SelI.fa -t blast4096_45000-48000.fa]
Hostname: [am75250]
C4 Alignment:
------------
Query: SPP00000015_1.0 # Protein # Selenoprotein I (SelI) # Homo sapiens # Complete
Target: Contig4096:subseq(45000,3000)
Model: protein2genome
Raw score: 125
Query range: 99 -> 156
Target range: 1101 -> 1272
100 : AlaTyrThrLeuAspGlyValAspGlyLysGlnAlaArgArgThrAsnSerSerThrProLe : 120
||| ! |||!!!||||||:!!||||||||||||||||||!:! !..!! !|||:!!|||||
AlaGlnThrPheAspGlyIleAspGlyLysGlnAlaArgLysLeuGlyMetSerSerProLe
1102 : GCGCAGACGTTCGACGGAATCGACGGGAAGCAGGCCAGAAAGCTGGGCATGAGCTCGCCCCT : 1162
121 : uGlyGluLeuPheAspHisGlyLeuAspSerTrpSerCysValTyrPheValValThrValT : 141
||||:!!|||.!.|||||||||||||||:!! ! |||!:!!:! ! !! ! !!
uGlyGlnLeuLeuAspHisGlyLeuAspAlaValValThrValPheTyrProTyrIleCysC
1163 : GGGGCAGCTGCTGGACCACGGCCTCGACGCCGTAGTGACGGTCTTCTACCCGTACATCTGCT : 1225
142 : yrSerIlePheGlyArgGlySerThrGlyValSerValPheValLeu : 156
!!!!:!!!:! ! !!|||! !||| !:!!!!!:!!.!!:!!:!!
ysThrLeuTyrProGlyGlyPheThrPheIleThrLeuLeuLeuVal
1226 : GCACGCTCTACCCCGGCGGCTTCACGTTCATCACCCTCCTCCTGGTG : 1272
vulgar: SPP00000015_1.0 99 156 . Contig4096:subseq(45000,3000) 1101 1272 + 125 M 57 171
-- completed exonerate analysis
Command line: [exonerate -m p2g -q SelR1.fa -t rblast4165_81000-84000.fa]
Hostname: [am75211]
-- completed exonerate analysis
Command line: [exonerate -m p2g -q SelR3.fa -t rblast1507_0-1000.fa]
Hostname: [am75211]
-- completed exonerate analysis
Command line: [exonerate -m p2g -q SelS.fa -t rblast4088_27500-28500.fa]
Hostname: [am75211]
-- completed exonerate analysis
Command line: [exonerate -m p2g -q SelU3.fa -t rblast1916_1000-4000.fa]
Hostname: [am75211]
-- completed exonerate analysis
Command line: [exonerate -m p2g -q SelV.fa -t blast4158_21000-24000.fa]
Hostname: [am75250]
-- completed exonerate analysis
Command line: [exonerate -m p2g -q SelW2.fa -t rblast4102_6000-9000.fa]
Hostname: [am75211]
-- completed exonerate analysis
Command line: [exonerate -m p2g -q TR1.fa -t blast4075_18500-21500.fa]
Hostname: [am75250]
C4 Alignment:
------------
Query: SPP00000034_1.0 # Protein # Thioredoxin reductase 1 (TR1) # Homo sapiens # Complete
Target: Contig4075:subseq(18500,3000)
Model: protein2genome
Raw score: 574
Query range: 12 -> 262
Target range: 471 -> 1259
13 : TyrAspLeuIleIleIleGlyGlyGlySerGlyGlyLeuAlaAlaAlaLysGluAlaAlaGl : 33
||||||!!! !:!!||||||||||||!!!||||||||||||||||||||||||||||||..
TyrAspPheAlaValIleGlyGlyGlyCysGlyGlyLeuAlaAlaAlaLysGluAlaAlaSe
472 : TACGACTTCGCCGTCATTGGCGGAGGTTGCGGTGGATTGGCCGCCGCTAAGGAGGCTGCCAG : 532
34 : nTyrGlyLysLysValMetValLeuAspPheValThrProThrProLeuGlyThrArgTrpG : 54
.! ||| !|||..!!!::!!!!!|||!:!||| !|||:!!|||! !||||||! !||||
rLeuGlyAlaLysThrIleLeuPheAspTyrValArgProSerProArgGlyThrThrTrpG
533 : CTTGGGAGCGAAGACGATTCTATTTGATTATGTGCGCCCTTCTCCGCGTGGCACAACTTGGG : 595
55 : lyLeuGlyGlyThrCysValAsnValGlyCysIleProLysLysLeuMetHisGlnAlaAla : 74
||||||||||||||||||||||||||||||||||||||||||||||||||||| ! |||!.!
lyLeuGlyGlyThrCysValAsnValGlyCysIleProLysLysLeuMetHisTyrAlaGly
596 : GTCTCGGCGGTACCTGTGTGAATGTCGGCTGTATCCCGAAGAAGTTGATGCACTATGCTGGC : 655
75 : LeuLeuGlyGlnAlaLeuGlnAspSerArgAsnTyrGlyTrpLysValGluGluThrValLy : 95
:!!||||||!!.:!!! !!!.||| !!..! !|||||| !.!.||| !
IleLeuGlyHisSerSerHisAspArgGluAlaLeuGlyTrpGlyAsnHisGluGly---Pr
656 : ATTTTGGGCCACTCGTCGCATGATCGTGAGGCGCTTGGTTGGGGCAACCATGAGGGT---CC : 715
96 : sHisAspTrpAspArgMetIleGlu >>>> Target Intron 1 >>>> AlaValGl : 106
!|||||||||! !|||:!::!!.!. 38 bp .!!|||||
oHisAspTrpGlyArgLeuValAsn++ ++ThrValGl
716 : GCACGATTGGGGCCGTCTTGTTAACgt.........................agACCGTGCA : 786
107 : nAsnHisIleGlySerLeuAsnTrpGlyTyrArgValAlaLeuArgGluLysLysValValT : 127
||||:!!||| !! !||||||!::..!|||||| !!.!|||! !..!|||!!.|||! !|
nAsnTyrIleLysMetLeuAsnPheSerTyrArgSerGlyLeuMetSerLysAsnValGluT
787 : GAACTACATCAAGATGCTCAACTTTTCCTACCGCTCGGGTCTCATGTCGAAAAACGTCGAAT : 849
128 : yrGluAsnAlaTyrGlyGlnPheIleGlyProHisArgIleLysAlaThrAsnAsnLysGly : 147
||! |||||| !.!....!! !! ! !:!!! !:!:! ! !|||:!!!!.||||||
yrValAsnAlaMetAlaSerLeuAlaAspLysAsnThrValThrTyrThrAspLysLysGly
850 : ACGTCAACGCCATGGCGAGCCTTGCCGATAAAAACACGGTGACATACACCGACAAGAAGGGG : 909
148 : LysGluLysIleTyrSerAlaGluArgPheLeuIleAlaThrGlyGluArgProArgTyr<- : 168
:!!:!! ! !!.!|||:!!!.! ! |||||||||! !|||! !||||||! ! !
GluLysHisGlnIleLysAlaLysAsnValLeuIleAlaIleGlyAlaArgProThrIlePr
910 : GAGAAGCATCAGATCAAGGCCAAGAATGTGCTGATCGCTATCGGAGCGCGTCCCACCATCCC : 972
169 : >LeuGlyIleProGlyAspLysGluTyrCysIleSerSerAspAspLeuPheSerLeuProT : 188
!! !:!: !|||! !!!:|||! |||!!!||||||||||||.!.||| ! !
oSerAspValLysGlyAlaTrpAspTyrSerIleThrSerAspAspLeuMetSerArgLysG
973 : CAGCGATGTGAAGGGTGCATGGGACTACTCAATCACCTCTGACGATTTAATGAGCAGGAAGG : 1035
189 : yrCysProGlyLysThrLeuValValGlyAlaSerTyrValAlaLeuGluCysAlaGlyPhe : 208
! ! !||||||||||||:!!||||||!.!|||!:!||||||||||||||||||||||||
luProValGlyLysThrLeuIleValGlyGlySerPheValAlaLeuGluCysAlaGlyPhe
1036 : AGCCCGTCGGCAAGACTCTTATCGTGGGCGGTTCGTTCGTTGCTCTGGAGTGCGCCGGCTTC : 1095
209 : LeuAlaGlyIleGlyLeuAspValThrValMetValArgSerIleLeuLeuArgGlyPheAs : 229
|||.!!..!!!:|||! !|||||||||||| |||||||||:!!:!!||||||||||||||
LeuThrSerMetGlyTyrAspValThrValAlaValArgSerLeuIleLeuArgGlyPheAs
1096 : CTCACCTCGATGGGTTACGACGTCACCGTCGCCGTCCGTTCTCTTATACTTAGAGGCTTCGA : 1158
230 : pGlnAspMetAlaAsnLysIleGlyGluHisMetGluGluHisGlyIleLysPheIleArgG : 250
|!::.!. |||:!!|||:!! !|||! !||| !! ! !|||! !||||||! !.!
pArgGlnCysAlaAspLysValGlnGluLeuMetLeuAlaThrGlyThrLysPheLysAsnG
1159 : CCGCCAGTGTGCGGACAAGGTCCAGGAGCTTATGTTGGCGACTGGCACCAAGTTCAAGAACG : 1221
251 : lnPheValProIleLysValGluGlnIleGluAlaGly : 262
! !!|||||| ! :!! :!!:!!||| !|||
lyValValProGlnAlaIleThrLysLeuGluAsnGly
1222 : GAGTTGTGCCGCAAGCCATCACCAAACTCGAGAATGGG : 1259
vulgar: SPP00000034_1.0 12 262 . Contig4075:subseq(18500,3000) 471 1259 + 574 M 81 243 G 1 0
M 9 27 5 0 2 I 0 34 3 0 2 M 64 192 G 0 3 M 95 285
C4 Alignment:
------------
Query: SPP00000034_1.0 # Protein # Thioredoxin reductase 1 (TR1) # Homo sapiens # Complete
Target: Contig4075:subseq(18500,3000)
Model: protein2genome
Raw score: 341
Query range: 282 -> 491
Target range: 1295 -> 1989
283 : GluTyrAsnThrValMetLeuAlaIleGlyArgAspAlaCysThrArg<->LysIleGlyLe : 302
|||!:!:!!|||:!!|||! !|||! !||||||..!!.!! !:!!||| :!!! ! !||
GluPheAspThrLeuMetTyrAlaThrGlyArgSerValSerSerArgMetGlnLysGluLe
1296 : GAATTTGACACGCTGATGTACGCTACCGGCCGTTCTGTGTCCTCAAGGATGCAAAAGGAGCT : 1356
303 : uGluThrValGlyValLysIleAsnGluLysThrGlyLysIleProValThrAspGluGluG : 323
|... !||||||:!!|||.!!!:!||| !||||||||| :!!.!!!!:!!:! :
uSerAspValGlyIleLysPheSerGlu---TyrGlyLysIle---IleAlaGluAspGlyL
1357 : CAGCGACGTGGGCATCAAATTCAGCGAG---TACGGCAAAATC---ATTGCCGAGGATGGCA : 1413
324 : lnThrAsnValProTyrIleTyrAlaIleGlyAspIleLeuGluAspLysValGluLeuThr : 343
!!|||::!||| ! !:!:||||||:!!||||||:!::!!|||! !!!. !! |||.!!
ysThrSerValGluGlyValTyrAlaValGlyAspValValGluGlyAsnProAlaLeuAla
1414 : AGACCTCTGTGGAGGGCGTGTACGCCGTCGGTGACGTGGTGGAGGGCAACCCTGCTCTCGCA : 1473
344 : ProValAlaIleGlnAlaGlyArgLeuLeuAlaGlnArgLeuTyrAlaGlySerThrValLy : 364
|||||||||:!!:!!! !|||..!|||||||||!::|||:!!!:!!.!..!||| ! !||
ProValAlaValLysAspGlyGluLeuLeuAlaArgArgIlePheGlyAsnSerAspLysLy
1474 : CCGGTGGCCGTTAAGGACGGTGAACTGCTTGCCCGCAGGATCTTCGGCAACTCAGACAAGAA : 1536
365 : s >>>> Target Intron 1 >>>> CysAspTyrGluAsnValProThrThrValPh : 375
| 37 bp !||| .!. !!:!!|||! ! !|||||
s++ ++LeuAspLeuAsnTyrIleProMetCysValPh
1537 : Ggt.........................agCTTGATCTGAACTACATTCCCATGTGTGTCTT : 1606
376 : eThrProLeuGluTyrGlyAlaCysGlyLeuSerGluGluLysAlaValGluLysPheGlyG : 396
|||||||! !||||||!.! !||||||:!!|||||||||! !||| !:!! !!:!||||
eThrProTyrGluTyrAlaArgCysGlyIleSerGluGluThrAlaSerLysLeuTyrGlyG
1607 : CACGCCGTACGAATACGCGCGGTGTGGAATTTCGGAGGAGACAGCCTCCAAATTGTACGGGG : 1669
397 : luGluAsnIleGluValTyrHisSerTyrPheTrpProLeuGluTrpThrIleProSerArg : 416
||!!: .!!!!:||||||! !.! ! !:! ! !!||||||!::.!! ! ! !|||
luAsp---PheAspValTyrLeuLysGluTyrThrThrLeuGluPheAlaAlaValHisArg
1670 : AAGAC---TTTGATGTCTACCTGAAGGAGTACACGACTCTGGAGTTTGCGGCTGTTCACCGC : 1726
417 : AspAsn >>>> Target Intron 2 >>>> AsnLysCysTyrAlaLysIleIleCys : 427
!!:!!. 42 bp !! ||| :!!|||!!:||||||
GluLys++ -+ProThrCysLeuSerLysMetIleCys
1727 : GAGAAGgt.........................cgCCTACTTGCCTGTCCAAAATGATCTGC : 1801
428 : AsnThrLysAspAsnGluArgValValGlyPheHisValLeuGlyProAsnAlaGlyGluVa : 448
||| ! :!!! ! !|||||||||!!.||| !!:!!||||||||||||||||||:!
------LysLysAspGlyThrValValGlyLeuHisPheValGlyProAsnAlaGlyGluIl
1802 : ------AAGAAGGACGGTACCGTCGTCGGCTTGCACTTCGTCGGCCCTAACGCCGGCGAAAT : 1858
449 : lThrGlnGlyPheAlaAlaAlaLeuLysCysGlyLeuThrLysLysGlnLeuAspSerThrI : 469
!! !||||||.!!..! !|||:!!::: |||:!!||||||:!!:!!:!!|||..!||||
eMetGlnGlyLeuCysMetAlaValArgLysGlyIleThrLysGluGluIleAspAspThrI
1859 : AATGCAGGGCCTCTGTATGGCTGTGCGCAAGGGAATCACCAAAGAGGAGATTGACGACACCA : 1921
470 : leGlyIleHisProValCysAlaGluValPheThrThrLeuSerValThrLysArgSerGly : 489
||||||||||||||..! !|||||| !|||..!!.!|||!!!|||||||||!:!||||||
leGlyIleHisProThrAspAlaGluSerPheValAsnLeuThrValThrLysLysSerGly
1922 : TCGGTATCCACCCGACTGACGCTGAGAGCTTCGTAAATCTCACTGTCACGAAGAAATCCGGC : 1981
490 : AlaSer : 491
! !|||
GluSer
1982 : GAGTCA : 1989
vulgar: SPP00000034_1.0 282 491 . Contig4075:subseq(18500,3000) 1295 1989 + 341 M 16 48 G 0 3 M 13 39
G 1 0 M 4 12 G 1 0 M 47 141 5 0 2 I 0 33 3 0 2 M 33 99 G 1 0 M 20 60 5 0 2 I 0 38 3 0 2 M 9 27 G 2 0
M 62 186
-- completed exonerate analysis
Command line: [exonerate -m p2g -q TR2.fa -t blast4075_18500-21500.fa]
Hostname: [am75250]
C4 Alignment:
------------
Query: SPP00000035_1.0 # Protein # Thioredoxin reductase 2 (TR2) # Homo sapiens # Complete
Target: Contig4075:subseq(18500,3000)
Model: protein2genome
Raw score: 660
Query range: 39 -> 311
Target range: 471 -> 1322
40 : TyrAspLeuLeuValValGlyGlyGlySerGlyGlyLeuAlaCysAlaLysGluAlaAlaGl : 60
||||||!!! !|||:!!|||||||||!!!||||||||||||..!|||||||||||||||..
TyrAspPheAlaValIleGlyGlyGlyCysGlyGlyLeuAlaAlaAlaLysGluAlaAlaSe
472 : TACGACTTCGCCGTCATTGGCGGAGGTTGCGGTGGATTGGCCGCCGCTAAGGAGGCTGCCAG : 532
61 : nLeuGlyArgLysValAlaValValAspTyrValGluProSerProGlnGlyThrArgTrpG : 81
.|||||| !|||..! !:!! !!|||||||||...|||||||||!::||||||! !||||
rLeuGlyAlaLysThrIleLeuPheAspTyrValArgProSerProArgGlyThrThrTrpG
533 : CTTGGGAGCGAAGACGATTCTATTTGATTATGTGCGCCCTTCTCCGCGTGGCACAACTTGGG : 595
82 : lyLeuGlyGlyThrCysValAsnValGlyCysIleProLysLysLeuMetHisGlnAlaAla : 101
||||||||||||||||||||||||||||||||||||||||||||||||||||| ! |||!.!
lyLeuGlyGlyThrCysValAsnValGlyCysIleProLysLysLeuMetHisTyrAlaGly
596 : GTCTCGGCGGTACCTGTGTGAATGTCGGCTGTATCCCGAAGAAGTTGATGCACTATGCTGGC : 655
102 : LeuLeuGlyGlyLeuIleGlnAspAlaProAsnTyrGlyTrpGluValAlaGlnProValPr : 122
:!!|||||| !! ! !!.||| ! ! !||||||! ! !:!! ! ||
IleLeuGlyHisSerSerHisAspArgGluAlaLeuGlyTrpGlyAsnHisGluGly---Pr
656 : ATTTTGGGCCACTCGTCGCATGATCGTGAGGCGCTTGGTTGGGGCAACCATGAGGGT---CC : 715
123 : oHisAspTrpArgLysMetAlaGlu >>>> Target Intron 1 >>>> AlaValGl : 133
|||||||||| !!::::!:!.!.!. 38 bp .!!|||||
oHisAspTrpGlyArgLeuValAsn++ ++ThrValGl
716 : GCACGATTGGGGCCGTCTTGTTAACgt.........................agACCGTGCA : 786
134 : nAsnHisValLysSerLeuAsnTrpGlyHisArgValGlnLeuGlnAspArgLysValLysT : 154
||||:!!:!!|||! !||||||!::..!:!!||| ! |||..!...!:!!!.|||:!!|
nAsnTyrIleLysMetLeuAsnPheSerTyrArgSerGlyLeuMetSerLysAsnValGluT
787 : GAACTACATCAAGATGCTCAACTTTTCCTACCGCTCGGGTCTCATGTCGAAAAACGTCGAAT : 849
155 : yrPheAsnIleLysAlaSerPheValAspGluHisThrValCysGlyValAlaLysGlyGly : 174
|| !!||| !! !||||||.!!!.!|||:!!:!!|||||| !..!! !||| !|||
yrValAsnAlaMetAlaSerLeuAlaAspLysAsnThrValThrTyrThrAspLysLysGly
850 : ACGTCAACGCCATGGCGAGCCTTGCCGATAAAAACACGGTGACATACACCGACAAGAAGGGG : 909
175 : LysGluIleLeuLeuSerAlaAspHisIleIleIleAlaThrGlyGlyArgProArgTyrPr : 195
:!!:!! !! !:!!!.!||| ! :!!:!::!:||||||! !|||!.!||||||! ! !||
GluLysHisGlnIleLysAlaLysAsnValLeuIleAlaIleGlyAlaArgProThrIlePr
910 : GAGAAGCATCAGATCAAGGCCAAGAATGTGCTGATCGCTATCGGAGCGCGTCCCACCATCCC : 972
196 : oThrHisIleGluGlyAlaLeuGluTyrGlyIleThrSerAspAspIlePheTrpLeuLysG : 216
|!:! !!:!::!!||||||! !!!:|||..!|||||||||||||||:!!.!. ! !||||
oSerAspValLysGlyAlaTrpAspTyrSerIleThrSerAspAspLeuMetSerArgLysG
973 : CAGCGATGTGAAGGGTGCATGGGACTACTCAATCACCTCTGACGATTTAATGAGCAGGAAGG : 1035
217 : luSerProGlyLysThrLeuValValGlyAlaSerTyrValAlaLeuGluCysAlaGlyPhe : 236
|| !! !||||||||||||:!!||||||!.!|||!:!||||||||||||||||||||||||
luProValGlyLysThrLeuIleValGlyGlySerPheValAlaLeuGluCysAlaGlyPhe
1036 : AGCCCGTCGGCAAGACTCTTATCGTGGGCGGTTCGTTCGTTGCTCTGGAGTGCGCCGGCTTC : 1095
237 : LeuThrGlyIleGlyLeuAspThrThrIleMetMetArgSerIleProLeuArgGlyPheAs : 257
||||||..!!!:|||! !|||..!|||:!! :!:||||||:!! !||||||||||||||
LeuThrSerMetGlyTyrAspValThrValAlaValArgSerLeuIleLeuArgGlyPheAs
1096 : CTCACCTCGATGGGTTACGACGTCACCGTCGCCGTCCGTTCTCTTATACTTAGAGGCTTCGA : 1158
258 : pGlnGlnMetSerSerMetValIleGluHisMetAlaSerHisGlyThrArgPheLeuArgG : 278
|!::||| :!!..!! !||| |||! !||| !:!! !||||||!:!||| !!.!|
pArgGlnCysAlaAspLysValGlnGluLeuMetLeuAlaThrGlyThrLysPheLysAsnG
1159 : CCGCCAGTGTGCGGACAAGGTCCAGGAGCTTATGTTGGCGACTGGCACCAAGTTCAAGAACG : 1221
279 : lyCysAlaProSerArgValArgArgLeuProAspGlyGlnLeuGlnValThrTrpGluAsp : 298
|| !!.!|||..! !:!!! !!:!||| !:!!|||::!||| ! :!! !!:: |||
lyValValProGlnAlaIleThrLysLeuGluAsnGlyArgLeuTyrIleGluPheThrAsp
1222 : GAGTTGTGCCGCAAGCCATCACCAAACTCGAGAATGGGAGACTATACATCGAATTCACCGAT : 1281
299 : SerThrThrGlyLysGluAspThrGlyThrPheAspThr : 311
.!!!:!:!!! !:!! |||||| !! !!:!! !|||
GlySerSerAspGluPheAspThrLeuMetTyrAlaThr
1282 : GGCAGCTCTGACGAATTTGACACGCTGATGTACGCTACC : 1322
vulgar: SPP00000035_1.0 39 311 . Contig4075:subseq(18500,3000) 471 1322 + 660 M 81 243 G 1 0 M 9 27 5
0 2 I 0 34 3 0 2 M 181 543
C4 Alignment:
------------
Query: SPP00000035_1.0 # Protein # Thioredoxin reductase 2 (TR2) # Homo sapiens # Complete
Target: Contig4075:subseq(18500,3000)
Model: protein2genome
Raw score: 153
Query range: 445 -> 514
Target range: 1782 -> 1983
446 : CysTyrValLysMetValCysLeuArgGluProProGlnLeuValLeuGlyLeuHisPheLe : 466
||| !||||||:!!||| !!:!!!: !|||:!!||||||||||||:!
CysLeuSerLysMetIleCysLysLysAsp------GlyThrValValGlyLeuHisPheVa
1783 : TGCCTGTCCAAAATGATCTGCAAGAAGGAC------GGTACCGTCGTCGGCTTGCACTTCGT : 1837
467 : uGlyProAsnAlaGlyGluValThrGlnGlyPheAlaLeuGlyIleLysCysGlyAlaSerT : 487
!||||||||||||||||||:!!! !||||||.!!..!:!!!.!:!:::: ||| !!!!
lGlyProAsnAlaGlyGluIleMetGlnGlyLeuCysMetAlaValArgLysGlyIleThrL
1838 : CGGCCCTAACGCCGGCGAAATAATGCAGGGCCTCTGTATGGCTGTGCGCAAGGGAATCACCA : 1900
488 : yrAlaGlnValMetArgThrValGlyIleHisProThrCysSerGluGluValValLysLeu : 507
! ! !:!!:!! !|||:!!||||||||||||||| !:!!|||... !!|||!!.|||
ysGluGluIleAspAspThrIleGlyIleHisProThrAspAlaGluSerPheValAsnLeu
1901 : AAGAGGAGATTGACGACACCATCGGTATCCACCCGACTGACGCTGAGAGCTTCGTAAATCTC : 1960
508 : ArgIleSerLysArgSerGly : 514
! !:!!!!!|||!:!||||||
ThrValThrLysLysSerGly
1961 : ACTGTCACGAAGAAATCCGGC : 1983
vulgar: SPP00000035_1.0 445 514 . Contig4075:subseq(18500,3000) 1782 1983 + 153 M 10 30 G 2 0 M 57 171
C4 Alignment:
------------
Query: SPP00000035_1.0 # Protein # Thioredoxin reductase 2 (TR2) # Homo sapiens # Complete
Target: Contig4075:subseq(18500,3000)
Model: protein2genome
Raw score: 132
Query range: 308 -> 386
Target range: 1298 -> 1529
309 : PheAspThrValLeuTrpAlaIleGlyArgValProAspThrArgSerLeuAsnLeuGluLy : 329
|||||||||:!!:!!!::|||! !|||||| ! !..!:!!||| :!!.!. !|||
PheAspThrLeuMetTyrAlaThrGlyArgSerValSerSerArg---MetGlnLysGluLe
1299 : TTTGACACGCTGATGTACGCTACCGGCCGTTCTGTGTCCTCAAGG---ATGCAAAAGGAGCT : 1356
330 : sAlaGlyValAspThrSerProAspThrGlnLysIleLeuValAspSerArgGluAlaThrS : 350
::!! !|||! !! !!.! !..! ! ! ! !:!!:!!! !..! !! !||||
uSerAspValGlyIleLysPheSerGluTyrGlyLysIleIleAlaGluAspGlyLysThrS
1357 : CAGCGACGTGGGCATCAAATTCAGCGAGTACGGCAAAATCATTGCCGAGGATGGCAAGACCT : 1419
351 : erValProHisIleTyrAlaIleGlyAspValValGluGlyArgProGluLeuThrProIle : 370
||||| ! !:!:||||||:!!||||||||||||||||||!.!|||! |||.!!|||:!:
erValGluGlyValTyrAlaValGlyAspValValGluGlyAsnProAlaLeuAlaProVal
1420 : CTGTGGAGGGCGTGTACGCCGTCGGTGACGTGGTGGAGGGCAACCCTGCTCTCGCACCGGTG : 1479
371 : AlaIleMetAlaGlyArgLeuLeuValGlnArgLeuPheGlyGlySer : 386
|||:!!! !! !|||..!||||||!.!!::|||:!!||||||..!|||
AlaValLysAspGlyGluLeuLeuAlaArgArgIlePheGlyAsnSer
1480 : GCCGTTAAGGACGGTGAACTGCTTGCCCGCAGGATCTTCGGCAACTCA : 1529
vulgar: SPP00000035_1.0 308 386 . Contig4075:subseq(18500,3000) 1298 1529 + 132 M 15 45 G 1 0 M 62 186
-- completed exonerate analysis
Command line: [exonerate -m p2g -q TR3.fa -t blast4075_18500-21500.fa]
Hostname: [am75250]
C4 Alignment:
------------
Query: SPP00000036_1.0 # Protein # Thioredoxin reductase 3 (TR3) # Homo sapiens # Complete
Target: Contig4075:subseq(18500,3000)
Model: protein2genome
Raw score: 584
Query range: 267 -> 517
Target range: 471 -> 1259
268 : TyrAspLeuIleIleIleGlyGlyGlySerGlyGlyLeuSerCysAlaLysGluAlaAlaIl : 288
||||||!!! !:!!||||||||||||!!!|||||||||:!!..!|||||||||||||||!
TyrAspPheAlaValIleGlyGlyGlyCysGlyGlyLeuAlaAlaAlaLysGluAlaAlaSe
472 : TACGACTTCGCCGTCATTGGCGGAGGTTGCGGTGGATTGGCCGCCGCTAAGGAGGCTGCCAG : 532
289 : eLeuGlyLysLysValMetValLeuAspPheValValProSerProGlnGlyThrSerTrpG : 309
!|||||| !|||..!!!::!!!!!|||!:!||| !|||||||||!::||||||!!!||||
rLeuGlyAlaLysThrIleLeuPheAspTyrValArgProSerProArgGlyThrThrTrpG
533 : CTTGGGAGCGAAGACGATTCTATTTGATTATGTGCGCCCTTCTCCGCGTGGCACAACTTGGG : 595
310 : lyLeuGlyGlyThrCysValAsnValGlyCysIleProLysLysLeuMetHisGlnAlaAla : 329
||||||||||||||||||||||||||||||||||||||||||||||||||||| ! |||!.!
lyLeuGlyGlyThrCysValAsnValGlyCysIleProLysLysLeuMetHisTyrAlaGly
596 : GTCTCGGCGGTACCTGTGTGAATGTCGGCTGTATCCCGAAGAAGTTGATGCACTATGCTGGC : 655
330 : LeuLeuGlyGlnAlaLeuCysAspSerArgLysPheGlyTrpGluTyrAsnGlnGlnValAr : 350
:!!||||||!!.:!!! ! !||| !!..! !.!!||||||! |||!!.:!!! !!
IleLeuGlyHisSerSerHisAspArgGluAlaLeuGlyTrpGly---AsnHisGluGlyPr
656 : ATTTTGGGCCACTCGTCGCATGATCGTGAGGCGCTTGGTTGGGGC---AACCATGAGGGTCC : 715
351 : gHisAsnTrpGluThrMetThrLys >>>> Target Intron 1 >>>> AlaIleGl : 361
!|||:!!|||! !:!:..!!!. 38 bp .!!:!:||
oHisAspTrpGlyArgLeuValAsn++ ++ThrValGl
716 : GCACGATTGGGGCCGTCTTGTTAACgt.........................agACCGTGCA : 786
362 : nAsnHisIleSerSerLeuAsnTrpGlyTyrArgLeuSerLeuArgGluLysAlaValAlaT : 382
||||:!!|||!.!! !||||||!::..!||||||! !.!!|||! !..!||| !|||! !|
nAsnTyrIleLysMetLeuAsnPheSerTyrArgSerGlyLeuMetSerLysAsnValGluT
787 : GAACTACATCAAGATGCTCAACTTTTCCTACCGCTCGGGTCTCATGTCGAAAAACGTCGAAT : 849
383 : yrValAsnSerTyrGlyGluPheValGluHisHisLysIleLysAlaThrAsnLysLysGly : 402
||||||||:!! !.!....!!!.!!!: ! :!!! !:!:! ! !|||:!!|||||||||
yrValAsnAlaMetAlaSerLeuAlaAspLysAsnThrValThrTyrThrAspLysLysGly
850 : ACGTCAACGCCATGGCGAGCCTTGCCGATAAAAACACGGTGACATACACCGACAAGAAGGGG : 909
403 : GlnGluThrTyrTyrThrAlaAlaGlnPheValIleAlaThrGlyGluArgProArgTyr<- : 423
:!!:!! ! ! !! !||| !.!. ! :!!||||||! !|||! !||||||! ! !
GluLysHisGlnIleLysAlaLysAsnValLeuIleAlaIleGlyAlaArgProThrIlePr
910 : GAGAAGCATCAGATCAAGGCCAAGAATGTGCTGATCGCTATCGGAGCGCGTCCCACCATCCC : 972
424 : >LeuGlyIleGlnGlyAspLysGluTyrCysIleThrSerAspAspLeuPheSerLeuProT : 443
!! !:!::!!|||! !!!:|||! ||||||||||||||||||.!.||| ! !
oSerAspValLysGlyAlaTrpAspTyrSerIleThrSerAspAspLeuMetSerArgLysG
973 : CAGCGATGTGAAGGGTGCATGGGACTACTCAATCACCTCTGACGATTTAATGAGCAGGAAGG : 1035
444 : yrCysProGlyLysThrLeuValValGlyAlaSerTyrValAlaLeuGluCysAlaGlyPhe : 463
! ! !||||||||||||:!!||||||!.!|||!:!||||||||||||||||||||||||
luProValGlyLysThrLeuIleValGlyGlySerPheValAlaLeuGluCysAlaGlyPhe
1036 : AGCCCGTCGGCAAGACTCTTATCGTGGGCGGTTCGTTCGTTGCTCTGGAGTGCGCCGGCTTC : 1095
464 : LeuAlaGlyPheGlyLeuAspValThrValMetValArgSerIleLeuLeuArgGlyPheAs : 484
|||.!!..!.!.|||! !|||||||||||| |||||||||:!!:!!||||||||||||||
LeuThrSerMetGlyTyrAspValThrValAlaValArgSerLeuIleLeuArgGlyPheAs
1096 : CTCACCTCGATGGGTTACGACGTCACCGTCGCCGTCCGTTCTCTTATACTTAGAGGCTTCGA : 1158
485 : pGlnGluMetAlaGluLysValGlySerTyrMetGluGlnHisGlyValLysPheLeuArgL : 505
|!:::!! |||!!:|||||| !..! !||| ! ! !|||..!|||||| !!.!
pArgGlnCysAlaAspLysValGlnGluLeuMetLeuAlaThrGlyThrLysPheLysAsnG
1159 : CCGCCAGTGTGCGGACAAGGTCCAGGAGCTTATGTTGGCGACTGGCACCAAGTTCAAGAACG : 1221
506 : ysPheIleProValMetValGlnGlnLeuGluLysGly : 517
! !!:!:||| ! :!! :!!||||||!!.|||
lyValValProGlnAlaIleThrLysLeuGluAsnGly
1222 : GAGTTGTGCCGCAAGCCATCACCAAACTCGAGAATGGG : 1259
vulgar: SPP00000036_1.0 267 517 . Contig4075:subseq(18500,3000) 471 1259 + 584 M 77 231 G 1 0
M 13 39 5 0 2 I 0 34 3 0 2 M 64 192 G 0 3 M 95 285
C4 Alignment:
------------
Query: SPP00000036_1.0 # Protein # Thioredoxin reductase 3 (TR3) # Homo sapiens # Complete
Target: Contig4075:subseq(18500,3000)
Model: protein2genome
Raw score: 345
Query range: 529 -> 744
Target range: 1271 -> 1983
530 : GluGlyThrGluThrIleGluGlyValTyrAsnThrValLeuLeuAlaIleGlyArgAspSe : 550
||| !|||!!: !! !...! !! !!:!:!!|||:!!:!!! !|||! !||||||..!
GluPheThrAspGlySerSerAspGluPheAspThrLeuMetTyrAlaThrGlyArgSerVa
1272 : GAATTCACCGATGGCAGCTCTGACGAATTTGACACGCTGATGTACGCTACCGGCCGTTCTGT : 1332
551 : rCysThrArg<->LysIleGlyLeuGluLysIleGlyValLysIleAsnGluLysSerGlyL : 570
!! !:!!||| :!!! ! !|||... ! :!:|||:!!|||.!!!:!||| ! !||||
lSerSerArgMetGlnLysGluLeuSerAspValGlyIleLysPheSerGlu---TyrGlyL
1333 : GTCCTCAAGGATGCAAAAGGAGCTCAGCGACGTGGGCATCAAATTCAGCGAG---TACGGCA : 1392
571 : ysIleProValAsnAspValGluGlnThrAsnValProTyrValTyrAlaValGlyAspIle : 590
||||| !!.!.!.||| ! :!!|||::!||| ! !||||||||||||||||||:!:
ysIleIleAlaGluAsp---GlyLysThrSerValGluGlyValTyrAlaValGlyAspVal
1393 : AAATCATTGCCGAGGAT---GGCAAGACCTCTGTGGAGGGCGTGTACGCCGTCGGTGACGTG : 1449
591 : LeuGluAspLysProGluLeuThrProValAlaIleGlnSerGlyLysLeuLeuAlaGlnAr : 611
:!!|||! !!!.|||! |||.!!|||||||||:!!:!!..!|||:!!|||||||||!::||
ValGluGlyAsnProAlaLeuAlaProValAlaValLysAspGlyGluLeuLeuAlaArgAr
1450 : GTGGAGGGCAACCCTGCTCTCGCACCGGTGGCCGTTAAGGACGGTGAACTGCTTGCCCGCAG : 1512
612 : gLeuPheGlyAlaSerLeuGluLys >>>> Target Intron 1 >>>> CysAspTy : 622
|:!!|||||| !||| !:!!||| 37 bp !|||
gIlePheGlyAsnSerAspLysLys++ ++LeuAspLe
1513 : GATCTTCGGCAACTCAGACAAGAAGgt.........................agCTTGATCT : 1582
623 : rIleAsnValProThrThrValPheThrProLeuGluTyrGlyCysCysGlyLeuSerGluG : 643
! ! !!:!!|||! ! !||||||||||||! !||||||!.! ! ||||||:!!|||||||
uAsnTyrIleProMetCysValPheThrProTyrGluTyrAlaArgCysGlyIleSerGluG
1583 : GAACTACATTCCCATGTGTGTCTTCACGCCGTACGAATACGCGCGGTGTGGAATTTCGGAGG : 1645
644 : luLysAlaIleGluValTyrLysLysGluAsnLeuGluIleTyrHisThrLeuPheTrpPro : 663
||! !||| !:!!:!!||| !|||:!!!!!!!::!!|||! ! ! !!:! ! !!
luThrAlaSerLysLeuTyr---GlyGluAspPheAspValTyrLeuLysGluTyrThrThr
1646 : AGACAGCCTCCAAATTGTAC---GGGGAAGACTTTGATGTCTACCTGAAGGAGTACACGACT : 1702
664 : LeuGluTrpThrValAlaGlyArgGluAsn >>>> Target Intron 2 >>>> Asn : 674
||||||!::.!!!.!!.! !||||||!!. 42 bp !
LeuGluPheAlaAlaValHisArgGluLys++ -+Pro
1703 : CTGGAGTTTGCGGCTGTTCACCGCGAGAAGgt.........................cgCCT : 1777
675 : ThrCysTyrAlaLysIleIleCysAsnLysPheAspHisAspArgValIleGlyPheHisIl : 695
|||||| :!!|||!!:||||||!!.||| !!! !! !|||:!!|||!!.|||.!
ThrCysLeuSerLysMetIleCysLysLys------AspGlyThrValValGlyLeuHisPh
1778 : ACTTGCCTGTCCAAAATGATCTGCAAGAAG------GACGGTACCGTCGTCGGCTTGCACTT : 1834
696 : eLeuGlyProAsnAlaGlyGluValThrGlnGlyPheAlaAlaAlaMetLysCysGlyLeuT : 716
!:!!||||||||||||||||||:!!! !||||||.!!..! !|||:!!::: |||:!!|
eValGlyProAsnAlaGlyGluIleMetGlnGlyLeuCysMetAlaValArgLysGlyIleT
1835 : CGTCGGCCCTAACGCCGGCGAAATAATGCAGGGCCTCTGTATGGCTGTGCGCAAGGGAATCA : 1897
717 : hrLysGlnLeuLeuAspAspThrIleGlyIleHisProThrCysGlyGluValPheThrThr : 736
|||||:!! !:!!||||||||||||||||||||||||||| !!.!||| !|||..!!.!
hrLysGluGluIleAspAspThrIleGlyIleHisProThrAspAlaGluSerPheValAsn
1898 : CCAAAGAGGAGATTGACGACACCATCGGTATCCACCCGACTGACGCTGAGAGCTTCGTAAAT : 1957
737 : LeuGluIleThrLysSerSerGly : 744
||| :!!||||||!.!||||||
LeuThrValThrLysLysSerGly
1958 : CTCACTGTCACGAAGAAATCCGGC : 1983
vulgar: SPP00000036_1.0 529 744 . Contig4075:subseq(18500,3000) 1271 1983 + 345 M 24 72 G 0 3 M 13 39
G 1 0 M 8 24 G 1 0 M 43 129 5 0 2 I 0 33 3 0 2 M 30 90 G 1 0 M 23 69 5 0 2 I 0 38 3 0 2 M 11 33 G 2 0
M 58 174
-- completed exonerate analysis