Show alternative alignments in bam or sam files

Hi All,

I want to check if there are alternative alignments in my results via BWA alignment.

I have both bam and sam files, but I don't want to download the complete sam files to my computer.

I know there is a way to just view the results with alternative alignments.

Should I use samtools or should I create a vcd file to view the alternatives?

Can I just use

vcftools view -A input.file > output.file

Any ideas?

Thanks

alignment

• 2,739 views • link

updated 4.1 years ago by Ram

• written 11.3 years ago by Crystal •

You can extract alternative (non-primary) alignments from the bam file using: samtools view -f 100 input.bam

• 0 views • link

updated 4.1 years ago by Ram

• written 11.3 years ago by Ashutosh Pandey

Thank you.

I did samtools view -f input.bam > output. It seems the output is a txt file. It provides information, but I'm confused with the format. So I'm wondering if I can get the output as a bam or sam file and then I can convert it to a readable txt file.

• 0 views • link

updated 4.1 years ago by Ram

• written 11.3 years ago by Crystal •

samtools view -f 100 -h > Output.sam

for the sam format with the header OR

samtools view -f 100 -bh -o Output.bam

for the bam file.

• 0 views • link

updated 4.1 years ago by Ram

• written 11.3 years ago by Ashutosh Pandey

0 answers

No answers yet.

More posts like this

bcftools view to failed reader data
written 3.0 years ago by Jihyun •

It's so confused for data editing with vcftools - bcftools fusion preprocessing. I want to view all results in bcftools, vcftool and zcat when I …
Any option in bwa mem aligner for headless sam output?
written 3.9 years ago by robjohn70000

Hi, In the bowtie2 implementation for read alignment to a reference, there is an option (**--no-sq** Suppress @SQ SAM header lines.) to obtain a headless …
how much a mapping quality will HISAT2 give to a unique alignment?
written 7.7 years ago by Zeason •

I just want to find some unique alignment from my sam file with samtools. so i just want to know which figure should be added …
Automating a SNP calling workflow
written 7.8 years ago by zhou_1228 •

I wrote a program to automate SNP calling process, but not sure if it works since I am new to shell scricpts. If I can …
How can I restore paired-end fastq files from a sorted bam file with fastq remain sorted?
written 8.7 years ago by lghust2011

I have two paired-end FASQ files named fq1.fastq and fq2.fastq. Then I use BWA MEM to map these two files with the paired-end mode, like …
How to quickly count the number of reads at a specific coordinate in a bam file
written 9.2 years ago by klin •

Hi all, I am trying to count the number of reads (or alignments) for specific genomic locations in a bam file. I don't want reads …
Weird issue with bam file generated by samtools
written 9.6 years ago by lakhujanivijay

I generated a bam file following the below steps: **indexed reference genome** bwa index reference.fasta **aligned using bwa mem** bwa mem ref.fasta f1.fq.gz f2.fq.gz > …
bwa mem algorithm filter secondary alignment
written 9.9 years ago by xudongyang •

Hello,everyone I want to use the bwa mem algorithm to analysis my PE 100 seqdata (BWA Version: 0.7.12-r1039). I used "-M" option to mark the …
Split BAM file into two subsamples
written 10.5 years ago by colin.kern

I want to create two random subsamples from a BAM file. I don't want to just do "samtools view -s 0.5" twice because I want …
Align metagenomic shotgun sequencing reads to primers
written 10.6 years ago by Crystal •

Hi All, I'm trying to use BWA to align metagenomics reads (~100bp) to six primers (~15bp). All the steps work fine but when I convert …