Skip to content

Commit

Permalink
BMC RNA-seq notebook
Browse files Browse the repository at this point in the history
  • Loading branch information
willbradshaw committed Dec 19, 2023
1 parent 92fbff1 commit 35b2f01
Show file tree
Hide file tree
Showing 32 changed files with 110,007 additions and 9 deletions.
1,851 changes: 1,851 additions & 0 deletions data/2023-12-19_rna-seq-workflow/adapter_stats.tsv

Large diffs are not rendered by default.

49 changes: 49 additions & 0 deletions data/2023-12-19_rna-seq-workflow/basic_stats.tsv
Original file line number Diff line number Diff line change
@@ -0,0 +1,49 @@
sample percent_gc mean_seq_len n_read_pairs percent_duplicates n_bases_approx per_base_sequence_quality per_tile_sequence_quality per_sequence_quality_scores per_base_sequence_content per_sequence_gc_content per_base_n_content sequence_length_distribution sequence_duplication_levels overrepresented_sequences adapter_content stage
1A 51 133.8680771459585 65476301 96.2554841175091 17400000000 pass/pass fail/fail warn/pass fail/fail warn/pass pass/pass warn/warn fail/fail warn/warn pass/pass cleaned
1C 51 132.7943818065906 40169799 95.3012205067167 10600000000 pass/pass fail/fail warn/pass fail/fail pass/pass pass/pass warn/warn fail/fail warn/warn pass/pass cleaned
2A 51 133.8210224632644 59990791 95.96230957773497 16000000000 pass/pass fail/fail warn/pass fail/fail pass/pass pass/pass warn/warn fail/fail warn/warn pass/pass cleaned
2C 51 135.9326951584075 40816343 95.22968500547375 11000000000 pass/pass fail/fail warn/warn fail/fail warn/pass pass/pass warn/warn fail/fail warn/warn pass/pass cleaned
6A 51 135.0995848516524 59731901 95.9628713760691 16000000000 pass/pass fail/fail warn/pass fail/fail warn/pass pass/pass warn/warn fail/fail warn/warn pass/pass cleaned
6C 51 133.0791847046676 44075930 95.46884644870812 11600000000 pass/pass fail/fail warn/warn fail/fail pass/pass pass/pass warn/warn fail/fail warn/warn pass/pass cleaned
SS1 51 132.6186958110153 50058144 93.48542280139144 13200000000 pass/pass fail/fail warn/pass fail/fail pass/pass pass/pass warn/warn fail/fail warn/warn pass/pass cleaned
SS2 50 130.4142505343276 32226953 92.96294155910216 8300000000 pass/pass fail/fail warn/warn fail/fail fail/pass pass/pass warn/warn fail/fail fail/warn pass/pass cleaned
1A 51 139.36127232383092 22856964 84.03839818359835 6200000000 pass/pass fail/fail warn/warn fail/fail warn/pass pass/pass warn/warn fail/fail warn/warn pass/pass dedup
1C 51 139.20760796494162 14350237 85.5109357195937 3900000000 pass/pass fail/fail warn/warn fail/fail warn/pass pass/pass warn/warn fail/fail warn/warn pass/pass dedup
2A 51 139.63674225763492 20427705 87.47183915601975 5600000000 pass/pass fail/fail warn/warn fail/fail pass/pass pass/pass warn/warn fail/fail warn/warn pass/pass dedup
2C 51 140.06660894623036 16008195 85.8902935536823 4400000000 pass/pass warn/fail warn/warn fail/fail warn/pass pass/pass warn/warn fail/fail warn/warn pass/pass dedup
6A 51 140.12831535289862 22072047 86.76742520720975 6100000000 pass/pass fail/fail warn/pass fail/fail warn/pass pass/pass warn/warn fail/fail warn/warn pass/pass dedup
6C 51 138.70188014663512 15982530 86.34695432901567 4300000000 pass/pass fail/fail warn/warn fail/fail warn/pass pass/pass warn/warn fail/fail warn/warn pass/pass dedup
SS1 51 138.43742527490832 19871170 76.04310701440284 5400000000 pass/pass fail/fail warn/warn fail/fail warn/pass pass/pass warn/warn fail/fail warn/warn pass/pass dedup
SS2 50 137.00414863944457 10599017 74.87558407316746 2800000000 warn/pass fail/fail warn/warn fail/fail fail/pass pass/pass warn/warn fail/fail warn/warn pass/pass dedup
1A 51.5 137.15324877221587 3889731 71.38358361496456 1066700000 pass/pass warn/fail warn/warn fail/warn warn/pass pass/pass warn/warn fail/fail warn/warn pass/pass host
1C 51.5 136.3974139190085 2403289 66.14407938418768 655400000 pass/pass warn/fail warn/pass fail/warn warn/pass pass/pass warn/warn fail/fail warn/warn pass/pass host
2A 51.5 137.1554242864068 3475895 68.43122611356736 953200000 pass/pass warn/fail warn/pass fail/warn warn/pass pass/pass warn/warn fail/fail warn/warn pass/pass host
2C 51 138.04065215417853 2753950 70.4192284764822 760100000 pass/pass warn/fail warn/warn fail/warn warn/pass pass/pass warn/warn fail/fail warn/warn pass/pass host
6A 51.5 137.85122092139602 3710804 75.93665530937844 1022800000 pass/pass fail/fail warn/pass fail/warn warn/pass pass/pass warn/warn fail/fail warn/warn pass/pass host
6C 51 136.21487199700704 2993680 69.91068399861464 815400000 pass/pass warn/fail warn/warn fail/warn warn/pass pass/pass warn/warn fail/fail fail/warn pass/pass host
SS1 51 136.87461348350567 4949336 64.81063504842454 1354600000 pass/pass warn/fail warn/warn fail/warn warn/pass pass/pass warn/warn fail/fail warn/warn pass/pass host
SS2 50 136.53833216842872 3036014 57.60979884422031 828900000 pass/pass fail/fail warn/warn fail/fail warn/warn pass/pass warn/warn fail/fail warn/warn pass/pass host
1A 51 150 78451849 94.93462712046542 23400000000 pass/pass fail/fail warn/warn fail/fail pass/pass pass/pass pass/pass fail/fail warn/warn fail/fail raw_concat
1C 51 150 47376584 94.16942468765848 14200000000 pass/pass fail/fail warn/pass fail/fail pass/pass pass/pass pass/pass fail/fail warn/warn fail/fail raw_concat
2A 51 150 71130999 94.95240423949721 21200000000 pass/pass fail/fail warn/pass fail/fail pass/pass pass/pass pass/pass fail/fail warn/warn fail/fail raw_concat
2C 51 150 49158246 93.42677984822313 14600000000 pass/pass fail/fail warn/warn fail/fail pass/pass pass/pass pass/pass fail/fail warn/warn fail/fail raw_concat
6A 51 150 71192267 94.817930410275 21200000000 pass/pass fail/fail warn/pass fail/fail pass/pass pass/pass pass/pass fail/fail warn/warn fail/fail raw_concat
6C 51 150 53655715 93.74559506275138 16000000000 pass/pass fail/fail warn/warn fail/fail pass/pass pass/pass pass/pass fail/fail warn/warn fail/fail raw_concat
SS1 51 150 61096105 91.61742981894301 18200000000 pass/pass fail/fail warn/warn fail/fail pass/pass pass/pass pass/pass fail/fail warn/warn fail/fail raw_concat
SS2 50 150 39146173 90.72005001240018 11600000000 warn/pass fail/fail warn/warn fail/fail warn/pass pass/pass pass/pass fail/fail fail/warn fail/fail raw_concat
1A 51.5 136.75124365577642 3910447 71.11941270165376 1069200000 pass/pass warn/fail warn/warn fail/warn warn/pass pass/pass warn/warn fail/fail warn/warn pass/pass ribo_initial
1C 51.5 135.99041018098131 2414227 65.9670095078155 656400000 pass/pass warn/fail warn/pass fail/warn warn/pass pass/pass warn/warn fail/fail warn/warn pass/pass ribo_initial
2A 51.5 136.78050495746933 3489482 68.34385523795663 954300000 pass/pass warn/fail warn/pass fail/warn warn/pass pass/pass warn/warn fail/fail warn/warn pass/pass ribo_initial
2C 51 137.62945690846857 2783656 69.58906768418473 766000000 pass/pass warn/fail warn/warn fail/warn warn/pass pass/pass warn/warn fail/fail warn/warn pass/pass ribo_initial
6A 51.5 137.53608104443185 3724130 75.83672415303403 1024100000 pass/pass warn/fail warn/pass fail/warn warn/pass pass/pass warn/warn fail/fail warn/warn pass/pass ribo_initial
6C 51 135.7960863042931 3007081 69.84557780252267 816400000 pass/pass warn/fail warn/warn fail/warn warn/pass pass/pass warn/warn fail/fail fail/warn pass/pass ribo_initial
SS1 51 136.47500906765998 4984748 64.34726643881876 1360300000 pass/pass warn/fail warn/warn fail/warn warn/pass pass/pass warn/warn fail/fail warn/warn pass/pass ribo_initial
SS2 50 136.18220473815117 3048362 57.36884524615661 830100000 pass/pass warn/fail warn/warn fail/fail warn/warn pass/pass warn/warn fail/fail warn/warn pass/pass ribo_initial
1A 52 130.88457029195476 522855 53.6722671779584 136800000 pass/pass warn/fail warn/pass fail/warn warn/warn pass/pass warn/warn fail/fail warn/warn pass/pass ribo_secondary
1C 51 130.53803565010548 364487 51.95283131410389 95000000 pass/pass fail/fail warn/pass fail/warn warn/warn pass/pass warn/warn fail/fail warn/warn pass/pass ribo_secondary
2A 51 131.40944459835464 498126 54.62744078695249 130800000 pass/pass fail/fail warn/pass fail/warn warn/warn pass/pass warn/warn fail/fail warn/warn pass/pass ribo_secondary
2C 51 133.09488399981325 406939 55.434787839904914 108200000 pass/pass warn/fail warn/pass fail/warn warn/warn pass/pass warn/warn fail/fail warn/warn pass/pass ribo_secondary
6A 51.5 132.6047124689108 496555 51.228388753494215 131600000 pass/pass warn/fail warn/pass fail/warn warn/warn pass/pass warn/warn warn/fail warn/warn pass/pass ribo_secondary
6C 51 130.40937472408305 441709 53.00244718657197 115100000 pass/pass warn/fail warn/warn fail/warn warn/warn pass/pass warn/warn fail/fail warn/warn pass/pass ribo_secondary
SS1 51 131.85058949359808 808745 45.528429335705326 213100000 pass/pass warn/fail warn/pass fail/warn warn/pass pass/pass warn/warn warn/warn warn/pass pass/pass ribo_secondary
SS2 50 134.13924738140605 726344 39.85374206117173 194800000 pass/pass fail/fail warn/warn fail/fail warn/warn pass/pass warn/warn warn/warn warn/pass pass/pass ribo_secondary
33 changes: 33 additions & 0 deletions data/2023-12-19_rna-seq-workflow/bracken_counts.tsv
Original file line number Diff line number Diff line change
@@ -0,0 +1,33 @@
name taxonomy_id taxonomy_lvl kraken_assigned_reads added_reads new_est_reads fraction_total_reads sample
Bacteria 2 D 253305 260 253565 0.97434 1A
Archaea 2157 D 3621 6 3627 0.01394 1A
Eukaryota 2759 D 656 15 671 0.00258 1A
Viruses 10239 D 2378 0 2378 0.00914 1A
Bacteria 2 D 193411 160 193571 0.97017 2C
Archaea 2157 D 2648 4 2652 0.01329 2C
Eukaryota 2759 D 1538 34 1572 0.00788 2C
Viruses 10239 D 1726 0 1726 0.00865 2C
Bacteria 2 D 234400 242 234642 0.97591 2A
Archaea 2157 D 3298 6 3304 0.01374 2A
Eukaryota 2759 D 200 5 205 8.5e-4 2A
Viruses 10239 D 2283 0 2283 0.0095 2A
Bacteria 2 D 165562 162 165724 0.97464 1C
Archaea 2157 D 2334 3 2337 0.01374 1C
Eukaryota 2759 D 224 5 229 0.00135 1C
Viruses 10239 D 1744 0 1744 0.01026 1C
Bacteria 2 D 361317 1704 363021 0.97672 SS1
Archaea 2157 D 5476 44 5520 0.01485 SS1
Eukaryota 2759 D 1597 213 1810 0.00487 SS1
Viruses 10239 D 1321 0 1321 0.00355 SS1
Bacteria 2 D 203068 204 203272 0.97502 6C
Archaea 2157 D 2547 4 2551 0.01224 6C
Eukaryota 2759 D 138 3 141 6.8e-4 6C
Viruses 10239 D 2514 0 2514 0.01206 6C
Bacteria 2 D 253688 2001 255689 0.9684 SS2
Archaea 2157 D 3455 42 3497 0.01324 SS2
Eukaryota 2759 D 505 122 627 0.00237 SS2
Viruses 10239 D 4218 0 4218 0.01598 SS2
Bacteria 2 D 241267 207 241474 0.97392 6A
Archaea 2157 D 3470 4 3474 0.01401 6A
Eukaryota 2759 D 345 6 351 0.00142 6A
Viruses 10239 D 2640 0 2640 0.01065 6A
Loading

0 comments on commit 35b2f01

Please sign in to comment.