Skip to content

Instantly share code, notes, and snippets.

@peterjc

peterjc/SAMv1_padded.fasta

Created May 14, 2012 10:00
Show Gist options
  • Save peterjc/2693037 to your computer and use it in GitHub Desktop.
Save peterjc/2693037 to your computer and use it in GitHub Desktop.
Download ZIP
SAM/BAM examples from the specification
Raw
SAMv1_padded.fasta
>ref
AGCATGTTAGATAA**GATAGCTGTGCTAGTAGGCAGTCAGCGCCAT
Raw
SAMv1_padded.sam
@HD VN:1.3 SO:coordinate
@SQ SN:ref LN:47
ref 516 ref 1 0 14M2D31M * 0 0 AGCATGTTAGATAAGATAGCTGTGCTAGTAGGCAGTCAGCGCCAT *
r001 163 ref 7 30 14M1D3M = 39 41 TTAGATAAAGGATACTG *
* 768 ref 8 30 1M * 0 0 * * CT:Z:.;Warning;Note=Ref wrong?
r002 0 ref 9 30 3S6M1D5M * 0 0 AAAAGATAAGGATA * PT:Z:1;4;+;homopolymer
r003 0 ref 9 30 5H6M * 0 0 AGCTAA * NM:i:1
r004 0 ref 18 30 6M14N5M * 0 0 ATAGCTTCAGC *
r003 16 ref 31 30 6H5M * 0 0 TAGGC * NM:i:0
r001 83 ref 39 30 9M = 7 -41 CAGCGCCAT
Raw
SAMv1_unpadded.fasta
>ref
AGCATGTTAGATAAGATAGCTGTGCTAGTAGGCAGTCAGCGCCAT
Raw
SAMv1_unpadded.sam
@HD VN:1.3 SO:coordinate
@SQ SN:ref LN:45
r001 163 ref 7 30 8M2I4M1D3M = 37 39 TTAGATAAAGGATACTG *
r002 0 ref 9 30 3S6M1P1I4M * 0 0 AAAAGATAAGGATA *
r003 0 ref 9 30 5H6M * 0 0 AGCTAA * NM:i:1
r004 0 ref 16 30 6M14N5M * 0 0 ATAGCTTCAGC *
r003 16 ref 29 30 6H5M * 0 0 TAGGC * NM:i:0
r001 83 ref 37 30 9M = 7 -39 CAGCGCCAT *
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment