Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	443,412	T→C	intergenic (+79995/‑128607)	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	443,412	0	T	C	100.0%	46.0 / NA	20	intergenic (+79995/‑128607)	unknown/unknown	tRNA‑SeC(p)‑TCA/SSU rRNA ## 16S rRNA, small subunit ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (8/12);  total (8/12)

CTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGTTCTACTTATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTCCGTTAAATTTGAGGGTTGTTAGTTGTTGGTGATGTGTAAAATTTACGTAGTGTTACTTAATTTAAAACCATATAAATTTAGAAAATTTATTTCTAGTTTTATATGGTTATTTTTCTAAAATCTCGCGGTCGGCTACACGCAAAACTGCTT  >  1/443188‑443538                                                                                                                                                                                                                                  |                                                                                                                               cTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGGAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTTTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAg                                                                                                      >  1:793873/1‑251 (MQ=255)         cTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttt                                                                                              <  1:126024/251‑5 (MQ=255)         cTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttt                                                                                              <  1:622206/251‑5 (MQ=255)              tgtTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttttatta                                                                                         <  1:821071/251‑10 (MQ=255)              tgtTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttttatta                                                                                         <  1:783339/251‑10 (MQ=255)                    aGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTTTTTATGTTTATATGATGATTCTGCGTTAAATTTTAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAgtgt                                                                                                   >  1:322279/1‑235 (MQ=255)                      caaacaaaGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttttattaaattattc                                                                                 <  1:543917/251‑18 (MQ=255)                       aaacaaaGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGGAATAC‑CCATCGTTTTTGTTTATGTTTATATGATGATTTTGCGTTAAATTTTAGGGGGGTTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttt                                                                                              >  1:158898/1‑233 (MQ=255)                          caaaGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttttattaaattattcatac                                                                             <  1:859192/251‑22 (MQ=255)                                      gctgAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttttattaaatta                                                                                    <  1:482745/232‑15 (MQ=255)                                         gAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTTTGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttt                                                                                              >  1:826309/1‑215 (MQ=255)                                                            tGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGGAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgtt                                                                                               >  1:775368/1‑196 (MQ=255)                                                                                     gCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttttatta                                                                                         >  1:440256/1‑171 (MQ=255)                                                                                                        aaaTAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttt                                                                                              <  1:258670/156‑5 (MQ=255)                                                                                                          aTAGTAGTTAAGCTCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttttattaaat                                                                                      <  1:284319/162‑13 (MQ=255)                                                                                                                     gctCTACTCATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttttattaaat                                                                                      <  1:366774/149‑13 (MQ=255)                                                                                                                                                      aaaCGGCGGGTTAATACTCCGCCGTTTTTCTTTTCATGTACATGATGATTCTGCTTTAAATTTGAGGGTAATTAGTCGCTGATGATGTGTAAAATTGACGTAGCGCTATTTGGCGTTAAACCACATAAGTTTAGAAATGTTGTTTCTAGGTTTATGTGGTTATTTTTCTAAAATCTCGCGGTCGGCAACACgtccggtaatgt  <  1:82718/203‑12 (MQ=255)                                                                                                                                                                           tcgtTTTTGTTTATGTTTATATGATGATTCTGCGTTAAATTTGAGGGCGATTAGTCGTTATTCATGTGTAAAATTTACGTAGTGttgttt                                                                                              <  1:72476/90‑5 (MQ=255)                                                                                                                                                                                                tgatgaTTCTGCTTTAAATTTGAGGGTAATTAGTCGCTGATGATGTGTAAAATTGACGTAGCGCTATTTGGCGTTAAACCACATAAGTTTAGAAATGTTGTTTCTAGGTTTATGTGGTTATTTTTCTAAAATCTCGCGGTCGGCAacac              >  1:770647/1‑149 (MQ=255)                                                                                                                                                                                                          gctttAAATTTGAGGGTAATTAGTCGCTGATGATGTGTAAAATTGACGTAGCGCTATTTGGCGTTAAACCACATAAGTTTAGAAATGTTGTTTCTAGGTTTATGTGGTTATTTTTCTAAAATCTCGCGGTCGGCAACACgtccggtaa     >  1:450714/4‑140 (MQ=255)                                                                                                                                                                                                                                  |                                                                                                                               CTGAAAGGCTGAGTGTTGAAGCAAACAAAGCTTAACTGCTGAATAAGCAGCTTAGAAAATGGTAGCTCTTACTGAAAGTTCACTGCCGAACAGGCAGCTTAGAAAATAGTAGTTAAGTTCTACTTATTCACTGCCGAACAGGCAGACAAAAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTCCGTTAAATTTGAGGGTTGTTAGTTGTTGGTGATGTGTAAAATTTACGTAGTGTTACTTAATTTAAAACCATATAAATTTAGAAAATTTATTTCTAGTTTTATATGGTTATTTTTCTAAAATCTCGCGGTCGGCTACACGCAAAACTGCTT  >  1/443188‑443538

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