Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • bwa aln fails

    I'm running bwa verision 0.5.9

    I have 6 fastq files each 33 GB.

    I'm running:

    bwa aln -l 28 -t 6 <ref_file> <in_file> -f <out_file>

    Four of the six fastq files align to the indexed genome correctly. All six files were generated from the same library, adapters have been clipped. However 2 fastq files run through:

    [bwa_aln] 17bp reads: max_diff = 2
    [bwa_aln] 38bp reads: max_diff = 3
    [bwa_aln] 64bp reads: max_diff = 4
    [bwa_aln] 93bp reads: max_diff = 5
    [bwa_aln] 124bp reads: max_diff = 6
    [bwa_aln] 157bp reads: max_diff = 7
    [bwa_aln] 190bp reads: max_diff = 8
    [bwa_aln] 225bp reads: max_diff = 9

    then stall here (for hours) like its frozen. I'm not getting a segmentation fault or any error. The program just ceases to continue.

    Has anybody seen this before or have a fix???

  • #2
    First reasonable check, are those two really look the same as the other four? If it's fastq.gz, unzip it first, then head or tail on it.

    Comment


    • #3
      Are you running them at the same time, if so try to rerun the 2 files separately?
      Thanks,

      Comment


      • #4
        How big is your machine?
        Have you checked your inputs?
        Bad input and/or not enough memory causes problems with BWA.

        See my response in this thread for an explanation:
        Discussion of next-gen sequencing related bioinformatics: resources, algorithms, open source efforts, etc

        Comment


        • #5
          Well, not really sure what I did. I changed my pbs script, dropped number of threads from 6 to 4 and it seems to be running....whatever works man.

          Comment

          Latest Articles

          Collapse

          • seqadmin
            Recent Advances in Sequencing Analysis Tools
            by seqadmin


            The sequencing world is rapidly changing due to declining costs, enhanced accuracies, and the advent of newer, cutting-edge instruments. Equally important to these developments are improvements in sequencing analysis, a process that converts vast amounts of raw data into a comprehensible and meaningful form. This complex task requires expertise and the right analysis tools. In this article, we highlight the progress and innovation in sequencing analysis by reviewing several of the...
            05-06-2024, 07:48 AM
          • 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...
            04-22-2024, 07:01 AM

          ad_right_rmr

          Collapse

          News

          Collapse

          Topics Statistics Last Post
          Started by seqadmin, 05-14-2024, 07:03 AM
          0 responses
          19 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 05-10-2024, 06:35 AM
          0 responses
          44 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 05-09-2024, 02:46 PM
          0 responses
          54 views
          0 likes
          Last Post seqadmin  
          Started by seqadmin, 05-07-2024, 06:57 AM
          0 responses
          42 views
          0 likes
          Last Post seqadmin  
          Working...
          X