Get unique alignment stats from SAM file with a combined genome
As the title states, I have sam/bam files from alignments of our reads against an index built from two fastA genomes concatenated. From these sam files, how would I go about finding out which reads aligned to chromosome1-X in species 1 and which reads aligned to chromsome1-X in species 2.
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Building Salmon and STAR index using rat GRCr8 reference genome
written by Ajina •
-
Truncating alignments in a SAM/BAM file
written by everj •
-
Sequence divergence of whole genomes
written by Colaptes
-
How to extract only the alignments to the chloroplast reference genome of a sam file
written by macielrodriguez2 •
-
Unmapped reads to fastq (hisat2 / samtools / bedtools) not matching pair-end reads
written by juan.crescente •
-
Making BAM or SAM files equivalent
written by keturner123 •
-
Filter sam/bam files
written by second_exon
-
merge BAM files from different lanes
written by e.amiri79 •
-
Align short reads to multifasta reference - cant see in IGV
written by kanacska •
-
Has anyone used Rhizoctonia solani AG1-IB genome for mapping using tophat2/bowtie2?
written by mirza
Use Qualimap, if you want visual results of alignments. You can also use
samtools idxstatsto get a quick recap of alignment stats.If you actually want to extract reads that are aligning to specific chromosomes then use
bazam(https://github.com/ssadedin/bazam ). You cal also usesamtools viewalong with an intervals file in BED format.Since they share the same chromosomes idxstats gives the total alignments to both, I'll check out bazam and Qualimap though and report back if they work for me, thank you.
In that case, there is no way to distinguish the alignments unless your reads have index sequences (distinct for individual species). Another thing to keep in mind. Not all aligners retain index sequence information (they may drop remainder of the string after a white space in fastq header) in the resulting SAM/BAM file. I know that
bbmap.shfrom BBMap suite retrains this information.