Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,193,141	G→A	G211G (GGC→GGT)	mnmA ←	tRNA‑specific 2‑thiouridylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,193,141	0	G	A	100.0%	70.2 / NA	22	G211G (GGC→GGT)	mnmA	tRNA‑specific 2‑thiouridylase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (9/13);  total (9/13)

CAACCCGACAGACATCAGCCGCGGGTGTTCATGGCCCTGAGCGACAACCAGAATGTTGTTTTCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGGCCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGC  >  NC_000913/1192913‑1193351                                                                                                                                                                                                                                     |                                                                                                                                                                                                                   cAACCCGGCAGACATCAGCCGCGGGTGGTCATGGCCCTGAGCGACAACCAGAATGTTGTTTTCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTa                                                                                                                                                                                               <  2:88524/250‑1 (MQ=255)             caTCAGCCGCGGGTGTTCATGGCCCTGAGCGACAACCAGAATGTTGTTTTCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGa                                                                                                                                                                                  <  1:88519/251‑1 (MQ=255)                     gcgGGTGTTCATGGCCCTGAGCGACAACCAGAATGTTGTTTTCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATg                                                                                                                                                                          <  1:197511/251‑1 (MQ=255)                           gtTCATGGCCCTGAGCGACAACCATGATGTTGTTTTCGACGTCTTTGTCCTCCACATACCACGGTTCTTTGGTACCTTCTTTGGGGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGGGATACATCAGCCCCTGGTGCTAGCCAATTTCATCGCCAGCGGCGGTAATGATTTGGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTGACGGTGTCCGATGATGCAAATGCCGGTa                                                                                                                                                                    <  2:120321/251‑1 (MQ=255)                                            aCAACCAGAATGTTGTTTTCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTcgc                                                                                                                                                   <  2:132305/251‑1 (MQ=255)                                                 cAGAATGTTGTTTTCGATGTCTTTGTCCCCCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGATCCACACACCTTTACACTGACCGCGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAAAGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTg                                                                                                                                               >  2:126329/1‑250 (MQ=255)                                                   gAATGTTGTTTTCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTgacc                                                                                                                                            >  2:1887/1‑251 (MQ=255)                                                                           ccaccaCATACCACTGTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTa                                                                                                                    <  1:234528/251‑1 (MQ=255)                                                                                                                                ccTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTg                                                               >  1:31396/1‑251 (MQ=255)                                                                                                                                                   tgcgATACCTCAGCCCATGGTGCTCACCAATTTCCTCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGCGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGCCCGGGAACAGGCTTTGCGCACTCTGCTCATGGCTGAGCGtata                                             >  2:54593/4‑250 (MQ=255)                                                                                                                                                      gATACAGCAGCCCCTGGTGCTCACCAATTTCAGCGCCATCGTCGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTGCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTGTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAgg                                         <  2:133236/251‑1 (MQ=255)                                                                                                                                                           aTCAGCCCCTGGTGCTTCCCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTa                                    <  2:140291/251‑1 (MQ=255)                                                                                                                                                             cAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAg                                   <  2:237948/250‑1 (MQ=255)                                                                                                                                                                 cccTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTg                              >  1:114582/1‑251 (MQ=255)                                                                                                                                                                 cccTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGAACGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTg                              <  2:61824/251‑1 (MQ=255)                                                                                                                                                                 cccTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATGTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACGCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTGTTTCCAGTTCGCCGACCGGGAACAGGCTGTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTg                              <  2:131953/251‑1 (MQ=255)                                                                                                                                                                       tGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATaaa                                           >  2:29409/1‑232 (MQ=255)                                                                                                                                                                       tGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATaaa                                           <  1:29409/232‑1 (MQ=255)                                                                                                                                                                                   tCATCGCCATCGACGGAAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGAAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGATCTGGACTTTATTGCTGTcc              >  2:206228/1‑249 (MQ=255)                                                                                                                                                                                              gACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACgcagc  >  2:86760/1‑250 (MQ=255)                                                                                                                                                                                                            tGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGctct                               <  2:210275/207‑1 (MQ=255)                                                                                                                                                                                                            tGCCCGGTTGCGCCGGGAGATAACGACCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGctct                               >  1:210275/1‑207 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                   CAACCCGACAGACATCAGCCGCGGGTGTTCATGGCCCTGAGCGACAACCAGAATGTTGTTTTCGACGTCTTTGTCCACCACATACCACGGTTCTTCGGTACCTTCTTTGGTGCCACCGATACCCAGACCTTTACGCTGACCGAGAGTGTGATACATCAGCCCCTGGTGCTCACCAATTTCATCGCCATCGACGGTAATGATTTTGCCCGGTTGCGCCGGGAGATAACGGCCCAGGAACTCGCGGAATTTACGTTCGCCAATGAAGCAAATGCCGGTAGAGTCTTTTTTCTTCGCGGTGACCAGACCAAGATCTTCAGCAATCTTACGCACCTGCGGTTTTTCCAGTTCGCCGACCGGGAACAGGCTTTGCGCAATCTGCTCATGGCTGAGCGTATAAAGGAAGTAGCTCTGGTCTTTATTGCTGTCCAGACCACGCAGC  >  NC_000913/1192913‑1193351

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