Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • Making bam files in Pysam

    I have made a script that should generate bam files with the help of Pysam. But the problem is that that code does not work.

    The mistake in the code is somewhere in the following code:
    =============================================

    inputSamFile = pysam.Samfile(inputBamFile, "rb")
    outputSamFile = pysam.Samfile(outputBamFile, "wb", template=inputSamFile)
    ============================================

    When I run this code, the following error message is generated:
    ===============================================

    [martijn@animal provisional]$ python DPBamFileMaker.py TCGA-A6-2672-10A-01W-0833-10_SOLiD.bam man.bam
    Traceback (most recent call last):
    File "DPBamFileMaker.py", line 117, in <module>
    main(sys.argv)
    File "DPBamFileMaker.py", line 114, in main
    convertInputBamFileToOutputBamFileWithOnlyDPs(inputBamFile, outputBamFile)
    File "DPBamFileMaker.py", line 49, in convertInputBamFileToOutputBamFileWithOnlyDPs
    outputSamFile = pysam.Samfile(outputBamFile, "wb", template=inputSamFile)
    File "csamtools.pyx", line 463, in csamtools.Samfile.__cinit__ (pysam/csamtools.c:4758)
    File "csamtools.pyx", line 609, in csamtools.Samfile._open (pysam/csamtools.c:6330)
    IOError: could not open file `man.bam`
    ==================================================

    I don't understan this error.

    What am I doing wrong in this code? Can anyone tell me?

  • #2
    Do you have write permissions it this directory?

    Comment

    Latest Articles

    Collapse

    • seqadmin
      Choosing Between NGS and qPCR
      by seqadmin



      Next-generation sequencing (NGS) and quantitative polymerase chain reaction (qPCR) are essential techniques for investigating the genome, transcriptome, and epigenome. In many cases, choosing the appropriate technique is straightforward, but in others, it can be more challenging to determine the most effective option. A simple distinction is that smaller, more focused projects are typically better suited for qPCR, while larger, more complex datasets benefit from NGS. However,...
      10-18-2024, 07:11 AM
    • seqadmin
      Non-Coding RNA Research and Technologies
      by seqadmin




      Non-coding RNAs (ncRNAs) do not code for proteins but play important roles in numerous cellular processes including gene silencing, developmental pathways, and more. There are numerous types including microRNA (miRNA), long ncRNA (lncRNA), circular RNA (circRNA), and more. In this article, we discuss innovative ncRNA research and explore recent technological advancements that improve the study of ncRNAs.

      Nobel Prize for MicroRNA Discovery
      This week,...
      10-07-2024, 08:07 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, 11-01-2024, 06:09 AM
    0 responses
    18 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 10-30-2024, 05:31 AM
    0 responses
    18 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 10-24-2024, 06:58 AM
    0 responses
    24 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 10-23-2024, 08:43 AM
    0 responses
    53 views
    0 likes
    Last Post seqadmin  
    Working...
    X