Hi Everyone,
I have a metagenomics(16s) sequence reads and they are about 200 million in number. I want to filter out some reads by quality score and I am not able to find out how to get it. My reads look somewhat like this.
@HWI-ST1035:115:C0RG7ACXX:5:1101:1481:2050 1:N:0:
TACGGAGGGTCCGAGCGTTATCCGGAATTATTGGGTTTAAAGGGTCCGCAGGCGGGCAATTGAGTCAGGGGTGAAATGGTGCGGCTCAACCGTAGCACTGCCCTTGATACTGGTTGTCTTGAGTCATTGTGAAGTGGCCGGAATATGTAGG
+
B@CFFFFFHCCFHIJJJIJIIJJJJFHIJJJJJJJ75CGIIJII6=CHFFDDDDD557(9+>>A3>9?B9;?CCACD>>C>BBDDDCDAB?(<(+:@CDD849A:CCC>CC44<5@CAC:>:4:>(++:3(:+:>?<<B59?@DCE###
@HWI-ST1035:115:C0RG7ACXX:5:1101:1363:2055 1:N:0:
TACAGAGGATGCAAGCGTTATCCGGAATGATTGGGCGTAAAGCGTCTGTAGGTNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNCNNNNNNNNNNNNNGGN
+
@@@DFFFF?DFH<FHIIGIIIIIIIHEHGIIIIIIIIFFHHIII<CFH@EHI###################################################################################################
@HWI-ST1035:115:C0RG7ACXX:5:1101:1497:2056 1:N:0:
TACGTAGGGTGCAAGCGTTGTCCGGAATTACTGGGCGTAAAGAGCTCGTAGGCGGTTTGTCGCGTCGGCTGTGAAAACCAGCAGCTCAACTGTTGGCTTGCAGGCGATACGGGAAGACTTGAGTATTTCAGGGGGGACTGGAATTCCTGGG
+
?@@D4ADDFA<DDFGIIGI6GHIIIIIIIICGIIIDFFFIII<FCCFF;EIFFFC>8;=@:8(85>7&08;?ADDB(9<2<?<@BBBBBB:A::AB:>?BBAA<.959>B<BB5&8+9ABBBB@###########################
Please can someone help me to calculate the quality score of each read. I have searched a lot but could not find any information. This is my first time to deal with such data.
Thanks for any help!!!!
I have a metagenomics(16s) sequence reads and they are about 200 million in number. I want to filter out some reads by quality score and I am not able to find out how to get it. My reads look somewhat like this.
@HWI-ST1035:115:C0RG7ACXX:5:1101:1481:2050 1:N:0:
TACGGAGGGTCCGAGCGTTATCCGGAATTATTGGGTTTAAAGGGTCCGCAGGCGGGCAATTGAGTCAGGGGTGAAATGGTGCGGCTCAACCGTAGCACTGCCCTTGATACTGGTTGTCTTGAGTCATTGTGAAGTGGCCGGAATATGTAGG
+
B@CFFFFFHCCFHIJJJIJIIJJJJFHIJJJJJJJ75CGIIJII6=CHFFDDDDD557(9+>>A3>9?B9;?CCACD>>C>BBDDDCDAB?(<(+:@CDD849A:CCC>CC44<5@CAC:>:4:>(++:3(:+:>?<<B59?@DCE###
@HWI-ST1035:115:C0RG7ACXX:5:1101:1363:2055 1:N:0:
TACAGAGGATGCAAGCGTTATCCGGAATGATTGGGCGTAAAGCGTCTGTAGGTNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNNCNNNNNNNNNNNNNGGN
+
@@@DFFFF?DFH<FHIIGIIIIIIIHEHGIIIIIIIIFFHHIII<CFH@EHI###################################################################################################
@HWI-ST1035:115:C0RG7ACXX:5:1101:1497:2056 1:N:0:
TACGTAGGGTGCAAGCGTTGTCCGGAATTACTGGGCGTAAAGAGCTCGTAGGCGGTTTGTCGCGTCGGCTGTGAAAACCAGCAGCTCAACTGTTGGCTTGCAGGCGATACGGGAAGACTTGAGTATTTCAGGGGGGACTGGAATTCCTGGG
+
?@@D4ADDFA<DDFGIIGI6GHIIIIIIIICGIIIDFFFIII<FCCFF;EIFFFC>8;=@:8(85>7&08;?ADDB(9<2<?<@BBBBBB:A::AB:>?BBAA<.959>B<BB5&8+9ABBBB@###########################
Please can someone help me to calculate the quality score of each read. I have searched a lot but could not find any information. This is my first time to deal with such data.
Thanks for any help!!!!
Comment