WGS sequencing of milk microbiome
raw reads were filtered using kneaddata (Huttenhower Lab) to eliminate human reads, reads shorter than 85 bp, and low quality reads (SLIDINGWINDOW: 4:20), and only paired reads were retained for analysis