Hi All,
I am thinking of using Mpileup from Samtools to get the information on the read quality and read depth from the BAM files. In that case which script should I use?
I have been recommended the script below:
ref=my_fasta.fasta samtools mpileup -f $ref -r chromosome_1:1018000-1018100 my_bam.bam my_bam2.bam
P.S. I dont have to do Variant calling, just need the read information.
Thanks in advance for your help.
1 answer
I would not use mpileup at all - if you need only alignment statistics. Use BedTools. For example here : A: How to get coverage for an alignment file The read quality you can get from FASTQC http://www.bioinformatics.babraham.ac.uk/projects/fastqc/.
But to be clear:
First what you should do is to check your reads by FASTQC (if the quality is bad do trimming), then align them to reference and then perform BedTools.
Best,
Agata
Log in to answer this question.
-
Missing short indels from vcf
written by kbranger •
-
Why does bcftools call from different max-depth (-d) value in mpileup give inconstant GT values?
written by magnolia •
-
No variants found ( have problem in the handbook)
written by quanyu •
-
Generating vcf files per individual instead of one big file from bam files
written by salman_96 •
-
merge individual runs after bcftools mpileup | bcftools call
written by Leonor •
-
Variant position showed up in mpileup but not in VCF call
written by leemhmark •
-
Why Does Mpileup Skip some variants ?
written by sami •
-
Variant Calling (Samtools Mpileup): Error
written by venuraherath •
-
Error in running a Samtools mpileup script for getting base frequency information
written by jigarnt •
-
Samtools mpileup quality scores all zero
written by fwuffy
