Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,494,499	C→T	intergenic (‑93/‑124)	idnD ← / → idnK	L‑idonate 5‑dehydrogenase/D‑gluconate kinase, thermosensitive

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,494,499	0	C	T	100.0%	85.7 / NA	26	intergenic (‑93/‑124)	idnD/idnK	L‑idonate 5‑dehydrogenase/D‑gluconate kinase, thermosensitive
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (12/14);  total (12/14)

TTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATTTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGTAGACACCTCGGAAGCATTTAAGCGGTTTTAACTGTCATTATTTGTGATGAAGATCACGTCAGAAAATTGTTACATTACTATGTTACGCATAACGTGATGTGCCTTGTAATTCTTATCAGTAGAAATAAAAAAACGTGAAATTATTATGCCGCCAGGCGTAGTATCGCAGCAGGTAAGATGATTCAGGAGATTTTAAATGGCGGGTGAAAGCTTTATTTTGATGGGCGTTTCAGGGAGTGGTAAAACATTAATTGGTAGCAAGGTTGCCGCGTTATTATCTGCTAAATTTATTGATGGTGACGATCTTCATCCAGC  >  NC_000913/4494255‑4494741                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   tttCCTTCCGGATAGTAATGTAAAGCGGAACCGCAAAATCCACCGCGTGTTATTTGTACTTATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCGCCCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCATGTAGACTACTCAGAAAGATTTAAGCGATGTTAACTGTCATGAGTTGGGATGAAGAGCACGTCAGAAAACTGTTATATTAc                                                                                                                                                                                                                                               <  2:148794/250‑1 (MQ=255)               ataatGTAAATCGGAACCGCCAATTCCACCTCGGGTCATCTGTACTAATGTTCCATTATTATTCCACTCTATCGTATGCTCGGTAACGCCAACAGTTTTCTTGCCTGCAACAACGCAGGTCTGTGTTTTCCCTTGCATAAATTCCTCACTGGTCAGGTAGACACCTCGGAAGCATTTAAGCGGTTTTAACTGTCATTATTTGTGATGAAGATCACGTCAGAAAATTGTTATATTACTATGTTACGCATAAc                                                                                                                                                                                                                                >  2:108035/1‑251 (MQ=255)                      aaaTCGGAACCGCAAATTCCACCTCGGGTTATTTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACACAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGTAGACACCTCGGAAGCATTTAAGCGGTTTTAACTGTCATTATTTGTGATGAAGATCACGTCAGAAAATTGTTATATTACTATGTTACGCATAACGTGAtg                                                                                                                                                                                                                          <  1:68230/250‑1 (MQ=255)                                               gggTTATTTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGTAGACACCTCGGAAGCATTTAAGCGGTTTTAACTGTCATTATTTGTGATGAAGATCACGTCAGAAAATTGTTATATTACTATGTTACGCATAACGTGATGTGCCTTGTAATTCTTATCAGTAGaa                                                                                                                                                                                                 >  2:12786/1‑250 (MQ=255)                                                      ttGTACTAAGGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATGCCTCACTGGTCAGGTAGACACCTCGGAAGCATTTAAGCGGTTTTAACTGTCATTATTTGTGATGAAGATCACGTCAGAAAATTGTTATATTACTATGTTACGCATAACGTGATGTGCCTTGTAATTCTTATCAGTAGAAATaaaaa                                                                                                                                                                                          <  2:132740/250‑1 (MQ=255)                                                                 ttCCATTATTATTCCAATCTCTCGTCTGCTCGGTAACGGCAACCGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGTAGACACCTCGGAAGCATTTACGCGGTTTTAACTGTCATTATTTGTGATGAAGATCACGTCAGAAAATTGTTATATTACTATGTTACGCATAACGTGATGTGCCTTGTAATTCTTATCAGTAGAAATAAAAAAACGTGAAatt                                                                                                                                                                               >  2:37605/1‑250 (MQ=255)                                                                 ttCCATTATTATTACAATCTATCGTCTGCTCGGTCACGGAAACAGTTTTCTTGCCCGCAACAACGCAGGAATGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGTAGACACCTCGGAAGCATTTAAGCGGTTTTAACTGTCATTATTTGTGATGACGATCACGTCAGAAAATTGTTATATTACTATGTTACGCATAACGTGATGTGCCTTGTAATTCTTATCAGTAGAAATAAAAAAACGTGAAatt                                                                                                                                                                               >  2:199812/1‑250 (MQ=255)                                                                     attattattCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGTAGACACCTCGGAAGCATTTAAGCGGTTTTAACTGTCATTATTTGTGATGAAGATCACGTCAGAAAATTGTTATATTACTATGTTACGCATAACGTGATGTGCCTTGTAATTCTTATCAGTAGAAATAAAAAAACGTGAAATtattat                                                                                                                                                                          <  2:136793/251‑1 (MQ=255)                                                                        attattCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGTAGACACCTCGGAAGCATTTAAGCGGTTTTAACTGTCATTATTTGTGATGAAGATCACGTCAGAAAATTGTTATATTACTATGTTACGCATAACGTGATGTGCCTTGTAATTCTTATCAGTAGAAATAAAAAAACGTGAAATTATTATgcc                                                                                                                                                                       >  1:76116/1‑251 (MQ=255)                                                                              ccAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGTAGACACCTCGGAAGCATTTAAGCGGTTTTAACTGTCATTATTTGTGATGAAGATCACGTCAGAAAATTGTTATATTACTATGTTACGCATAACGTGATGTGCCTTGTAATTCTTATCAGTAGAAATAAAAAAACGTGAAATTATTATGCCGCCAg                                                                                                                                                                  <  1:194131/250‑1 (MQ=255)                                                                                 aTCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGTAGACACCTCGGAAGCATTTAAGCGGTTTTAACTGTCATTATTTGTGATGAAGATCACGTCAGAAAATTGTTATATTACTATGTTACGCATAACGTGATGTGCCTTGTAATTCTTATCAGTAGAAATAAAAAAACGTGAAATTATTATGCCGCCAGGCgt                                                                                                                                                              >  1:260440/1‑251 (MQ=255)                                                                                     agcGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCATGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGTAGACACCTCGGAAGCATTTAAGCGGTTTTAACTGTCATTATTTGTGATAAAGATCACGTCAGAAAATTGTTATATTACTATGTTACGCATAACGTGATGTGCCTTGTAATTCTTATCAGTAGAAATAAAAAAACGTGAAATTATTATGCCGCCAGGCgtaaaa                                                                                                                                                          <  2:231271/249‑4 (MQ=255)                                                                                     aTCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGTAGACACCTCGGAAGCATTTAAGCGGTTTTAACTGTCATTATTTGTGATGAAGATCACGTCAGAAAATTGTTATATTACTATGTTACGCATAACGTGATGTGCCTTGTAATTCTTATCAGTa                                                                                                                                                                                                    <  1:287022/209‑1 (MQ=255)                                                                                     aTCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGTAGACACCTCGGAAGCATTTAAGCGGTTTTAACTGTCATTATTTGTGATGAAGATCACGTCAGAAAATTGTTATATTACTATGTTACGCATAACGTGATGTGCCTTGTAATTCTTATCAGTa                                                                                                                                                                                                    >  2:287022/1‑209 (MQ=255)                                                                                       cGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGTAGACACCTCGGAAGCATTTAAGCGGTTTTAACTGTCATTATTTGTGATGAAGATCACGTCAGAAAATTGTTATATTACTATGTTACGCATAACGTGATGTGCCTTGTAATTCTTATCAGTAGAAATAAAAAAACGTGAAATTATTATGCCGCCAGGCGTAGTATc                                                                                                                                                        >  1:228098/1‑251 (MQ=255)                                                                                                         aaCAGTTTTCTTTCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTTGTCAGGTATACACCTCGGAAGCATTTAAGCGGTTTTAACTGTCATTATTTGTGATGAAGATCACGTCAGAAAATTGTTATATTACTATGTTACGCATAACGTGATGTGCCTTTTAATTCTTATCAGTAGAAATAAAAAAACGTGAAATTATTATGCCGCCAGGCGTAGTATCGCAGCAGTTAAGATGAtt                                                                                                                                      <  1:51471/251‑1 (MQ=255)                                                                                                                            acaacGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGTAGACACCTCGGAAGCATTTAAGCGGTTTTAACTGTCATTATTTGTGATGAAGATCACGTCAGAAAATTGTTATATTACTATGTTACGCATAACGTGATGTGCCTTGTAATTCTTATCAGTAGAAATAAAAAAACGTGAAATTATTATGCCGCCAGGCGTAGTATCGCAGCAGGTAAGATGATTCAGGAGATTTTAAATGGCg                                                                                                                   <  1:253630/251‑1 (MQ=255)                                                                                                                                                                 cACTGGTCAGGTAGACACCTCGGAAGCATTTAAGCGGTTTTAACTGTCATTATTTGTGATGAAGATCACGTCAGAAAATTGTTATATTACTATGTTACGCATAACGTGATGTGCCTTGTAATTCTTATCAGTAGAAATAAAAAAACGTGAAATTATTATGCCGCCAGGCGTAGTATCGCAGCAGGTAAGATGATTCAGGAGATTTTAAATGGCGGGTGAAAGCTTTATTTTGATGGGCGTTTCAGGGAg                                                                                <  2:225847/249‑1 (MQ=255)                                                                                                                                                                 cACTGGTCAGGTAGACACCTCGGAAGCATTTAAGCGGTTTTAACTGTCATTATTTGTGATGAAGATCACGTCAGAAAATTGTTATATTACTATGTTACGCATAACGTGATGTGCCTTGTAATTCTTATCAGTAGAAATAAAAAAACGTGAAATTATTATGCCGCCAGGCGTAGTATCGCAGCAGGTAAGATGATTCAGGAGATTTTAAATGGCGGGTGAAAGCTTTATTTTGATGGGCGTTTCAGGGAGTg                                                                              <  2:179600/251‑1 (MQ=255)                                                                                                                                                                 cACTGGTCAGGTAGACACCTCGGAAGCATTTAAGCGGTTTTAACTGTCATTATTTGTGATGAAGATCACGTCAGAAAATTGTTATATTACTATGTTACGCATAACGTGATGTGCCTTGTAATTCTTATCAGTAGAAATAAAAAAACGTGAAATTATTATGCCGCCAGGCGTAGTATCGCAGCAGGTAAGATGATTCAGGAGATTTTAAATGGCGGGTGAAAGCTTTATTTTGATGGGCGTTTCAGGGAGTg                                                                              >  2:58417/1‑251 (MQ=255)                                                                                                                                                                        cAGGTAGACACCTCGGAAGCATTTAAGCGGTTTTAACTGTCATTATTTGTGATGAAGATCACGTCAGAAAATTGTTATATTACTATGTTACGCATAACGTGATGTGCCTTGTAATTCTTATCAGTAGAAATAAAAAAACGTGAAATTATTATGCCGCCAGGCGTAGTATCGCAGCAGGTAAGATGATTCAGGAGATTTTAAATGGCGGGTGAAAGCTTTATTTTGATGGGCGTTTCAGGGAGTGGTAAAAc                                                                       <  2:107935/251‑1 (MQ=255)                                                                                                                                                                                                aaGCGGGTTTAACTGTCATTATTTGTGATTAAGATCACGTCAGAAAATTGTTATATTACTATGTTACGCATAACGTGATGTGCCTTGTAATTCTTATCAGTAGAAATAAAAAAACGTGAAATTATTATGCCGCCAGGCGTAGTATCGCAGCAGGTCAGATGATTCAGGAGATTTTAAATGGCGGGTGAAAGCGTTATTTTGATGGGCGTTTCATGGAGTTGTAAAACATTAAATTGTAGCAATGTTGCcgc                                               >  2:221045/1‑251 (MQ=255)                                                                                                                                                                                                aaGCGGGTTTAACTGTCATTATTTGGGATGAAGATCACGTCAGAAAATTGTTATATTACTATGTTACGCATAACGTGATGTGCCTTGTAATTCTTATCAGTAGAAATAAAAAAACGTGAAATTATTATGCCGCCAGGCGTAGTATCGCAGCAGGTAAGATGATTCAGGAGATTTTAAATGGCGGGGGAAAGCTTTATTTTGATGGGCGTTTCAGGGAGTGGTAAAACATTAATTGGTAGAAAGGTTGCcgc                                               >  1:13039/1‑251 (MQ=255)                                                                                                                                                                                                              gTCATTATTTGTGATGAAGATCACGTCAGAAAATTGTTATATTACTATGTTACGCATAACGTGATGTGCCTTGTAATTCTTATCAGTAGAAATAAAAAAACGTGAAATTATTATGCCGCCAGGCGTAGTATCGCAGCAGGTAAGATGATTCAGGAGATTTTAAATGGCGGGTGAAAGCGTTATTTTGATGGGCGTTTCAGGGAGTGGTAAAACAGTAATTGGTAGCAAGGTTGCCGCGTTATTATCTGc                                   <  1:28532/249‑1 (MQ=255)                                                                                                                                                                                                                                         aGAAAATTGTTATATTACTATGTTACGCATAACGTGATGTGCCTTGTAATTCTTATCAGTAGAAATAAAAAAACGTGAAATTATTATGCCGCCAGGCGTAGTATCGCAGCAGGTAAGATGATTCAGGAGATTTTAAATGGCGGGTGAAAGCTTTATTTTGATGGGCGTTTCAGGGAGTGGTAAAACATTAATTGGTAGCAAGGTTGCCGCGTTATTATCTGCTAAATTTATTGATGGTGACGATCTTCATc      >  2:58293/1‑251 (MQ=255)                                                                                                                                                                                                                                             aaTTGTTATATTACTATGTTACGCATAACGTGATGTGCCTTGTAATTCTTATCAGTAGAAATAAAAAAACGTGAAATTATTATGCCGCCAGGCGTAGTATCGCAGCAGGTAAGATGATTCAGGAGATTTTAAATGGCGGGTGAAAGCTTTATTTTGATGGGCGTTTCAGGGAGTGGTAAAACATTAATTGGTAGCAAGGTTGCCGCGTTATTATCTGCTAAATTTATTGATGGTGACGATCTTCATCCAGc  <  1:7180/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                   TTTCCTTCCTGATAATAATGTAAATCGGAACCGCAAATTCCACCTCGGGTTATTTGTACTAATGTTCCATTATTATTCCAATCTATCGTCTGCTCGGTAACGGCAACAGTTTTCTTGCCCGCAACAACGCAGGACTGTGTTTTCACTTGCATAAATTCCTCACTGGTCAGGTAGACACCTCGGAAGCATTTAAGCGGTTTTAACTGTCATTATTTGTGATGAAGATCACGTCAGAAAATTGTTACATTACTATGTTACGCATAACGTGATGTGCCTTGTAATTCTTATCAGTAGAAATAAAAAAACGTGAAATTATTATGCCGCCAGGCGTAGTATCGCAGCAGGTAAGATGATTCAGGAGATTTTAAATGGCGGGTGAAAGCTTTATTTTGATGGGCGTTTCAGGGAGTGGTAAAACATTAATTGGTAGCAAGGTTGCCGCGTTATTATCTGCTAAATTTATTGATGGTGACGATCTTCATCCAGC  >  NC_000913/4494255‑4494741

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