# NOT RUN {
# Show the counts for each length
seq_lengths(mammals_16S)
# Plot the distribution of lengths
seqLengths <- seq_lengths(mammals_16S)
barplot(seqLengths,
main = "Frequency of sequence lengths for 16S mammal primers",
xlab="Amplicon length (in bp)",
ylab=("Frequency"))
# }
Run the code above in your browser using DataLab