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

-F option yielding incongruous output #119

Open
cxr5298 opened this issue Feb 25, 2019 · 0 comments
Open

-F option yielding incongruous output #119

cxr5298 opened this issue Feb 25, 2019 · 0 comments

Comments

@cxr5298
Copy link

cxr5298 commented Feb 25, 2019

While using SNAP with an Illumina ampliseq dataset I get kind of surprising output when using SNAP's -F filter function.

As it is I run the data I'm splitting it out into two buckets one for aligned reads and unaligned reads, now conventional wisdom would suggest that the read totals across both the aligned and unaligned sets would add up to the read count of the run as a whole.

However it doesn't, it adds up as being greater than the original read count. Moreover there is overlap between the unaligned and aligned reads when I dive into the individual .sam files. I also looked at the supposed unaligned reads and they in fact align to my original reference, its just SNAP that's giving them mapq scores of 0. For the life of me I cannot tell why this is happening. Unfortunately I am not able to share any data so I understand if any input you can offer is limited. I will share what I can below:

The call to SNAP:
~/project/snap-aligner paired ~/data/indx/ ~/data/sample1_R1.fastq ~/data/sample1_R2.fastq -F a -o sample1_aln.sam ~/project/snap-aligner paired ~/data/indx/ ~/data/sample1_R1.fastq ~/data/sample1_R2.fastq -F u -o sample1_unaln.sam

The indexed reference I'm using is a small database of only 48 sequences ranging from 95 to 225 bases in length.

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