Tag: FreeBayes

A Benchmark of Genetic Variant Calling Pipelines Using Metagenomic Short-Read Sequencing

Introduction Short-read metagenomic sequencing is the technique most widely used to explore the natural habitat of millions of bacteria. In comparison with 16S rRNA sequencing, shotgun metagenomic sequencing (MGS) provides sequence information of the whole genomes, which can be used to identify different genes present in an individual bacterium and…

Continue Reading A Benchmark of Genetic Variant Calling Pipelines Using Metagenomic Short-Read Sequencing

An FGFR2 mutation as the potential cause of a new phenotype including early-onset osteoporosis and bone fractures: a case report | BMC Medical Genomics

Anamnesis vitae A 13 year old male born was as result of the VII pregnancy, from unrelated parents. Other pregnancies resulted in: I-II silent miscarriage in the second trimester; III – female, born in 2003 (III-3 Fig. 1) that has the following phenotypic features: genu valgum, hip dysplasia, combined thoracolumbar scoliosis,…

Continue Reading An FGFR2 mutation as the potential cause of a new phenotype including early-onset osteoporosis and bone fractures: a case report | BMC Medical Genomics

r – Fst calculation from VCF files

I have four vcf files, SNPs_s1.vcf, SNPs_s2.vcf, SNPs_s3.vcf, and SNPs_s4.vcf, which contain information about SNPs. These vcf files were obtained by using the following methods: the initial input files were short-paired reads I did mapping with minimap2 ./minimap2 -ax sr ref.fa read1.fq.gz read2.fq.gz > aln.sam converted to bam file samtools…

Continue Reading r – Fst calculation from VCF files

Building reference dbSNP file using WGS samples

Building reference dbSNP file using WGS samples 2 Dear scientific community, I have to call variants from WGS samples of citrus. I used GATK pipeline for post processing of aligned reads but reference dbSNP file is not available for citrus sinensis. I am using bootstraping method. Removed duplicates and called…

Continue Reading Building reference dbSNP file using WGS samples

Tool for structural variation detection of plant WGS data

Tool for structural variation detection of plant WGS data 1 Dear all, I have been calling SNPs and indels from WGS data but I have never performed structural variation calling. Please suggest any tool to identify structural variations from plants WGS data. Does freebayes facilitate this? Looking for your valuable…

Continue Reading Tool for structural variation detection of plant WGS data

Filtering for primary and secondary reads using sam flags (0 properly paired reads in alignment step)

Hey everybody, I have just performed the alignment of paired-end reads to a reference using bwa mem with the -M flag, ran samtools markdup and flagstat. Flagstat produced the following output: 4985084 + 0 in total (QC-passed reads + QC-failed reads) 1806492 + 0 primary 3178592 + 0 secondary 0…

Continue Reading Filtering for primary and secondary reads using sam flags (0 properly paired reads in alignment step)

Challenges in Variant Calling and Genotyping with Short-Read Data Mapped to a Pangenome Graph: Seeking Guidance

Challenges in Variant Calling and Genotyping with Short-Read Data Mapped to a Pangenome Graph: Seeking Guidance 0 Hello all, We are reaching out since we have some practical questions regarding variant calling and analyses of short-read data mapped to a pangenome graph. We are working on a project aimed to…

Continue Reading Challenges in Variant Calling and Genotyping with Short-Read Data Mapped to a Pangenome Graph: Seeking Guidance

Building reference dbSNP file for citrus sinensis using 80 WGS samples

Building reference dbSNP file for citrus sinensis using 80 WGS samples 1 Dear scientific community, I have to call variants from 80 WGS samples of citrus. I used GATK pipeline for post processing of aligned reads but reference dbSNP file is not available for citrus sinensis. I am using bootstraping…

Continue Reading Building reference dbSNP file for citrus sinensis using 80 WGS samples

Bootstrapping for BQSR of 80 WGS samples

Bootstrapping for BQSR of 80 WGS samples 0 Dear scientific community, I have to call variants from 80 WGS samples of citrus. I used GATK pipeline for post processing of aligned reads but reference dbSNP file is not available for citrus sinensis. I am using bootstraping method. Removed duplicates and…

Continue Reading Bootstrapping for BQSR of 80 WGS samples

Tracking SARS-CoV-2 Omicron lineages using real-time reverse transcriptase PCR assays and prospective comparison with genome sequencing

Population and clinical samples The province of Alberta has a population of approximately 4.4 million. During July 19 to December 31, 2022 (the period of this study), molecular testing for SARS-CoV-2 was performed for individuals at risk for severe illness who would benefit from early treatment with antivirals, healthcare workers,…

Continue Reading Tracking SARS-CoV-2 Omicron lineages using real-time reverse transcriptase PCR assays and prospective comparison with genome sequencing

GenomicStar hiring Bioinformatics Data Analyst in Chennai, Tamil Nadu, India

Position: Bioinformatics Data Analyst Company: GenomicStar – A SaaS based Precision Medicine Platform Job Purpose: The Bioinformatics Data Analyst will play a critical role in the development, validation, and implementation of data analysis pipelines for various genetic test panels, incorporating state-of-the-art AI techniques. This role is integral to our commitment…

Continue Reading GenomicStar hiring Bioinformatics Data Analyst in Chennai, Tamil Nadu, India

Solved We are now going to call variants with two different

We are now going to call variants with two different approaches from the files we have been working with all course. Please use the following files, parameters, and listed versions of the software for this assignment. We will use the reference Ebola genome: /data/compres/refs/AF086833.2.fasta And this set of paired-end sequences:…

Continue Reading Solved We are now going to call variants with two different

Mismatch repair deficiency is not sufficient to elicit tumor immunogenicity

Mice All animal use was approved by the Department of Comparative Medicine at the Massachusetts Institute of Technology (MIT) and the Institutional Animal Care and Use Committee under protocol no. 0714-076-17. Mice were housed with a 12-h light/12-h dark cycle with temperatures in the range 20–22 °C and 30–70% humidity. KrasLSL-G12D…

Continue Reading Mismatch repair deficiency is not sufficient to elicit tumor immunogenicity

sarek: Introduction

Introduction nf-core/sarek is a workflow designed to detect variants on whole genome or targeted sequencing data. Initially designed for Human, and Mouse, it can work on any species with a reference genome. Sarek can also handle tumour / normal pairs and could include additional relapses. The pipeline is built using…

Continue Reading sarek: Introduction

Virtual Variant Detection Workshop: September 11-14, 2023

News:Virtual Variant Detection Workshop: September 11-14, 2023 0 Virtual Variant Detection Workshop Dates: September 11-14, 2023 (4 days) Time: 9.00am – 12.00pm Location: Online Cost: $400 (UConn affiliates including UConn Health) $500(External Participants) The workshop will cover an introduction to linux and high performance computing, an introduction to variant detection,…

Continue Reading Virtual Variant Detection Workshop: September 11-14, 2023

Analysis of drug-resistance characteristics of MDR-TB

Introduction Tuberculosis (TB), caused by Mycobacterium tuberculosis (MTB), continues to pose a significant global public health challenge.1–3 In 2021, there were 10.6 million new TB cases worldwide, with China accounting for 780,000 cases, making it the third highest burden country for TB globally.4 The World Health Organization (WHO) estimates that…

Continue Reading Analysis of drug-resistance characteristics of MDR-TB

Problems with bcf concat

Problems with bcf concat 1 Hello! I have multiple VCFs files for the same sample. I used the commad bcftools concat -a -Oz -o concat.vcf.gz 242487_vs_242482.bcftools.vcf.gz 242487_vs_242482.freebayes.vcf.gz But i have the next error Checking the headers and starting positions of 2 files [W::hts_idx_load3] The index file is older than the…

Continue Reading Problems with bcf concat

The Biostar Herald for Thursday, August 24, 2023

The Biostar Herald publishes user submitted links of bioinformatics relevance. It aims to provide a summary of interesting and relevant information you may have missed. You too can submit links here. This edition of the Herald was brought to you by contribution from Istvan Albert, and was edited by Istvan…

Continue Reading The Biostar Herald for Thursday, August 24, 2023

VCF merge or concatenate?

I would like to double check whether to use VCF concatenate or VCF merge on my chromosome files. I have done SNP calling using FreeBayes, but split this by chromosomes in order to call SNPs in parallel. I also split some particularly large chromosomes by chromosome position with no overlap,…

Continue Reading VCF merge or concatenate?

Bioinformatics Analyst Job in Minnesota

The Research Informatics (RI) Bioinformatics group within the University of Minnesota Supercomputing Institute (MSI) is hiring a full-time Bioinformatics Analyst to support basic biomedical and applied clinical genetic research for the Department of Laboratory Medicine and Pathology (LM&P) at the University of Minnesota. The analyst in this position will…

Continue Reading Bioinformatics Analyst Job in Minnesota

Proper Bowtie2 settings for snp calling applications

Proper Bowtie2 settings for snp calling applications 1 Hi, I’d like to map illumina paired-end reads to the reference genome by Bowtie2 for downstream snp calling with Freebayes. The reads are already adapters and quality trimmed. I plan to specify the following Bowtie2 options: –end-to-end –very-sensitive –no-mixed –no-discordant I’m wondering…

Continue Reading Proper Bowtie2 settings for snp calling applications

26 All paired end Illumina libraries were inspected with FastQC v0115

26All paired-end Illumina libraries were inspected with FastQC v0.11.5() and trimmed and filteredwith Trimmomatic v0.36 (66). Read mapping was performed with sppIDer pipeline(81) using a combined reference including the genome assembly of CBS9595 (asrepresentative of P1 lineage, Supplementary File 1 for more details) and CBS14141(as representative of P2 lineage, Supplementary…

Continue Reading 26 All paired end Illumina libraries were inspected with FastQC v0115

vcftools – Fst statistics

vcftools – Fst statistics 1 I have .vcf file with over 78 000 variations in 95 experiment subjects. I used vcftools to exctract fixation index score for general differentiation between Case and Control groups. ./vcftools –vcf SCH_freebayes_sort_ostE_SnpEff_SnpSift.vcf –weir-fst-pop Case.txt –weir-fst-pop Control.txt –out FstExon_Case_vs_Control After filtering, kept 95 out of 95…

Continue Reading vcftools – Fst statistics

Unable to find FASTA index entry for ‘/mountpoint/fastQ/Homo_sapiens_assembly38-3.sorted.bed.gz’

Unable to find FASTA index entry for ‘/mountpoint/fastQ/Homo_sapiens_assembly38-3.sorted.bed.gz’ 1 I’m trying to use freebayes to call variants through the following parameters : freebayes -f $REF -r $REGIONS -b 03.align/$FQ.align.sort.marked.bam > 04.freebayes/$FQ.var.vcf The error reports: unable to find FASTA index for bed file. I have both .fasta and .fasta.fai file and…

Continue Reading Unable to find FASTA index entry for ‘/mountpoint/fastQ/Homo_sapiens_assembly38-3.sorted.bed.gz’

SNP analysis with an assembly

Hi there, I am new in SNP analyses so before starting doing anything I would like to check if my pipeline is correct. What I have now is : RNA-seq samples (.fq.gz) + Trinity assembly (from those reads). My model organism has not an assembled genome, that’s the way I…

Continue Reading SNP analysis with an assembly

Vancomycin intermediate-resistant Staphylococcus haemolyticu | IDR

Wanyang Dong,1,* Qi Peng,1,* Xiaohua Tang,1,2,* Tian Zhong,1 Shunan Lin,1 Ziling Zhi,1 Jingyi Ye,1 Bixia Yang,1 Ning Sun,1,3 Wenchang Yuan1 1Guangzhou Key Laboratory for Clinical Rapid Diagnosis and Early Warning of Infectious Diseases, KingMed School of Laboratory Medicine, Guangzhou Medical University, Guangzhou, 510180, People’s Republic of China; 2Third Affiliated Hospital…

Continue Reading Vancomycin intermediate-resistant Staphylococcus haemolyticu | IDR

Finding more effective ways to combat leishmaniasis by working with parasites

Schematic representation of the experimental design. a Experimental conditions. Three experimental conditions were tested: WT and HR strains growing without drug challenge, and HR strains growing under drug challenge. Four types of samples were evaluated per experimental condition (input, monosome, light polysomes, and heavy polysomes). The experiment was done in…

Continue Reading Finding more effective ways to combat leishmaniasis by working with parasites

Tinybio Slackbot – a large language model slackbot for bioinformatics

How is it different from typing the same into ChatGPT? And it does not take 2-3 minutes … Can you give me two variant callers? Sure, variant callers are bioinformatics software tools used to identify variants (like SNPs and indels) from sequence data. Here are two commonly used variant callers:…

Continue Reading Tinybio Slackbot – a large language model slackbot for bioinformatics

Could not get first alignment from target

Can you share some of the image as text for easier understanding? It seems like there might be an issue with your BAM file or the region you are trying to call variants on. To help diagnose the issue, please follow these steps: 1. Check if your BAM file is…

Continue Reading Could not get first alignment from target

python – Samtools shared library libcrypto.so.1.0.0 not found

I am trying to run a snakemake pipeline with packages installed from a main.yml file. One of the dependencies of the package is samtools, Previously, having the dependency listed as samtools did not seem to run into any problems. The actual yml file is: name: yevo_pipeline_env channels: – defaults -…

Continue Reading python – Samtools shared library libcrypto.so.1.0.0 not found

Variant Calling with Multiple Individuals in Freebayes

Variant Calling with Multiple Individuals in Freebayes 1 Hello, I am currently a new student in Bioinformatics and was trying to do variant calling with Freebayes software. I have 10 individuals and also have a reference genome. I have already sorted and indexed all my .bam files ready for variant…

Continue Reading Variant Calling with Multiple Individuals in Freebayes

How to fix FreeBayes error “unable to find FASTA index entry for bed files”

How to fix FreeBayes error “unable to find FASTA index entry for bed files” 0 I’m trying to use freebayes to call variants through the following parameters using RNA-Seq data: $freebayes –fasta-reference $ref \ –region $BED \ $bam \ ${outdir}/${sample}.freebayes.vcf and my bed files like this: chr1:11868-12227 chr1:12612-12721 chr1:12974-13052 chr1:13220-14501…

Continue Reading How to fix FreeBayes error “unable to find FASTA index entry for bed files”

Outgroup for plastome NJ-tree (from pseudoalignment based on SNPs)

Outgroup for plastome NJ-tree (from pseudoalignment based on SNPs) 0 Good morning, I have whole-genome sequencing data for several samples of one plant species. I performed mapping of all reads on the merged nuclear+plastid+mitochondrial reference genome of this species (bwa-mem2), extract part of the alignment, that was refered to plastome,…

Continue Reading Outgroup for plastome NJ-tree (from pseudoalignment based on SNPs)

Why 0/0 genotype refers to SNP and ALT != REF (in vcf file after freebayes)

Why 0/0 genotype refers to SNP and ALT != REF (in vcf file after freebayes) 1 Good morning, Could you please explain me, why 0/0 genotypes in vcf are SNPs? According to this description of GT: 0/0 the sample is a homozygous reference 0/1 the sample is heterozygous (carries both…

Continue Reading Why 0/0 genotype refers to SNP and ALT != REF (in vcf file after freebayes)

reference for freebayes or samtools mpileup after extracting chromosome from alignment

reference for freebayes or samtools mpileup after extracting chromosome from alignment 1 Good evening, I have extracted one chromosome from alignment map (.bam), using samtools view: samtools view -b map.bam chr1 > map_chr1.bam Now I would like to perform SNP calling using freebayes. Is it correct to use chr1.fasta as…

Continue Reading reference for freebayes or samtools mpileup after extracting chromosome from alignment

Navigating the Bioinformatics Workflow for Whole Exome Sequencing: A Step-by-Step Guide

Next-generation sequencing (NGS), which makes millions to billions of sequence reads at a fast rate, has greatly sped up genomics research. At the moment, Illumina, Ion Torrent/Life Technologies, 454/Roche, Pacific Bioscience, Nanopore, and GenapSys are all NGS platforms that can be used. They can produce reads of 100–10,000 bp in…

Continue Reading Navigating the Bioinformatics Workflow for Whole Exome Sequencing: A Step-by-Step Guide

SV tool recommandation for detecting 50bp> deletion sequence

SV tool recommandation for detecting 50bp> deletion sequence 0 Hello. I am a 1st year grad student in genetics lab. I am new to bioinformatics. computer: MS-Linux ubuntu-miniconda3 data: fish tissue gDNA NGS(MGISEQ,, 150 paried end read) picard bam, bai file My goals are 1. compare sample seq with NCBI…

Continue Reading SV tool recommandation for detecting 50bp> deletion sequence

Please supply a reference FASTA/GBK/EMBL file with –reference

Error: Please supply a reference FASTA/GBK/EMBL file with –reference 0 Hi, i am trying to run snippy on multiple genomes, however it gives following error like Please supply a reference FASTA/GBK/EMBL file with –reference even after providing the reference file. i don’t understand why it is happen and here is…

Continue Reading Please supply a reference FASTA/GBK/EMBL file with –reference

Bioinformatics Analyst in Minneapolis, MN for University of Minnesota, Twin Cities

Details Posted: 13-Aug-22 Location: Minneapolis, Minnesota Salary: 43944.32 – 123022.07 Categories: Research Support – Laboratory/Non-Laboratory Staff/Administrative Additional Information: 2 openings available. The Research Informatics Solutions (RIS) group within the University of Minnesota Supercomputing Institute (MSI) is hiring two full-time Bioinformatics Analysts to support research at the University of Minnesota. Analysts…

Continue Reading Bioinformatics Analyst in Minneapolis, MN for University of Minnesota, Twin Cities

Freebayes-parallel with large bam file – individual threads running for >6 days

Context: I’m trying to call variants on a sequencing project using pooled genotyping-by-sequencing. Pools consist of 94 samples each, alongside a number of individuals. Sequence data was demultiplexed and then aligned to a reference genome using hisat2, and the resultant bams were merged with samtools merge. The problem bam is…

Continue Reading Freebayes-parallel with large bam file – individual threads running for >6 days

snp – Reference variant detected as altered one in bam file

I received (from manufacturer) several .bam files and I used four callers (samtools, freebayes, haplotypecaller, deepvariant) to find some sequence variants. In obtained .vcf files, I took a closer look to some calls. I found interesting, homozygous one rs477033 (C/G Ref/Alt) with flag ‘COMMON=0’ and very low MAF. I also…

Continue Reading snp – Reference variant detected as altered one in bam file

[moiexpositoalonsolab/grenepipe] freebayes causes early error about number of threads

Hi Lucas, got a weird one for you. If I change the caller from hapotypecaller to freebayes, I get the error below. It’s doubly strange because it seems to occur well before freebayes would be used in the pipeline. [Sat Dec 11 11:13:02 2021] rule samtools_stats: input: dedup/111D03-1.bam output: qc/samtools-stats/111D03-1.txt…

Continue Reading [moiexpositoalonsolab/grenepipe] freebayes causes early error about number of threads

Phylogeographic reconstruction of the marbled crayfish origin

Procambarus fallax collections and PCR genotyping Animals were collected from various wild populations (Table S1) in compliance with state and local regulations (Georgia department of natural resources scientific collection permit 115621108, state of Florida collection permits S-19-10 and S-20-04). DNA was isolated from abdominal muscle tissue using SDS-based extraction and precipitation…

Continue Reading Phylogeographic reconstruction of the marbled crayfish origin

High frequency of an otherwise rare phenotype in a small and isolated tiger population

Significance Small and isolated populations have low genetic variation due to founding bottlenecks and genetic drift. Few empirical studies demonstrate visible phenotypic change associated with drift using genetic data in endangered species. We used genomic analyses of a captive tiger pedigree to identify the genetic basis for a rare trait,…

Continue Reading High frequency of an otherwise rare phenotype in a small and isolated tiger population

Exctracting amino acid substitutions

Exctracting amino acid substitutions 0 Good day, I’m trying to develop a pipeline to determine mutations which are responsible for amino acid changes in genes associated with antibiotic resistance. I have roughly 300 bacrtial isolates. My approach so far has not been fruitful, in short this is what i tried:…

Continue Reading Exctracting amino acid substitutions

FreeBayes VCF output with FORMAT unknown

Hey, I am looking for a way to add samples ID names to the FORMAT in my vcf file. I have 10 sorted Bam files. I used Freebayes to create vcf files and my next step is merging all 10 files for VcfSampleCompare. And for that I need to define…

Continue Reading FreeBayes VCF output with FORMAT unknown

FreeBayes Error

FreeBayes Error 0 Hi- I am trying to use FreeBayes to identify polymorphic sites and it will not run. It keeps throwing this error “unable to find fasta index entry for GL630349”. Does anyone know what to do in this situation? Thank you Galaxy Coursera FreeBayes Polymorphisms • 16 views…

Continue Reading FreeBayes Error

VCF Filter On Small Genomes

VCF Filter On Small Genomes 0 Hi guys, I am working on a yeast species (Candida glabrata) NGS data to find any mutations related to drug resistance. I am new in bioinformatics so I am using Galaxy.eu to get use to algorithms. There is literature about some genes that mutations…

Continue Reading VCF Filter On Small Genomes