different bam from same NGS data
1
0
Entering edit mode
5.9 years ago
gab ▴ 20

I was given some bam files I used to call somatic variants. The high number of SNP and mutations i found seemed suspicious, so I asked for the original fastq files to redo the allignment. While waiting for them I tried to get fastq from those bam files using bedtools bamtofastq . I used bwa mem to check allignment and the new bam file yeld less mutations and SNPs. Is there a procedure I didn't understand in the bedtools function that might have given me a fastq file different than the original?

bam bwa • 1.1k views
ADD COMMENT
1
Entering edit mode

are you using the same parameters/software/ remove duplicates etc... ?

ADD REPLY
0
Entering edit mode

I guess so, although I'm still waiting for the guy who did the first analysis to reply

ADD REPLY
1
Entering edit mode

And that you aligned to the same reference?

ADD REPLY
0
Entering edit mode

Yes, same exact file fasta

ADD REPLY
1
Entering edit mode
5.9 years ago
jsporter ▴ 60

The BAM file will not necessarily include the whole read from the original FASTQ file since it can be clipped, and this could lead to different alignment results.

ADD COMMENT

Login before adding your answer.

Traffic: 2819 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