In regard to 7 yrs old discussion about aligning linear to circular sequences - is anyone aware of aligners that can handle such alignment properly? For alignment purposes, for example, CLC-Bio Genomic Workbench does not seem to be aware if sequence is circular, although it can show rCRS mtDNA as circular. Yet it cannot properly align to it a linear sequence that is should align between pos. 16022-origin-1280 to circular rCRS. Nor can Geneious, for another example.
•
written
by
yaximik •
2 answers
I know this is an old post - but I wanted to clarify that CLC Genomics Workbench does indeed support mapping reads to circular references. The issue is that (currently) our track viewer is only linear - something we plan to improve on in 2020 as a number of customers have requested CIRCOS like viewers (for a variety of use cases).
Thus - you can indeed map (correctly) to a circular reference in CLC, and reads that align to the origin will be correctly mapped. But when you export them in SAM/BAM format (for example to visualize in a tool that supports circularized coverage maps) you run into issues with the SAM format not supporting circular references.
We have a FAQ on this specific topic actually (mapping in CLC, exporting to BAM, and reimporting into CLC you'll see that there's a loss of information due to the SAM/BAM format not supporting these origin spanning reads correctly. however, if the original CLC data in the same view is correct). Link is here: https://secure.clcbio.com/helpspot/index.php?pg=kb.printer.friendly&id=11#p419
Happy to answer any other questions on this.
Checkout my referral link
•
written
by
6mrsbrun •
Log in to answer this question.
-
Best Tool for Fast RNA-Seq Transcriptome Alignment to Extract Unmapped Reads
written by aboll •
-
Multiple sequence alignment of mtDNA HV1-3 in R
written by irieljoerin •
-
Aligning multiple overlapping DNA sequence reads to predicted sequence not working
written by samuelhaysom •
-
Circular Genome Alignment
written by isaac.rstefan •
-
Efficiently aligning a lot of reads on the same small reference sequence
written by cedric.stroobandt •
-
Any alignment programs that accept regular expressions within input sequences?
written by Ghoti •
-
Splice aware aligner - what does it mean?
written by jonasmst
-
Find genome given accession number or gene symbol
written by likethewind89 •
-
Clustalw Ignoring Custom Weight Matrix
written by interfect •
-
Good Multiple Sequence Alignment Tool For Ordered/Disordered Proteins
written by miquelduranfrigola
Since both of those are commercial programs have you contacted their respective tech supports about this?
Yes, I guess they are scratching pumpkins... I just wanted to know if anyone in the community is aware of other possibilities, like open source...
Is it NGS data ? Because i did it on times on NGS data with bwa and to make this possible i just merge 2 references together ....
No, this is just alignment of a contig previously assembled separately. CLC Bio developers confirmed that their Genomics Workbench aligner cannot handle circularity. As workaround they suggested either to use their Map to a Reference tool, which can handle circularity, or to align to a duplicated circular sequence.
Please use
ADD COMMENT/ADD REPLYwhen responding to existing posts to keep threads logically organized.BTW, have you looked at this: Aligning Circular Sequences (old thread but has ideas/tools).
Have you considered Circlator?