Number of reads mapped with each sequence in a multiple reference sequences
Hi all,
I have single-end reads and my aim is to map all the reads with multiple reference sequences in one FASTA file. These reference sequences are on the same gene location but they contain different variants. Does anyone know how I can enumerate the number of reads mapped with each sequence in my reference sequence list? Many thanks.
• 1,988 views
•
link
updated
by
natasha.sernova
•
written
by
KT •
•
written
by
KT •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
Which Reference Hg38 or T2T?? Both??
written by Jon •
-
RNA Seq analysis
written by Micro_bioinfo •
-
reference fasta creation
written by geneart$$ •
-
trouble with 'too short' reads using STAR. each mate maps fine, but the paired end mapping does not
written by james.alan.gregory •
-
How minimap2 mapping works
written by SUDOsundu •
-
Quick question: samtool idxstats and multiple mapping
written by Aurelie MLB
-
Counting single end reads
written by valopes •
-
Re-merge split BAM files
written by KT •
-
Align short reads to multiple reference sequence
written by Sandeep
-
How To Delete The All Fastq Reads Which Includes A Potential 50Bp Illumina Single End Pcr Primer 1
written by Tonyzeng •
Is it what you need?
Extracting reads mapped to any of multiple references
Depending on how many differences there are in the reads and the reference and the fact that most reads are going to multi-map (which you would need to adjust settings for with an NGS aligner) this may require some trial and error to come up with settings that will work.
See this thread for ideas: bwa mem perfectly align reads over their entire length
This sounds like an ideal use-case for Salmon or Kallisto.