Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • FASTQC help

    I'm not a java programmer so I need some help doing something I hope there's someone here who can assist!

    According to the fastqc documentation, for calculating duplication levels it considers only the first 200,000 reads (to be memory efficient).... Memory is not an issue for me. I want a comprehensive output on my seq duplication level. Can anyone help me in modifying the module which does the above said task to accept and process all reads?

    Reason: I want to compare the results to those from picard's ELC.

    Warning: Is FASTQC open source? Can it be modified without permission? It's GPL v3 released so i think it can be modified openly

    Thanks in advance!

  • #2
    Originally posted by arkal View Post
    I'm not a java programmer so I need some help doing something I hope there's someone here who can assist!

    According to the fastqc documentation, for calculating duplication levels it considers only the first 200,000 reads (to be memory efficient).... Memory is not an issue for me. I want a comprehensive output on my seq duplication level. Can anyone help me in modifying the module which does the above said task to accept and process all reads?
    If you want to change this limit the value you need to change would be the line:

    Code:
    private final int OBSERVATION_CUTOFF = 200000;
    ..in OverRepresentedSequences.java. Or you could just remove the check all together by deleting the later lines:

    Code:
    if (uniqueSequenceCount == OBSERVATION_CUTOFF) {
      frozen = true;
    }
    You'll probably also need to change the default memory allocation in the launcher since you'll probably run out of memory in the 250MB allocation which is provided by default (edit the fastqc wrapper and increase the 250 to a much larger number).


    Originally posted by arkal View Post
    Reason: I want to compare the results to those from picard's ELC.
    Even with these changes you probably won't get the same answer as ELC as they're not measuring the same thing. The headline value in FastQC is the percentage of sequences which are not unique, which isn't quite the same thing as the percentage of different sequences. I'd prefer to have the latter value in FastQC but inferring that from the sample which we take has proved to be tricky (though I think we now know how to do this for a future release).

    Originally posted by arkal View Post
    Warning: Is FASTQC open source? Can it be modified without permission? It's GPL v3 released so i think it can be modified openly
    Yes, FastQC is open source and you can modify it for your own uses with no restrictions at all. The only restriction from the GPL is that if you distribute a modified version of the program you need to release the source for that under GPLv3 as well.

    Hope this helps

    Simon.

    Comment

    Latest Articles

    Collapse

    • seqadmin
      Strategies for Sequencing Challenging Samples
      by seqadmin


      Despite advancements in sequencing platforms and related sample preparation technologies, certain sample types continue to present significant challenges that can compromise sequencing results. Pedro Echave, Senior Manager of the Global Business Segment at Revvity, explained that the success of a sequencing experiment ultimately depends on the amount and integrity of the nucleic acid template (RNA or DNA) obtained from a sample. “The better the quality of the nucleic acid isolated...
      03-22-2024, 06:39 AM
    • seqadmin
      Techniques and Challenges in Conservation Genomics
      by seqadmin



      The field of conservation genomics centers on applying genomics technologies in support of conservation efforts and the preservation of biodiversity. This article features interviews with two researchers who showcase their innovative work and highlight the current state and future of conservation genomics.

      Avian Conservation
      Matthew DeSaix, a recent doctoral graduate from Kristen Ruegg’s lab at The University of Colorado, shared that most of his research...
      03-08-2024, 10:41 AM

    ad_right_rmr

    Collapse

    News

    Collapse

    Topics Statistics Last Post
    Started by seqadmin, Yesterday, 06:37 PM
    0 responses
    10 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, Yesterday, 06:07 PM
    0 responses
    10 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-22-2024, 10:03 AM
    0 responses
    51 views
    0 likes
    Last Post seqadmin  
    Started by seqadmin, 03-21-2024, 07:32 AM
    0 responses
    67 views
    0 likes
    Last Post seqadmin  
    Working...
    X