SEQanswers

Go Back   SEQanswers > Applications Forums > Genomic Resequencing



Similar Threads
Thread Thread Starter Forum Replies Last Post
problems understanding pileup format pi101 Bioinformatics 2 11-14-2012 02:47 PM
samtools pileup BioSlayer Bioinformatics 9 01-16-2012 03:10 AM
samtools pileup output file charlie_sequencing Bioinformatics 1 02-10-2011 12:51 PM
SAMtool's pileup format - reference base question mfischer Bioinformatics 11 10-19-2010 04:12 PM
SAMtools - * character in pileup file JosieR Bioinformatics 9 09-08-2010 11:14 AM

Reply
 
Thread Tools
Old 01-22-2013, 10:10 AM   #21
jgibbons1
Senior Member
 
Location: Worcester, MA

Join Date: Oct 2009
Posts: 133
Default

Using pileup with the -f argument allows you to supply the faidx indexed reference sequence file. I used this option and it fixed my problem.
jgibbons1 is offline   Reply With Quote
Old 03-05-2013, 12:49 PM   #22
nupurgupta
Member
 
Location: New Jersey

Join Date: Aug 2010
Posts: 29
Default Same problem

Did you find a solution to the null problem please?

Quote:
Originally Posted by ericpante View Post
Hi everybody,

I am have similar problems with samtools 0.1.18. I would like to have reference characters listed in a pileup files, but I have problems with headers.

samtools faidx AGSbrut.fasta
samtools view -q 20 -buh -t AGSbrut.fasta.fai A.sam | samtools sort - A
samtools view -q 20 -buh -t AGSbrut.fasta.fai S.sam | samtools sort - S
samtools mpileup -B -f AGSbrut_index.fai A.bam S.bam > AS.mpileup

[fai_build_core] different line length in sequence 'null'.
Segmentation fault

I hypothesized that this 'null' sequence may be a blank line; so I looked for it manually and with sed, with no luck. I also looked for other potential problems based on what was previously reported (no extra spaces, characters, etc in reference sequence names in fai and sam files). I also tried to re-head the file, with no success:

samtools view -HS -t AGSbrut.fasta.fai A.sam > Aheader.sam
samtools reheader Aheader.sam A.bam > Aheaded.bam

[bam_header_read] EOF marker is absent. The input is probably truncated.

All insights are welcome!
thank you, eric
nupurgupta is offline   Reply With Quote
Reply

Tags
mosaik, pileup, reference, samtools

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 02:40 PM.


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