Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Question: graphmap align pacbio data #87

Open
zwm0502 opened this issue Mar 9, 2018 · 0 comments
Open

Question: graphmap align pacbio data #87

zwm0502 opened this issue Mar 9, 2018 · 0 comments

Comments

@zwm0502
Copy link

zwm0502 commented Mar 9, 2018

Hello,
I want to ask some questions about using graphmap to align Pacbio fastq data.
I use this command ,
graphmap align --out-fmt m5 --threads 8 -C --ordered --evalue 1e-5 --error-rate 0.15 -r ecoliMG1655.fasta -d query.fq -o out.m5
I found the detailed imformation of m5 format on this website https://github.com/PacificBiosciences/blasr/wiki/Blasr-Output-Format.

qName qLength qStart qEnd qStrand tName tLength tStart tEnd tStrand score numMatch numMismatch numIns numDel mapQV qAlignedSeq matchPattern tAlignedSeq

I write script to get the number of insertion,deletion and substitution in the alignment result from qAlignedSeq and tAlignedSeq.
But I found that the number of substitution is different from the numMismatch. In addition, match pattern is not consistent with the detailed information in qAlignedSeq and tAlignedSeq sometimes.

In the example of attachment I found 2888 substitutions from dealing with qAlignedSeq and tAlignedSeq. But there are 45 numMismatch written in the line.
Could you please tell me why? How do you define the substitution?
Is the condition that a base is substituted by another one will be counted?
And what's the meaning of '' in the matchPattern part? The sum of numMismatch,numIns and numDel is not equal to the number of '' in some lines but in somelines it equal to.
So is it suitable for graphmap to align pacbio data? Or is it customized for nanopore data?
Thanks very much for all your kind help in advance.
1.m5.txt

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant