MuTect short read preprocessing
Hi, I am programming to do short read preprocessing for BAM files. I have no idea about how to know if paired-end reads are overlapping and the overlapping region, and the mismatch bases for a read (I know the field MD gives the mismatching position, but I did not understand it. Why it is a string not a numeric array?) Thanks very much for your time.
Best Jing
• 1,218 views
•
link
written
by
jing.mengrabbit •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Determine fragment alignment overlaps by parsing CIGAR string in the MC tag
written by abcoxyzide •
-
Clipping bases from non-overlapping regions of mate reads in BAM file
written by lechu •
-
MD tags with diff ref sequence
written by radwa.raed •
-
Adding information tags with types, positions and counts of matches and mismatches in every read in…
written by lechu •
-
samtools pileup format: insertions and deletions' quality
written by jing.mengrabbit •
-
Question on SAMTOOLS MPILEUP output
written by bioplanet •
-
Question about MD tag in BAM file
written by abascalfederico
-
In CWL, how to make arrays optional input
written by jeltje.van.baren
-
Pysam fetching reads with indels: Returns wrong positions of bases
written by Lylthera •
-
Is it possible to reconstruct alignment from CIGAR and MD strings alone?
written by Lynxoid
