Question: Chimeric alignment with minimap2
0
gravatar for Nicolas Rosewick
2.2 years ago by
Belgium, Brussels
Nicolas Rosewick8.6k wrote:

Hi,

On minimap2 github (https://github.com/lh3/minimap2) it says :

For long RNA-seq reads, minimap2 may produce chimeric alignments potentially caused by gene fusions/structural variations or by an intron longer than the max intron length -G (200k by default). For now, it is not recommended to apply an excessively large -G as this slows down minimap2 and sometimes leads to false alignments.

How can I extract such chimeric reads ? Are these only cis-fusions (on the same chromosome) or also trans fusions (different chromosomes) ? Is there a way to tune the chimeric detection ?

Thanks

nanopore chimeric minimap2 • 1.6k views
ADD COMMENTlink modified 2.2 years ago by Macspider3.0k • written 2.2 years ago by Nicolas Rosewick8.6k
0
gravatar for Macspider
2.2 years ago by
Macspider3.0k
Vienna - BOKU
Macspider3.0k wrote:

If the output is a SAM-formatted file, you can retrieve them by including only reads that carry the 0x800 flag, using samtools view -f 0x800.

ADD COMMENTlink written 2.2 years ago by Macspider3.0k
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: 838 users visited in the last hour