Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	689,437	G→T	intergenic (+198522/‑7271)	unknown → / → unknown	tRNA‑SeC(p)‑TCA/SSU rRNA ## 16S rRNA, small subunit ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	689,437	0	G	T	100.0%	15.3 / NA	6	intergenic (+198522/‑7271)	unknown/unknown	tRNA‑SeC(p)‑TCA/SSU rRNA ## 16S rRNA, small subunit ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (4/2);  total (4/2)

CGGTAGGGCTAAATACAAAGACAGTGGAGGAGTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGGATATGTGGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCTTTATAGAGTAGCCTAAAAACTGCACGTTTTACGTGCGGATTGTTATTTTATGATTTTTTTGTACAGGGTGTAATTTTACCTGAA  >  1/689411‑689675                           |                                                                                                                                                                                                                                               cGGTAGGGCTAAATACAAAGGTAGTGTAGGAGTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCTTTATAGAGTAGCCTAAAAACTGCACGTTTTACGTGCGGATTGTTATTTTATGATTTTTTTGTACAGGgtg                >  1:901509/1‑251 (MQ=255)   gTAGGGCTAAATACAAAGGTAGTGTAGGAGTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCTTTATAGAGTAGCCTAAAAACTGCACGTTTTACGTGCGGATTGTTATTTTATGATTTTTTTGTACAGGGTGTa              >  1:239399/1‑251 (MQ=255)       ggCTAAATACAAAGGTAGTGTAGGAGTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCTTTATAGAGTAGCCTAAAAACTGCACGTTTTACGTGCGGATTGTTATTTTATGATTTTTTTGTACAGGGTGTAAttt          >  1:88441/1‑251 (MQ=255)       ggCTAAATACAAAGGTAGTGTAGGAGTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCTTTATAGAGTAGCCTAAAAACTGCACGTTTTACGTGCGGATTGTTATTTTATGATTTTTTTGTACAGGGTGTAAt            <  1:173696/249‑1 (MQ=255)             aTACAAAGGTAGTGTAGGAGTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCTTTATAGAGTAGCCTAAAAACTGCACGTTTTACGTGCGGATTGTTATTTTATGATTTTTTTGTACAGGGTGTAATTTTAcc      <  1:91189/249‑1 (MQ=255)                caaaggtaGTGTAGGAGTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGAATATGTGGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCTTTATAGAGTAGCCTAAAAACTGCACGTTTTACGTGCGGATTGTTATTTTATGATTTTTTTGTACAGGGTGTAATTTTACCTGaa  >  1:555476/8‑250 (MQ=255)                           |                                                                                                                                                                                                                                               CGGTAGGGCTAAATACAAAGACAGTGGAGGAGTATATAAGGAATCAGGAAAAGGAAGACATGATTCAGGATAATTTATCGAAGAAAGGATATGTGGACCCCTTTAAGGGGTAGGTCATAGAGTTATTTCACGGTTGTTGGTTTTCATATGCCCTTTGAAGGGCATGTGAAGTACTAGGTCTTTATAGAGTAGCCTAAAAACTGCACGTTTTACGTGCGGATTGTTATTTTATGATTTTTTTGTACAGGGTGTAATTTTACCTGAA  >  1/689411‑689675

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑