BSMAP-methatio-methykit workflow

In [4]:
#Setting Variables
#file ID
fid="CgLarv_T3D3_sept"
#where is bsmap
bsmap="/Users/Shared/Apps/bsmap-2.73/"
#fastq files location R1 location
R1="/Volumes/web/trilobite/Crassostrea_gigas_HTSdata/batterbox/BiGo_Larvae/filtered_Bs_CgLarve_T3D3_GCCAAT_L007_R1.fastq.gz"
#genome file 
genome="/Volumes/web/whale/ensembl/ftp.ensemblgenomes.org/pub/release-21/metazoa/fasta/crassostrea_gigas/dna/Crassostrea_gigas.GCA_000297895.1.21.dna_sm.genome.fa"
In [5]:
cd /Volumes/web/Mollusk/bs_larvae_exp/
/Volumes/web/Mollusk/bs_larvae_exp

In [6]:
mkdir {fid}
In [7]:
cd {fid}
/Volumes/web/Mollusk/bs_larvae_exp/CgLarv_T3D3_sept

In [8]:
!{bsmap}bsmap -a {R1} -d {genome} -o bsmap_out.sam -p 1

BSMAP v2.73
Start at:  Wed Feb 12 17:36:53 2014

Input reference file: /Volumes/web/whale/ensembl/ftp.ensemblgenomes.org/pub/release-21/metazoa/fasta/crassostrea_gigas/dna/Crassostrea_gigas.GCA_000297895.1.21.dna_sm.genome.fa 	(format: FASTA)
Load in 7658 db seqs, total size 557717710 bp. 10 secs passed
total_kmers: 43046721
Create seed table. 35 secs passed
max number of mismatches: read_length * 8% 	max gap size: 0
kmer cut-off ratio:5e-07
max multi-hits: 100	max Ns: 5	seed size: 16	index interval: 4
quality cutoff: 0	base quality char: '!'
min fragment size:28	max fragemt size:500
start from read #1	end at read #4294967295
additional alignment: T in reads => C in reference
mapping strand: ++,-+
Single-end alignment(1 threads)
Input read file: /Volumes/web/trilobite/Crassostrea_gigas_HTSdata/batterbox/BiGo_Larvae/filtered_Bs_CgLarve_T3D3_GCCAAT_L007_R1.fastq.gz 	(format: gzipped FASTQ)
Output file: bsmap_out.sam	 (format: SAM)
Thread #0: 	50000 reads finished. 41 secs passed
Thread #0: 	100000 reads finished. 46 secs passed
Thread #0: 	150000 reads finished. 52 secs passed
Thread #0: 	200000 reads finished. 57 secs passed
Thread #0: 	250000 reads finished. 63 secs passed
Thread #0: 	300000 reads finished. 68 secs passed
Thread #0: 	350000 reads finished. 73 secs passed
Thread #0: 	400000 reads finished. 79 secs passed
Thread #0: 	450000 reads finished. 85 secs passed
Thread #0: 	500000 reads finished. 90 secs passed
Thread #0: 	550000 reads finished. 95 secs passed
Thread #0: 	600000 reads finished. 101 secs passed
Thread #0: 	650000 reads finished. 107 secs passed
Thread #0: 	700000 reads finished. 112 secs passed
Thread #0: 	750000 reads finished. 117 secs passed
Thread #0: 	800000 reads finished. 122 secs passed
Thread #0: 	850000 reads finished. 127 secs passed
Thread #0: 	900000 reads finished. 132 secs passed
Thread #0: 	950000 reads finished. 138 secs passed
Thread #0: 	1000000 reads finished. 143 secs passed
Thread #0: 	1050000 reads finished. 148 secs passed
Thread #0: 	1100000 reads finished. 154 secs passed
Thread #0: 	1150000 reads finished. 159 secs passed
Thread #0: 	1200000 reads finished. 164 secs passed
Thread #0: 	1250000 reads finished. 169 secs passed
Thread #0: 	1300000 reads finished. 174 secs passed
Thread #0: 	1350000 reads finished. 179 secs passed
Thread #0: 	1400000 reads finished. 185 secs passed
Thread #0: 	1450000 reads finished. 190 secs passed
Thread #0: 	1500000 reads finished. 195 secs passed
Thread #0: 	1550000 reads finished. 201 secs passed
Thread #0: 	1600000 reads finished. 206 secs passed
Thread #0: 	1650000 reads finished. 211 secs passed
Thread #0: 	1700000 reads finished. 216 secs passed
Thread #0: 	1750000 reads finished. 221 secs passed
Thread #0: 	1800000 reads finished. 227 secs passed
Thread #0: 	1850000 reads finished. 233 secs passed
Thread #0: 	1900000 reads finished. 240 secs passed
Thread #0: 	1950000 reads finished. 246 secs passed
Thread #0: 	2000000 reads finished. 252 secs passed
Thread #0: 	2050000 reads finished. 258 secs passed
Thread #0: 	2100000 reads finished. 263 secs passed
Thread #0: 	2150000 reads finished. 268 secs passed
Thread #0: 	2200000 reads finished. 273 secs passed
Thread #0: 	2250000 reads finished. 278 secs passed
Thread #0: 	2300000 reads finished. 282 secs passed
Thread #0: 	2350000 reads finished. 287 secs passed
Thread #0: 	2400000 reads finished. 292 secs passed
Thread #0: 	2450000 reads finished. 297 secs passed
Thread #0: 	2500000 reads finished. 301 secs passed
Thread #0: 	2550000 reads finished. 306 secs passed
Thread #0: 	2600000 reads finished. 311 secs passed
Thread #0: 	2650000 reads finished. 316 secs passed
Thread #0: 	2700000 reads finished. 320 secs passed
Thread #0: 	2750000 reads finished. 325 secs passed
Thread #0: 	2800000 reads finished. 330 secs passed
Thread #0: 	2850000 reads finished. 335 secs passed
Thread #0: 	2900000 reads finished. 339 secs passed
Thread #0: 	2950000 reads finished. 344 secs passed
Thread #0: 	3000000 reads finished. 349 secs passed
Thread #0: 	3050000 reads finished. 354 secs passed
Thread #0: 	3100000 reads finished. 359 secs passed
Thread #0: 	3150000 reads finished. 364 secs passed
Thread #0: 	3200000 reads finished. 368 secs passed
Thread #0: 	3250000 reads finished. 373 secs passed
Thread #0: 	3300000 reads finished. 378 secs passed
Thread #0: 	3350000 reads finished. 383 secs passed
Thread #0: 	3400000 reads finished. 387 secs passed
Thread #0: 	3450000 reads finished. 392 secs passed
Thread #0: 	3500000 reads finished. 397 secs passed
Thread #0: 	3550000 reads finished. 402 secs passed
Thread #0: 	3600000 reads finished. 407 secs passed
Thread #0: 	3650000 reads finished. 411 secs passed
Thread #0: 	3700000 reads finished. 416 secs passed
Thread #0: 	3750000 reads finished. 421 secs passed
Thread #0: 	3800000 reads finished. 426 secs passed
Thread #0: 	3850000 reads finished. 431 secs passed
Thread #0: 	3900000 reads finished. 436 secs passed
Thread #0: 	3950000 reads finished. 440 secs passed
Thread #0: 	4000000 reads finished. 445 secs passed
Thread #0: 	4050000 reads finished. 450 secs passed
Thread #0: 	4100000 reads finished. 455 secs passed
Thread #0: 	4150000 reads finished. 459 secs passed
Thread #0: 	4200000 reads finished. 464 secs passed
Thread #0: 	4250000 reads finished. 469 secs passed
Thread #0: 	4300000 reads finished. 474 secs passed
Thread #0: 	4350000 reads finished. 479 secs passed
Thread #0: 	4400000 reads finished. 483 secs passed
Thread #0: 	4450000 reads finished. 488 secs passed
Thread #0: 	4500000 reads finished. 493 secs passed
Thread #0: 	4550000 reads finished. 498 secs passed
Thread #0: 	4600000 reads finished. 503 secs passed
Thread #0: 	4650000 reads finished. 507 secs passed
Thread #0: 	4700000 reads finished. 512 secs passed
Thread #0: 	4750000 reads finished. 519 secs passed
Thread #0: 	4800000 reads finished. 525 secs passed
Thread #0: 	4850000 reads finished. 532 secs passed
Thread #0: 	4900000 reads finished. 538 secs passed
Thread #0: 	4950000 reads finished. 545 secs passed
Thread #0: 	5000000 reads finished. 551 secs passed
Thread #0: 	5050000 reads finished. 557 secs passed
Thread #0: 	5100000 reads finished. 564 secs passed
Thread #0: 	5150000 reads finished. 570 secs passed
Thread #0: 	5200000 reads finished. 576 secs passed
Thread #0: 	5250000 reads finished. 583 secs passed
Thread #0: 	5300000 reads finished. 589 secs passed
Thread #0: 	5350000 reads finished. 596 secs passed
Thread #0: 	5400000 reads finished. 602 secs passed
Thread #0: 	5450000 reads finished. 608 secs passed
Thread #0: 	5500000 reads finished. 613 secs passed
Thread #0: 	5550000 reads finished. 618 secs passed
Thread #0: 	5600000 reads finished. 622 secs passed
Thread #0: 	5650000 reads finished. 627 secs passed
Thread #0: 	5700000 reads finished. 632 secs passed
Thread #0: 	5750000 reads finished. 637 secs passed
Thread #0: 	5800000 reads finished. 641 secs passed
Thread #0: 	5837666 reads finished. 645 secs passed
Total number of aligned reads: 4263808 (73%)
Done.
Finished at Wed Feb 12 17:47:38 2014
Total time consumed:  645 secs

In [9]:
!python {bsmap}methratio.py -d {genome} -u -z -g -o methratio_out.txt -s {bsmap}samtools bsmap_out.sam
@ Wed Feb 12 17:47:46 2014: reading reference /Volumes/web/whale/ensembl/ftp.ensemblgenomes.org/pub/release-21/metazoa/fasta/crassostrea_gigas/dna/Crassostrea_gigas.GCA_000297895.1.21.dna_sm.genome.fa ...
@ Wed Feb 12 17:52:41 2014: reading bsmap_out.sam ...
[samopen] SAM header is present: 7658 sequences.
[sam_read1] reference 'NM:i:0' is recognized as '*'.
Parse error at line 1765919: unmatched CIGAR operation
@ Wed Feb 12 17:53:54 2014: combining CpG methylation from both strands ...
@ Wed Feb 12 17:54:23 2014: writing methratio_out.txt ...
@ Wed Feb 12 17:58:22 2014: done.
total 1359950 valid mappings, 12989572 covered cytosines, average coverage: 1.15 fold.

In []: