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

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)