How do I figure out the percent of reads that went into each of the bins that mapped to ? In the beginning, I concatenated my .fa files for my index file to mapp against all of my paired fastq files. Now i just want to know the number of reads from the fastq files make up each one of the bins.
How to view read mapping from Bowtie2, Samtools, etc
I am wanting to be able to view the percentage of reads from metagenomic bins to the paired fastq files. I have already completed making .bam files, and .sam files and sorted, and indexed the sam files with samtools.
I just need the percent of reads from my bins that are found in any of the seven paired end reads that I have. Assume that I have all of the correct files and am just looking to view the output.
Thanks for your help in advance.
1 answer
samtools idxstats <aln.bam>
Retrieve and print stats in the index file. The output is TAB-delimited with each line consisting of reference sequence name, sequence length, # mapped reads and # unmapped reads.
That previous command shows the reference sequence name (NODE_###_##...) and then all of the contig mapped numbers. How does adding up the reads mapped to each sequece show me the percent of reads from each bin when I don't know what contigs from each fastq went into which bin ?
Log in to answer this question.
More posts like this
-
BAM 2 FASTQ using Picard and Samtools
written by kcarey •
-
Maxbin2 usage for paired end reads
written by v.berriosfarias •
-
How do I find the percentage of paired-end fastq files containing a given string
written by lkianmehr
-
Can featureCounts use sortedbycoordinate files?
written by tryseq •
-
Determine mapped reads per scaffold in bam file: alternatives to samtools idxstat
written by cssulliv •
-
Stop BWA from storing unmapped reads
written by CC •
-
samtools covert to bam, sort and index all gz.sam
written by bioguy24
-
Best way of Merging Multiple Lanes generated by HiSeq 2000, index BAM or VCF?
written by venuraherath •
-
How can I restore paired-end fastq files from a sorted bam file with fastq remain sorted?
written by lghust2011
-
samtools mpileup Aborted (Core dumped)
written by rafa.rios.50 •
