Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • SplitNCigarReads exception

    Hi all,

    I'm trying to use the GATK pipeline for analyzing RNA-seq data and to detect variants. According to the GATK Best Practices, I need to allow reads for N cigars.
    So I tried
    java -jar GenomeAnalysisTK.jar -T SplitNCigarReads -R /human_g1k_v37.fasta -I Sample.dedupped.bam -o split.bam -U ALLOW_N_CIGAR_READS

    I started working with version 3.2.2. This works fine for a while, until it comes to a certain line in the BAM-file and spits out this error:
    Bad input: Cannot split this read (might be an empty section between Ns, for example 1N1D1N): 77M54369N2D79N24M

    I looked around on the internet and found that it helps sometimes to add -fixNDN to the command. But this results in the following error:
    Message: Code exception.

    Futhermore, it could be a bug, so I installed the latest version of GATK (3.3.0) but I get the exact same errors (both with and without -fixNDN).

    Any suggestions are most welcomed!
    Thanks,
    Lien

  • #2
    Hi,

    I think I run into this issue the other day, too. Try the latest nightly build of GATK. That fixed it for me. This is a bug that has been fixed if I recall it right.

    Cheers

    Comment


    • #3
      Indeed, the latest build fixed it.

      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...
        Yesterday, 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, Today, 06:58 AM
      0 responses
      6 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, Yesterday, 08:18 AM
      0 responses
      15 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, Yesterday, 08:04 AM
      0 responses
      14 views
      0 likes
      Last Post seqadmin  
      Started by seqadmin, 06-03-2024, 06:55 AM
      0 responses
      13 views
      0 likes
      Last Post seqadmin  
      Working...
      X