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

Whether multiple mapping affect ribocode result #42

Open
Hilarialy opened this issue May 12, 2022 · 1 comment
Open

Whether multiple mapping affect ribocode result #42

Hilarialy opened this issue May 12, 2022 · 1 comment

Comments

@Hilarialy
Copy link

Hilarialy commented May 12, 2022

Sorry to bother you.
I followed your workflow use STAR arguments --quantMode TranscriptomeSAM and --outFilterMultimapNmax 1,this is the command looks like

STAR --outFilterType BySJout --runThreadN 10 --outFilterMismatchNmax 2 \
--genomeDir /reference/GRCm38/STAR \
--readFilesIn /my/rmrRNA/${sample}_trim_norrna.fq  \
--outFileNamePrefix ${sample} \
--outSAMtype BAM SortedByCoordinate \
--quantMode TranscriptomeSAM GeneCounts \
--outFilterMultimapNmax 1 --outFilterMatchNmin 16 --alignEndsType EndToEnd

the output Aligned.toTranscriptome.out.bam file still have multiple mapping sequence ( NH:i > 1 ) , because STAR wiil output all records in this toTranscriptome bam file , Does these multiple mapping records affect the Ribocode result to find uORF ?

@Hilarialy Hilarialy changed the title Does multiple mapping affect ribocode result ? Whether multiple mapping affect ribocode result May 12, 2022
@zhengtaoxiao
Copy link
Member

STAR excludes those reads aligned to multiple genome regions, not the transcriptome. It is expected to see "NH:i>1" in transcritome bam file for those genes having multiple overlapped isoforms. It has no influence on RiboCode result.

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

2 participants