SEQanswers

Go Back   SEQanswers > Sequencing Technologies/Companies > Illumina/Solexa



Similar Threads
Thread Thread Starter Forum Replies Last Post
Samtools "is recognized as '*'" "truncated file" error axiom7 Bioinformatics 3 11-26-2014 03:53 AM
Samtools error: "different target sequence name" evansdan Bioinformatics 2 12-02-2013 03:00 AM
DEXSeq error in estimateDispersions: match.arg(start.method, c("log(y)", "mean")) fpadilla Bioinformatics 14 07-03-2013 03:11 PM
samtools installation error "cannot find -lcurses" ihelbig Bioinformatics 5 02-01-2013 12:48 PM
bowtie2 failing with "Could not parse passthrough output:" error ManuSetty Bioinformatics 0 04-13-2012 10:59 PM

Reply
 
Thread Tools
Old 04-11-2016, 01:46 AM   #1
znasim09
Member
 
Location: Seoul, Korea

Join Date: Sep 2015
Posts: 22
Default Samtools mpileup error "[W::sam_read1] parse error at line 1"

Hello everyone,

I am analyzing WGS data of my mutant (Raw data in fastq), I performed the alignment by using Bowtie and BWA as:

Using Bowtie:
bowtie -S reference.fa file.sam
then converted the sam file to bam file
samtools view -bS -o file.bam file.sam

Using BWA:
bwa aln reference.fa file.fastq > file.sai
bwa samse reference.fa file.sai file.fastq gzip > file.sam.gz

then converted sam.gz to bam
samtools view -bt reference.fa file.sam.gz | samtools sort -o file.bam

Then I tried to use mpileup on both bam files but got similar errors:

samtools mpileup -v reference.fa file.bam > outfile

[mpileup] 2 samples in 2 input files
<mpileup> Set max per-file depth to 4000
[W::sam_read1] parse error at line 1

Can anyone explain whats wrong and how it can be corrected?
Thanks
znasim09 is offline   Reply With Quote
Old 04-11-2016, 01:56 AM   #2
dschika
Member
 
Location: Zurich

Join Date: Mar 2010
Posts: 56
Default

Check samtools mpileup man page.
You may want to try -f instead of -v ?
dschika is offline   Reply With Quote
Old 04-11-2016, 04:24 AM   #3
GenoMax
Senior Member
 
Location: East Coast USA

Join Date: Feb 2008
Posts: 6,576
Default

Are you using the newest versions of bwa/samtools?
GenoMax is offline   Reply With Quote
Old 04-11-2016, 05:25 AM   #4
znasim09
Member
 
Location: Seoul, Korea

Join Date: Sep 2015
Posts: 22
Default

@dschika will it then give vcf format output?

@GenoMax,

I am using bwa-0.7.13 and
samtools-1.3
znasim09 is offline   Reply With Quote
Old 04-11-2016, 05:39 AM   #5
dschika
Member
 
Location: Zurich

Join Date: Mar 2010
Posts: 56
Default

Ok, probably you want to use -v and -f, my bad.

As you read the manual which you can find here, you know what -v and -f do:

Quote:
-f, --fasta-ref FILE The faidx-indexed reference file in the FASTA format. The file can be optionally compressed by bgzip.
-v, --VCF Compute genotype likelihoods and output them in the variant call format (VCF). Output is bgzip-compressed VCF unless -u option is set.
dschika 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 01:54 AM.


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