Bash Script for a specific SNP
Beginner here. Trying to look for a specific allele presence/absence in a group of fasta files for SARS-CoV-2 genomes. Also not sure if I should use the first/original sequence out of Wuhan as the reference genome or the first genome sequenced with the allele of interest to call the SNP. I am writing this just with bash, I don't really know R or Python enough to use it. Any help, tips, tricks much appreciated!
1 answer
Not sure what kind of data you are starting with but there are well established pipelines for this analysis. You should use one of those if you can.
SARS Cov2 NGS data analysis pipeline
COVID-19 analysis pipelines
Log in to answer this question.
More posts like this
-
What do solid red/green reads in IGV mean (looking at SARS-CoV-2 BAM)
written by Cornelius •
-
Get cleaved SARS-CoV-2 protein sequences (all major variants)?
written by ngarber •
-
vg for gene presence-absence analysis
written by liorglic
-
Binding site mutation
written by Mohamed Elsherbieny •
-
Forum: GISAID must be called to release EpiCoV data
written by diego1530 •
-
how to understand this "Jaccard distance" ?
written by 2001linana •
-
How to understand single nucleotide polymorphism (SNP)?
written by 2001linana •
-
Selection of Multiple Sequence Alignment Tool
written by pfee418 •
-
If I've identified the gene via BLAST, can I translate it to protein and assume its correct?
written by christinejgu •
-
large deletion detection in a haploid microbial eukaryote
written by thomas.welch •
