# nhmmscan :: search DNA sequence(s) against a DNA profile database # HMMER 3.3 (Nov 2019); http://hmmer.org/ # Copyright (C) 2019 Howard Hughes Medical Institute. # Freely distributed under the BSD open source license. # - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - # query sequence file: /tmp/tmpqydies2m/files/f/a/5/dataset_fa559c91-0436-48cb-a47d-20062d4af284.dat # target HMM database: localref.hmm # per-seq hits tabular output: /tmp/tmpqydies2m/files/4/7/2/dataset_4723128d-8711-4edb-8afe-a68476a6c127.dat # hits output in Dfam format: /tmp/tmpqydies2m/files/1/a/e/dataset_1aeea50b-1263-4750-8237-4fc3f10c0999.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 hit: E-value score bias Model start end Description ------- ------ ----- -------- ----- ----- ----------- 4e-11 41.3 7.5 MADE1 302390 302466 MADE1 (MAriner Derived Element 1), a TcMar-Mariner DNA transposon 1.9e-08 32.8 8.3 MADE1 174456 174498 MADE1 (MAriner Derived Element 1), a TcMar-Mariner DNA transposon 6.3e-08 31.0 6.7 MADE1 302466 302389 MADE1 (MAriner Derived Element 1), a TcMar-Mariner DNA transposon 4.9e-06 25.0 7.0 MADE1 174493 174456 MADE1 (MAriner Derived Element 1), a TcMar-Mariner DNA transposon ------ inclusion threshold ------ 2.2 6.9 7.2 MADE1 304073 304103 MADE1 (MAriner Derived Element 1), a TcMar-Mariner DNA transposon Annotation for each hit (and alignments): >> MADE1 MADE1 (MAriner Derived Element 1), a TcMar-Mariner DNA transposon score bias Evalue hmmfrom hmm to alifrom ali to envfrom env to mod len acc ------ ----- --------- ------- ------- --------- --------- --------- --------- --------- ---- ! 41.3 7.5 4e-11 4 80 .] 302390 302466 .. 302387 302466 .. 80 0.88 Alignment: score: 41.3 bits xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx....xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx RF MADE1 4 ggttggtgcaaaagtaattgcggtttttgccattacttttaatggc....aaaaaccgcaattacttttgcaccaacctaa 80 ggt ggtgcaaaa aattg ggtttttgccatt cttttaat gc aaaa g a t ctttt caccaa ctaa humanchr1/239220001-239550000 302390 GGTCGGTGCAAAATCAATTGTGGTTTTTGCCATTGCTTTTAATTGCttttAAAA---GT-AATGCTTTTACACCAATCTAA 302466 89*******************************************966644554...34.4578**************997 PP >> MADE1 MADE1 (MAriner Derived Element 1), a TcMar-Mariner DNA transposon score bias Evalue hmmfrom hmm to alifrom ali to envfrom env to mod len acc ------ ----- --------- ------- ------- --------- --------- --------- --------- --------- ---- ! 32.8 8.3 1.9e-08 1 43 [. 174456 174498 .. 174456 174518 .. 80 0.91 Alignment: score: 32.8 bits xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx RF MADE1 1 ttaggttggtgcaaaagtaattgcggtttttgccattactttt 43 ttaggtt gtgcaaaagtaattg ggtttttg cattactttt humanchr1/239220001-239550000 174456 TTAGGTTAGTGCAAAAGTAATTGTGGTTTTTGTCATTACTTTT 174498 589************************************9986 PP >> MADE1 MADE1 (MAriner Derived Element 1), a TcMar-Mariner DNA transposon score bias Evalue hmmfrom hmm to alifrom ali to envfrom env to mod len acc ------ ----- --------- ------- ------- --------- --------- --------- --------- --------- ---- ! 31.0 6.7 6.3e-08 1 78 [. 302466 302389 .. 302466 302387 .. 80 0.80 Alignment: score: 31.0 bits xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx....xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx RF MADE1 1 ttaggttggtgcaaaagtaattgcggttttt....gccattacttttaatggcaaaaaccgcaattacttttgcaccaacct 78 ttag ttggtg aaaag a t tttt gc atta +aatggcaaaaacc caatt ttttgcacc acc humanchr1/239220001-239550000 302466 TTAGATTGGTGTAAAAGCATT-A---CTTTTaaaaGCAATTAAAAGCAATGGCAAAAACCACAATTGATTTTGCACCGACCA 302389 6899************97543.2...23333455566666666666799*****************************9985 PP >> MADE1 MADE1 (MAriner Derived Element 1), a TcMar-Mariner DNA transposon score bias Evalue hmmfrom hmm to alifrom ali to envfrom env to mod len acc ------ ----- --------- ------- ------- --------- --------- --------- --------- --------- ---- ! 25.0 7.0 4.9e-06 43 80 .] 174493 174456 .. 174513 174456 .. 80 0.94 Alignment: score: 25.0 bits xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx RF MADE1 43 taatggcaaaaaccgcaattacttttgcaccaacctaa 80 taatg caaaaacc caattacttttgcac aacctaa humanchr1/239220001-239550000 174493 TAATGACAAAAACCACAATTACTTTTGCACTAACCTAA 174456 5899*******************************986 PP >> MADE1 MADE1 (MAriner Derived Element 1), a TcMar-Mariner DNA transposon score bias Evalue hmmfrom hmm to alifrom ali to envfrom env to mod len acc ------ ----- --------- ------- ------- --------- --------- --------- --------- --------- ---- ? 6.9 7.2 2.2 41 71 .. 304073 304103 .. 304053 304109 .. 80 0.85 Alignment: score: 6.9 bits xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx RF MADE1 41 tttaatggcaaaaaccgcaattacttttgca 71 tt a tgg aaaaa ca tta ttttgca humanchr1/239220001-239550000 304073 TTAAGTGGGAAAAAATACACTTATTTTTGCA 304103 456789************************8 PP Internal pipeline statistics summary: ------------------------------------- Query sequence(s): 1 (660000 residues searched) Target model(s): 1 (80 nodes) Residues passing SSV filter: 60770 (0.0921); expected (0.02) Residues passing bias filter: 35792 (0.0542); expected (0.02) Residues passing Vit filter: 1612 (0.00244); expected (0.001) Residues passing Fwd filter: 1194 (0.00181); expected (1e-05) Total number of hits: 5 (0.000405) # CPU time: 0.02u 0.00s 00:00:00.02 Elapsed: 00:00:00.02 # Mc/sec: 2289.06 // [ok]