Skip to content

Instantly share code, notes, and snippets.

Show Gist options
  • Save nh13/3240167 to your computer and use it in GitHub Desktop.
Save nh13/3240167 to your computer and use it in GitHub Desktop.
TMAP aln ROC Single Ends
# tmap mapall -f ${REF} -r r0.fq -v stage1 map1 map2 map3 > r0.fq.map1map2map3.sam
40 77347 0
39 77817 0
38 78306 0
37 78962 0
36 79452 0
35 79794 0
34 80128 0
33 80742 0
32 81118 0
31 81435 0
30 81803 0
29 82515 0
28 82939 0
27 83297 0
26 83648 0
25 84347 0
24 84766 0
23 85110 0
22 85418 0
21 86202 0
20 86506 0
19 86806 0
18 87087 0
17 87980 0
16 88249 0
15 88583 0
14 88839 0
13 89765 0
12 90014 0
11 90400 1
10 90662 1
9 91712 1
8 92777 27
7 93178 28
6 93418 30
5 94652 44
4 94943 59
3 95493 79
2 95957 119
1 96397 242
# tmap mapall -f ${REF} -r r0.fq -v stage1 map4 > r0.fq.map4.sam
40 84212 0
39 84470 0
38 84724 0
37 85033 0
36 85271 0
35 85438 0
34 85632 0
33 85988 0
32 86200 0
31 86356 0
30 86538 0
29 86929 0
28 87130 0
27 87326 0
26 87512 0
25 87951 0
24 88187 0
23 88367 0
22 88545 0
21 89067 0
20 89283 0
19 89472 0
18 89658 0
17 90241 0
16 90433 0
15 90620 0
14 90764 0
13 91431 2
12 91588 3
11 91798 5
10 91930 5
9 92727 11
8 93559 33
7 93831 41
6 93986 52
5 95126 105
4 95275 126
3 95741 185
2 96120 265
1 96380 366
@nh13
Copy link
Author

nh13 commented Aug 2, 2012

Working on a single CPU:
"map1 map2 map3" command took 207.73 seconds (reference load time was 4.5 seconds).
"map4" command took 41.15 seconds (reference load time was 4.5 seconds).

Note that I specify the reference load time into memory here because the dataset is so small that it may be a dominant factor. Additionally, if the reference is pre-loaded into memory (with "tmap server") or the OS caches it in memory, the load time is reduced.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment