BOL: Related items

kSNP3.0: SNP detection and phylogenetic analysis of genomes without genome alignment or reference genome

Jit — Fri, 08 Dec 2017 16:48:40 -0600

Sept. 20, 2017 Version 3.1 released. Major upgrade. Version 3.1 fixes the problems with SNP annotation that arose when NCBI discontinued use of GI numbers. Please read carefully the Preface (page 3) and the File of annotated genomes section (pages 9-10) in the version 3.1 User Guide. Thanks to Tom Slezak for revsing the get_genbank_file3 script and to Tod Stuber (USDA) for testing version 3.1 even though he doesn't need the annotation feature. All users are encouraged to upgrade to version 3.1.

Address of the bookmark: https://sourceforge.net/projects/ksnp/files/

Snakemake workflow: dna-seq-gatk-variant-calling

Jit — Thu, 25 Jul 2019 12:55:07 -0500

This Snakemake pipeline implements the GATK best-practices workflow for calling small genomic variants.

Address of the bookmark: https://github.com/snakemake-workflows/dna-seq-gatk-variant-calling

Best Practices for Variant Calling with the GATK

biogeek — Sat, 22 Feb 2020 03:07:31 -0600

The presentations below were filmed during the March 2015 GATK Workshop, part of the BroadE Workshop series. At the time of this workshop, the current version of Broad’s Genome Analysis Toolkit (GATK) was version 3.3.

Genome Analysis Toolkit

03/19/15	Introduction to High-Throughput Sequencing data formats and methods	Joel Thibault	PDF	Video
03/19/15	Introduction to the GATK	Geraldine Van der Auwera	PDF	Video
03/19/15	Mapping, processing, and duplicate marking with Picard tools	Matt Sooknah	PDF	Video
03/19/15	Mapping and processing RNAseq	Ami Levy-Moonshine	PDF	Video
03/19/15	Indel realignment	Mark Fleharty	PDF	Video
03/19/15	Base quality score recalibration	David Roazen	PDF	Video
03/19/15	Introduction to variant discovery: calling cohorts	Louis Bergelson	PDF	Video
03/19/15	Variant calling and joint genotyping	Sheila Chandran	PDF	Video
03/19/15	Variant quality score recalibration	Bertrand Haas	PDF	Video
03/19/15	Introduction to working with variants	Yossi Farjoun	PDF	Video
03/19/15	Genotype refinement	Laura Gauthier	PDF	Video
03/19/15	Annotation and variant evaluation	David Benjamin	PDF	Video

Address of the bookmark: https://www.broadinstitute.org/partnerships/education/broade/best-practices-variant-calling-gatk-1

AIRVF: a filtering toolbox for precise variant calling in Ion Torrent sequencing

BioStar — Fri, 22 Dec 2017 00:31:06 -0600

AIRVF that works on flowgram, raw and mapped reads and called variants to reduce artifact-driven false variant calls. Tests on sequencing data of standard reference material showed up to ∼98% reduction of false variants when combined to conventional public pipelines and ∼48% to the in-house commercial solution, with a minimal loss of sensitivity

The program with a detailed manual is available at https://sourceforge.net/projects/airvf/

Address of the bookmark: https://sourceforge.net/projects/airvf/

The new corona variant has 23 mutations in all, which is unusually huge !

Shruti Paniwala — Wed, 23 Dec 2020 03:50:50 -0600

The new SARS-CoV-2 version, B.1.1.7, which was first seen in the third week of September in Kent and Greater London, has since spread to other locations in the UK. According to the COVID-19 Genomics UK Consortium (COG-UK Consortium) that analysed the genome data of the virus and identified the variant, the new variant has been spreading "rapidly" over the last four weeks and has now been detected in other locations in the UK, suggesting further spread of the variant in the region.

According to a preliminary report posted on December 19 by the COG-UK Consortium scientists, as of December 15, 1,623 variant genomes have been sequenced. In a December 21 tweet, COG-UK Consortium said that it added 2,963 more genome sequences of SARS-CoV-2, of which 942 (32%) belong to the new variant. The Consortium intends to sequence 20,000 more SARS-CoV-2 genomes in the next two weeks to further ascertain the spread of the variant.

There is no clear proof, at least not yet, that it does cause severe pandemic. But there is a justification for seriously taking the possibility. Another coronavirus lineage in South Africa has acquired one specific mutation that is also present in B.1.1.7. This variant is increasingly spreading across South Africa's coastal regions. And doctors have observed in preliminary research that individuals infected with this variant bear a higher viral load-a higher concentration of the virus in their upper respiratory tract. In many viral diseases, this is associated with more severe symptoms.

cisMuton: caller that detects SNVs/indels by comparing target (foreground) and control (background) samples

Jit — Tue, 11 Feb 2020 05:55:00 -0600

cisMuton is a caller that detects SNVs/indels by comparing target (foreground) and control (background) samples.

cisMuton calls mutations from target capture regions, which are defined by the overlapping regions of ${GROUP}.target.bed and ${GROUP}.gene.bed.

Please read the Quick Start Guide for cisCall first for an outline of how to run cisMuton.

Address of the bookmark: https://static.ciscall.org/cisCall5_doc/index.html

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

Rahul Nayak — Mon, 13 Nov 2017 05:10:23 -0600

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, Nextseq 500/550, MiniSeq...and other Illumina 1.8 or newer formats

Address of the bookmark: https://github.com/OpenGene/AfterQC

MinION_GC: An R script to do some QC on MinION data

Radha Agarkar — Sun, 03 Dec 2017 15:19:18 -0600

Other tools focus on getting data out of the fastq or fast5 files, which is slow and computationally intensive. The benefit of this approach is that it works on a single, small, .txt summary file. So it's a lot quicker than most other things out there: it takes about a minute to analyse a 4GB flowcell on my laptop.

https://github.com/roblanf/minion_qc

Address of the bookmark: https://github.com/roblanf/minion_qc

NextSV: a meta-caller for structural variants from low-coverage long-read sequencing data

Jit — Mon, 06 Aug 2018 17:24:53 -0500

NextSV, a meta SV caller and a computational pipeline to perform SV calling from low coverage long-read sequencing data. NextSV integrates three aligners and three SV callers and generates two integrated call sets (sensitive/stringent) for different analysis purpose. The output of NextSV is in ANNOVAR-compatible bed format. Users can easily perform downstream annotation using ANNOVAR and disease gene discovery using Phenolyzer.

Address of the bookmark: https://github.com/Nextomics/NextSV

BamView: a free interactive display of read alignments in BAM data files

Neel — Fri, 09 Nov 2018 13:43:22 -0600

To run the application on UNIX from the downloaded jar file run the UNIX:

java -mx512m -jar BamView.jar

and extra command line options are given when '-h' is used:

java -jar BamView.jar -h

BAM files can be specified on the command line with the '-a' option:

java -mx512m -jar BamView.jar -a pathToFile/sorted.bam

If a BAM filename is not given on the command line BamView will prompt for a file to be entered. The BAM index file should have the same name as the BAM file but with a '.bai' suffix. Multiple BAM files can be loaded and overlaid in the viewer. To make this easier BamView will read in files that contain a list of filenames.

Address of the bookmark: http://bamview.sourceforge.net/