Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,183,713	T→C	intergenic (‑62/+3)	yegX ← / ← thiD	putative hydrolase/bifunctional hydroxymethylpyrimidine kinase/phosphomethylpyrimidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,183,713	0	T	C	100.0%	123.4 / NA	36	intergenic (‑62/+3)	yegX/thiD	putative hydrolase/bifunctional hydroxymethylpyrimidine kinase/phosphomethylpyrimidine kinase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (22/14);  total (22/14)

ATGGAAATTAGCCCAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTTTCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGTGCTCACCACCAGGCGTGGAAGTGGTGAACCGGACCAATACCGTGACCAACTTCCAGCGTGTCGGCCTGGGCTAACGCCGATGAAAGCCAGCTTTTTGCCTCCTGTACGGTGTCAGCCCAGTTTGTATGGCGCGGGC  >  NC_000913/2183584‑2183849                                                                                                                                  |                                                                                                                                         aTGGAAATTAGCCCATGAGCGCATACTAGGGAGGTGAATTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTTTCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGCGCTCaccacc                                                                                                                                <  1:386396/140‑1 (MQ=255) aTGGAAATTAGCCCAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTTTCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGCGCTCaccacc                                                                                                                                <  1:491332/140‑1 (MQ=255)  tGGAAATTAGCCCAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTTTCATCCCTGAATGTCAGTGCCAAAGGCTAACAATAAACAAAGCAGCTATAGAACGGCGCTCAccacca                                                                                                                               >  2:386396/1‑140 (MQ=255)         tAGCCCAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTTTCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGCGCTCACCACCAGGCGTgg                                                                                                                        <  2:25400/140‑1 (MQ=255)            cccAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTTTCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGCGCTCACCACCAGGCGTGGAAg                                                                                                                     >  1:676571/1‑140 (MQ=255)            cccAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTTTCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGCGCTCACCACCAGGCGTGGAAg                                                                                                                     >  1:547591/1‑140 (MQ=255)                      cATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTTTCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGCGCTCACCACCAGGCGTGGAAGTGGTGAACCg                                                                                                           <  2:676571/140‑1 (MQ=255)                              agagGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTTTCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGCGCTCACCACCAGGCGTGGAAGTGGTGAACCGGACCAATa                                                                                                   >  1:72961/1‑140 (MQ=255)                                agGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTTTCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGCGCTCACCACCAGGCGTGGAAGTGGTGAACCGGACCAATAcc                                                                                                 >  2:76307/1‑140 (MQ=255)                                agGTGAGTTATTTTCGACTGGTAATTCTTAGCTGCATTGGTTTCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGCGCTCACCACCATGCGAGGCAGTGGTTAACCGGACCAATAcc                                                                                                 >  1:385901/1‑140 (MQ=255)                                 ggTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTTTCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAA‑‑‑AGCTATAGTACGGCGCTCACCACCAGGCGTGGAAGTGGTGAACCGGACCAATACCGTGa                                                                                             >  2:476354/1‑140 (MQ=255)                                   tGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTTTCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGCGCTCACCACCAGGCGTGGAAGTGGTGAACCGGACCAATACCGTg                                                                                              >  2:628443/1‑140 (MQ=255)                                   tGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTTTCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGCGCTCACCACCAGGCGTGGAAGTGGTGAACCGGACCAATACCGTg                                                                                              >  2:66664/1‑140 (MQ=255)                                     aGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTTTCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGCGCTCACCACCAGGCGTGGAAGTGGTGAACCGGACCAATACCGTGAc                                                                                            >  1:517719/1‑140 (MQ=255)                                     aGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTTTCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGCGCTCACCACCAGGCGTGGAAGTGGTGAACCGGACCAATACCGTGAc                                                                                            >  1:67786/1‑140 (MQ=255)                                        ttttttCGACTGGTAATTCTTAGCTGCATTGGTTTCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGCGCTCACCACCAGGCGTGGAAGTGGCGAACCGGACCAATACCGTGACCaa                                                                                         <  2:547591/140‑1 (MQ=255)                                                   ggTAATTCTTAGCTGCATTGGTTTCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGCGCTCACCACCAGGCGTGGAAGTGGTGAACCGGACCAATACCGTGACCAACTTCCAGCgtg                                                                              <  2:700574/140‑1 (MQ=255)                                                   ggTAATTCTTAGCTGCATTGGTATCATCCTTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGCGCTCACCACCAGGCGTGGAAGTGGTGAACCGGACCAATACCGTGACCAACTTCCAGCgtg                                                                              <  1:366113/140‑1 (MQ=255)                                                                  cATTGGTTTCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGCGCTCACCACCAGGCGTGGAAGTGGTGAACCGGACCTATAACGTTACCAACTTTCAGCGTGTCGGCCTGGGCTAAc                                                               >  2:258431/1‑140 (MQ=255)                                                                  cATTGGTTTCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGCGCTCACCACCAGGCGTGGAAGTGGTGAACCGGACCAATAGCGTAACCAACTTGCAGCCTGTCGGCCTGGGCTAAc                                                               >  2:238958/1‑140 (MQ=255)                                                                   aTTGGTTTCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGCGCTCACCACCAGGCGTGGAAGTGGTGAACCGGACCAATACCGTGACCAACTTCCAGCGTGTCGGCCTGGGCTAACg                                                              >  1:746609/1‑140 (MQ=255)                                                                            aTCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGCGCTCACCACCAGGCGTGGAAGTGGTGAACCGGACCAATACCGTGACCAACTTCCAGCGTGTCGGCCTGGGCTAACGCCGATGaaa                                                     >  2:357480/1‑140 (MQ=255)                                                                                    aTGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGCGCTCACCACCAGGCGTGGAAGTGGTGAACCGGACCAATACCGTGACCAACTTCCAGCGTGTCGGCCTGGGCTAACGCCGATGAAAGCCAGCtt                                             <  2:72961/140‑1 (MQ=255)                                                                                      gTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGCGCTCACCACCAGGCGTGGAAGTGGTGAACCGGACCAATACCGTGACCAACTTCCAGCGTGTCGGCGTGGCCTAACGCCGTTGAAAGCCAGCtttt                                           >  1:263783/1‑140 (MQ=255)                                                                                           tGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGCGCTCACCACCAGGCGTGGAAGTGGTGAACCGGACCAATACCGTGACCAACTTCCAGCGTGTCGGCCTGGGCTAACGCCGATGAAAGCCAGCTTTTTGc                                        <  1:704936/138‑1 (MQ=255)                                                                                           tGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGCGCTCACCACCAGGCGTGGAAGTGGTGAACCGGACCAATACCGTGACCAACTTCCAGCGTGTCGGCCTGGGCTAACGCCGATGAAAGCCAGCTTTTTGc                                        >  2:704936/1‑138 (MQ=255)                                                                                                  ggCTGACAATAACCAAAACAGCTATAGTACGTCGCTCACAACCAGGCGAGGAAGTGGTGAACCGGACCAATACCGTGACCAACTTCCAGCGTGTCGGCCTGGGCTAACGCCGATGAAAGCCAGCTTTTTGCCTCCTGTAc                               <  2:263783/140‑1 (MQ=255)                                                                                                                aaaGCAGCTATAGTAAGGCGCTTACCACCTGGCGTGGCAGTGGTGAACCGGACCAATACCGTGACCAACTTCCAGCGTGTCGGCCTGGGCTAACGCCGATGAAAGCCAGCTTTTTGCCTCCTGTACGGTGTCAGCCCAGt                 <  1:262203/140‑1 (MQ=255)                                                                                                                       cTATAGTACGGCGCTCACCACCAGGCGTGGAAGTGGTGAACCGGACCAATACCGTGACCAACTTCCAGCGTGTCGGCCTGGGCTAACGCCGATGAAAGCCAGTTTTTTGCATCCGGTTCGGTGTCAGCCCAGTTTGTATg          >  1:322795/1‑140 (MQ=255)                                                                                                                       cTATAGTACGGCGCTCACCACCAGGCGTGGAAGTGGTGAACCGGACCAATACCGTGACCAACTTCCAGCGTGTCGGCCTGGGCTAACGCCGATGAAAGCCAGCTTTTTGCCTCCTGTACGGTGTCAGCCCAGTTTGTATg          >  1:779289/1‑140 (MQ=255)                                                                                                                       cTATAGTACGGCGCTCACCACCAGGCGTGGAAGTGGTGAACCGGACCAATACCGTGACCAACTTCCAGCGTGTCGGCCTGGGCTAACGCCGATGAAAGCCAGCTTTTTGCCTCCTGTACGGTGTCAGCCCAGTTTGTATg          >  2:819533/1‑140 (MQ=255)                                                                                                                       cTATAGTACGGCGCTCACCACCAGGCGTGGAAGTGGTGAACCGGACCAATACCGTGACCAACTTCCAGCGTGTCGGCCTGGGCTAACGCCGATGAAAGCCAGATTTTTGCGTCCGGTACGGTGTGAGCCCAGTTTGTATg          >  1:231757/1‑140 (MQ=255)                                                                                                                          taGTACGGCGCTCACTACCAGGTGTTGAACTGGTGAACCGGACCAATACCGTGACCAACTTCCAGCGTGTCGGCCTGGGCTAACGCCGATGAAAGCCAGCTTTTTGCCTCCTGTACGGTGTCAGCCCAGTTTGTATGgcg       <  1:238958/140‑1 (MQ=255)                                                                                                                          taGTACGGCGCTCACAACCAGGCGTTGAACTGGTGAAACGGACCAATACCGTGACCAACTTCCAGCGTGTCGGCCTGGGCTAACGCCGATGAAAGCCAGCTTTTTGCCTCCTGTACGGTGTCAGCCCAGTTTGTATGgcg       <  1:258431/140‑1 (MQ=255)                                                                                                                            gTACGGCGCTCACCACCAGGCGTGGAAGTGGTGAACCGGACCAATACCGTGACCAACTTCCAGCGTGTCGGCCTGGGCTAACGCCGATGAAAGCCAGCTTTTTGCCTCCTGTACGGTGTCAGCCCAGTTTGTATGGcgcg     >  2:593043/1‑140 (MQ=255)                                                                                                                               cGGCGCTCACCACCAGGCGTGGAAGTGGTGAACCGGACCAATACCGTGACCAACTTCCAGCGTGTCGGCCTGGGCTAACGCCGATGAAAGCCAGCTTTTTGCCTCCTGTACGGTGTCAGCCCAGTTTGTATGGCGCGGGc  <  1:66664/140‑1 (MQ=255)                                                                                                                                  |                                                                                                                                         ATGGAAATTAGCCCAAGAGCGCATAATAGAGAGGTGAGTTTTTTTCGACTGGTAATTCTTAGCTGCATTGGTTTCATCCCTGAATGTCAGTGCCAAAGGCTGACAATAACCAAAGCAGCTATAGTACGGTGCTCACCACCAGGCGTGGAAGTGGTGAACCGGACCAATACCGTGACCAACTTCCAGCGTGTCGGCCTGGGCTAACGCCGATGAAAGCCAGCTTTTTGCCTCCTGTACGGTGTCAGCCCAGTTTGTATGGCGCGGGC  >  NC_000913/2183584‑2183849

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