001、
root@PC1:/home/test# ls
a.fastq test.py
root@PC1:/home/test# cat a.fastq ## 测试fastq文件
@DJB775P1:248:D0MDGACXX:7:1202:12362:49613
TGCTTACTCTGCGTTGATACCACTGCTTAGATCGGAAGAGCACACGTCTGAA
+
JJJJJIIJJJJJJHIHHHGHFFFFFFCEEEEEDBD?DDDDDDBDDDABDDCA
@DJB775P1:248:D0MDGACXX:7:1202:12782:49716
CTCTGCGTTGATACCACTGCTTACTCTGCGTTGATACCACTGCTTAGATCGG
+
IIIIIIIIIIIIIIIHHHHHHFFFFFFEECCCCBCECCCCCCCCCCCCCCCC
root@PC1:/home/test# cat test.py ## 测试程序
#!/usr/bin/python
in_file = open("a.fastq", "r")
dict1 = {}
idx = 0
for i in in_file:
idx += 1
i = i.strip()
if idx % 4 == 1:
key = i
dict1[key] = ""
elif idx % 4 == 2:
dict1[key] = len(i)
for i in dict1:
print(i + "-------" + str(dict1[i]))
in_file.close()
root@PC1:/home/test# python test.py ## 执行程序
@DJB775P1:248:D0MDGACXX:7:1202:12362:49613-------52
@DJB775P1:248:D0MDGACXX:7:1202:12782:49716-------52
参考：https://www.jianshu.com/p/5ee54bea4cb0
标签：DJB775P1,dict1,python,fastq,碱基,1202,test,home From： https://www.cnblogs.com/liujiaxin2018/p/16588622.html