To call variants can I use my aligned WGS data as a reference genome
Hi,
I have several yeast WGS data for different time data points. I wanna call variants depending on the first time point WGS data(I already map my samples with bowtie). But I could not figure out how to use my first time point WGS data as a reference rather than S288C fasta. I would be glad if someone can give me some ideas.
Thank you
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Deseq2 5 level condition - building contrast
written by annaA •
-
Different blast results between NR or local genome
written by pablo
-
Labels in Enhanced Plot
written by nishant.asw3 •
-
Strand for GTEx WGS data?
written by spitzer •
-
Recommendations for Whole Genome Sequencing Workflow
written by jrleary
-
Minet works for time-course expression data
written by joseph.houjue •
-
sciClone + clonEvol on WES
written by Kasthuri
-
combining fasta files
written by br.tania •
-
Finding Repeat Genes In Unused Reads
written by Lee Katz
-
List Of Genes With Ires' In Yeast?
written by Jason
the bam are mapped on S288C fasta isn't it ? why would you need another reference ?
Yes they are but I wanna compare one of my samples with others , rather S288C.
Then create an alternate Reference with data for the sample: https://gatk.broadinstitute.org/hc/en-us/articles/360037594571-FastaAlternateReferenceMaker