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 -t $threads ./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
ADD COMMENT
2
Entering edit mode

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

ADD REPLY
0
Entering edit mode

The output is:

>chr1
ADD REPLY
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?

ADD COMMENT
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!

ADD REPLY

Login before adding your answer.

Traffic: 1031 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6