This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Find mismatch and INDEL position

HI everyone,

I have a short sequence, which has been mapped to the genome with known numbers of mismatches and INDELs. Below shows an example:

target sequence: ATGGATCGTACGTAACGCTAGCATGACT, query sequence: ATGCGTCCGCGTAACG, mismatch No.: 6, INDEL No.: 3.

based on the above information, how to generate the following pattern, in other words, how to output the mismatch and INDEL position information?

ATG-GATCGTACGTAACGCTAGCATGACT
||| | ||    || |||
ATGCG-TCCGCG-TAACG

Thank you very much.

mismatch script indel
• 4,563 views • link
updated by Ram Ram • written by biolab biolab

Ok, so when you say how to output do you mean 'how to represent them properly' or 'how to write a code to represent them the way you displayed'? I think I got a bit confused there.

• 0 views • link
written by Jordan Jordan

Hi Jordan,

I mean "how to represent them properly?" I know little about CIGAR string. Any tools that helps align my short sequence to the genome and output mismatch and INDEL information will be useful to me. But one important thing is that I would allow fuzzy match. Thanks!

• 0 views • link
updated by Ram Ram • written by biolab biolab

I see. You can use BLAST and choose the option Align two or more sequences.

You can put in query and target sequence in fasta format. For e.g., in your case:

>query_seq
ATGCGTCCGCGTAACG

>target_seq
ATGGATCGTACGTAACGCTAGCATGACT

You need to give your match, mismatch and gap scores. Any tool which does alignment needs them.

This will display the type of alignment you mentioned in your question. Not a CIGAR string though.

• 0 views • link
updated by Ram Ram • written by Jordan Jordan

1 answer

You mean something like a CIGAR string?

3M1I1M1D2M4X1D5M
  • M - Match
  • X - Mismatch
  • I - Insertion
  • D - Deletion
• 0 views • link
updated by Ram Ram • written by Jordan Jordan

Log in to answer this question.

More posts like this