1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
|
# LAST version 1243
#
# a=11 b=2 A=11 B=2 F=12 e=40 d=29 x=39 y=31 z=39 D=1e+06 E=8.45245e+08
# R=01 u=0 s=2 S=0 M=0 T=0 m=10 l=1 n=10 k=1 w=1000 t=3.08611 j=3 Q=0
# /home/mcfrith/tmp/x
# Reference sequences=1 normal letters=491
# lambda=0.327048 K=0.105437
#
# A R N D C Q E G H I L K M F P S T W Y V B J Z X *
# A 4 -1 -2 -2 0 -1 -1 0 -2 -1 -1 -1 -1 -2 -1 1 0 -3 -2 0 -2 -1 -1 -1 -4
# R -1 5 0 -2 -3 1 0 -2 0 -3 -2 2 -1 -3 -2 -1 -1 -3 -2 -3 -1 -2 0 -1 -4
# N -2 0 6 1 -3 0 0 0 1 -3 -3 0 -2 -3 -2 1 0 -4 -2 -3 4 -3 0 -1 -4
# D -2 -2 1 6 -3 0 2 -1 -1 -3 -4 -1 -3 -3 -1 0 -1 -4 -3 -3 4 -3 1 -1 -4
# C 0 -3 -3 -3 9 -3 -4 -3 -3 -1 -1 -3 -1 -2 -3 -1 -1 -2 -2 -1 -3 -1 -3 -1 -4
# Q -1 1 0 0 -3 5 2 -2 0 -3 -2 1 0 -3 -1 0 -1 -2 -1 -2 0 -2 4 -1 -4
# E -1 0 0 2 -4 2 5 -2 0 -3 -3 1 -2 -3 -1 0 -1 -3 -2 -2 1 -3 4 -1 -4
# G 0 -2 0 -1 -3 -2 -2 6 -2 -4 -4 -2 -3 -3 -2 0 -2 -2 -3 -3 -1 -4 -2 -1 -4
# H -2 0 1 -1 -3 0 0 -2 8 -3 -3 -1 -2 -1 -2 -1 -2 -2 2 -3 0 -3 0 -1 -4
# I -1 -3 -3 -3 -1 -3 -3 -4 -3 4 2 -3 1 0 -3 -2 -1 -3 -1 3 -3 3 -3 -1 -4
# L -1 -2 -3 -4 -1 -2 -3 -4 -3 2 4 -2 2 0 -3 -2 -1 -2 -1 1 -4 3 -3 -1 -4
# K -1 2 0 -1 -3 1 1 -2 -1 -3 -2 5 -1 -3 -1 0 -1 -3 -2 -2 0 -3 1 -1 -4
# M -1 -1 -2 -3 -1 0 -2 -3 -2 1 2 -1 5 0 -2 -1 -1 -1 -1 1 -3 2 -1 -1 -4
# F -2 -3 -3 -3 -2 -3 -3 -3 -1 0 0 -3 0 6 -4 -2 -2 1 3 -1 -3 0 -3 -1 -4
# P -1 -2 -2 -1 -3 -1 -1 -2 -2 -3 -3 -1 -2 -4 7 -1 -1 -4 -3 -2 -2 -3 -1 -1 -4
# S 1 -1 1 0 -1 0 0 0 -1 -2 -2 0 -1 -2 -1 4 1 -3 -2 -2 0 -2 0 -1 -4
# T 0 -1 0 -1 -1 -1 -1 -2 -2 -1 -1 -1 -1 -2 -1 1 5 -2 -2 0 -1 -1 -1 -1 -4
# W -3 -3 -4 -4 -2 -2 -3 -2 -2 -3 -2 -3 -1 1 -4 -3 -2 11 2 -3 -4 -2 -2 -1 -4
# Y -2 -2 -2 -3 -2 -1 -2 -3 2 -1 -1 -2 -1 3 -3 -2 -2 2 7 -1 -3 -1 -2 -1 -4
# V 0 -3 -3 -3 -1 -2 -2 -3 -3 3 1 -2 1 -1 -2 -2 0 -3 -1 4 -3 2 -2 -1 -4
# B -2 -1 4 4 -3 0 1 -1 0 -3 -4 0 -3 -3 -2 0 -1 -4 -3 -3 4 -3 0 -1 -4
# J -1 -2 -3 -3 -1 -2 -3 -4 -3 3 3 -3 2 0 -3 -2 -1 -2 -1 2 -3 3 -3 -1 -4
# Z -1 0 0 1 -3 4 4 -2 0 -3 -3 1 -1 -3 -1 0 -1 -2 -2 -2 0 -3 4 -1 -4
# X -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -4
# * -4 -4 -4 -4 -4 -4 -4 -4 -4 -4 -4 -4 -4 -4 -4 -4 -4 -4 -4 -4 -4 -4 -4 -4 1
#
# Coordinates are 0-based. For - strand matches, coordinates
# in the reverse complement of the 2nd sequence are used.
#
# name start alnSize strand seqSize alignment
#
# batch 0
a score=118 EG2=1.8 E=2.4e-11
s Q2LCP8 110 232 + 491 LILIYSSIIGMLISMEAHNLITLFLSLEISSICFYILALNKNSRKGIEGGLKYYIIGGIATTILLLGIVSIYKSTGSLMYTDLLVIVMENTDDYRIQMGIALIVLGLILKLGIAPFHGWLIDVYEGAGMLMTFYLTITQKlvtlmvlinlyqnliiytnaimftngliILILVTLVVGTIGSLRQQKLIRFIAYSAIVNSALLILMLAGSMSEELMVNSVYYLINYIIGLTV
s chrM 5258 639 + 16775 LICTVSLIMGTSITISSNH*ILA*TGLEINTLAIIPLISKSHHPRAIEATIKYFLTQSTASALILFSSITNA*STGQ*DITQL------NHPTSCLILTIAIAI-----KLGLVPFHF*FPEVLQGSSLITALLLSTLIKLPPITLLLLTSQSL----Nttllt---llaISSTLIGG*MG-LNQTQTRKILAFSSISHLG*IIIIISYNPQLTILTFILYTIMTSTVFLSL
a score=72 EG2=6.3e+06 E=9e-05
s Q2LCP8 379 52 + 491 ILVLMYLAGLPPMTNFISKILILLPYMVLGRVYLTMLAFFLSVGVMIYYMNL
s chrM 5972 156 + 16775 ILTLLSLAGLPPLTGFMPK*LIIQELTKQEITPiatiitilsllslfFYLRL
a score=55 EG2=1.6e+09 E=0.024
s Q2LCP8 101 87 + 491 KKTKITDEGLILIYSSIIGMLISMEAHNLITLFLSLEISSICFYILALN-KNSRKGIEGGLKYYIIGGIATTILLLGIVSIYKSTGSL
s chrM 11752 255 + 16775 KRKRIFISTLIIIQPFII---LAFSATELILFYISFEATLIPTLILITR*GNQPERLSAGIYLLFYTLISSLPLLVSILYLHTNTGTL
a score=51 EG2=6e+09 E=0.09
s Q2LCP8 170 49 + 491 TTILLLGIVS-IYKSTGSLMYTDLLVIVMENTDDYRIQMGIALIVLGLIL
s chrM 13892 141 + 16775 TTCLCLGALSTLFAATCALTQNDIKKIIAFSTSS---QLGLIIVTIGLDL
a score=48 EG2=1.6e+10 E=0.24
s Q2LCP8 394 27 + 491 FISKILILLPYMVLGRVYLTMLAFFLS
s chrM 11767 81 + 16775 FISTLIIIQPFIILAFSATELILFYIS
a score=46 EG2=3.1e+10 E=0.47
s Q2LCP8 16 38 + 491 YILPLII-LIVLSISIKKDVQRVHIIIQSLKLtlivimi
s chrM 14080 115 + 16775 yiipmlr\LIIHSLNGEQDIRKIGCLQKTLPITTSCLTI
a score=44 EG2=5.9e+10 E=0.91
s Q2LCP8 27 57 + 491 SISIK-KDVQRVHIIIQSLKLtlivimiviGIEESIYVKLNGHLIKTELIKLFEYLLL
s chrM 13941 170 + 16775 ALSLK/NDIKKIIAFSTSSQLGLIIVTIGLDLPQLAFLHISTHAFFKAILFLCSGLII
a score=42 EG2=1.1e+11 E=1.8
s Q2LCP8 166 24 + 491 GGIATTILLLGIVSIYKSTGSLMY
s chrM 54 72 + 16775 GGILCIIVHTFIYHIYYGTGNIYY
a score=42 EG2=1.1e+11 E=1.8
s Q2LCP8 292 22 + 491 LRQQKLIRFIAYSAIVNSALLI
s chrM 12313 66 + 16775 LRQTDLKSLIAYSSVSHMGLVI
a score=42 EG2=1.1e+11 E=1.8
s Q2LCP8 388 15 + 491 LPPMTNFISKILILL
s chrM 12589 45 + 16775 LPPTTNLIAELTIIV
a score=41 EG2=1.6e+11 E=2.4
s Q2LCP8 10 66 + 491 MIKYSIYILPLIILIVLSISIKKDVQRVHIIIQSLKLtlivimiviGIEESI--YVKLNGHLIKTELI
s chrM 11677 204 + 16775 LLVLSCWFLPLIIIASQGHLQHEPHKRKRIFISTLIIIQPFIILAFSATELILFYISFEATLIPTLIL
a score=41 EG2=1.6e+11 E=2.4
s Q2LCP8 15 44 + 491 IYILPLIILIVLSISIKKDVQRVHIIIQSLKLtlivimiviGIE
s chrM 8444 129 + 16775 LYLLTLILIEKLS-SNTVDAQEVELI*TILPAIVlvllalpslq
a score=41 EG2=1.6e+11 E=2.4
s Q2LCP8 117 200 + 491 IIGMLISMEAHNLITLFLSLE-ISSICFYILALNKNSRKGIEGGLKYYIIGGIATTILLLGIVSIYKSTGSLMYTDLLVIVMENTDDYRIQMGIALIVLGLILKLGIAPFHGWLIDVYEGAGMLMTFYLTITQKlvtlmvlinlyqnliiytnaimftngliILILVTLVVGTIGSLRQQKLIRFIAYSAIVNSALLILML
s chrM 13445 576 + 16775 liailTLTIANNIFLLFVG*EGVGIMSFLLIG**QGRAEANTAALQAIIYNRIGDIGLILSIA*L---ASSLNT*EIQQITHPNQTPTLPLLGLILAATGKSAQFGL---HP*LPAAIEGPTPVSALLHSSTIVVAGIFLLIRTHP--FLSSNKTALTTCLCLGALSTLFAATC-ALTQNDIKKIIAFSTSSQLGLIIVTI
a score=41 EG2=1.6e+11 E=2.4
s Q2LCP8 375 42 + 491 AIVYILVLMYLAGLPPMTNFISKILILLPYMVLGRVYLTMLA
s chrM 5606 108 + 16775 SLITALLLSTLIKLPPIT------LLLLTSQSLNttlltlla
a score=40 EG2=2.2e+11 E=3.4
s Q2LCP8 374 30 + 491 GAIVYILVLMYLAGLPPM---TNFISKILILLP
s chrM 7064 99 + 16775 gasvdlaifhylagVSSILGAINFITTIINIKP
a score=43 EG2=8.2e+10 E=1.3
s Q2LCP8 202 20 + 491 DYRIQMGIALIVLGLILKLG
s chrM 9413 57 - 16775 DRRIrigvg-iglglllqlG
a score=41 EG2=1.6e+11 E=2.4
s Q2LCP8 0 34 + 491 MVCLFENLMNMIKYSIYILPLIILIVLSISIKKD
s chrM 196 102 - 16775 LVSLGVSFVSLALFLVYLGGMLVVFVYSVSLAAD
a score=40 EG2=2.2e+11 E=3.4
s Q2LCP8 406 20 + 491 VLGRVYLTMLAFFLSVGVMI
s chrM 2209 60 - 16775 IIPrvvtmmaavfvmvgivi
# batch 1
# Query sequences=2
|
