# 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.3.fa
# target HMM database:             /N/dc2/projects/marinovg/genomes/PFam/PFam27.0/Pfam-A.hmm
# output directed to file:         H1.3.PFam-27-A
# per-seq hits tabular output:     H1.3.hmmscan.tblout
# number of worker threads:        16
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

Query:       histone_H1.3_[Homo_sapiens]gi|4885377|ref|NP_005311.1  [L=221]
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
    ------- ------ -----    ------- ------ -----   ---- --  --------       -----------
      3e-22   78.5   2.3      3e-22   78.5   1.6    1.8  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 !   78.5   1.6     2e-26     3e-22       3      77 .]      39     110 ..      37     110 .. 0.95
   2 ?   -2.6   0.2      0.39   5.7e+03      39      39 ..     160     160 ..     136     184 .. 0.59

  Alignments for each domain:
  == domain 1    score: 78.5 bits;  conditional E-value: 2e-26
                                                            --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.3_[Homo_sapiens]gi|4885377|ref|NP_005311.1  39 PPVSELITKAVAASKERSGVSLAALKKALAAaGYDVEKNN----SRIKLGLKSLVS 90 
                                                            9****************************99689988888....************ PP

                                                            CTSEE-TT..TTTSEEEESS CS
                                         Linker_histone  58 sGeLvqvkgtGasGsfklak 77 
                                                            +G+Lvq+kgtGasGsfkl+k
  histone_H1.3_[Homo_sapiens]gi|4885377|ref|NP_005311.1  91 KGTLVQTKGTGASGSFKLNK 110
                                                            ******************97 PP

  == domain 2    score: -2.6 bits;  conditional E-value: 0.39
                                                            C CS
                                         Linker_histone  39 i 39 
                                                            +
  histone_H1.3_[Homo_sapiens]gi|4885377|ref|NP_005311.1 160 K 160
                                                            2 PP



Internal pipeline statistics summary:
-------------------------------------
Query sequence(s):                         1  (221 residues)
Target model(s):                       14831  (2610332 nodes)
Passed MSV filter:                       886  (0.0597397); expected 296.6 (0.02)
Passed bias filter:                      307  (0.0206999); expected 296.6 (0.02)
Passed Vit filter:                        45  (0.00303419); expected 14.8 (0.001)
Passed Fwd filter:                         7  (0.000471984); 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.21u 0.52s 00:00:00.73 Elapsed: 00:00:01.99
# Mc/sec: 289.89
//