[nilshomer]

Thanks for nilshomer's help!
Now the error message is disappeared.
But I only can see a black screen when I typed the command:
samtools tview aln.sorted.bam ref.fasta
I want to see how many short reads are aligned on certain one position, I don't know if the command can show the information.

When you enter tview, you are put at the first contig, first position. Try navigating to the region of interest (type "g" etc.). I think the help screen is the "?" key.

[iloveneworleans]

Thanks for nilshomer's help!
Now the error message is disappeared.
But I only can see a black screen when I typed the command:
samtools tview aln.sorted.bam ref.fasta
I want to see how many short reads are aligned on certain one position, I don't know if the command can show the information.

[nilshomer]

I have a problem about converting sam format file into bam format file with samtools.
I downloaded the latest version:v0.1.5c
and my fasta file is about 3G; I met a "Bus error" when I run the below command,
samtools view -bt h_sapiens.fa -o output.bam h_sapiens.sam

The error message is:
[sam_header_read2] 23631565 sequence loaded.
Bus error

Does anyone meet the kind of problem before?
Could you tell me how to solve this problem?
thanks a lot

I think the problem is the "-t" option:

         -t FILE  list of reference names and lengths (force -S) [null]

It requires the names and lengths of the contigs in the reference. To create this file, use "samtools faidx" on your reference. A ".fai" should be created, which is what you should use instead for the "-t" above.

[iloveneworleans]

I have a problem about converting sam format file into bam format file with samtools.
I downloaded the latest version:v0.1.5c
and my fasta file is about 3G; I met a "Bus error" when I run the below command,
samtools view -bt h_sapiens.fa -o output.bam h_sapiens.sam

The error message is:
[sam_header_read2] 23631565 sequence loaded.
Bus error

Does anyone meet the kind of problem before?
Could you tell me how to solve this problem?
thanks a lot

[jess]

Does anyone know of any way SAM generates a .snp file?? This is needed as the input when i use the SNPfilter option in samtools.pl...so i was wondering if i am not aware of any option which generated .snp file. Let me knw ASAP

[lparsons]

I have noticed the same behavior with samtools tview. The biggest issue I have is with very deep alignments since I can find no way to scroll up or down in large chunks.

I also have a question about how people generally use tview. I have an alignment from a virus that is very deep (2000x) since it is such a small genome. I suspect this may be why there are MANY indels reported by BWA. However, most of those indels do not pass varFilter. Is there a way to regenerate the SAM/BAM file with the "false" variations removed? In other words, to view the SAM/BAM as filtered by samtools varFilter? Thanks.

[pparg]

I tried samtools tview to browse the alignment results in bam format. However, several keys do not work, including:
Arrows Small scroll movement
J, K Large scroll movement
backspace Scroll back one screen
..
I haven’t checked other keys. But would like to report what I saw. BTW, is there any other ways besides tview to view the alignment results in sam/bam format? Thanks a lot!

[zee]

I was looking for the same thing, however the -s also gave me the individual bases. I settled with a perl one-liner:

samtools pileup  -cf genome.fasta  file.bam |perl -lane 'print join"\t",@F[0..7]'

For more complicated things I am using the Bio:: DB::Sam pileup methods.
Hope this helps.

Bit of a newbie question. I've been trying to use the pileup analysis on a BWA dataset. Is there any way to switch of the read bases, read quality and alignment quality information in the output file and get a summarized format instead?

I'm looking at a small number of sequences that have a coverage of 50.000X upwards, and as a result the pileup output becomes almost unmanageable.

Thanks!

[ohofmann]

Bit of a newbie question. I've been trying to use the pileup analysis on a BWA dataset. Is there any way to switch of the read bases, read quality and alignment quality information in the output file and get a summarized format instead?

I'm looking at a small number of sequences that have a coverage of 50.000X upwards, and as a result the pileup output becomes almost unmanageable.

Thanks!

[gcrdb]

Bio:B::Sam perl APIs need to start from BAM files (-bam) , not SAM files(no "-sam" at all). I only have SAM files which from bwa, all I need is to convert SAM to BAM.
I am stuck with SAM files.....
samtools import ref .fai in.sam out.bam
got error:
[sam_header_read2] 22 sequences loaded.
[sam_read1] reference '-143963499' is recognized as '*'.
Parse error at line 1: invalid CIGAR operation
Aborted

thanks,

[lh3]

Lincoln has released SAM/BAM perl APIs a few days (weeks?) ago. It is here:

http://search.cpan.org/~lds/Bio-SamTools-1.00/lib/Bio/DB/Sam.pm

Compiling this module requires samtools C source codes. Bio:B::Sam is known to work with samtools-0.1.4 and 0.1.5 (released today).

BTW, the latest samtools supports opening BAM files over FTP. For example:

samtools tview ftp://ftp.ncbi.nih.gov/1000genomes/f...32.2009_06.bam

[gcrdb]

hi, I have trouble conveting sam to bam.. I tried both:

samtools import ref .fai in.sam out.bam
got error:
[sam_header_read2] 22 sequences loaded.
[sam_read1] reference '-143963499' is recognized as '*'.
Parse error at line 1: invalid CIGAR operation
Aborted

samtools view -bt ref .fai -o in.sam out.bam
and got similar error:
[sam_header_read2] 22 sequences loaded.
[sam_read1] reference '' is recognized as '*'.
[main_samview] truncated file.

thanks,

[ElMichael]

Hi,
could anybody, please, explain the output format of the wgsim_eval.pl script?
I used this script to evaluate aln.sam file after making alignment with BWA.

06x 1654169 / 3308330 3308330 5.000e-01
05x 31765 / 63530 3371860 5.000e-01
04x 4938 / 9872 3381732 5.000e-01
03x 163891 / 327252 3708984 5.001e-01
02x 65120 / 129918 3838902 5.001e-01
01x 2669 / 5090 3843992 5.001e-01
00x 113748 / 141416 3985408 5.109e-01

BTW, in the BWA-man is written that " These reads are mapped with bowtie, bwa, maq and soap... The resultant alignments were then evaluated with wgsim_eval.pl script. "
How could I use this script for alignments from other programs such as bowtie, soap?
thanks,
Mike.

[krawitzp]

genome likelihood format

Hi,
where can I find further documentation on the genome likelihood format 3.0 ?
thanks,
peter

[corthay]

I appreciate that you immediately replied to my question.
I would like to handle the sam format files.