How to perform in R the variant calling on .sam and. bam input done with BBMap?
How to perform in R the variant calling on .sam and. bam input done with BBMap? Is it a proper way or it is better to use Freebayes, BALSA, GATK or SAMTools for that?
alignment
variant
r
• 2,198 views
•
link
written
by
usgenome2018 •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
allele matching between identical genotypes
written by rj.rezwan •Hi, I have the variant calling output in the form of `vcf` files. Variant calling was done with `GATK haplotypecaller`. I am interested in performing …
-
reference for freebayes or samtools mpileup after extracting chromosome from alignment
written by poecile.pal •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 …
-
samtools variant calling workflow
written by mirzaHi, I am working on variant calling on fungal genomes. I have Illumina HiSeq reads. I am new to this and am following this workflow …
-
Is it necessery to perform GATK-BQSR before performing variant calling using freebayes.
written by sktbanerjee1 •I am trying to call variants from whole exome paired end sequences aligned to HG38p7 using the Freebayes variant caller. As it is a general …
-
workflow for variant calling using GATK in multiple paired end reads
written by rezaI have four fastq file that correspond two paired end lane (lane1: L1R1.fastq, L1R2.fastq and lane2:L2R1,L2R2). Mapping to reference performed using below command via “bwa …
-
What R packages and functions can perform variant calling on .sam and .bam files?
written by elvissober •<p>What R packages and their functions can perform variant calling on .sam and .bam files? Need a commented clear code example. Found RSamtools. No concrete …
-
How to do variant calling on .sam and .bam files using 'bio-samtools'?
written by elvissober •How to do variant calling on .sam files using 'bio-samtools'? I have the code: ``` #!/usr/bin/ruby require 'bio-samtools' bam = Bio::DB::Sam.new(:bam=>"my_sorted.bam", :fasta=>'ref.fasta') bam.open ``` What …
-
Filter/Annotate VCF with dbSNP during variant calling workflow
written by alonsHi all, I would like to filter / annotate my VCF and/or BAM files using dbSNP as part of the input during the steps of …
-
Issue of using freebayes to do the variant calling
written by jz6002 •Hi, this is my first time to use freebayes. I met a problem when I tried to use it to do the variant calling for …
-
Snp Detection In Viruses: Does A Haploid Caller Like Freebayes Need To Be Used Instead Of Something…
written by JacobS<p>I've done variant calling for many eukaryotic organisms, and always followed a simple, highly vetted pipeline of mapping my reads to a reference, generating pileup …
Hello usgenome2018!
Questions similar to yours can already be found at:
We have closed your question to allow us to keep similar content in the same thread.
If you disagree with this please tell us why in a reply below. We'll be happy to talk about it.
Cheers!