Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • bedtools2-24.0 error

    Not sure what this error is. Would the bed file being used cause this? Thank you

    Code:
    cmccabe@DTV-A5211QLM:~/Desktop/NGS$ coverageBed -hist -d -abam /home/cmccabe/Desktop/NGS/pool_I_090215/IonXpress_008_150902_newheader.bam -b /home/cmccabe/Desktop/NGS/bed/sorted_unix_5column_xgen_targets.bed > /home/cmccabe/Desktop/NGS/pool_I_090215/output.bam.hist.txt
    
    ***** ERROR: -counts, -d, -mean, and -hist are all mutually exclusive options. *
    Bed file
    Code:
    chr1	955542	955763	+	AGRN:exon.1
    chr1	957570	957852	+	AGRN:exon.2
    chr1	976034	976270	+	AGRN:exon.2;AGRN:exon.3;AGRN:exon.4

  • #2
    As the error says, you are trying to use options that are incompatible (can't be used) with each other.

    Comment


    • #3
      Strange that same command used to work in v. 2.20. Thank you .

      Basically, I am just trying to get the counts per base in the bam.

      So maybe I do not need the -hist option?

      Comment


      • #4
        I did the below command to output the per base coverage per target::

        Code:
        coverageBed -hist -abam /home/cmccabe/Desktop/NGS/pool_I_090215/IonXpress_008_150902_newheader.bam -b /home/cmccabe/Desktop/NGS/bed/sorted_unix_5column_xgen_targets.bed > /home/cmccabe/Desktop/NGS/pool_I_090215/output.bam.hist.txt
        output.bam.hist.txt

        Code:
        chr1    14251    14393    X28LU:04862:12482    0    +    14251    14393    0,0,0    1    142,    0,    0    142    142    1.0000000
        chr1    16224    16400    X28LU:08504:06628    0    -    16224    16400    0,0,0    1    176,    0,    0    176    176    1.0000000
        chr1    16324    16500    X28LU:06201:09146    3    -    16324    16500    0,0,0    1    176,    0,    0    176    176    1.0000000
        this output looks very different then 2.20.0 and I noticed that there is a new release 2.25.0, maybe I should use that? Thank you

        The target bed looks like this, could that be part of the problem?

        Code:
        chr1    955542    955763    +    AGRN:exon.1
        chr1    957570    957852    +    AGRN:exon.2
        chr1    976034    976270    +    AGRN:exon.2;AGRN:exon.3;AGRN:exon.4
        chr1    976542    976787    +    AGRN:exon.3;AGRN:exon.5
        Thank you .

        Comment

        Latest Articles

        Collapse

        • seqadmin
          Essential Discoveries and Tools in Epitranscriptomics
          by seqadmin


          The field of epigenetics has traditionally concentrated more on DNA and how changes like methylation and phosphorylation of histones impact gene expression and regulation. However, our increased understanding of RNA modifications and their importance in cellular processes has led to a rise in epitranscriptomics research. “Epitranscriptomics brings together the concepts of epigenetics and gene expression,” explained Adrien Leger, PhD, Principal Research Scientist on Modified Bases...
          Yesterday, 07:01 AM
        • seqadmin
          Current Approaches to Protein Sequencing
          by seqadmin


          Proteins are often described as the workhorses of the cell, and identifying their sequences is key to understanding their role in biological processes and disease. Currently, the most common technique used to determine protein sequences is mass spectrometry. While still a valuable tool, mass spectrometry faces several limitations and requires a highly experienced scientist familiar with the equipment to operate it. Additionally, other proteomic methods, like affinity assays, are constrained...
          04-04-2024, 04:25 PM

        ad_right_rmr

        Collapse

        News

        Collapse

        Topics Statistics Last Post
        Started by seqadmin, 04-11-2024, 12:08 PM
        0 responses
        39 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 04-10-2024, 10:19 PM
        0 responses
        41 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 04-10-2024, 09:21 AM
        0 responses
        35 views
        0 likes
        Last Post seqadmin  
        Started by seqadmin, 04-04-2024, 09:00 AM
        0 responses
        55 views
        0 likes
        Last Post seqadmin  
        Working...
        X