Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	438,112	A→G	S370G (AGC→GGC)	thiI →	tRNA 4‑thiouridine(8) synthase ThiI

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	438,112	0	A	G	100.0%	91.1 / NA	25	S370G (AGC→GGC)	thiI	tRNA 4‑thiouridine(8) synthase ThiI
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (14/11);  total (14/11)

GGCCCGCCAGATTGGCACCGAAGACTTTGCTCGCACGATGCCGGAATATTGTGGTGTGATCTCCAAAAGCCCGACGGTGAAAGCAGTTAAATCGAAGATTGAAGCGGAAGAAGAGAAGTTCGACTTCAGCATTCTCGATAAAGTGGTTGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCCAGCAGACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGGCTTCGGCCCGAACGACGTGATCCTCGATAT  >  NZ_CP009273/437985‑438242                                                                                                                                |                                                                                                                                   ggCCCGCCAGATTGGCACCGAAGACTTTGCTCGCACGATGCCGGAATATTGTGGTGTGATCTCCAAAAGCCCGACGGTGAAAGCAGTTAAATCGAAGATTGAAGCGGAAGAAGAGAAGTTCGACTTCGGCATTCTCGa                                                                                                                          >  2:346471/1‑138 (MQ=255)  gCCCGCCAGATTGGCACCGAAGACTTTGCTCGCACGATGCCGGAATATTGTGGTGTGATCTCCAAAAGCCCGACGGTGAAAGCAGTTAAATCGAAGATTGAAGCGGAAGAAGAGAAGTTCGACTTCGGCATTCTCGAt                                                                                                                         >  2:530896/1‑138 (MQ=255)                cACCGAAGACTTTGCTCGCACGATGCCGGAATATTGTGGTGTGATCTCCAAAAGCCCGACGGTGAAAGCAGTTAAATCGAAGATTGAAGCGGAAGAAGAGAAGTTCGACTTCGGCATTCTCGATAAAGTGGTTGAGGaa                                                                                                          <  1:590388/139‑1 (MQ=255)                     aaGACTTTGCTCGCACGATGCCGGAATATTGTGGTGTGATCTCCAAAAGCCCGACGGTGAAAGCAGTTAAATCGAAGATTGAAGCGGAAGAAGAGAAGTTCGACTTCGGCATTCTCGATAAAGTGGTTGAGGAAGCGaa                                                                                                     <  2:263241/139‑1 (MQ=255)                                        gCCGGAATATTGTGGTGTGATCTCCAAAAGCCCGACGGTGAAAGCAGTTAAATCGAAGATTGAAGCGGAAGAAGAGAAGTTCGACTTCGGCATTCTCGATAAAGTGGTTGAGGAAGCGAATAACGTTGATATCCGCGa                                                                                   >  1:727435/1‑138 (MQ=255)                                                 ttGTGGTGTGATCTCCAAAAGCCCGACGGTGAAAGCAGTTAAATCGAAGATTGAAGCGGAAGAAGAGAAGTTCGACTTCGGCATTCTCGATAAAGTGGTTGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCcag                                                                         >  1:714118/1‑139 (MQ=255)                                                 ttGTGGTGTGATCTCCAAAAGCCCGACGGTGAAAGCAGTTAAATCGAAGATTGAAGCGGAAGAAGAGAAGTTCGACTTCGGCATTCTCGATAAAGTGGTTGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCcag                                                                         >  1:563906/1‑139 (MQ=255)                                                       tgtgATCTCCAAAAGCCCGACGGTGAAAGCAGTTAAATCGAAGATTGAAGCGGAAGAAGAGAAGTTCGACTTCGGCATTCTCGATAAAGTGGTTGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCCAGCAGAcc                                                                   >  1:259805/1‑139 (MQ=255)                                                         tgATCTCCAAAAGCCCGACGGTGAAAGCAGTTAAATCGAAGATTGAAGCGGAAGAAGAGAAGTTCGACTTCGGCATTCTCGATAAAGTGGTTGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCCAGCAGAcc                                                                   >  2:22633/1‑137 (MQ=255)                                                           aTCTCCAAAAGCCCGACGGTGAAAGCAGTTAAATCGAAGATTGAAGCGGAAGAAGAGAAGTTCGACTTCGGCATTCTCGATAAAGTGGTTGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCCAGCAGACCGAg                                                                >  2:328462/1‑138 (MQ=255)                                                           aTCTCCAAAAGCCCGACGGTGAAAGCAGTTAAATCGAAGATTGAAGCGGAAGAAGAGAAGTTCGACTTCGGCATTCTCGATAAAGTGGTTGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCCAGCAGACCGAg                                                                >  1:541463/1‑138 (MQ=255)                                                           aTCTCCAAAAGCCCGACGGTGAAAGCAGTTAAATCGAAGATTGAAGCGGAAGAAGAGAAGTTCGACTTCGGCATTCTCGATAAAGTGGTTGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCCAGCAGACCGAg                                                                <  1:126446/138‑1 (MQ=255)                                                                cAAAAGCCCGACGGTGAAAGCAGTTAAATCGAAGATTGAAGCGGAAGAAGAGAAGTTCGACTTCGGCATTCTCGATAAAGTGGTTGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCCAGCAGACCGAGCAGGa                                                           <  1:488004/138‑1 (MQ=255)                                                                cAAAAGCCCGACGGTGAAAGCAGTTAAATCGAAGATTGAAGCGGAAGAAGAGAAGTTCGACTTCGGCATTCTCGATAAAGTGGTTGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCCAGCAGACCGAGCAGGa                                                           <  2:27054/138‑1 (MQ=255)                                                                            ggTGAAAGCAGTTAAATCGAAGATTGAAGCGGAAGAAGAGAAGTTCGACTTCGGCATTCTCGATAAAGTGGTTGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCCAGCAGACCGAGCAGGAAGTGgtggaagtg                                              <  1:530896/139‑1 (MQ=255)                                                                                       ttAAATCGAAGATTGAAGCGGAAGAAGAGAAGTTCGACTTCGGCATTCTCGATAAAGTGGTTGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCCAGCAGACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCaa                                   >  1:73461/1‑139 (MQ=255)                                                                                               aaGATTGAAGCGGAAGAAGAGAAGTTCGACTTCGGCATTCTCGATAAAGAGGTTGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCCAGCAGACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGGCTTCg                           >  1:410003/1‑139 (MQ=255)                                                                                                  aTTGAAGCGGAAGAAGAGAAGTTCGACTTCGGCATTCTCGATAAAGTGGTTGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCCAGCAGACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGGCTTCGGcc                        <  2:563906/139‑1 (MQ=255)                                                                                                     gAAGCGGAAGAAGAGAAGTTCGACTTCGGCATTCTCGATAAAGTGGTTGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCCAGCAGACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGGCTTCGGCCCg                      <  1:150522/138‑1 (MQ=255)                                                                                                            aagaagAGAAGTTCGACTTCGGCATTCTCGATAAAGTGGTTGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCCAGCAGACCGAGCAGGAAgtggtg                                                    >  1:647833/1‑101 (MQ=255)                                                                                                            aagaagAGAAGTTCGACTTCGGCATTCTCGATAAAGTGGTTGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCCAGCAGACCGAGCAGGAAgtggtg                                                    <  2:647833/101‑1 (MQ=255)                                                                                                                      gTTCGACTTCGGCATTCTCGATAAAGTGGTTGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCCAGCAGACCGAGCAGGa                                                           >  1:762665/1‑84 (MQ=255)                                                                                                                      gTTCGACTTCGGCATTCTCGATAAAGTGGTTGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCCAGCAGACCGAGCAGGa                                                           <  2:762665/84‑1 (MQ=255)                                                                                                                      gTTCGACTTCGGCATTCTCGATAAAGTGGTTGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCCAGCAGACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGGCTTCGGCCCGAACGACGTGATCCTCGa     >  1:198814/1‑138 (MQ=255)                                                                                                                         cGACTTCGGCATTCTCGATAAAGTGGTTGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCCAGCAGACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGGCTTCGGCCCGAACGACGTGATCCTCGatat  <  2:407829/138‑1 (MQ=255)                                                                                                                                |                                                                                                                                   GGCCCGCCAGATTGGCACCGAAGACTTTGCTCGCACGATGCCGGAATATTGTGGTGTGATCTCCAAAAGCCCGACGGTGAAAGCAGTTAAATCGAAGATTGAAGCGGAAGAAGAGAAGTTCGACTTCAGCATTCTCGATAAAGTGGTTGAGGAAGCGAATAACGTTGATATCCGCGAAATCGCCCAGCAGACCGAGCAGGAAGTGGTGGAAGTGGAAACCGTCAATGGCTTCGGCCCGAACGACGTGATCCTCGATAT  >  NZ_CP009273/437985‑438242

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