# hmmscan :: search sequence(s) against a profile database
# HMMER 3.2 (June 2018); http://hmmer.org/
# Copyright (C) 2018 Howard Hughes Medical Institute.
# Freely distributed under the BSD open source license.
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
# query sequence file:             /tmp/tmpp4O0Ju/files/000/dataset_20.dat
# target HMM database:             localref.hmm
# per-seq hits tabular output:     /tmp/tmpp4O0Ju/files/000/dataset_22.dat
# per-dom hits tabular output:     /tmp/tmpp4O0Ju/files/000/dataset_23.dat
# pfam-style tabular hit output:   /tmp/tmpp4O0Ju/files/000/dataset_24.dat
# max ASCII text line length:      unlimited
# Vit filter P threshold:       <= 0.001
# Fwd filter P threshold:       <= 1e-05
# random number seed set to:       4
# number of worker threads:        1
# - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -

Query:       humanchr1/239220001-239550000  [L=330000]
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
    ------- ------ -----    ------- ------ -----   ---- --  -------- -----------
    5.7e-18   51.8  27.9      2e-12   34.1   0.7    9.6  5  MADE1     MADE1 (MAriner Derived Element 1), a TcMar-Mariner DNA transposon


Domain annotation for each model (and alignments):
>> MADE1  MADE1 (MAriner Derived Element 1), a TcMar-Mariner DNA transposon
   #    score  bias  c-Evalue  i-Evalue hmmfrom  hmm to    alifrom  ali to    envfrom  env to     acc
 ---   ------ ----- --------- --------- ------- -------    ------- -------    ------- -------    ----
   1 ?   -4.5   0.0         1         1      30      54 ..   80044   80068 ..   80033   80072 .. 0.81
   2 ?   -6.6   3.3         1         1      13      71 ..  154012  154072 ..  154011  154076 .. 0.75
   3 !   27.4   0.7   2.4e-10   2.4e-10       1      43 [.  174456  174498 ..  174456  174577 .. 0.66
   4 !   34.1   0.7     2e-12     2e-12       2      80 .]  302388  302466 ..  302387  302466 .. 0.86
   5 ?    2.8   0.7     0.011     0.011      27      75 ..  304060  304107 ..  304022  304109 .. 0.62

  Alignments for each domain:
  == domain 1  score: -4.5 bits;  conditional E-value: 1
                                      xxxxxxxxxxxxxxxxxxxxxxxxx RF
                          MADE1    30 ttgccattacttttaatggcaaaaa 54   
                                      t g catt  ttt aatggcaaa a
  humanchr1/239220001-239550000 80044 TAGTCATTCATTTCAATGGCAAATA 80068
                                      457899***************9865 PP

  == domain 2  score: -6.6 bits;  conditional E-value: 1
                                       xxxxxxxxxxxxxxxxxxxxxxxxx......xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx RF
                          MADE1     13 aaaagtaattgcggtttttgccatt......acttttaatggcaaaaaccgcaattacttttgca 71    
                                       aaaagta tt +   ttttgc att      a  tttaa  gcaaa a +    tta  tttgca
  humanchr1/239220001-239550000 154012 AAAAGTAGTTTTCAATTTTGCAATTtgaccaATATTTAAATGCAAATATT----TTATATTTGCA 154072
                                       78999999999999999999999984444444457777777899988765....67777778776 PP

  == domain 3  score: 27.4 bits;  conditional E-value: 2.4e-10
                                       xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx RF
                          MADE1      1 ttaggttggtgcaaaagtaattgcggtttttgccattactttt 43    
                                       ttaggtt gtgcaaaagtaattg+ggtttttg cattactttt
  humanchr1/239220001-239550000 174456 TTAGGTTAGTGCAAAAGTAATTGTGGTTTTTGTCATTACTTTT 174498
                                       79**************************************996 PP

  == domain 4  score: 34.1 bits;  conditional E-value: 2e-12
                                       xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx...xxxxxxxxxxxxxxxxxxxxxxxxxxxxxx RF
                          MADE1      2 taggttggtgcaaaagtaattgcggtttttgccattacttttaatggca...aaaaccgcaattacttttgcaccaacctaa 80    
                                       t ggt ggtgcaaaa  aattg+ggtttttgccatt cttttaat gc    aaaa  g  a t ctttt caccaa ctaa
  humanchr1/239220001-239550000 302388 TTGGTCGGTGCAAAATCAATTGTGGTTTTTGCCATTGCTTTTAATTGCTtttAAAA--G-TAATGCTTTTACACCAATCTAA 302466
                                       56899******************************************962223333..3.45578**************997 PP

  == domain 5  score: 2.8 bits;  conditional E-value: 0.011
                                       xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx RF
                          MADE1     27 tttttgccattacttttaatggcaaaaaccgcaattacttttgcaccaa 75    
                                       tttt g c  ta  tt a tgg aaaaa ++ca tta ttttgca  aa
  humanchr1/239220001-239550000 304060 TTTTAGACTATA-GTTAAGTGGGAAAAAATACACTTATTTTTGCATTAA 304107
                                       222223222222.3556779************************98765 PP



Internal pipeline statistics summary:
-------------------------------------
Query sequence(s):                         1  (330000 residues searched)
Target model(s):                           1  (80 nodes)
Passed MSV filter:                         1  (1); expected 0.0 (0.02)
Passed bias filter:                        1  (1); expected 0.0 (0.02)
Passed Vit filter:                         1  (1); expected 0.0 (0.001)
Passed Fwd filter:                         1  (1); expected 0.0 (1e-05)
Initial search space (Z):                  1  [actual number of targets]
Domain search space  (domZ):               1  [number of targets reported over threshold]
# CPU time: 0.14u 0.01s 00:00:00.15 Elapsed: 00:00:00.14
# Mc/sec: 177.58
//
[ok]