BWA aln, according to GDC DNA Seq Analysis Pipeline
Hi everyone,
I am following the DNA Seq Analysis Pipeline by GDC (https://docs.gdc.cancer.gov/Data/Bioinformatics_Pipelines/DNA_Seq_Variant_Calling_Pipeline/).
They give the following hint for bwa aln: 'If the quality scores are encoded as Illumina 1.3 or 1.5, use BWA aln with the "-l" flag.'
How do I get this information from my bam files?
Thank you!
Best
Lynn
• 1,376 views
•
link
updated
by
Zhenyu Zhang
•
written
by
loy_loy •
•
written
by
loy_loy •
0 answers
No answers yet.
Log in to answer this question.
More posts like this
-
How to choose BWA-MEM2 settings
written by mgranada3 •
-
WES data analysis: the rationale behind merging BAM files and marking duplicates
written by Thanh •
-
When to merge BAM or SAM files?
written by Thanh •
-
Ploidy in WGS samples
written by loy_loy •
-
Why does TCGA RNA-seq pipeline starts from a BAM file?
written by シン •
-
How to align incomplete pairs and singled end reads?
written by loy_loy •
-
SortSam Error after BWA alignment
written by loy_loy •
-
Bwa Mem: How To Specify The Fastq Phred Format?
written by Giovanni M Dall'Olio
-
What Is The Default Quality Encoding Expected By Bwa?
written by Panos
-
Fastq Format Problem When Using Bwa Aln -Q Option (Read Trim)
written by Bioscientist
If you have the fastq files you can determine it (more info in here ) because version 1.5 doesn't use "@" and "A" for the phred scores
Right. Basically either you know from others or you need to guess from the phred score enums. Btw, you are not expected to encounter these 1.3 or 1.5 encoding unless you are working on some super old data (like a small fraction of TCGA).