BLASR: How to subsample alignments from a region

Hello!

I have an alignment of ~20x PacBio long reads to my reference genome and I am trying to extract alignments within only a certain region, say 2mb-10mb (sort of like samtools view ) but I am not sure how to do this.

I have alignments in the m4 format.

I think I should be specifying a region table but I am not sure.

Can anybody help?

pacificbiosciences blasr alignment pacbio

• 2,827 views • link

updated 4.1 years ago by Ram

• written 10.7 years ago by shwethacm

Thanks, rhall.

I did eventually parse it myself but I was looking to check if I could be lazy about it and do it quickly using a preexisting tool :)

Well, now I know!

Thanks a ton anyway.

Shwetha

• 0 views • link

updated 6.7 years ago by Ram

• written 10.7 years ago by shwethacm

1 answer

I'm not aware of any software to process m4 files, but they are straightforward to parse with simple scripts, https://github.com/PacificBiosciences/blasr/wiki/Blasr-Output-Format. If you are more comfortable with sam format, blasr can output sam instead of m4, the -sam parameter.

• 0 views • link

updated 6.7 years ago by Ram

• written 10.7 years ago by rhall

More posts like this

How to extract a multiple sequence alignment file from a CRAM file?
written 1.4 years ago by Charles Plessy

Hello everyone, I have a multiple query genomes aligned to a target genome stored in a CRAM file where each query genome is represented by …
Samtools: exclude secondary and supplementary alignments
written 1.9 years ago by Alexandra •

Hello! I have a question regarding how samtools -F handles certain flags. I have bam files from RNA-seq and wish to exclude supplementary and secondary …
Samtools view -e not work
written 3.6 years ago by octpus616 •

Hi, I recall that samtools provided a powerful expression based filtering feature `samtools view -e 'filter_expression'`, which seemed to work fine in the past, but …
Extract paired-end reads with one end falls within a given region from a bam file
written 3.9 years ago by Wang Cong •

I have paired-end sequencing reads and two ends of a read may not come from the same region. I want to extract those that have …
export alignment from BAM file
written 5.1 years ago by luca •

Dear all, I'd like to export an alignment in a fasta (or similar) format from a bam file. I know there are tools to extract …
How to trim alignments in a BAM file to a restricted genomic window?
written 5.3 years ago by francois •

Say I have long reads which I have aligned to hg38. Now I want to trim my alignments to a specific region of chr20, say …
Why Samtools view don't extract sequence out of the given region?
written 6.3 years ago by kamanovae •

Hi! I have nanopore sequencing results in fastq format. I need to know the coverage of specific genes. To do this, I aligned the sequence …
stampy "stuck" aligning
written 10.3 years ago by AP

Hi all, I am using stampy for aligning against a reference genome but it seems to get stuck during the alignment process. For a single …
Samtools And Flagstat
written 13.0 years ago by nbvasani

Hello All Biostar community user, I am trying to use flagstat command option in samtools. Usually people get below output when they type `samtools.pl` in …
Why Import Command In Samtools Is Not Shown In The Short Help Of Samtools?
written 15.1 years ago by Pablo Marin-Garcia

<p>I am converting sam to bam using 'samtools import' command.</p> <p>I am curious why if I do 'samtools' I can not see the 'import' command?</p> …