Hello all,
I performed a read mapping of one species against a reference genome. I converted the file_sorted.bam to a file.mpileup. Now, I'd like to create a file containing shared (identical) positions between the reference genome and the mapped species. Also, these shared positions must have not only a good mapping quality but also a read depth > 10. Does anyone know the best way to do this?
Thank you!
•
written
by
samocarp •
1 answer
I would suggest (1) using VarScan to identify differences and similarities between reference and your sample or (2) writing a Python script if you know Python. If you use VarScan you could do something like this:
java -Xmx8g -jar VarScan.v2.3.9.jar mpileup2snp file.mpileup --min-coverage 10
If you write a Python script you could just take a look at the pileup file format spec here: http://samtools.sourceforge.net/pileup.shtml. It's pretty easy to identify SNPs (differences between sample and reference) and read depth in a pileup file. In your case since you want sites that are the same between reference and sample you will want straight periods or commas (a match to the reference genome on the forward or reverse strands, respectively) in your fifth column.
Log in to answer this question.
-
reference-guided assembly tools for short read data?
written by Hani •
-
Mean samples coverage and mapping quality in Qualimap
written by Eli •
-
% of reads unmapped: too short - STAR 2.5.2b
written by quan.clement •
-
ChIP Seq poor mapping quality and strange fast QC result
written by e0223167 •
-
How can I Annotate my miRNA data against MirBase?
written by k.kathirvel93
-
How reliable are quality scores?
written by harshitsahay.iitb •
-
Get identical mapped nucleotide length from bowtie2 generated SAM files
written by hdy
-
low number mapped reads in ribosome profiling
written by laura.fancello •
-
Calculating Mapping Coverage
written by Hatem Elshazly •
-
The Need Of Repeatmask For Read-Depth Based Cnv Detection Algorithm
written by Bioscientist
