Minimap2: supplementary alignment PAF
Hi!
I am trying to parse PAF file generated from minimap2.1 to get primary and supplementary alignments of reads but excluding secundary alignments. In SAM format I would search for flags: 0, 16, 2048 and 2064. Unfortunately I don't currently find a nice way to get the same information out of PAF file.
Does anyone knows how to retrieve supplementary/primary but not secundary alignments from PAF file?
1 answer
Hi,
According to minimap2 man page, the "tp" tag in the thirteenth column indicates the type of alignment: either primary, secondary or inversion.
Therefore, if you want only primary or supplementary, you can simply filter as follow :
awk '($13~"tp:A:P") {print $0}' file.paf > filtered_file.paf
Log in to answer this question.
More posts like this
-
Samtools: exclude secondary and supplementary alignments
written by Alexandra •
-
Extracting complete supplementary alignments for long reads from a bam file
written by aacomeau •
-
Reads with highest MAPQ values from SAM files are showing mismatches to reference sequence and IGV …
written by Mo •
-
Very high coverage Nanopore alignment Hg38
written by pablo
-
question regarding sequencing alignment
written by Lin •
-
[pysam] How to find both both primary and supplementary read？
written by octpus616 •
-
minimap2 aligning sequencing reads to barcodes
written by matthew.tarnowski1 •
-
Multiple mapping with longreads - Get primary mapping
written by ricardoguerreiro2121
-
Difference between chimeric alignments and multiple mapping
written by Vanilla
-
Solution To Find Out Read Orientation As Well As Strand To Which Reads Align From Sam File
written by k2bhide
