SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
BWA sam and Samtools sam->bam conversion problem maasha Bioinformatics 6 06-05-2013 07:39 AM
BWA sampe problem patel Bioinformatics 9 10-24-2011 05:11 AM
bwa sampe problem sszelinger Bioinformatics 3 06-15-2011 05:34 AM
samtools-0.1.12a and bwa-0.5.9rc1 installing problem on fedora14 ketan_bnf Bioinformatics 2 01-14-2011 06:39 PM
BWA - Samtools problem joa_ds Bioinformatics 18 11-08-2009 06:55 PM

Reply
 
Thread Tools
Old 06-30-2010, 11:24 PM   #1
genelab
Member
 
Location: honston

Join Date: Nov 2009
Posts: 27
Default BWA - samtools problem

Hello everyone

There are some topics concerning BWA and Samtools, but i still can't find the reasons of the trouble.

Here is what I count

bwa index -a bwtsw human.fasta
bwa aln human.fasta read.fq > aln_sa.sai
bwa samse human.fasta aln_sa.sai read.fq > aln.sam
samtools view -b -h -S -T human.fasta aln.sam > aln.bam
samtools sort aln.bam aln.sort
samtools index aln.sort.bam

when i finished these steps and get the sorted bam file, i then run

java -jar ViewSam.jar I=aln.sort.bam

it got the following error message:


Exception in thread "main" java.lang.RuntimeException: SAM validation error: ERROR: Record 1, Read name 1_R9DycrMpuN1, MRNM should not be set for unpaired read.
at net.sf.samtools.SAMUtils.processValidationErrors(SAMUtils.java:334)
at net.sf.samtools.BAMFileReader$BAMFileIterator.advance(BAMFileReader.java:444)
at net.sf.samtools.BAMFileReader$BAMFileIterator.<init>(BAMFileReader.java:408)
at net.sf.samtools.BAMFileReader$BAMFileIterator.<init>(BAMFileReader.java:398)
at net.sf.samtools.BAMFileReader.getIterator(BAMFileReader.java:201)
at net.sf.samtools.SAMFileReader.iterator(SAMFileReader.java:283)
at net.sf.samtools.SAMFileReader.iterator(SAMFileReader.java:38)
at net.sf.picard.sam.ViewSam.doWork(ViewSam.java:68)
at net.sf.picard.cmdline.CommandLineProgram.instanceMain(CommandLineProgram.java:150)
at net.sf.picard.sam.ViewSam.main(ViewSam.java:58)

how can i fix this?
Thanks!!!
genelab is offline   Reply With Quote
Old 07-02-2010, 05:53 AM   #2
lh3
Senior Member
 
Location: Boston

Join Date: Feb 2008
Posts: 693
Default

Give the first line of your alignment, by samtools.
lh3 is offline   Reply With Quote
Old 07-02-2010, 11:13 PM   #3
genelab
Member
 
Location: honston

Join Date: Nov 2009
Posts: 27
Default

I have found one way to figure out the problem

I used the function "setValidationStringency" in samtools


thx!
genelab is offline   Reply With Quote
Reply

Thread Tools

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off




All times are GMT -8. The time now is 12:01 PM.


Powered by vBulletin® Version 3.8.9
Copyright ©2000 - 2020, vBulletin Solutions, Inc.
Single Sign On provided by vBSSO