Hi there, I'm working on an allele skew project where we are utilizing two mouse strains. In order to perform some QC, we would like to extract strain specific SNP info and we would like to know the SNP position based on the strain of interest, not the reference mm10. Therefore, we were thinking about performing a liftover to "convert" the coordinates. Anyways, I identified GATK LiftoverVCF as a possible tool but it requires a chain file that's seemingly not available. Would anyone have any suggestions for either making a strain specific chain file or a different tool that doesn't require a chain file? Thank you!
The specific question that we want to answer: How many uniquely mapped reads overlap with a SNP position of each strain?
1 answer
The method I decided to use is the Coordinate remapping tool available through NCBI. It's very powerful and removes the need for a chain file. I may still end up using flo but here's the those two links:
Coordinate Remapping: https://www.ncbi.nlm.nih.gov/genome/tools/remap
flo GitHub: https://github.com/wurmlab/flo#using-flo
Log in to answer this question.
-
Liftover from mm10 to CAST/PWK – Very Low Mapping Rate & Missing Chain File
written by daffodil •
-
Comparing and assessing the quality of gene annotations between two mouse strain genomes
written by TheWonderer •
-
Converting Genome Reference Consortium to UCSC, i.e. GRCh* -> hg*
written by DavidStreid •
-
Annovar: Mouse database download
written by v.chimankar •
-
LiftOver a VCF file
written by Susmita Mandal
-
Using GATK4 Mutect2 on mouse data, need a genome snp reference
written by jc.szamosi •
-
assembly of mouse transcriptome from BAM file
written by regulatinggeneexpression •
-
Best tools for lifting over genome coordinates for non model organism (custom) reference genomes
written by William
-
Feasibility of functional Impact prediction of mouse variants based on liftedOver mm10 to hg19 co-o…
written by Parthav Jailwala •
-
Identify Virus Subpopulations with 454 Sequencing
written by pld
I'm going to try and use https://github.com/wurmlab/flo#using-flo as a method for creating the required chain file. Let me know if there are other suggestions.