I am working with a data set of ONT MinION metagenome data aligning it against a fasta file where each entry is a different species using a series of aligners (BWA, Bowtie2, SNAP, and Minimap2) to see which aligner yields the best results. However for each alinger I am getting more alignments to the specific species than there are reads in the input file.
For example a fastq containing 100,000 reads of Mouse and 50,000 reads of Wood Mouse Herpes virus aligned against my database will return 300,000 alignments for Mouse and 120,000 alignments for Wood Mouse Herpes Virus.
I understand that this is in some part due to the fact that some of these aligners report secondary and chimeric alignments but I was wondering if there wasnt a known aligner or aligner configuration wherein I could report only a single best alignment for each read in my input file?
1 answer
Log in to answer this question.
-
Question about RNA-Seq data alignment
written by mohammedtoufiq91
-
Recommended mapper
written by biobiu
-
Removing paralogs from multiple sequence alignments
written by mglasena •
-
Big differences between two alignment tools outputs: BWA vs NanoPipe
written by alexandru.bologa.marian •
-
Which NGS aligner to find inexact matches
written by Anand Rao
-
Is HISAT2 able to output only the best alignment?
written by RoseString •
-
Differences between two SAM/BAM files with samtools
written by dshulgin
-
Aligner for super-short sequences
written by IV
-
report all possible alignments
written by cmo
-
reporting the best alignment(s) only in pairwise2
written by sumithrasank75
Have you looked at the alignments to see those being produced by
bwaandbowtie2seem reasonable/logical? Look at the CIGAR strings and lengths of alignments.minimap2is the only bonafide aligner for long reads (not sure about SNAP) and should produce results that you should compare the others to.As far as I can tell the alignments appear valid just that there's more than there should be. This phenomenon is across the board whether I'm using
minmap2,bowtie2,snap, orBWA. I'm not too familiar with CIGAR strings but from what I can tell the alignments all appear to be good. If I wanted to double check the CIGAR strings how would I do that?how would you treat the reads/alignments that are truly double in your reference?
In the case where the read aligns twice to the same species I'd consider it a single 'hit' for the species in question. In the case where a read aligns to multiple species it would be something I'd need to investigate further to determine it's origin. I'm more preoccupied with identifying origin of read than I am the specifics of its alignment to the species in question as the metagenome Im working is going to be consisting of only a handful of potential organisms.