Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,490,141	C→T	*333* (TAG→TAA)	idnR ←	transcriptional repressor, 5‑gluconate‑binding

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,490,141	0	C	T	100.0%	110.7 / NA	35	*333* (TAG→TAA)	idnR	transcriptional repressor, 5‑gluconate‑binding
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (18/17);  total (18/17)

AGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGTTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATCAAAACGGGGAGTAATTAC  >  NC_000913/4490020‑4490269                                                                                                                          |                                                                                                                                 gcagtaaaaacagttctgtgtccggcgtgcgggcaattaacaggggttcactcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTg                                                                                                                 <  2:801717‑M2/32‑1 (MQ=255)     taaaaacagttctgtgtccggcgtgcgggcaattaacaggggttcactcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAAGAACTTAAAGCGTGTTGCCGTGATaa                                                                                                             <  2:318636‑M2/36‑1 (MQ=255)           cagttctgtgtccggcgtgcgggcaattaacaggggttcactcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGa                                                                                                       <  2:518356‑M2/42‑1 (MQ=255)                        ggcgtgcgggcaattaacaggggttcactcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCaa                                                                                          <  1:145728‑M2/55‑1 (MQ=255)                        ggcgtgcgggcaattaacaggggttcactcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCaa                                                                                          <  2:697066‑M2/55‑1 (MQ=255)                        ggcgtgcgggcaattaacaggggttcactcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCaa                                                                                          <  1:532446‑M2/55‑1 (MQ=255)                        ggcgtgcgggcaattaacaggggttcactcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCaa                                                                                          <  2:313352‑M2/55‑1 (MQ=255)                                                   ctcatgagatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTAtttt                                                               >  2:383734‑M2/58‑139 (MQ=255)                                                          gatatcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCC‑‑GATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGc                                                      <  1:35158‑M2/89‑1 (MQ=255)                                                              tcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGc                                                                                                                     <  1:103906‑M2/28‑1 (MQ=255)                                                              tcctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGc                                                                                                                     >  2:103906‑M2/47‑74 (MQ=255)                                                                ctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAg                                                  >  2:569124‑M2/45‑139 (MQ=255)                                                                ctccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAg                                                  >  1:505835‑M2/45‑139 (MQ=255)                                                                  ccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGc                                                                                                                     <  2:138542‑M2/28‑1 (MQ=255)                                                                  ccagtttcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGc                                                                                                                     >  1:138542‑M2/43‑70 (MQ=255)                                                                        tcctgcctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACAttt                                          >  2:252386‑M2/37‑139 (MQ=255)                                                                             cctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGggcg                                     >  1:522254‑M2/32‑139 (MQ=255)                                                                             cctggagtataggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGggcg                                     <  2:749253‑M2/108‑1 (MQ=255)                                                                                       aggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAg                                                                                       >  2:236471‑M2/22‑79 (MQ=255)                                                                                       aggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAg                                                                                       <  1:236471‑M2/58‑1 (MQ=255)                                                                                        ggcaacgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGAcc                          <  1:623060‑M2/119‑1 (MQ=255)                                                                                            acgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAata                      >  1:761469‑M2/17‑139 (MQ=255)                                                                                             cgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAatat                     >  1:413426‑M2/16‑139 (MQ=255)                                                                                             cgtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAatat                     >  2:278182‑M2/16‑139 (MQ=255)                                                                                              gtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCAt                                <  2:378362‑M2/113‑1 (MQ=255)                                                                                              gtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCAt                                >  1:378362‑M2/15‑127 (MQ=255)                                                                                              gtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATc                    >  1:79146‑M2/15‑139 (MQ=255)                                                                                              gtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATc                    >  2:293324‑M2/15‑139 (MQ=255)                                                                                              gtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATc                    >  2:72380‑M2/15‑139 (MQ=255)                                                                                              gtccggcaaagagaGAAGAGAAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATc                    >  1:86591‑M2/15‑139 (MQ=255)                                                                                                                agaAAATAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATCAAAACGGGGAGTAATTAc  <  2:761469/139‑1 (MQ=255)                                                                                                                  aaaaTAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTAtttt                                                               >  2:235801/1‑76 (MQ=255)                                                                                                                  aaaaTAACTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTAtttt                                                               <  1:235801/76‑1 (MQ=255)                                                                                                                       aaCTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATCaaaa                >  1:591062/1‑118 (MQ=255)                                                                                                                       aaCTTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATCaaaa                <  2:591062/118‑1 (MQ=255)                                                                                                                          |                                                                                                                                 AGTTCTGTGTCCGGCGTGCGGGCAATTAACAGGGGTTCACTCATGAGATATCCTCCAGTTTCCTGCCTGGAGTATAGGCAACGTCCGGCAAAGAGAATACGTGTAATGAAGAGAAAATAACCTAAAGCGTGTTGCCGTGATAAATCTGATACCCGAGGTCAACAGTGTTGTGGTTGTGATCGTTATTTTTAATTTTGCTCAGCAACATTTGGGCGGCCATCCGACCAATATCAAAACGGGGAGTAATTAC  >  NC_000913/4490020‑4490269

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