Hello all.
I have created BAM files for 10 samples . My goal is to identify SNPs that are shared among samples (or by a subset of my sample) , as well as those that are sample-specific . Ultimately, I'd like to identify the genotype (haplotype) of each sample, and reconstruct their relationship.
The main issue I am having is that coverage is highly uneven among samples. For example, a region X may be covered in only 4 samples, while a region Z may be covered only in 3).
I am looking at extracting all regions that are covered (1X and over, or from 2X to 50X) in all samples (or a subset of samples), and then use these to call SNPs and extract genotypes.
Any help on this matter will be greatly appreciated. Thank you very much!
Best, Salmo
0 answers
No answers yet.
Log in to answer this question.
-
Filter multisample vcf by group
written by Theo •
-
When identifying DEGs, should all the libraries be normalized together?
written by Estefania •
-
How to create a phylogenetic tree from 30 VCF files
written by lgspeight •
-
How to Identify set of associated SNPs with single SNP using 1001 genomes
written by ES.67 •
-
How to report coverage at overlapping regions of paired-end reads?
written by francois •
-
How to evaluate "control" genes in target sequencing data?
written by jinxinhao1988
-
Heterozygous genome assembly software: actions on regions with copy number more than 2X?
written by dmathog •
-
pairwise calculation of relationship based on selected SNPs with MAF>0.3
written by J.F.Jiang
-
Haplotype frequencies from 1000 genomes
written by whittlemr •
-
Haplotype Frequency Calculator For Vcf Files
written by diviya.smith •
