Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,034,982	3 bp→GTA	intergenic (‑28526/‑235)	unknown ← / → unknown	tRNA‑Gly‑GCC/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	1,034,982	0	A	G	100.0%	55.3 / NA	25	intergenic (‑28526/‑237)	unknown/unknown	tRNA‑Gly‑GCC/SSU rRNA ## 16S rRNA, small subunit ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (12/13);  total (12/13)
*	1	1,034,983	0	C	T	100.0%	53.9 / NA	25	intergenic (‑28527/‑236)	unknown/unknown	tRNA‑Gly‑GCC/SSU rRNA ## 16S rRNA, small subunit ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/13);  total (12/13)
*	1	1,034,984	0	T	A	100.0%	54.4 / NA	25	intergenic (‑28528/‑235)	unknown/unknown	tRNA‑Gly‑GCC/SSU rRNA ## 16S rRNA, small subunit ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (12/13);  total (12/13)

GTGGTGCGGCACTTGGCGGTTTGATTGGAAGCCAAGTGGGCAAATAATTTATTTTACTTAACTCGGAAAGCCAAAAGTGCGGTCAGAAATTCAAGCGTTTTTCTGACCGCATTTTTTTGATCTCCATAATTTTTTATCTCAAAGGAAAACACAATTTGGCGAGTTTATGAGCAAATAATAAAAAAGTTCAATTTTTTGTAAAAAAACTCTTGCAAAG‑AATTAGTAAATACCTATAATACGCGGGCACAACGACGCGCCGTTGTAGGGTTTTATAGACAGCGCGTCGTTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATAATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCG  >  1/1034777‑1035162                                                                                                                                                                                                              |||                                                                                                                                                                                    gtggtgCGGCACTTGGCGGTTTGATTGGAAGCCAAGTGGGCAAATAATTTATTTTACTTAACTCGGAAAGCCAAAAGTGCGGTCAGAAATTCAAGCGTTTTTCTGACCGCATTTTTTTGATCTCCATAATTTTTTATCTCAAAGGAAAACACAATTTGGCGAGTTTATGAGCAAATAATAAAAAAGTTCAATTTTTTGTAAAAAAGTACTTGCAAGGTAAGGA‑TAAAAGCCTATAATACGCGGGCACAacg                                                                                                                                         >  1:198120/1‑251 (MQ=38)    gtgCGGCACTTAGCGGTTTGATTGGAAGCCAAGTGGGCAAATAATTTATTTTACTTAACTCGGAAAGCCAAAAGTGCGGTCAGAAATTCAAGCGTTTTTCTGACCGCATTTTTTTGATCTCCATAATTTTTTATCTCAAAGGAAAACACAATTTGGCGAGTTTATGAGCAAATAATAAAAAAGTTCAATTTTTTGTAAAAAAGTACTTGCAAGGTAAGGA‑TAAAAGCCTATAATACGCGGGCACAacgacg                                                                                                                                      >  1:421433/1‑251 (MQ=38)       cGGCACTTGGCGGTTTGATTGGAAGCCAAGTGGGCAAATAATTTATTTTACTTAACTCGGAAAGCCAAAAGTGCGGTCAGAAATTCAAGCGTTTTTCTGACCGCATTTTTTTGATCTCCATAATTTTTTATCTCAAAGGAAAACACAATTTGGCGAGTTTATGAGCAAATAATAAAAAAGTTCAATTTTTTGTAAAAAAGTACTTGCAAGGTAAGGA‑TAAAAGCCTATAATACGCGGGCACAACGACgcgc                                                                                                                                   <  1:169811/251‑1 (MQ=38)                 cGGTTTGATTGGAAGCCAAGTGGGCAAATAATTTATTTTACTTAACTCGGAAAGCCAAAAGTGCGGTCAGAAATTCAAGCGTTTTTCTGACCGCATTTTTTTGATCTCCATAATTTTTTATCTCAAAGGAAAACACAATTTGGCGAGTTTATGAGCAAATAATAAAAAAGTTCAATTTTTTGTAAAAAAGTACTTGCAAGGTAAGGA‑TAAAAGCCTATAATACGCGGGCACAACGACGCGCCGTTGTAggg                                                                                                                         <  1:77441/251‑1 (MQ=37)                 cGGTTTGATTGGAAGCCAAGTGGGCAAATAATTTATTTTACTTAACTCGGAAAGCCAAAAGTGCGGTCAGAAATTCAAGCGTTTTTCTGACCGCA‑TTTTTTGATCTCCATAATTTTTTATCTCAAAGGAAAACACAATTTGGCGAGTTTATGAGCAAATAATAAAAAAGTTCAATTTTTTGTAAAAAAGTACTTGCAAGGTAAGGA‑TAAAAGCCTATAATACGCGGGCACAACGACGCGCCGTTGTAggg                                                                                                                         <  1:554485/250‑1 (MQ=37)                         ttGGAAGCCAAGTGGGCAAATAATTTATTTTACTTAACTCGGAAAGCCAAAAGTGCGGTCAGAAATTCAAGCGTTTTTCTGACCGCATTTTTTTGATCTCCATAATTTTTTATCTCAAAGGAAAACACAATTTGGCGAGTTTATGAGCAAATAATAAAAAAGTTCAATTTTTTGTAAAAAAGTACTTGCAAGGTAAGGA‑TAAAAGCCTATAATACGCGGGCACAACGACGCGCCGTTGTAGGGTTTTATAg                                                                                                                 <  1:107148/251‑1 (MQ=37)                            gAAGCCAAGTGGGCAAATAATTTATTTTACTTAACTCGGAAAGCCAAAAGTGCGGTCAGAAATTCAAGCGTTTTTCTGACCGCATTTTTTTGATCTCCATAATTTTTTATCTCAAAGGAAAACACAATTTGGCGAGTTTATGAGCAAATAATAAAAAAGTTCAATTTTTTGTAAAAAAGTACTTGCAAGGTAAGGA‑TAAAAGCCTATAATACGCGGGCACAACGACGCGCCGTTGTAggg                                                                                                                         >  1:806376/1‑240 (MQ=37)                            gAAGCCAAGTGGGCAAATAATTTATTTTACTTAACTCGGAAAGCCAAAAGCGCGGTCAGAAATTCAAGCGTTTTTCTGACCGCATTTTTTTGATCTCCATAATTTTTTATCTCAAAGGAAAACACAATTTGGCGAGTTTATGAGCAAATAATAAAAAAGTTCAATTTTTTGTAAAAAAGTACTTGCAAGGTAAGGA‑TAAAAGCCTATAATACGCGGGCACAACGACGCGCCGTTGTAGGGTTTTATAGACa                                                                                                              >  1:532233/1‑251 (MQ=25)                                          aaaTAATTTATTTTACTTAACTCGGAAAGCCAAAAGTGCGGTCAGAAATTCAAGCGTTTTTCTGACCGCATTTTTTTGATCTCCATAATTTTTTATCTCAAAGGAAAACACAATTTGGCGAGTTTATGAGCAAATAATAAAAAAGTTCAATTTTTTGTAAAAAAGTACTTGCAAGGTAAGGA‑TAAAAGCCTATAATACGCGGGCACAACGACGCGCCGTTGTAGGGTTTtat                                                                                                                   >  1:867558/1‑232 (MQ=25)                                          aaaTAATTTATTTTACTTAACTCGGAAAGCCAAAAGTGCGGTCAGAAATTCAAGCGTTTTTCTGACCGCATTTTTTTGATCTCCATAATTTTTTATCTCAAAGGAAAACACAATTTGGCGAGTTTATGAGCAAATAATAAAAAAGTTCAATTTTTTGTAAAAAAGTACTTGCAAGGTAAGGA‑TAAAAGCCTATAATACGCGGGCACAACGACGCGCCGTTGTAGGGTTTTATAGACAGCGCGTCGtttttt                                                                                                >  1:708186/1‑251 (MQ=25)                                                   atttTACTTAACTCGGAAAGCCAAAAGTGCGGTCAGAAATTCAAGCGTTTTTCTGACCGCATTTTTTTGATCTCCATAATTTTTTATCTCAAAGGAAAACACAATTTGGCGAGTTTATGAGCAAATAATAAAAAAGTTCAATTTTTTGTAAAAAAGTACTTGCAAGGTAAGGA‑TAAAAGCCTATAATACGCGGGCACAACGACGCGCCGTTGTAGGGTTTTATAGAc                                                                                                               >  1:797658/1‑227 (MQ=25)                                                   atttTACTTAACTCGGAAAGCCAAAAGTGCGGTCAGAAATTCAAGCGTTTTTCTGACCGCA‑TTTTTTGATCTCCATAATTTTTTATCTCAAAGGAAAACACAATTTGGCGAGTTTATGAGCAAATAATAAAAAAGTTCAATTTTTTGTAAAAAAGTACTTGCAAGGTAAGGA‑TAAAAGCCTATAATACGCGGGCACAACGACGCGCCGTTGTAGGGTTTTATAGACAGCGCGTCGTTTTTTGTTCTTTaa                                                                                       >  1:68217/1‑250 (MQ=21)                                                                   aaaGCCAAAAGTGCGGTCAGAAATTCAAGCGTTTTTCTGACCGCATTTTTTTGATCTCCATAATTTTTTATCTCAAAGGAAAACACAATTTGGCGAGTTTATGAGCAAATAATAAAAAAGTTCAATTTTTTGTAAAAAAGTACTTGCAAGGTAAGGA‑TAAAAGCCTATAATACGCGGGCACAACGACGCGCCGTTGTAGGGTTTTATAGACAGCGCGTCGTTTTTTGTTCTTTAACAATCAATCATACaa                                                                        <  1:746276/250‑1 (MQ=21)                                                                    aaGCCAAAAGTGCGGTCAGAAATTCAAGCGTTTTTCTGACCGCATTTTTTTGATCTCCATAATTTTTTATCTCAAAGGAAAACACAATTTGGCGAGTTTATGAGCAAATAATAAAAAAGTTCAATTTTTTGTAAAAAAGTACTTGCAAGGTAAGGA‑TAAAAGCCTATAATACGCGGGCACAACGACGCGCCGTTGTAGGGTTTTATAGACAGCGCGTCGTTTTTTGTTCTTTAACAATCAATCATACAAAc                                                                      <  1:896755/251‑1 (MQ=21)                                                                       ccAAAAGTGCGGTCAGAAATTCAAGCGTTTTTCTGACCGCATTTTTTTGATCTCCATAATTTTTTATCTCAAAGGAAAACACAATTTGGCGAGTTTATGAGCAAATAATAAAAAAGTTCAATTTTTTGTAAAAAAGTACTTGCAAGGTAAGGA‑TAAAAGCCTATAATACGCGGGCACAACGACGCGCCGTTGTAGGGTTTTATAGACAGCGCGTCGTTTTTTGTTCTTTAACAATCAATCATACAAACtg                                                                    <  1:314764/250‑1 (MQ=18)                                                                        cAAAAGTGCGGTCAGAAATTCAAGCGTTTTTCTGACCGCATTTTTTTGATCTCCATAATTTTTTATCTCAAAGGAAAACACAATTTGGCGAGTTTATGAGCAAATAATAAAAAAGTTCAATTTTTTGTAAAAAAGTACTTGCAAGGTAAGGA‑TAAAAGCCTATAATACGCGGGCACAACGACGCGCCGTTGTAGGGTTTTATAGACAGCGCGTCGTTTTTTGTTCTTTAACAATCAATCATACAAACtgtg                                                                  <  1:476513/251‑1 (MQ=18)                                                                                    cAGAAATTCAAGCGTTTTTCTGACCGCATTTTTTTGATCTCCATAATTTTTTATCTCAAAGGAAAACACAATTTGGCGAGTTTATGAGCAAATAATAAAAAAGTTCAATTTTTTGTAAAAAAGTACTTGCAAGGTAAGGA‑TAAAAGCCTATAATACGCGGGCACAACGACGCGCCGTTGTAGGGTTTTATAGACAGCGCGTCGTTTTTTGTTCTTTAACAATCAATCATACaa                                                                        >  1:550647/1‑233 (MQ=18)                                                                                       aaaTTCAAGCGTTTTTCTGACCGCATTTTTTTGATCTCCATAATTTTTTATCTCAAAGGAAAACACAATTTGGCGAGTTTATGAGCAAATAATAAAAAAGTTCAATTTTTTGTAAAAAAGTACTTGCAAGGTAAGGA‑TAAAAGCCTATAATACGCGGGCACAACGACGCGCCGTTGTAGGGTTTTATAGACAGCGCGTCGTTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGAt                                                   <  1:43455/251‑1 (MQ=18)                                                                                       aaaTTCAAGCGTTTTTCTGACCGCATTTTTTTGATCTCCATAATTTTTTATCTCAAAGGAAAACACAATTTGGCGAGTTTATGAGCAAATAATAAAAAAGTTCAATTTTTTGTAAAAAAGTACTTGCAAGGTAAGGA‑TAAAAGCCTATAATACGCGGGCACAACGACGCGCCGTTGTAGGGTTTTATAGACAGCGCGTCGTTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGAt                                                   <  1:30888/251‑1 (MQ=18)                                                                                                             gCATTTTTTTGATCTCCATAATTTTTTATCTCAAAGGAAAACACAATTTGGCGAGTTTATGAGCAAATAATAAAAAAGTTCAATTTTTTGTAAAAAAGTACTTGCAAGGTAAGGA‑TAAAAGCCTATAATACGCGGGCACAACGACGCGCCGTTGTAGGGTTTTATAGACAGCGCGTCGTTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAaatg                                       <  1:793014/241‑2 (MQ=14)                                                                                                              cATTTTTTTGATCTCCATAATTTTTTATCTCAAAGGAAAACACAATTTGGCGAGTTTATGAGCAAATAATAAAAAAGTTCAATTTTTTGTAAAAAAGTACTTGCAAGGTAAGGA‑TAAAAGCCTATAATACGCGGGCACAACGACGCGCCGTTGTAGGGTTTTATAGACAGCGCGTCGTTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGAt                                                   >  1:205842/1‑228 (MQ=17)                                                                                                                  tttttGATCTCCATAATTTTTTATCTCAAAGGAAAACACAATTTGGCGAGTTTATGAGCAAATAATAAAAAAGTTCAATTTTTTGTAAAAAAGTACTTGCAAGGTAAGGA‑TAAAAGCCTATAATACGCGGGCACAACGACGCGCCGTTGTAGGGTTTTATAGACAGCGCGTCGTTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGaa                         >  1:803828/1‑250 (MQ=14)                                                                                                                                        aTCTCAAAGGAAAACACAATTTGGCGAGTTTATGAGCAAATAATAAAAAAGTTCAATTTTTTGTAAAAAAGTACTTGCAAGGTAAGGA‑TAAAAGCCTATAATACGCGGGCACAACGACGCGCCGTTGTAGGGTTTTATAGACAGCGCGTCGTTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCg  >  1:306897/1‑251 (MQ=11)                                                                                                                                         tctcAAAGGAAAACACAATTTGGCGAGTTTATGAGCAAATAATAAAAAAGTTCAATTTTTTGTAAAAAAGTACTTGCAAGGTAAGGA‑TAAAAGCCTATAATACGCGGGCACAACGACGCGCCGTTGTAGGGTTTtat                                                                                                                   <  1:716553/137‑1 (MQ=11)                                                                                                                                         tctcAAAGGAAAACACAATTTGGCGAGTTTATGAGCAAATAATAAAAAAGTTCAATTTTTTGTAAAAAAGTACTTGCAAGGTAAGGA‑TAAAAGCCTATAATACGCGGGCACAACGACGCGCCGTTGTAGGGTTTTATAGACAGCGCGTCGTTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATGATTTTAAGAt                             <  1:764852/223‑1 (MQ=11)                                                                                                                                                                                                              |||                                                                                                                                                                                    GTGGTGCGGCACTTGGCGGTTTGATTGGAAGCCAAGTGGGCAAATAATTTATTTTACTTAACTCGGAAAGCCAAAAGTGCGGTCAGAAATTCAAGCGTTTTTCTGACCGCATTTTTTTGATCTCCATAATTTTTTATCTCAAAGGAAAACACAATTTGGCGAGTTTATGAGCAAATAATAAAAAAGTTCAATTTTTTGTAAAAAAACTCTTGCAAAG‑AATTAGTAAATACCTATAATACGCGGGCACAACGACGCGCCGTTGTAGGGTTTTATAGACAGCGCGTCGTTTTTTGTTCTTTAACAATCAATCATACAAACTGTGTGGGCACTTATAGATTCGTTTTAAATAATTTTAAGATTGAACTTAATAGGTGCCGGAAGATTCG  >  1/1034777‑1035162

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