diff --git a/ProcessFASTQ.ipynb b/ProcessFASTQ.ipynb index 3943225353d510c17aed04324529aa70099979b4..7a78e1cd3aa39a30cf8c354a00fd7a7360713fb3 100644 --- a/ProcessFASTQ.ipynb +++ b/ProcessFASTQ.ipynb @@ -683,7 +683,7 @@ }, { "cell_type": "code", - "execution_count": 240, + "execution_count": 242, "metadata": { "collapsed": false }, @@ -722,7 +722,11 @@ " #sub_iter = (rec for rec in filt_iter \\\n", " # if matches.readline().split('\\t')[1] == '4')\n", " \n", - " \n", + " for record in filt_iter:\n", + " arr = matches.readline().split('\\t')\n", + " seq_id = arr[0]\n", + " flag_field = arr[1]\n", + " \n", "\n", " # Write back fastq\n", " SeqIO.write(sub_iter,substracted,\"fastq\")"