Seqanswers Leaderboard Ad

Collapse

Announcement

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

  • AfterQC: Automatic Filtering, Trimming, Error Removing and Quality Control for fastq

    AfterQC

    Code on github: https://github.com/OpenGene/AfterQC

    Automatic Filtering, Trimming, Error Removing and Quality Control for fastq data

    AfterQC can simply go through all fastq files in a folder and then output three folders: good, bad and QC folders, which contains good reads, bad reads and the QC results of each fastq file/pair.
    Currently it supports processing data from HiSeq 2000/2500/3000/4000, X10, X5, Nextseq 500/550, MiniSeq...


    Features:
    AfterQC does following tasks automatically:

    * Filters reads with too low quality, too short length or too many N
    * Filters reads with abnormal PolyA/PolyT/PolyC/PolyG sequences
    * Does per-base quality control and plots the figures
    * Trims reads at front and tail, according to QC results
    * For pair-end sequencing data, AfterQC automatically corrects low quality wrong bases in overlapped area of read1/read2
    * Detects and eliminates bubble artifact caused by sequencer due to fluid dynamics issues
    * Single molecule barcode sequencing support: if all reads have a single molecule barcode (see duplex sequencing), AfterQC shifts the barcodes from the reads to the fastq query names
    * Support single-end sequencing or pair-end sequencing

    Dependency:
    AfterQC uses editdistance module, run following before using AfterQC:

    pip install editdistance


    Simple usage:
    1, Prepare your fastq files in a folder
    2, For single-end sequencing, the filenames in the folder should be *R1*
    For pair-end sequencing, the filenames in the folder should be *R1* and *R2*


    cd /path/to/fastq/folder
    python path/to/AfterQC/after.py


    Two folders will be automatically generated, a folder 'good' stores the good reads and a folder 'bad' stores the bad reads
    AfterQC will print some statistical information after it is done, such how many good reads, how many bad reads, and how many reads are corrected.

    Quality Control only
    If you only want to get quality control statistics, run:
    python after.py --qc_only
    Last edited by [email protected]; 05-07-2016, 05:31 AM.
    OpenGene(Libraries and tools for NGS data analysis),AfterQC(Fastq Filtering and QC)
    FusionDirect.jl( Detect gene fusion), SeqMaker.jl(Next Generation Sequencing simulation)

  • #2
    Hi Chen,

    thanks, this looks interesting. What do you mean by this "* Detects and eliminates bubble artifact caused by sequencer due to fluid dynamics issues"?
    How would you recognize such problems? Are the reads filtered out?

    Comment


    • #3
      Hi luc,

      When a sequencer is working, it may raise some bubbles from its liquid lanes. Bubbles can affect DNA synthesis reactions and wash processes. According to our investigation, Illumina NextSeq sequencers are much more probable to have bubbles. NextSeq sequencers are widely used in clinical environment, like cancer testing. For these quality critical applications, the reads in the bubble areas must be discarded due to they are usually wrong representations of the original DNA templates and may cause fake mutations. Since these reads usually have high quality scores, they cannot be filtered by regular quality filters. We found that there were more and longer polyG reads (a polyG read is a read has long sub-sequence of Guanine) in the bubble areas.

      Based on this phenomenon, this debubble tool can visualise the bubbles, furthermore detect their contours, and then filter out the reads in bubbles.
      OpenGene(Libraries and tools for NGS data analysis),AfterQC(Fastq Filtering and QC)
      FusionDirect.jl( Detect gene fusion), SeqMaker.jl(Next Generation Sequencing simulation)

      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
      11 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
      68 views
      0 likes
      Last Post seqadmin  
      Working...
      X