Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	441,907	G→A	S120N (AGC→AAC)	thiI →	tRNA uridine 4‑sulfurtransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	441,907	0	G	A	100.0%	39.6 / NA	14	S120N (AGC→AAC)	thiI	tRNA uridine 4‑sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (4/10);  total (4/10)

CGCAAAAGATGAAAACCAGCGTCTGGCTATTCGCGACGCTCTGACCCGTATTCCGGGTATCCACCATATTCTCGAAGTCGAAGACGTGCCGTTTACCGACATGCACGATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCCTGCTGATTAAAGGCCGCTACGAAGGTATTGGCGGTTTCCCGATCGGCACCCAGGAAGATGTGCTGTCGCTCATTTCCGGTGGTTTCGACTCCGGTGTTTCCAGTTATATGTTGATGCGTCGCG  >  NC_000913/441707‑442149                                                                                                                                                                                                         |                                                                                                                                                                                                                                                   cgcAAAAGATGAAAACCAGCGTCTGGCTATTCGCGACGCTCTGTCCCGTATTCCGGGTATCCACCATATTCTCGAAGTCGAAGACGTGCCGTTTACCGACATGCACGATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAACTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATAttn                                                                                                                                                                                                  <  1:300370/251‑2 (MQ=255)                                           gACCCGTATTCCGGGTATCCACCATATTCTCGAAGTCGAAGACGTGCCGTTTACCGACATGCACGATATTTGCGAGAAAGCGGTGGTTCAGTAGCGCGATCAGCTGGAAGTCAAAACCTTCTGCGGTCGCGTGAAGCGGCGTGGCAAACATGATTTTAACTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGc                                                                                                                                                                                                        <  1:504038/203‑1 (MQ=255)                                              ccGTATTCCGGGTATCCACCATATTCTCGAAGTCGAAGACGTGCCGTTTACCGACATGCACGATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAACGCGATTGATGTGGAACGTTACGTcggcggcg                                                                                                                                                                                                                    <  1:140184/188‑1 (MQ=255)                                                                                          cGTTTACCGACATGCACGATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAACTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTCAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGCAGACGATCGTCTCCTGCTGATTAAAGGCCGCTAc                                                                                                         >  1:277058/1‑251 (MQ=255)                                                                                           gTTTACCGACATGCACAATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGAGAAGCGCCGTGACAAACATGATTTTAACTCGATTGATGTGGACCGTTACGTCGGCGGCGGTTTAAATCATCATATTGAATc                                                                                                                                                                                              <  1:495457/165‑1 (MQ=255)                                                                                                      tGCACGATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCGTCTGCGTACGCGTGAAGCGCCGTGTCAAACATGATTTGAACTCGATTGAGGTTGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGaa                                                                                                                                                                                  <  1:252751/166‑1 (MQ=255)                                                                                                                            gTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAACTCGATTGAtgtgn                                                                                                                                                                                                                                      <  1:301047/92‑2 (MQ=255)                                                                                                                            gTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAACTCGATTGATGTgg                                                                                                                                                                                                                                      >  2:301047/1‑92 (MQ=255)                                                                                                                                                 gcgggAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAACTAGATTGATGTGGAACGTGACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGTTGAGCAAACCGGATGTGACTGTCCAACTGTAAGTGGAAGACGATCGTCTCCTGCTGATTAAAGGCCGCTACGAAGGTATTGGc                                                                                             <  1:287837/205‑1 (MQ=255)                                                                                                                                                                      cGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAACTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGCAGCTGACCAATCCGGATGTGACTGTCCATCtgng                                                                                                                                                 >  2:317894/1‑133 (MQ=255)                                                                                                                                                                       gTACGCGTGAAGCGCCGTGGCAAACATGATTTTANCTCANTNNATNTNGANCGTTa                                                                                                                                                                                                                               <  2:226549/56‑1 (MQ=255)                                                                                                                                                                       gTACGCGTGAAGCGCCGTGGCAAACATGATTTTAACTCGATTGATGTGGAACGTTa                                                                                                                                                                                                                               >  1:226549/1‑56 (MQ=255)                                                                                                                                                                                  gcgcCGTGGCAAACATGATTTTAACTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCCTGCTGATTAAAGGCCGCTACGAAGGTATTGGCGGTTTCCCGATCGGCACCCAGGAAGATGTGCTGTCGCTCATTTCCGGTGGTTTCGACTCCgg                               <  1:262699/237‑1 (MQ=255)                                                                                                                                                                                           cAAACATGATTTTAACTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCAGATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCCTGCTGATTAAAGGCCGCTACGAAGGTATTGGc                                                                                             <  1:488597/166‑1 (MQ=255)                                                                                                                                                                                                  gATTTTAACTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCCTGCTGATTAAAGGCCGCTACGAAGGTATTGGCGGTTTCCCGATCGGCACCCAGGAAGATGTGCTGTCGCTCATTTCCGGTGGTTTCGACTCCGGTGTTTCCAGTTATATGTTGATGCGTcgcn  <  1:317894/250‑2 (MQ=255)                                                                                                                                                                                                         |                                                                                                                                                                                                                                                   CGCAAAAGATGAAAACCAGCGTCTGGCTATTCGCGACGCTCTGACCCGTATTCCGGGTATCCACCATATTCTCGAAGTCGAAGACGTGCCGTTTACCGACATGCACGATATTTTCGAGAAAGCGTTGGTTCAGTATCGCGATCAGCTGGAAGGCAAAACCTTCTGCGTACGCGTGAAGCGCCGTGGCAAACATGATTTTAGCTCGATTGATGTGGAACGTTACGTCGGCGGCGGTTTAAATCAGCATATTGAATCCGCGCGCGTGAAGCTGACCAATCCGGATGTGACTGTCCATCTGGAAGTGGAAGACGATCGTCTCCTGCTGATTAAAGGCCGCTACGAAGGTATTGGCGGTTTCCCGATCGGCACCCAGGAAGATGTGCTGTCGCTCATTTCCGGTGGTTTCGACTCCGGTGTTTCCAGTTATATGTTGATGCGTCGCG  >  NC_000913/441707‑442149

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