Skip to content

Instantly share code, notes, and snippets.

@amwenger
Created January 6, 2015 01:20
Show Gist options
  • Save amwenger/a7b94361c9ed50c62bc7 to your computer and use it in GitHub Desktop.
Save amwenger/a7b94361c9ed50c62bc7 to your computer and use it in GitHub Desktop.
Test case for samtools merge sort order validation
@HD VN:1.4 SO:coordinate
@SQ SN:1 LN:249250621
@RG ID:TEST PL:LS454 SM:TESTSM
A1 147 1 69584 0 101M = 69330 -355 ATGAGGTCGATAGTTTTTATTGTGACCTTCCTAGGGTAATCAAACTTGCCTGTACAGATACCTACAGGCTAGATATTATGGTCATTGCTAACAGTGGTGTG CA<DBDDBECAAB>DFDDFEHHHA@DGGGEGIIHHAEGIIIIHGGGGGECGGGFFGIF?CGGHDIJJIIHAIGIJIGJIJIIGEIIGIHHFGHFFFFFCCC NM:i:0 MD:Z:101 AS:i:101 XS:i:101 RG:Z:TEST
B1 99 1 69585 0 101M = 69759 275 TGAGGTCGATAGTTTTTATTGTGACCTTCCTAGGGTAATCAAACTTGCCTGTACAGATACCTACAGGCTAGATATTATGGTCATTGCTAACAGTGGTGTGC ?@@DD;DAACC@<CGHI=E4A:?492+AEF@FGDH9CA9*:DG3DA9*BDGI@DD8*....8BD=>@D@7@@=7ACC>ED7?)..;6;CC>@>>;A>@### NM:i:0 MD:Z:101 AS:i:101 XS:i:101 RG:Z:TEST
C1 83 1 69588 0 101M = 69571 -118 GGTCGATAGTTTTTATTGTGACCTTCCTAGGGTAATCAAACTTGCCTGTACAGATACCTACAGGCTAGATATTATGGTCATTGCTAACAGTGGTGTGCTCA ?;5=A>A:BCCCCCDB;DEC?;>E@;D9EC@HGC@GCIGGGGAF;CCFD9GHGHF?HFDGIHFC9GFC<F@CCHCGHDFC@IHE?<:FDFBHDDDDDD@?= NM:i:0 MD:Z:101 AS:i:101 XS:i:101 RG:Z:TEST
D1 163 1 69585 0 101M = 69629 145 TGAGGTCGATAGTTTTTATTGTGACCTTCCTAGGGTAATCAAACTTGCCTGTACAGATACCTACAGGCTAGATATTATGGTCATTGCTAACAGTGGTGTGC CCCFFDFFHHGGHHJJJIIJGHIJJJJJJJJJJJJ?DHIJJJIJIJJJJIIJJJJJJJGIJIGIIJJJJIIGHFGHHFFDEFFFFFCCEDEEDCCC>?AB: NM:i:0 MD:Z:101 AS:i:101 XS:i:101 RG:Z:TEST
@HD VN:1.4 SO:coordinate
@SQ SN:1 LN:249250621
@RG ID:TEST2 PL:LS454 SM:TESTSM2
A2 147 1 69584 0 101M = 69330 -355 ATGAGGTCGATAGTTTTTATTGTGACCTTCCTAGGGTAATCAAACTTGCCTGTACAGATACCTACAGGCTAGATATTATGGTCATTGCTAACAGTGGTGTG CA<DBDDBECAAB>DFDDFEHHHA@DGGGEGIIHHAEGIIIIHGGGGGECGGGFFGIF?CGGHDIJJIIHAIGIJIGJIJIIGEIIGIHHFGHFFFFFCCC NM:i:0 MD:Z:101 AS:i:101 XS:i:101 RG:Z:TEST2
B2 99 1 69585 0 101M = 69759 275 TGAGGTCGATAGTTTTTATTGTGACCTTCCTAGGGTAATCAAACTTGCCTGTACAGATACCTACAGGCTAGATATTATGGTCATTGCTAACAGTGGTGTGC ?@@DD;DAACC@<CGHI=E4A:?492+AEF@FGDH9CA9*:DG3DA9*BDGI@DD8*....8BD=>@D@7@@=7ACC>ED7?)..;6;CC>@>>;A>@### NM:i:0 MD:Z:101 AS:i:101 XS:i:101 RG:Z:TEST2
C2 83 1 69588 0 101M = 69571 -118 GGTCGATAGTTTTTATTGTGACCTTCCTAGGGTAATCAAACTTGCCTGTACAGATACCTACAGGCTAGATATTATGGTCATTGCTAACAGTGGTGTGCTCA ?;5=A>A:BCCCCCDB;DEC?;>E@;D9EC@HGC@GCIGGGGAF;CCFD9GHGHF?HFDGIHFC9GFC<F@CCHCGHDFC@IHE?<:FDFBHDDDDDD@?= NM:i:0 MD:Z:101 AS:i:101 XS:i:101 RG:Z:TEST2
D2 163 1 69585 0 101M = 69629 145 TGAGGTCGATAGTTTTTATTGTGACCTTCCTAGGGTAATCAAACTTGCCTGTACAGATACCTACAGGCTAGATATTATGGTCATTGCTAACAGTGGTGTGC CCCFFDFFHHGGHHJJJIIJGHIJJJJJJJJJJJJ?DHIJJJIJIJJJJIIJJJJJJJGIJIGIIJJJJIIGHFGHHFFDEFFFFFCCEDEEDCCC>?AB: NM:i:0 MD:Z:101 AS:i:101 XS:i:101 RG:Z:TEST2
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment