Seqanswers Leaderboard Ad

Collapse

Announcement

Collapse
No announcement yet.
X
 
  • Filter
  • Time
  • Show
Clear All
new posts

  • Forum:Problem with homer in findMotifs.pl when using input and bg fasta



    I am using the following command:

    Code:
    findMotifs.pl input.fa fasta ./Output -fastaBg bg.fa -len 8,10,12 -norevopp
    The input and bg fasta have this structure and the bg fasta include all sequences in input plus many others:

    Code:
    >ENSMUST00000027125_Coq10b_mmu_chr1_55071635_55072702_+_utr_55071803_55072702(+)
    ATTTCTTTTGAATTCCGCTCCCTTCTGCACTCTCAGCTCGCTACTCTGTTCTTCGATGAAGTTGTGAAACAAATGGTAGC AGCCTTTGAAAGAAGAGCCTGTAAACTGTATGGTCCAGAGACAAACATACCTCGGGAATTAATGCTTCATGAAATTCACC ACACCTAAGAGGAAAATATTAGCTGCCTCCACCTACTCTTGGCTAGTTTGTTCACTTCTAGGAAGTCCTTTTACCATCTG` TTGAGAAGTCAGAAAGCATTTGTTAAACCTGCCTTGATTCTAAGCCCGTGCTGTTGAAAATTTGCACATTGAACATGGAC CCACTTGTACATAGAATTATTTCTTCAATCAAGTGTGACTCTAAGTATCATGTACATTTGCAGGCTCCGACCACCTTTGT AATAACGGATGTCATCACTGTTGCTAGGATACCACATTCCTCGTTTGAGTGTACAGATGAACAAGTCTTTTAATTCTCAC CTTACATGAAAAGGTTAGCTGAGATACAATGTGTGTTATATTAACCATATCATGTTTAAGTTATTAGGTTCAGAGTATTT GTAACTTATTGTTATTCGGCATGCCATATGGCTTAGGGTATTTGAATAATCATATATTTACCATTAAAACTGTGATTTAA AGTATTGCTAATGAAGTCTTAGCACTTTGGGTATTTTAATTGTTCTTATGGGTAGCAGTAGATGATTCAGTGTTGTTGGG

    But I get the following error:

    Code:
    Selected Options:
    
    Input file = /media/sequentia/synology_office2/Projects/554-Belloc/RIPseq/APAdetection/cpeb4_target.fa
    
    Promoter Set = fasta
    
    Output Directory = ./CPEB4
    
    Will use FASTA files for motif finding
    
    Target Sequences = /media/sequentia/synology_office2/Projects/554-Belloc/RIPseq/APAdetection/cpeb4_target.fa
    
    Background Sequences = /media/sequentia/synology_office2/Projects/554-Belloc/RIPseq/APAdetection/most_used_cpeb4.fa
    
    Motif length set at 8, 10, 12,
    
    Will not search the reverse strand
    
    Using custom gene IDs for GO analysis
    
    Parsing FASTA format files...
    
    
    
    
    Progress: Step4 - removing redundant promoters
    
    
    
    
    Progress: Step5 - adjusting background sequences for GC/CpG content...
    
    
    
    
    Sequences processed:
    
    0 total
    
    
    
    
    Frequency Bins: 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.6 0.7 0.8
    
    Freq Bin Count
    
    Illegal division by zero at /software/HOMER/bin/assignGeneWeights.pl line 63.
    
    
    
    
    Normalizing lower order oligos using homer2
    
    
    
    
    Reading input files...
    
    0 total sequences read
    
    Autonormalization: 1-mers (4 total)
    
    A inf% inf% -nan
    
    C inf% inf% -nan
    
    G inf% inf% -nan
    
    T inf% inf% -nan
    
    Autonormalization: 2-mers (16 total)
    
    AA inf% inf% -nan
    
    CA inf% inf% -nan
    
    GA inf% inf% -nan
    
    TA inf% inf% -nan
    
    AC inf% inf% -nan
    
    CC inf% inf% -nan
    
    GC inf% inf% -nan
    
    TC inf% inf% -nan
    
    AG inf% inf% -nan
    
    CG inf% inf% -nan
    
    GG inf% inf% -nan
    
    TG inf% inf% -nan
    
    AT inf% inf% -nan
    
    CT inf% inf% -nan
    
    GT inf% inf% -nan
    
    TT inf% inf% -nan
    
    Autonormalization: 3-mers (64 total)
    
    Normalization weights can be found in file: ./CPEB4/seq.autonorm.tsv
    
    Converging on autonormalization solution:
    
    ...............................................................................
    
    Final normalization: Autonormalization: 1-mers (4 total)
    
    A inf% inf% -nan
    
    C inf% inf% -nan
    
    (base) [B]idetoma@sequentia[/B]:[B]/media/sequentia/synology_office2/Projects/554-Belloc/RIPseq/motifs[/B]$ more nohup.out
    
    
    
    
    Selected Options:
    
    Input file = /media/sequentia/synology_office2/Projects/554-Belloc/RIPseq/APAdetection/cpeb4_target.fa
    
    Promoter Set = fasta
    
    Output Directory = ./CPEB4
    
    Will use FASTA files for motif finding
    
    Target Sequences = /media/sequentia/synology_office2/Projects/554-Belloc/RIPseq/APAdetection/cpeb4_target.fa
    
    Background Sequences = /media/sequentia/synology_office2/Projects/554-Belloc/RIPseq/APAdetection/most_used_cpeb4.fa
    
    Motif length set at 8, 10, 12,
    
    Will not search the reverse strand
    
    Using custom gene IDs for GO analysis
    
    Parsing FASTA format files...
    
    
    
    
    Progress: Step4 - removing redundant promoters
    
    
    
    
    Progress: Step5 - adjusting background sequences for GC/CpG content...
    
    
    
    
    Sequences processed:
    
    0 total
    
    
    
    
    Frequency Bins: 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.6 0.7 0.8
    
    Freq Bin Count
    
    Illegal division by zero at /software/HOMER/bin/assignGeneWeights.pl line 63.
    
    
    
    
    Normalizing lower order oligos using homer2
    
    
    
    
    Reading input files...
    
    0 total sequences read
    
    Autonormalization: 1-mers (4 total)
    
    A inf% inf% -nan
    
    C inf% inf% -nan
    
    G inf% inf% -nan
    
    T inf% inf% -nan
    
    Autonormalization: 2-mers (16 total)
    
    AA inf% inf% -nan
    
    CA inf% inf% -nan
    
    GA inf% inf% -nan
    
    TA inf% inf% -nan
    
    AC inf% inf% -nan
    
    CC inf% inf% -nan
    
    GC inf% inf% -nan
    
    TC inf% inf% -nan
    
    AG inf% inf% -nan
    
    CG inf% inf% -nan
    
    GG inf% inf% -nan
    
    TG inf% inf% -nan
    
    AT inf% inf% -nan
    
    CT inf% inf% -nan
    
    GT inf% inf% -nan
    
    TT inf% inf% -nan
    
    Autonormalization: 3-mers (64 total)
    
    Normalization weights can be found in file: ./CPEB4/seq.autonorm.tsv
    
    Converging on autonormalization solution:
    
    ...............................................................................
    
    Final normalization: Autonormalization: 1-mers (4 total)
    
    A inf% inf% -nan
    
    C inf% inf% -nan
    
    G inf% inf% -nan
    
    T inf% inf% -nan
    
    Autonormalization: 2-mers (16 total)
    
    AA inf% inf% -nan
    
    CA inf% inf% -nan
    
    GA inf% inf% -nan
    
    TA inf% inf% -nan
    
    AC inf% inf% -nan
    
    CC inf% inf% -nan
    
    GC inf% inf% -nan
    
    TC inf% inf% -nan
    
    AG inf% inf% -nan
    
    CG inf% inf% -nan
    
    GG inf% inf% -nan
    
    TG inf% inf% -nan
    
    AT inf% inf% -nan
    
    CT inf% inf% -nan
    
    GT inf% inf% -nan
    
    TT inf% inf% -nan
    
    Autonormalization: 3-mers (64 total)
    
    
    
    
    Progress: Step6 - Gene Ontology Enrichment Analysis
    
    Skipping...
    
    
    
    
    Progress: Step7 - Known motif enrichment
    
    
    
    
    Reading input files...
    
    0 total sequences read
    
    1006 motifs loaded
    
    Cache length = 11180
    
    Using hypergeometric scoring
    
    Checking enrichment of 1006 motif(s)
    
    |0% 50% 100%|
    
    =================================================================================
    
    Illegal division by zero at /software/HOMER/bin/findKnownMotifs.pl line 152.
    
    
    
    
    Progress: Step8 - De novo motif finding (HOMER)
    
    
    
    
    Scanning input files...
    
    !!! Something is wrong... are you sure you chose the right length for motif finding?
    
    !!! i.e. also check your sequence file!!!
    
    
    
    
    Scanning input files...
    
    !!! Something is wrong... are you sure you chose the right length for motif finding?
    
    !!! i.e. also check your sequence file!!!
    
    
    
    
    -blen automatically set to 2
    
    Scanning input files...
    
    !!! Something is wrong... are you sure you chose the right length for motif finding?
    
    !!! i.e. also check your sequence file!!!
    
    Use of uninitialized value in numeric gt (>) at /software/HOMER/bin/compareMotifs.pl line 1394.
    
    !!! Filtered out all motifs!!!
    
    Job finished

    What could be the problem?

Latest Articles

Collapse

  • seqadmin
    Strategies for Sequencing Challenging Samples
    by seqadmin


    Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
    03-22-2024, 06:39 AM
  • seqadmin
    Techniques and Challenges in Conservation Genomics
    by seqadmin



    The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

    Avian Conservation
    Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
    03-08-2024, 10:41 AM

ad_right_rmr

Collapse

News

Collapse

Topics Statistics Last Post
Started by seqadmin, Yesterday, 06:37 PM
0 responses
10 views
0 likes
Last Post seqadmin  
Started by seqadmin, Yesterday, 06:07 PM
0 responses
9 views
0 likes
Last Post seqadmin  
Started by seqadmin, 03-22-2024, 10:03 AM
0 responses
50 views
0 likes
Last Post seqadmin  
Started by seqadmin, 03-21-2024, 07:32 AM
0 responses
67 views
0 likes
Last Post seqadmin  
Working...
X