Category: beagle

Hardy-Weinberg Equilibrium and MAF filtering post imputation

Hardy-Weinberg Equilibrium and MAF filtering post imputation 0 Hi, I have 5 million variants after imputation using Beagle 5.2 after filering low confidence imputation ( DR2 < 0.9). Pre-imputation QC was that I removed rare variants MAF < 0.05 and use only variants with no missing genotype data (–geno 0…

Continue Reading Hardy-Weinberg Equilibrium and MAF filtering post imputation

Using a phased reference assembly to phase sample genotypes

Using a phased reference assembly to phase sample genotypes 0 I’ve been trying without success to phase the variants of about 30 samples in a 5Mb non-recombining region. I have a phased genome assembly of the region but so far I have not been successful in using that to inform…

Continue Reading Using a phased reference assembly to phase sample genotypes

map file for BEAGLE imputation

Forum:map file for BEAGLE imputation 0 Hello Community, I am having an issue with BEAGLE V5.2 throwing an error stating all my loci have the same genetic position when they most certainly do not. I have posted the error and a sample of the top of the map file I…

Continue Reading map file for BEAGLE imputation

No VCF records found in the specified interval

Beagle 5 error: No VCF records found in the specified interval 0 Hi, I am running into an issue while doing Imputation with Beagle 5 and not sure what is causing the error. I have vcf files converted from PLINK by the following command ./plink –bfile qcd_in–chr 20 –recode vcf-iid…

Continue Reading No VCF records found in the specified interval

The genome of Shorea leprosula (Dipterocarpaceae) highlights the ecological relevance of drought in aseasonal tropical rainforests

Sequencing of Shorea leprosula genome Sample collection Leaf samples of S. leprosula were obtained from a reproductively mature (diameter at breast height, 50 cm) diploid tree B1_19 (DNA ID 214) grown in the Dipterocarp Arboretum, Forest Research Institute Malaysia (FRIM). DNA extraction Genomic DNA was extracted from leaf samples using the…

Continue Reading The genome of Shorea leprosula (Dipterocarpaceae) highlights the ecological relevance of drought in aseasonal tropical rainforests

Cenozoic climatic changes drive evolution and dispersal of coastal benthic foraminifera in the Southern Ocean

1. Thomson, M. R. A. Geological and palaeoenvironmental history of the Scotia Sea region as a basis for biological interpretation. Deep Sea Res. II 51, 1467–1487 (2004). Article  ADS  Google Scholar  2. Maldonado, A. et al. A model of oceanic development by ridge jumping: Opening of the Scotia Sea. Glob….

Continue Reading Cenozoic climatic changes drive evolution and dispersal of coastal benthic foraminifera in the Southern Ocean

Cool Bioinformatics Scripts With Python

You can use this script in two ways read tons of millions of P values from stdin # python zcat pval.txt.gz | qqplot.py -out test -title “QQ plot on the fly” # julia zcat pval.txt.gz | qqplot.jl –out test –title “QQ plot on the fly” warning : If you have…

Continue Reading Cool Bioinformatics Scripts With Python

GitHub – Zilong-Li/BioScripts: Cool Bioinformatics Scripts

You can use make a QQ plot in the following ways. one-liner for reading tons of millions of P values from the pipe # python zcat pval.txt.gz | qqplot.py -out test -title “QQ plot on the fly” # julia (recommand to run it in the REPL) zcat pval.txt.gz | qqplot.jl…

Continue Reading GitHub – Zilong-Li/BioScripts: Cool Bioinformatics Scripts

Cool Bioinformatics Scripts | PythonRepo

You can use this script in two ways read tons of millions of P values from stdin # python zcat pval.txt.gz | qqplot.py -out test -title “QQ plot on the fly” # julia zcat pval.txt.gz | qqplot.jl –out test –title “QQ plot on the fly” warning : If you…

Continue Reading Cool Bioinformatics Scripts | PythonRepo

I can’t get a dossage file using PLINK

Hi, I have been trying to get a dosage file from vcf, map and fam files. For that, I have written this bash script : plink –fam plink.fam –map plink.map –dosage one.vcf –write-dosage However, I got this error: –dosage: Reading from one.vcf. Error: Line 1 of one.vcf has fewer tokens…

Continue Reading I can’t get a dossage file using PLINK

Beagle 5.2 Imputation Issue, understanding output

Hey everyone, I’m having some trouble with the imputation step in my GBS pipeline. I’m following the FastGBS pipeline basically as written up until this point. Before running Beagle, i’ve performed a filtering step using vcftools (though I get the same error whether or not I do this; I just…

Continue Reading Beagle 5.2 Imputation Issue, understanding output

Bioinformatics Scientist in Frederick, MD

Job DescriptionBioinformatics ScientistFull Time Direct Hire Remote positionAre you looking for bioinformatics work? Are you interested in joining a team of talented bioinformaticians dedicated to understanding the genetics of cancer? In this role you will:* Function as a scientific thought leader within for all aspects of GWAS and population genetics….

Continue Reading Bioinformatics Scientist in Frederick, MD

Beagle 5.2 for imputation, outofmemory or ‘bus error’ any help?

Beagle 5.2 for imputation, outofmemory or ‘bus error’ any help? 0 Hi, I was trying to use beagle 5.2 for imputation of genotype data (N=62k). I was using chr21 with chunks of 3MB (total 12 chunks running in parallel). This tool is very fast. The problem is every time I…

Continue Reading Beagle 5.2 for imputation, outofmemory or ‘bus error’ any help?

phase_trio.sh | searchcode

phase_trio.sh | searchcode PageRenderTime 24ms CodeModel.GetById 16ms app.highlight 5ms RepoModel.GetById 1ms app.codeStats 0ms /Phase/phase_trio.sh github.com/BioinformaticsArchive/fCNV Shell |…

Continue Reading phase_trio.sh | searchcode

Beagple 5.2 phasing error

Beagple 5.2 phasing error 0 Hi everyone, I’m trying to phase a multi-sample (12 samples) vcf file with the first chromosome. I got this vcf after pruning with plink and recode it back to vcf. The file looks like this: 1 112 . C T . . PR GT ./….

Continue Reading Beagple 5.2 phasing error

Error when Phasing with Beagle 5.2

Error when Phasing with Beagle 5.2 0 I’m having trouble phasing a multi-sample (9-samples) vcf file produced by gatk HaplotypeCaller with Beagle 5.2. I do not have a genetic map or reference panel. I am working with a very heterozygous group of organisms (sea urchins). When I run beagle with…

Continue Reading Error when Phasing with Beagle 5.2

linkage disequilibrium and haplotype analysis of GWAS .

  linkage disequilibrium and haplotype analysis of GWAS . 0   Hi all, I have GWAS data. I have my data in 22 chromosome files in plink format. I have imputed genotype with Sanger imputation server. I use plink for my analysis but because plink 1.9 no more supports –hap…

Continue Reading linkage disequilibrium and haplotype analysis of GWAS .