analyze


In [1]:
! ls


00_fastq_QC.ipynb  learning_python.ipynb

In [2]:
! pwd


/media/storage/Notebooks2/mariya

In [3]:
! head /media/storage/FASTQs/K562_HindIII_1.fastq


@NS500645:59:HCL32BGXY:1:11101:14163:1054 1:N:0:GCCAAT
CATTCNTAAAGAAAAGAATTTTCAACNCAGAATTTCATATCCAGCCAACTAAGCTAGCTTCAAGGAAATACATTT
+
AAAAA#EEAEEEEEEEEEEAE/EEEE#EEEEEEEEEEEEEEEEAE</EE/EEEEAEEEEEEEEEEEAAEE</EEE
@NS500645:59:HCL32BGXY:1:11101:4416:1054 1:N:0:GCCAAT
CGAGTNAGGGAAGCTAGCTTCCATTCNTTTGCTTCTGTTGTGTGTTTTTTCTTTTGTTTTTTTTTTGTTTTGGTT
+
AAAAA#EEEEEEEEEEEEEEAEEEE/#EEEEEEEEEEEEEEEEEAEEEEEEEEEEEEEEEEEEEEE/<EEEEEAE
@NS500645:59:HCL32BGXY:1:11101:14193:1055 1:N:0:GCCAAT
AAAGTNCCCTGCATGATAGCATTTGTATAAAGTTCAAAACAGACCAAATGGATCTCTAATTTGTAGAAGGTCAGG

In [5]:
! wc -l /media/storage/FASTQs/K562_HindIII_1.fastq


40000000 /media/storage/FASTQs/K562_HindIII_1.fastq

FASQ Quaulity check using TADbit


In [6]:
from pytadbit.utils.fastq_utils import quality_plot

In [8]:
quality_plot('/media/storage/FASTQs/K562_HindIII_1.fastq', r_enz='HindIII', nreads=1000000)


Out[8]:
(0.8355, 28.7747)

In [9]:
quality_plot('/media/storage/FASTQs/K562_MboI_1.fastq', r_enz='MboI', nreads=1000000)


Out[9]:
(2.7736, 35.1285)