Question: Chimeric alignment with minimap2
0
gravatar for Nicolas Rosewick
16 months ago by
Belgium, Brussels
Nicolas Rosewick7.5k 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.1k views
ADD COMMENTlink modified 16 months ago by Macspider2.8k • written 16 months ago by Nicolas Rosewick7.5k
0
gravatar for Macspider
16 months ago by
Macspider2.8k
Vienna - BOKU
Macspider2.8k 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 16 months ago by Macspider2.8k
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: 1391 users visited in the last hour