SEQanswers

Go Back   SEQanswers > Bioinformatics > Bioinformatics



Similar Threads
Thread Thread Starter Forum Replies Last Post
Problem calling consensus sequence with Samtools mpileup? Lspoor Bioinformatics 3 10-24-2011 08:43 AM
Editing fasta , reference base in snp calling samtools moriah Bioinformatics 2 08-10-2011 12:11 AM
samtools/mpileup heterozygous SNPs calling combiochem Bioinformatics 4 08-02-2011 08:05 AM
calling Heterozygous SNPs with samtools mpileup egatti Bioinformatics 1 07-21-2011 09:16 AM
Samtools mpileup. Which reads are excluded during SNP calling? wwmm933 Bioinformatics 5 05-27-2011 09:12 AM

Reply
 
Thread Tools
Old 07-15-2012, 09:25 PM   #1
ferris
Junior Member
 
Location: China

Join Date: Jun 2012
Posts: 2
Default Calling for help in using samtools mpileup: all Ns in reference bases

I'm mapping 36bp single ended reads to C.elegans cDNA reference by bwa, trying to make variant calling after getting mpileup files. The problem is the samtools mpileup doesn't work well with my reference sequence and show all Ns in the third column.
Like this:
AC8.4 3 N 1 ^!G I
AC8.4 4 N 1 A I
AC8.4 5 N 1 A I
AC8.4 6 N 1 T I
AC8.4 7 N 1 G I
AC8.4 8 N 1 C I
Do you know what's wrong?
P.S. my reference consists of 55783 chromosomes (transcripts), each one with about 20~10k bp length.
My commanding lines are like this:
bwa -n 1 -t 7 -l 10000 -o 0 XXX.fa XXX.fastq > XXX.sai
bwa samse -n 10 XXX.fa XXX.sai XXXfastq > XXX.sam
samtools view -S -b XXX.sam > XXX.bam
samtools sort XXX.bam XXX_sorted
samtools mpileup -f XXX.fa XXX_sorted.bam > result.mpileup

Plz help me, thank you very much!
ferris is offline   Reply With Quote
Old 07-16-2012, 06:10 AM   #2
ferris
Junior Member
 
Location: China

Join Date: Jun 2012
Posts: 2
Default

Solved by myself already!
Somebody told me that samtool mpileup doesn't work with references with "|" or "." symbol in their chromosome names. So I tried first by simplifying chromosome (transcript) names to refrain "|" and ".". However it didn't work.
Later I found empty lines and Ns in some of sequence of my transcripts. After I removed all the bad transcripts, samtools worked quite well!
Hope anybody has similar problems will get help!
ferris is offline   Reply With Quote
Reply

Tags
bwa, reference base, samtools mpileup

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 06:52 AM.


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