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

-m setting not outputting a sam file with an -NH tag #133

Open
elayton13 opened this issue May 20, 2022 · 0 comments
Open

-m setting not outputting a sam file with an -NH tag #133

elayton13 opened this issue May 20, 2022 · 0 comments

Comments

@elayton13
Copy link

Hello,

I am using bowtie version 1.1.1 for small RNA sequencing analysis. When using the -m setting to allow multi mapping, I am not getting an -NH tag in the resulting sam file and I think this is then causing me issues when I am trying to run featurecounts with the --fraction option (which relies on the -NH tag to fraction the counts). The reason I think this is because my featurecounts raw counts are not decimal numbers, and because taking out the --fraction parameter doesn't seem to change the number of counts.

My original bowtie command was:
bowtie -v 1 -m 10 -S -a --best --strata

and subsequent featurecounts command:
featureCounts -T 16 -M -t miRNA --fraction -g ID -s 1 -F GTF -a

I have tried editing the bowtie command to take out the -a, the --best --strata and change the -m but I am still not getting an -NH tag outputted in my sam files. I do get the following types of tags:
XA:i:1 MD:Z:5T12 NM:i:1

Any insight would be really appreciated!
Thank you
Emma

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