Question: Extracting Indels from two aligned DNA sequences
0
gravatar for manishagupta293
10 months ago by
manishagupta2930 wrote:

How i can extract Indels out from two aligned DNA sequences using python or biopython?

alignment • 341 views
ADD COMMENTlink modified 10 months ago by Hussain Ather890 • written 10 months ago by manishagupta2930
0
gravatar for Hussain Ather
10 months ago by
Hussain Ather890
National Institutes of Health, Bethesda, MD
Hussain Ather890 wrote:

SeqFIRE does what you want https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3394284/

Otherwise you could use samtools view to view them.

samtools view -h in.bam  \
| awk '$1 ~ "^@" || $6 ~ "I|D"' \
| samtools view -b - > out.bam
ADD COMMENTlink modified 10 months ago • written 10 months ago by Hussain Ather890

SeqFIRE use protein alignments. But we want to use DNA alignment only. Is their a way to extract indels using BioPython?

ADD REPLYlink written 10 months ago by manishagupta2930
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 804 users visited in the last hour