# hmmscan :: search sequence(s) against a profile database
# HMMER 3.0 (March 2010); http://hmmer.org/
# Copyright (C) 2010 Howard Hughes Medical Institute.
# Freely distributed under the GNU General Public License (GPLv3).
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
# query sequence file:             H1.5.fa
# target HMM database:             /N/dc2/projects/marinovg/genomes/PFam/PFam27.0/Pfam-A.hmm
# output directed to file:         H1.5.PFam-27-A
# per-seq hits tabular output:     H1.5.hmmscan.tblout
# number of worker threads:        16
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

Query:       histone_H1.5_[Homo_sapiens]gi|4885381|ref|NP_005313.1  [L=226]
Scores for complete sequence (score includes all domains):
   --- full sequence ---   --- best 1 domain ---    -#dom-
    E-value  score  bias    E-value  score  bias    exp  N  Model          Description
    ------- ------ -----    ------- ------ -----   ---- --  --------       -----------
    1.1e-22   79.8   1.8    1.1e-22   79.8   1.3    2.1  2  Linker_histone linker histone H1 and H5 family


Domain annotation for each model (and alignments):
>> Linker_histone  linker histone H1 and H5 family
   #    score  bias  c-Evalue  i-Evalue hmmfrom  hmm to    alifrom  ali to    envfrom  env to     acc
 ---   ------ ----- --------- --------- ------- -------    ------- -------    ------- -------    ----
   1 !   79.8   1.3   7.6e-27   1.1e-22       3      77 .]      41     112 ..      39     112 .. 0.96
   2 ?   -3.3   1.5      0.64   9.4e+03      34      39 ..     149     154 ..     123     185 .. 0.58

  Alignments for each domain:
  == domain 1    score: 79.8 bits;  conditional E-value: 7.6e-27
                                                            --HHHHHHHHCCCCCCCTSEEHHHHHHHHHC.CTCCCCCSTTHHHHHHHHHHHHHH CS
                                         Linker_histone   3 ppyseMiveaikalkerkGsSvqaikKyiea.nykveinasnlasllklaLkkgvk 57 
                                                            pp+se+i++a++a+ker+G S +a+kK ++a  y+ve+n+    s++kl Lk++v+
  histone_H1.5_[Homo_sapiens]gi|4885381|ref|NP_005313.1  41 PPVSELITKAVAASKERNGLSLAALKKALAAgGYDVEKNN----SRIKLGLKSLVS 92 
                                                            9******************************889988888....************ PP

                                                            CTSEE-TT..TTTSEEEESS CS
                                         Linker_histone  58 sGeLvqvkgtGasGsfklak 77 
                                                            +G+Lvq+kgtGasGsfkl+k
  histone_H1.5_[Homo_sapiens]gi|4885381|ref|NP_005313.1  93 KGTLVQTKGTGASGSFKLNK 112
                                                            ******************97 PP

  == domain 2    score: -3.3 bits;  conditional E-value: 0.64
                                                            CTCCCC CS
                                         Linker_histone  34 nykvei 39 
                                                            +  v++
  histone_H1.5_[Homo_sapiens]gi|4885381|ref|NP_005313.1 149 KKAVKK 154
                                                            322222 PP



Internal pipeline statistics summary:
-------------------------------------
Query sequence(s):                         1  (226 residues)
Target model(s):                       14831  (2610332 nodes)
Passed MSV filter:                      1050  (0.0707977); expected 296.6 (0.02)
Passed bias filter:                      290  (0.0195536); expected 296.6 (0.02)
Passed Vit filter:                        59  (0.00397815); expected 14.8 (0.001)
Passed Fwd filter:                        24  (0.00161823); expected 0.1 (1e-05)
Initial search space (Z):              14831  [actual number of targets]
Domain search space  (domZ):               1  [number of targets reported over threshold]
# CPU time: 0.29u 0.45s 00:00:00.74 Elapsed: 00:00:01.94
# Mc/sec: 304.09
//