BWA-mem alignment
1
0
Entering edit mode
5 weeks ago
ICfc97 • 0

Hi!

I'm doing an alignment to a reference genome by using the bwa-mem algorithm but I'm obtaining a bam file in which all the data seems to appear in the Chr1 as I see in IGV and the bam file obtained. The code that I'm using is:

threads=SLURM_JOB_CPUS_PER_NODE bwa index ./GRCh37.p13_alignCNVs.fasta bwa mem -tthreads ./GRCh37.p13_alignCNVs.fasta ./NA12878_1.fastq ./NA12878_2.fastq |
samtools view -S -b - |
samtools sort - -o sample_sorted.bam
samtools index sample_sorted.bam

WGS Data Genomics • 176 views
2
Entering edit mode

What is the output of grep '>' GRCh37.p13_alignCNVs.fasta?

0
Entering edit mode

The output is:

>chr1

2
Entering edit mode
5 weeks ago
ATpoint 50k

So, your reference file only has one chromosome "chr1", therefore it is not a surprise that all reads are aligned to chr1. So what exactly did you expect that would happen here? Could it be you downloaded an incomplete or wrong reference? What is the analysis goal?

0
Entering edit mode

That's true! The goal is to detect CNVs in the whole genome. I guess I downloaded the wrong reference genome, that was a silly mistake. Thanks a lot!