The "Ifs" of NGS QC and Trimming

1. Ensures Data Integrity
   If you want to minimize errors in downstream analyses, QC and trimming remove low-quality reads and bases, ensuring high-confidence data. This step is essential for reliable variant calling, assembly, and other applications.

2. Removes Contaminants
   If adapter sequences or contaminants are present in the raw reads, trimming can eliminate them. This prevents issues like misalignment or incorrect biological interpretations, ensuring cleaner data for analysis.

3. Improves Mapping and Assembly
   If your goal is better alignment to a reference genome or improved de novo assembly, trimming low-quality bases and adapters is critical. High-quality reads map more efficiently and generate more accurate assemblies.

4. Reduces Computational Load
   If you want to save computational resources, trimming reduces the dataset size, which speeds up processing and analysis. Clean datasets mean less computational time spent on processing low-quality data.

5. Prepares for Standardized Analyses
   If your project involves multiple datasets, QC and trimming ensure uniformity across them. This standardization makes comparisons valid and reproducible, particularly in large collaborative studies.

The "Buts" of NGS QC and Trimming

1. Risk of Over-Trimming
   But excessive trimming can lead to the loss of informative sequences, reducing read depth and potentially discarding biologically relevant data. This is especially critical in studies with limited sequencing depth.

2. Bias Introduction
   But trimming algorithms might introduce biases, especially if they inadvertently remove sequences with specific biological patterns. This can skew results and compromise biological insights.

3. Loss of Context in Paired-End Reads
   But trimming one read in a pair more than the other can lead to loss of pairing information. This complicates downstream analyses that rely on paired-end data, such as structural variant detection.

4. Time and Resource Intensive
   But running QC and trimming for large datasets can be computationally expensive and time-consuming. As sequencing depth increases, preprocessing becomes a bottleneck in the analysis pipeline.

5. Variable Standards
   But the criteria for trimming (e.g., quality threshold, minimum read length) can vary between tools and datasets. This variability may affect reproducibility and comparability of results across studies.

Balancing the "Ifs" and "Buts"

To maximize the benefits of QC and trimming while mitigating the challenges, consider the following best practices:

• Use QC Tools Wisely: Start with tools like FastQC to identify quality issues in your raw data. Visualizing quality metrics helps tailor your trimming parameters.

• Choose Reliable Trimming Tools: Tools like Trimmomatic, Cutadapt, and BBduk offer adaptive and customizable trimming options. Select one that aligns with your dataset and project goals.

• Set Reasonable Parameters: Avoid over-trimming by setting quality thresholds and minimum read lengths that balance data retention and quality improvement.

• Test Downstream Effects: Validate the impact of QC and trimming on downstream analyses, such as alignment efficiency, variant calling accuracy, or assembly quality.

• Document Your Workflow: Maintain detailed records of the parameters and tools used for QC and trimming. This ensures reproducibility and enables better troubleshooting.

Conclusion

NGS quality control and trimming are essential steps to ensure reliable and accurate data for analysis. While the "ifs" highlight the clear benefits of these steps, the "buts" remind us of the potential pitfalls. By adopting best practices and carefully balancing these considerations, you can optimize your preprocessing workflow and unlock the full potential of your sequencing data.

minimap -t 24 assembly.fasta long_reads.fastq.gz | racon -t 24 long_reads.fastq.gz - assembly.fasta racon_assembly.fasta
nucmer -p nucmer assembly.fasta racon_assembly.fasta
show-snps -C -T -r nucmer.delta

This reports Racon's changes in a table. You can exclude indels with the -I option in show-snps. 

This process (Racon -> MUMmer -> SNP table) solves the problem I originally raised in this issue. So as far as I'm concerned, you can close this issue (or keep it open if you still want to implement some kind of variant table).

We accomplish this by first calculating genotype likelihoods for every site. For a given site in a diploid genome, there are 10 possible genotypes (AA, AC, AG, AT, CC, CG, CT, GG, GT, TT). Referee takes as input the genotype likelihoods calculated for all 10 genotypes given the called reference base at each position.

Referee is a program to calculate a quality score for every position in a genome assembly. This allows for easy filtering of low quality sites for any downstream analysis.

https://github.com/gwct/referee

Address of the bookmark: https://gwct.github.io/referee/#

Computational methods used by the Shasta assembler include:

• Using a run-length representation of the read sequence. This makes the assembly process more resilient to errors in homopolymer repeat counts, which are the most common type of errors in Oxford Nanopore reads.
• Using in some phases of the computation a representation of the read sequence based on markers, a fixed subset of short k-mers (k ≈ 10).

More at https://chanzuckerberg.github.io/shasta/index.html

Address of the bookmark: https://github.com/chanzuckerberg/shasta

More at https://github.com/ekg/mutatrix

./mutatrix -S sample -P test/ -p 2 -n 10 reference.fasta

Address of the bookmark: https://github.com/ekg/mutatrix

Yau Group develops statistical and computational methods for the analysis of genomic datasets with a particular interest in cancer sequencing applications and the use of Bayesian Statistics.

Yau Group are currently have projects in somatic mutation analysis of heterogeneous cancers, data fusion or integration techniques and single cell genomics.

More @ http://www.well.ox.ac.uk/~cyau/index.html

The Department is at the forefront of developments in machine learning, data science, and computational methods for understanding complex biological systems, integrating these approaches into both its research and teaching. This role offers significant scope to shape the future of bioinformatics and quantitative biology at Oxford and beyond.

You will be expected to develop and lead a world-class research programme in bioinformatics, helping to strengthen Oxford's position in areas such as statistical genomics, computational biology, machine learning for biological data, engineering biology, systems biology, and biomedical data science. We value the breadth of approaches that contribute to modern bioinformatics. Whatever your particular focus, you will understand how quantitative and computational methods can be applied to important biological and biomedical questions, and how collaboration across disciplines can generate new insights and opportunities. Alongside your own research, you will support the development of doctoral students, postdoctoral researchers, and early-career academics. You will contribute to our teaching at undergraduate and postgraduate level, delivering lectures, supervising student projects, and supporting graduate research students.

We are looking for someone who combines scholarly distinction with a collaborative and ambitious approach to research. You will be a researcher with an international reputation in bioinformatics, with a strong record of producing influential research at the interface of statistics, computational biology, genomics, and data science, and of securing research funding. Your work should advance the field, whether through methodological innovation, biological discovery, or the development of new approaches to understanding complex biological and biomedical data. You will bring experience of shaping research strategy, building productive collaborations, and exercising academic leadership. Above all, we wish to appoint someone who shares our values of intellectual curiosity, academic excellence, scholarly integrity, and a commitment to fostering a positive, supportive, and inclusive environment for all.

If you would like to discuss this post, please contact the Head of the Department of Statistics, Professor Frank Windmeijer, frank.windmeijer@stats.ox.ac.uk, or for information about University College please contact the Senior Tutor, Dr Andrew Bell, andrew.bell@univ.ox.ac.uk. All enquiries will be treated in strict confidence and will not form part of the selection decision.

The closing date for applications is 12:00 noon UK time on Monday 21 September 2026. Interviews are expected to be held in November 2026.

Apply here https://my.corehr.com/pls/uoxrecruit/erq_jobspec_version_4.display_form

Abstract: Strand NGS includes comprehensive workflows for DNA-Seq, RNA-Seq, Small RNA-Seq, ChIP-Seq, MeDIP-Seq, and Methyl-Seq analysis. Each workflow includes a quality assessment and filter section, followed by a workflow-specific analysis section. The pipeline functionality in Strand NGS allows users to execute a sequence of analysis steps with specific parameters - all without any manual intervention. This simplifies the analysis in large scale sequencing projects where every sample needs to be processed identically.

In this webinar we will discuss the pre-packaged pipelines present in Strand NGS. The packaged pipelines have well-chosen default parameters and are suitable for users analyzing data for the first time in the tool. We will also show how advanced users can customize pipelines and share them with other Strand NGS users. Finally, we will show a brief glimpse of an elaborate pipeline that aligns reads, filters poor-quality matches, computes coverage metrics, identifies variants, checks for sample cross-contamination, and emails quality reports - all from within Strand NGS.

Speaker: Dr. Vamsi Veeramachaneni, Vice President - Bioinformatics, Strand Life Sciences

Details: Session 1: 2:30 PM IST, Session 2 : 10:30 PM IST
Register here: http://www.strand-ngs.com/webinar_registration

by Dr. Pandurang Kolekar, Bioinformatics Engineer, Strand Life Sciences

Abstract:

Unique Molecular Identifiers (UMIs) are short random nucleotide sequences that are increasingly being used in high-throughput sequencing experiments. In this webinar, we will highlight the UMI-friendly features of Strand NGS v3.1 including support for handling well known and customised UMI libraries, QC metrics, consensus alignment, UMI-based family size filters for read list, genome browser enabled with UMI-specific features and filters, UMI-aware variant calling parameters, and exporting UMI-tagged aligned samples. These all features together empower users to harness the potential of UMI-tagged NGS data for deeper insights. A case study demonstrating application of these UMI-based features in Strand NGS for low frequency variant calling in cfDNA sample will be presented.

UMI-tagged NGS libraries allow, ultra-sensitive detection of low frequency variants from liquid biopsy samples using DNA-Seq and accurate quantification of transcript-level expression using RNA-Seq. The recent release of Strand NGS v3.1, is equipped with the necessary features to efficiently analyse UMI-tagged NGS data helping researchers and labs involved in rare variant calling like in cfDNA based cancer diagnostics, and accurate transcript quantification with RNA-Seq.

Webinar Details:

Session 1: 13 Dec 2017, 2:30 PM IST
Session 2: 13 Dec 2017, 9:30 PM IST

Register here: http://www.strand-ngs.com/webinar_registration

Details:
Session 1: 28 Feb 2018, 9 AM CET
Session 2: 28 Feb 2018, 8 AM PST
Register here: http://www.strand-ngs.com/webinar_registration

About Speaker:

Dr. Suman Kapoor, Manager- Application Science at Strand Life Sciences, has over 11 years experience in molecular biology, next-generation sequencing based testing, clinical genomics, and personalized medicine for disease management and prenatal testing. Dr. Suman holds a Ph.D in Molecular and Cell Biology from Indian Institute of Science, Bangalore. Prior to joining Strand NGS team, Suman has worked extensively on protein synthesis in eubacteria and has experience working in CAP and NABL accredited lab validating and interpreting NGS based diagnostic tests.