Thanks! I will try it
Hi,
I have a list of mutations and their genomic positions.
for example:
chr1:1573209-1573209 T>A
and I have fasta sequences of the consensus sequence around the sites of these mutations.
for example:
>hg19_ct_UserTrack_3545_CDK11A range=chr1:1573157-1573247 5'pad=10 3'pad=10 strand=+ repeatMasking=lower
GAGGATGGTGTTGATCTCCCTCAGCGACGTGATCGGGAAGCCCTCCTTCT
CCTTCTCCATCTTCAGCCGCTTTAGAGCCACAATTTCATCT
Is there any tool that can insert the mutation into the fasta?
and more complicated - a tool that can combine several different mutations in the same fasta?
I wrote a perl script that creates a fasta for each mutation, but now I need to create all possible combinations of mutations (i.e. if there are 3 different mutations in a single fasta sequence, I need all the combinations: mut1 + mut2, mut2 + mut3, mut1 + mut2 + mut3), and I wonder if there is an existing tool that can do that.
Thanks in advance for any help
2 answers
Just convert the custom variant format to VCF and then use GATK's FastaAlternateReferenceMaker.
Also, you could use SimNGS: Using Simngs To Generate Reads
Log in to answer this question.
Hello lilachfr!
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!