Info
La pregunta está cerrada. Vuélvala a abrir para editarla o responderla.
Reconstruct Alignment changes (read & reference) from CIGAR and MD:Z tag in SAM files
1 visualización (últimos 30 días)
Mostrar comentarios más antiguos
Hi all,
Firstly, I'm not a bioinformatician but have bioinformatics toolbox for work.
I'm writing my own function to interpret a combination of the CIGAR & MD:Z tag from sam files, so I can see the resulting alignment pattern.
Line1: Some READ/REFERENCE OF TGAC'S
Line2: Some alignment pattern of "|| | |||" for example
Line3: Some READ/REFERENCE OF TGAC'S
What I'm wondering is if there is a function that will take the MD and CIGAR and output the identical character positions (maybe even take the read from column 10 in the SAM and give a predicted reference (although insertions would be unknown bases).
Anyone know of such a function for Matlab, either in Bio toolbox or on file-exchange?
Thanks
0 comentarios
Respuestas (0)
La pregunta está cerrada.
Ver también
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!