Info
Cette question est clôturée. Rouvrir pour modifier ou répondre.
Reconstruct Alignment changes (read & reference) from CIGAR and MD:Z tag in SAM files
1 vue (au cours des 30 derniers jours)
Afficher commentaires plus anciens
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 commentaires
Réponses (0)
Cette question est clôturée.
Voir également
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!