Iterative mapping
In [2]:
r_enz = 'HindIII'
In [4]:
from pytadbit.mapping.full_mapper import full_mapping
In [5]:
full_mapping('media/storage/db/reference_genome/Homo_sapiens/hg38/hg38.gem' ,
reads=2, clean=True, r_enz=r_enz,
frag_map=False,
out_map_dir='results/{0}/01_mapping/map{0}_r1' .format(r_enz),
fastq_path='/media/storage/FASTQs/K562_{0}_1.fastq'.format(r_enz),
windows=(1, 25), (1, 30), (1, 35), (1, 40), (1, 45), (1, 50), (1, 75), (25, 75))
In [6]:
full_mapping('media/storage/db/reference_genome/Homo_sapiens/hg38/hg38.gem' ,
reads=2, clean=True, r_enz=r_enz,
frag_map=False,
out_map_dir='results/{0}/01_mapping/map{0}_r2' .format(r_enz),
fastq_path='/media/storage/FASTQs/K562_{0}_2.fastq'.format(r_enz),
windows=(1, 25), (1, 30), (1, 35), (1, 40), (1, 45), (1, 50), (1, 75), (25, 75))
In [8]:
! mkdir -p results/
In [14]:
! mkdir -p results/HindIII/01_mapping
In [10]:
! mkdir -p results/MboI/01_mapping
In [18]:
! ls results/HindIII/
In [12]:
full_mapping('media/storage/db/reference_genome/Homo_sapiens/hg38/hg38.gem' ,
nthreads=2, clean=True, r_enz=r_enz,
frag_map=True,
out_map_dir='results/{0}/01_mapping/map{0}_r1' .format(r_enz),
fastq_path='/media/storage/FASTQs/K562_{0}_1.fastq'.format(r_enz))
In [27]:
ls results/HindIII/01_mapping/mapHindIII_r1/
In [21]:
full_mapping('media/storage/db/reference_genome/Homo_sapiens/hg38/hg38.gem' ,
nthreads=2, clean=True, r_enz=r_enz,
frag_map=True,
out_map_dir='results/{0}/01_mapping/map{0}_r1' .format(r_enz),
fastq_path='/media/storage/FASTQs/K562_{0}_1.fastq'.format(r_enz))
In [ ]:
full_mapping('media/storage/db/reference_genome/Homo_sapiens/hg38/hg38.gem' ,
nthreads=2, clean=True, r_enz=r_enz,
frag_map=True,
out_map_dir='results/{0}/01_mapping/map{0}_r2' .format(r_enz),
fastq_path='/media/storage/FASTQs/K562_{0}_2.fastq'.format(r_enz))
In [28]:
ls
In [29]:
! head results/HindIII/01__mapping/mapHindIII_r2/K562_HindIII_2_full_1-end.map
In [ ]: