Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	443,473	2 bp→TG	intergenic (+80056/‑128545)	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,473	0	A	T	100.0%	66.1 / NA	24	intergenic (+80056/‑128546)	unknown/unknown	tRNA‑SeC(p)‑TCA/SSU rRNA ## 16S rRNA, small subunit ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (15/9);  total (15/9)
*	1	443,474	0	T	G	100.0%	69.8 / NA	24	intergenic (+80057/‑128545)	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 G (15/9);  total (15/9)

AAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTCCGTTAAATTTGAGGGTTGTTAGTTGTTGGTGATGTGTAAAATTTACGTAGTGTTACTTAATTTAAAACCATATAAATTTAGAAAATTTATTTCTAGTTTTATATGGTTATTTTTCTAAAATCTCGCGGTCGGCTACACGCAAAACTGCTTTTGGTCGGCACGAATTAACTCGCGTAAAGAACCAAAATGTTGCAATACATTTTTCGATAACTGCATCACCGGGCAGTCCTTAATTCCCGTACGCAAAAAGATCGCCAATAGTTCATCGTCCGTCAAAGCGGCAGCACCGTATTTTAATAGCTTTTCCCTTGGCATTAATGCTTCAT  >  1/443337‑443714                                                                                                                                          ||                                                                                                                                                                                                                                                 aaaCGGCGGGTTAATACTCCGCCGTTTTTCTTTTCATGTACATGATGATTCTGCTTTAAATTTGAGGGTAATTAGTCGCTGATGATGTGTAAAATTGACGTAGCGCTATTTGGCGTTAAACCACATAAGTTTAGAAATGTTGTTTCTAGGTTTATGTGGTTATTTTTCTAAAATCTCGCGGTCGGCAACACgtccggtaatgt                                                                                                                                                                                  <  1:82718/203‑12 (MQ=255)                                           tgatgaTTCTGCTTTAAATTTGAGGGTAATTAGTCGCTGATGATGTGTAAAATTGACGTAGCGCTATTTGGCGTTAAACCACATAAGTTTAGAAATGTTGTTTCTAGGTTTATGTGGTTATTTTTCTAAAATCTCGCGGTCGGCAacac                                                                                                                                                                                              >  1:770647/1‑149 (MQ=255)                                                     gctttAAATTTGAGGGTAATTAGTCGCTGATGATGTGTAAAATTGACGTAGCGCTATTTGGCGTTAAACCACATAAGTTTAGAAATGTTGTTTCTAGGTTTATGTGGTTATTTTTCTAAAATCTCGCGGTCGGCAACACgtccggtaa                                                                                                                                                                                     >  1:450714/4‑140 (MQ=255)                                                                                                       gcgctatttggcgttaaaccacataagtttagaAATGTTGTTTCTAGGTTTATGTGGTTATTTTTCTAAAATCTCGCGGTCGGCAACACGtccggtaatgttgatttatgcaaaggttacataaagccactggcaagagttattcctgagcaagaagagacggctagtcagggagagtataccggtgagaatcgtgttgcgccgaacccgatgaaatctcgtagtcaatcttggtttttccgcccgggtta                            >  1:842229‑M1/16‑90 (MQ=255)                                                                                                          ctatttggcgttaaaccacataagtttAGAAATGTTGTTTCTAGGTTTATGTGGTTATTTTTCTAAAATCTCGCGGTCGGCAACACGtccggtaatgttgatttatgcaaaggttacataaagccactggcaagagttattcctgagcaagaagagacggctagtcagggagagtataccggtgagaatcgtgttgcgccgaacccgatgaaatctcgtagtcaatcttggtttttccgcccgggttatca                         >  1:419408‑M1/13‑87 (MQ=255)                                                                                                          ctatttggcgttaaaccacataagtttAGAAATGTTGTTTCTAGGTTTATGTGGTCATTTTTCTAAAATCTCGCGGTCGGCAACACGtccggtaatgttgatttatgcaaaggttacataaagccactggcaagagttattcctgagcaagaagagacggctagtcagggagagtataccggtgagaatcgtgttgcgccgaacccgatgaaatctcgtagtcaatcttggtttttccgcccgggttatca                         >  1:174329‑M1/13‑87 (MQ=255)                                                                                                            atttggcgttaaaccacataagtTTAGAAATGTTGTTTCTAGGTTTATGTGGTTATTTTTCTAAAATCTCGCGGTCGGCAACACGtccggtaatgttgatttatgcaaaggttacataaagccactggcaagagttattcctgagcaagaagagacggctagtcagggagagtataccggtgagaatcgtgttgcgccgaacccgatgaaatctcgtagtcaatcttggtttttccgcccgggttatca                         >  1:111497‑M1/11‑85 (MQ=255)                                                                                                            atttggcgttaaaccacataagtTTAGAAATGTTGTTTCTAGGTTTATGTGGTTATTTTTCTAAAATCTCGCGGTCGGCAACACGtccggtaatgttgatttatgcaaaggttacataaagccactggcaagagttattcctgagcaagaagagacggctagtcagggagagtataccggtgagaatcgtgttgcgccgaacccgatgaaatctcgtagtcaatcttggtttttccgcccgggttatca                         >  1:235911‑M1/11‑85 (MQ=255)                                                                                                            atttggcgttaaaccacataagtTTAGAAATGTTGTTTCTAGGTTTATGTGGTTATTTTTCTAAAATCTCGCGGTCGGCAACACGtccggtaatgttgatttatgcaaaggttacataaagccactggcaagagttattcctgagcaagaagagacggctagtcagggagagtataccggtgagaatcgtgttgcgc                                                                             >  1:783122‑M1/11‑85 (MQ=255)                                                                                                               tggcgttaaaccacataAGTTTAGAAATGTTGTTTCTAGGTTTATGTGGTTATTTTTCTAAAATCTCGCGGTCGGCAACACGtccggtaatgttgatttatgcaaaggttacataaagccactggcaagagttattctgagcaagaagagacggctagtcagggagagtataccggtgagaatcgtgttgcgccgaacccgatgaaatctcgtagtcaatcttggtttttccgcccgggttatcacttttc                    <  1:599060‑M1/244‑170 (MQ=255)                                                                                                                 gcgttaaaccacaTAAGTTTAGAAATGTTGTTTCTAGGTTTATGTGGTTATTTTTCTAAAATCTCGCGGTCGGCAACACGtccggtaatgttgatttatgcaaaggttacataaagccactggcaagagttattcctgagcaagaagagacggctagtcagggagagtataccggtgagaatcgtgttgcgccgaacccgatgaaatctcgtagtcaatcttggtttttccgcccgggttatcacttttc                   >  1:88256‑M1/6‑80 (MQ=255)                                                                                                                 gcgttaaaccacaTAAGTTTAGAAATGTTGTTTCTAGGTTTATGTGGTTATTTTTCTAAAATCTCGCGGTCGGCAACACGtccggtaatgttgatttatgcaaaggttacataaagccactggcaagagttattcctgagcaagaagagacggctagtcagggagagtataccggtgagaatcgtgttgcgccgaacccgatgaaatctcgtagtcaatcttggtttttccgcccgggttatcacttttc                   >  1:880404‑M1/6‑80 (MQ=255)                                                                                                                 gcgttaaaccacaTAAGTTTAGAAATGTTGTTTCTAGGTTTATGTGGTTATTTTTCTAAAATCTCGCGGTCGGCAACACGtccggtaatgttgatttatgcaaaggttacataaagccactggcaagagttattcctgagcaagaagagacggctagtcagggagagtataccggtgagaatcgtgttgcgccgaacccgatgaaatctcgtagtcaatcttggtttttccgcccgggttatcacttttc                   >  1:268599‑M1/6‑80 (MQ=255)                                                                                                                 gcgttaaaccacaTAAGTTTAGAAATGTTGTTTCTAGGTTTATGTGGTTATTTTTCTAAAATCTCGCGGTCGGCAACACGtccggtaatgttgatttatgcaaaggttacataaagccactggcaagagttattcctgagcaagaagagacggctagtcagggagagtataccggtgagaatcgtgttgcgccgaacccgatgaaatctcgtagtcaatcttggtttttccgcccgggttatcacttttc                   >  1:57300‑M1/6‑80 (MQ=255)                                                                                                                         ccacataagttTAGAAATGTTGTTTCTAGGTTTATGTGGTTATTTTTCTAAAATCTCGCGGTCGGCAACACGtccggtaatgttgatttatgcaaaggttacataaagccactggcaagagttattcctgagcaagaagagacggctagtcagggagagtataccggtgggaatcgtgttgcgccgaacccgatgaaatctcgtagtcaatcttggtttttccgcccgggttatcacttttccgagcaaca          >  1:450349‑M1/5‑72 (MQ=255)                                                                                                                            cataaGTTTAGAAATGTTGTTTCTAGGTTTATGTGGTTATTTTTCTAAAATCTCGCGGTCGGCAACACGtccggtaatgttgatttatgcaaaggttacataaagccactggcaagagttattcctgagcaagaagagacggctagtcagggagagtataccggtgagaatcgtgttgcgccgaacccgatgaaatctcgtagtcaatcttggtttttccgcccgggttatcacttttccgagcaacatta       <  1:346124‑M1/250‑183 (MQ=255)                                                                                                                               aagtttagaAATGTTGTTTCTAGGTTTATGTGGTTATTTTTCTAAAATCTCGCGGTCGGCAACACGtccggtaatgttgatttatgcaaaggttacataaagccactggcaagagttattcctgagcaagaagagacggctagtcagggagagtataccggtgagaatcgtgttgcgccgaacccgatgaaatctcgtagtcaatcttggtttttccgcccgggttatcacttttccgagcaacattatct    <  1:643011‑M1/248‑186 (MQ=255)                                                                                                                               aagtttagaAATGTTGTTTCTAGGTTTATGTGGTTATTTTTCTAAAATCTCGCGGTCGGCAACACGtccggtaatgttgatttatgcaaaggttacataaagccactggcaagagttattcctgagcaagaagagacggctagtcagggagagtataccggtgagaatcgtgttgcgccgaacccgatgaaatctcgtagtcaatcttggtttttccgcccgggttatcacttttccgagcaacattatct    <  1:849248‑M1/248‑186 (MQ=255)                                                                                                                               aagtttagaAATGTTGTTTCTAGGTTTATGTGGTTATTTTTCTAAAATCTCGCGGTCGGCAACACGtccggtaatgttgatttatgcaaaggttacataaagccactggcaagagttattcctgagcaagaagagacggctagtcagggagagtataccggtgagaatcgtgttgcgccgaacccgatgaaatctcgtagtcaatcttggtttttccgcccgggttatcacttttccgagcaacattatct    <  1:743567‑M1/248‑186 (MQ=255)                                                                                                                               aagtttagaAATGTTGTTTCTAGGTTTATGTGGTTATTTTTCTAAAATCTCGCGGTCGGCAACACGtccggtaatgttgatttatgcaaaggttacataaagccactggcaagagttattcctgagcaagaagagacggctagtcagggagagtataccggtgagaatcgtgttgcgccgaacccgatgaaatctcgtagtcaatcttggtttttccgcccgggttatcacttttccgagcaacattatct    <  1:410578‑M1/248‑186 (MQ=255)                                                                                                                                 gtttaGAAATGTTGTTTCTAGGTTTATGTGGTTATTTTTCTAAAATCTCGCGGTCGGCAACACGtccggtaatgttgatttatgcaaaggttacataaagccactggcaagagttattcctgagcaagaagagacggctagtcagggagagtataccggtgagaatcgtgttgcgccgaacccgatgaaatctcgtagtcaatcttggtttttccgcccgggttatcacttttccgagcaacattatctcg  >  1:189643‑M1/2‑64 (MQ=255)                                                                                                                                 gtttaGAAATGTTGTTTCTAGGTTTATGTGGTTATTTTTCTAAAATCTCGCGGTCGGCAACACGtccggtaatgttgatttatgcaaaggttacataaagccactggcaagagttattcctgagcaagaagagacggctagtcagggagagtataccggtgagaat                                                                                       >  1:401336‑M1/2‑64 (MQ=255)                                                                                                                                  tttAGAAATGTTGTTTCTAGGTTTATGTGGTTATTTTTCTAAAATCTCGCGGTCGGCAACACGtccggtaatgttgatttatgcaaaggttacataaagccactggcaagagttattcctgagcaagaagagacggctagtcagggagagtataccggtgagaatcgtgttgcgccgaacccgatgaaatctcgtagtcaatcttggtttttccgcccgggttatcacttttccgagcaacattatctc   <  1:42191‑M1/249‑187 (MQ=255)                                                                                                                                  tttAGAAATGTTGTTTCTAGGTTTATGTGGTTATTTTTCTAAAATCTCGCGGTCGGCAACACGtccggtaatgttgatttatgcaaaggttacataaagccactggcaagagttattcctgagcaagaagagacggctagtcagggagagtataccggtgagaatcgtgttgcgccgaacccgatgaaatctcgtagtcaatcttggtttttccgcccgggttatcacttttccgagcaacattatctc   <  1:274764‑M1/249‑187 (MQ=255)                                                                                                                                          ||                                                                                                                                                                                                                                                 AAACGGCGGGGTAATAC‑CCGTCGTTTTTGTTTATGTTTATATGATGATTCTCCGTTAAATTTGAGGGTTGTTAGTTGTTGGTGATGTGTAAAATTTACGTAGTGTTACTTAATTTAAAACCATATAAATTTAGAAAATTTATTTCTAGTTTTATATGGTTATTTTTCTAAAATCTCGCGGTCGGCTACACGCAAAACTGCTTTTGGTCGGCACGAATTAACTCGCGTAAAGAACCAAAATGTTGCAATACATTTTTCGATAACTGCATCACCGGGCAGTCCTTAATTCCCGTACGCAAAAAGATCGCCAATAGTTCATCGTCCGTCAAAGCGGCAGCACCGTATTTTAATAGCTTTTCCCTTGGCATTAATGCTTCAT  >  1/443337‑443714

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