Question: Problem with alignment
0
gravatar for nayshool
4.4 years ago by
nayshool20
Israel
nayshool20 wrote:

Hello everyone!

We have received from BGI china a few samples we sent for exome sequence + Basic analysis. We got the fastQ files and BAM files. I have made my own BAM files using BWA and Picard tools using the following commends:

BWA:
bwa mem human_g1k_v37.fasta file1.fq.gz file2.fq.gz file3.fq.gz file4.fq.gz > result.sam

PICARD:
java -jar SortSam.jar I=result.sam O=result.bam SORT_ORDER=coordinate

The size of the files that I got (2Gb) , was about a half from the BAM files BGI (4.2 Gb) sent me.
Here are the samtools flagstat results for both files:

my BAM file:

[gen-biorep@gen-biorep-3 VA6]$ samtools flagstat result.bam
28842836 + 0 in total (QC-passed reads + QC-failed reads)
0 + 0 duplicates
28824861 + 0 mapped (99.94%:-nan%)
28842836 + 0 paired in sequencing
14423043 + 0 read1
14419793 + 0 read2
28593748 + 0 properly paired (99.14%:-nan%)
28814004 + 0 with itself and mate mapped
10857 + 0 singletons (0.04%:-nan%)
165631 + 0 with mate mapped to a different chr
131466 + 0 with mate mapped to a different chr (mapQ>=5)




BGI BAM File:

[gen-biorep@gen-biorep-3 result_alignment]$ samtools flagstat VA6.rmdup.bam
57668392 + 0 in total (QC-passed reads + QC-failed reads)
3039053 + 0 duplicates
57369744 + 0 mapped (99.48%:-nan%)
57668392 + 0 paired in sequencing
28834196 + 0 read1
28834196 + 0 read2
56991680 + 0 properly paired (98.83%:-nan%)
57273460 + 0 with itself and mate mapped
96284 + 0 singletons (0.17%:-nan%)
198072 + 0 with mate mapped to a different chr
181774 + 0 with mate mapped to a different chr (mapQ>=5


why does half of my reads disappear? How can I solve this problem?

Thank you,

Omri

alignment next-gen • 1.5k views
ADD COMMENTlink modified 4.1 years ago by naruto1nrt20 • written 4.4 years ago by nayshool20

Hello nayshool!

It appears that your post has been cross-posted to another site: SeqAnswers (where I happen to have already replied).

This is typically not recommended as it runs the risk of annoying people in both communities.

ADD REPLYlink modified 4.4 years ago • written 4.4 years ago by Devon Ryan86k

i.e. Duplicate of http://seqanswers.com/forums/showthread.php?t=45597

ADD REPLYlink written 4.4 years ago by Peter5.7k
0
gravatar for naruto1nrt
4.1 years ago by
naruto1nrt20
Australia
naruto1nrt20 wrote:

In your command, file3 and file4 will be ignore. Your bam file didn't have these reads!

ADD COMMENTlink written 4.1 years ago by naruto1nrt20
Please log in to add an answer.

Help
Access

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1275 users visited in the last hour