Hi, I have two sequences and I would like to generate indels for them. As the alignment just gives the score, to represent the the alignment of the two sequence I would like to generate CIGAR .
Hi, I have two sequences and I would like to generate indels for them. As the alignment just gives the score, to represent the the alignment of the two sequence I would like to generate CIGAR .
Alternatively, just add some in manually if you just have two sequences and you don't need many indels. You can just randomly pick some positions.
So in other words, you want to align two sequences and get the resulting CIGAR string? Are you aligning them manually or with a program. In the latter case, just use a program that will generate CIGAR strings. In the former case, you can just generate one yourself (CIGAR strings aren't exactly complicated things).
Then use a local alignment program that will produce a CIGAR string (or at least something that you can then convert into one). Alternatively, if the programs output a formatted pairwise alignment, then just parse that (the only hard parts are the ends, which can be soft-clipped).
Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Just open the sequences with a text editor and insert/delete some bases? Maybe I'm missing something in your question...