%matplotlib inline
from sequana import *
f = FastQ(sequana_data('test.fastq', 'testing'))
f.count_reads()
250
f.identifier
Identifier (Illumina_1.8+)
f.next()
{'identifier': b'@HISEQ:426:C5T65ACXX:5:2302:1943:2127 1:N:0:GTGAAA', 'quality': b'BCCFFFFFHHHHHIIJJJJJJIIJJJJJJJJFHIHIJJJIJIIIIGHFFFFFFEEEEEEEDDDDDFDDDDDDDDD>CDDEDEEDDDDDDCCDDDDDDDDCD', 'sequence': b'CTTACCTTCGCATCAAGAGGTACCAACATGGAGGATGAGGCGGACCAATACTTTTCACATGATGATCCAATTAGTAGTGATCAATCCAGGTTCGGATGGTT'}
fq = FastQC(sequana_data('test.fastq', 'testing'))
fq.boxplot_quality()