Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • #16
    Originally posted by cmccabe View Post
    Here is the error:

    Code:
    Microsoft Windows [Version 6.1.7601]
    Copyright (c) 2009 Microsoft Corporation.  All rights reserved.
    
    C:\Users\cmccabe>cd "C:\Users\cmccabe\Desktop\Ion_Torrent_ Software\bedtools-2.1
    7.0"  (files here)
    
    C:\Users\cmccabe\Desktop\Ion_Torrent_ Software\bedtools-2.17.0>[COLOR="Red"][SIZE="2"]bin\bedtools[/SIZE][/COLOR] bedToBam -i Ampli
    cons.bed -g hg19_genome.txt > Amplicons.bam
    'bedToBam' is not recognized as an internal or external command,
    operable program or batch file.
    Thanks.
    I think you left out the executable program name "bedtools" when issuing the command. See the highlight in red.

    Comment


    • #17
      $ cd /cygdrive/c/Users/cmccabe/Desktop/Ion_Torrent_ Software/bedtools-2.17.0

      $ bin/bamToBed -i Amplicons.bed -g hg19_genome.txt > Amplicons.bam

      The command runs in perl however the input Amplicons.bed in 44kb and the output amplicons.bam is 0kb. That doesent seem right and below is a couple lines from the BED file. Thanks.

      chr1 2160290 2160475 AMPL2688828352 . SKI
      chr1 2160480 2160634 AMPL2688829530 . SKI

      Comment


      • #18
        Should you not be running "bedToBam" if you are generating a BAM file?

        Comment


        • #19
          What I am trying to do is to calculate predicted off-target reads based on a BED file?

          One sugesstion was to use bedToBam (from bedtools) followed by Picard's CalculateHSMetrics
          Thanks.

          Comment


          • #20
            Originally posted by cmccabe View Post
            What I am trying to do is to calculate predicted off-target reads based on a BED file?

            One sugesstion was to use bedToBam (from bedtools) followed by Picard's CalculateHSMetrics
            Thanks.
            In post #17 you had "bamToBed" listed instead of "bedToBam". Perhaps a typo or a wrong command?

            Comment


            • #21
              bed2bam

              I switched the command to bedToBam and got the same result. Maybe there is something wrong with the format of the BED file?

              Comment


              • #22
                Since you are trying to do this on windows there could be some issue with file formats etc. It would be tough to diagnose them through the forum.

                Do you think you would be up to trying this through the cygwin terminal (instead of windows terminal)?
                Last edited by GenoMax; 07-10-2013, 06:29 AM.

                Comment


                • #23
                  bed2bam

                  I am using the Cygwin Terminal to run the command. Thank you.

                  Comment


                  • #24
                    Can you try the following in cygwin terminal?
                    Code:
                    $ dos2unix Amplicons.bed
                    $ dos2unix hg19_genome.txt
                    This should reformat the files to unix file format.

                    Try running your bedToBam command again.

                    Comment


                    • #25
                      bed2bam

                      The file is is still 0kb, should I re-install cygwi and re-compile bedtools?

                      Comment


                      • #26
                        Is your bed file tab-delimited? The two lines you posted don't appear to be although tabs may have been lost by copying/pasting.

                        Comment


                        • #27
                          bed2bam

                          Yes, the format of the BED is Tab-deliminated.

                          Comment


                          • #28
                            Once I explicitly changed the format in the example you pasted I was able to create the bam file.

                            What about your hg19_genome file?

                            Comment


                            • #29
                              bed2bam

                              The hg19_genome is also tab-deliminated.

                              Comment


                              • #30
                                bed2Bam

                                I re-compiled in cygwin

                                did dos2unix Amplicons.bed
                                dos2unix hg19_genome.txt

                                both of those files are tab- deliminated

                                and the commands I used in Cygwin:

                                $ cd "/cygdrive/c/Users/cmccabe/Desktop/Ion_Torrent_ Software/bedtools-2.17.0"

                                $ bin/bamToBed -i Amplicons.bed -g hg19_genome.txt > Amplicons.bam

                                Amplicons.bed (went from 44kb to 43kb)
                                chr1 2160290 2160475 SKI . AMPL2688828352
                                chr1 2160480 2160634 SKI . AMPL2688829530
                                chr1 2160632 2160738 SKI . AMPL2688828058

                                hg19_genome.txt (1kb) - acc chromosomes and lenghths are in this file
                                chr1 249,250,621
                                chr2 243,199,373

                                Amplicons.bam is still 0kb and not working.

                                Comment

                                Latest Articles

                                Collapse

                                • seqadmin
                                  Best Practices for Single-Cell Sequencing Analysis
                                  by seqadmin



                                  While isolating and preparing single cells for sequencing was historically the bottleneck, recent technological advancements have shifted the challenge to data analysis. This highlights the rapidly evolving nature of single-cell sequencing. The inherent complexity of single-cell analysis has intensified with the surge in data volume and the incorporation of diverse and more complex datasets. This article explores the challenges in analysis, examines common pitfalls, offers...
                                  06-06-2024, 07:15 AM
                                • seqadmin
                                  Latest Developments in Precision Medicine
                                  by seqadmin



                                  Technological advances have led to drastic improvements in the field of precision medicine, enabling more personalized approaches to treatment. This article explores four leading groups that are overcoming many of the challenges of genomic profiling and precision medicine through their innovative platforms and technologies.

                                  Somatic Genomics
                                  “We have such a tremendous amount of genetic diversity that exists within each of us, and not just between us as individuals,”...
                                  05-24-2024, 01:16 PM

                                ad_right_rmr

                                Collapse

                                News

                                Collapse

                                Topics Statistics Last Post
                                Started by seqadmin, Yesterday, 07:49 AM
                                0 responses
                                12 views
                                0 likes
                                Last Post seqadmin  
                                Started by seqadmin, 06-20-2024, 07:23 AM
                                0 responses
                                14 views
                                0 likes
                                Last Post seqadmin  
                                Started by seqadmin, 06-17-2024, 06:54 AM
                                0 responses
                                16 views
                                0 likes
                                Last Post seqadmin  
                                Started by seqadmin, 06-14-2024, 07:24 AM
                                0 responses
                                25 views
                                0 likes
                                Last Post seqadmin  
                                Working...
                                X