Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	251,629	C→T	G244G (GGC→GGT)	dinB →	DNA polymerase IV

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	251,629	0	C	T	100.0%	114.4 / NA	32	G244G (GGC→GGT)	dinB	DNA polymerase IV
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (13/19);  total (13/19)

TTGTGATTACGCCGGCAGAAGTTCCGGCATTTTTACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAAGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTCGACGATTTTCAGCAAACCACCCAGGAGCACGTCTGGCCGCGGCTGAATAAAGCTGATCTAATCGCCACCGCGCGTAAAACCTGGGATGAACG  >  NC_000913/251382‑251865                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                             ttGTGATTACGCCGGCAGAAGTTCCGGCATTTTTACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAAGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGtgtg                                                                                                                                                                                                                                           >  1:194087/1‑251 (MQ=255)                       tCCGGCATTTTTACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAAGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGat                                                                                                                                                                                                                      <  1:72351/250‑1 (MQ=255)                        ccGGCATTTTTACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAAGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGatat                                                                                                                                                                                                                    <  1:154628/251‑1 (MQ=255)                                tttACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAAGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCAc                                                                                                                                                                                                             <  1:352497/250‑1 (MQ=255)                                            aCCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAAGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCAGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAAtgt                                                                                                                                                                                                 <  1:59848/250‑1 (MQ=255)                                                                     ggCAAAGTCTCAGCGGCAAAACTGGAAGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATc                                                                                                                                                                                     >  1:180063/1‑237 (MQ=255)                                                                     ggCAAAGTCTCAGCGGCAAAACTGGAAGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATc                                                                                                                                                                                     <  2:180063/237‑1 (MQ=255)                                                                        aaaGTCTCAGCGGCAAAACTGGAAGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCgg                                                                                                                                                                     >  1:78245/1‑250 (MQ=255)                                                                         aaGTCTCAGCGGCAAAACTGGAAGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTTGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGCTTATCGAGCGGCTGTATCCGGa                                                                                                                                                                    >  1:284487/1‑250 (MQ=255)                                                                               cAGCGGCAAAACTGGAAGCGATGGGGCTGCGGACCTGCGGTGATTTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAAGGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGaa                                                                                                                                                             <  2:342932/251‑1 (MQ=255)                                                                                              aaGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAg                                                                                                                                              <  2:202335/251‑1 (MQ=255)                                                                                                   aTGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTaa                                                                                                                                          >  1:308823/1‑250 (MQ=255)                                                                                                     ggggCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAAc                                                                                                                                       >  1:50145/1‑251 (MQ=255)                                                                                                       ggCTGCGGACCTGCGGAGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAACCCt                                                                                                                                     >  2:326577/1‑251 (MQ=255)                                                                                                                        gATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCg                                                                                                                    <  1:271061/251‑1 (MQ=255)                                                                                                                               aaaaGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAgggg                                                                                                             >  1:12879/1‑251 (MQ=255)                                                                                                                                   gtgtGATCTGGTGATGCTGCTTAAACGCTTTGGCATATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGTTTATCGAGCGGCTGTATCCGGAGCTTGAACGCCGTCTGGCAATGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGa                                                                                                         <  1:341270/251‑1 (MQ=255)                                                                                                                                        aTCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAAtt                                                                                                     >  2:92291/1‑250 (MQ=255)                                                                                                                                            ggTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAAtt                                                                                                     <  1:186912/246‑1 (MQ=255)                                                                                                                                               gATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTc                                                                                              <  1:87012/250‑1 (MQ=255)                                                                                                                                                    tgctTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTCGACGAt                                                                                        >  2:320186/1‑251 (MQ=255)                                                                                                                                                     gctTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTGAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTCGACGAtt                                                                                       <  1:59210/251‑1 (MQ=255)                                                                                                                                                       ttAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTCGACGAtttt                                                                                     >  2:3280/1‑251 (MQ=255)                                                                                                                                                        tAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTCGACGAtttt                                                                                     <  1:210666/250‑1 (MQ=255)                                                                                                                                                                ttGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGAAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTCGACGATTTTCAGCAAacc                                                                            >  1:6487/1‑251 (MQ=255)                                                                                                                                                                    cAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTCGACGATTTTCAGCAAACCAccc                                                                        <  2:39753/251‑1 (MQ=255)                                                                                                                                                                             ccGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTCGACGATTTTCAGCAAACCACCCAGGAGCACg                                                               >  2:298006/1‑251 (MQ=255)                                                                                                                                                                                    ttGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTCGACGATTTTCAGCAAACCACCCAGGAGCACGTCTGGc                                                         <  1:249653/250‑1 (MQ=255)                                                                                                                                                                                                     ggggTTTGACGAACGAGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGGTGGCTGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCGGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTCGACGATTTTCAGCAAACCACCCAGGAGCACGTCTGGCCGCGGCTGAATAAAGCTg                                       <  2:304354/251‑1 (MQ=255)                                                                                                                                                                                                      gggATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTCGACGATTTGCAGCAAACCACCCAGGAGCACGTCTGGCCGCGGCTGAATAAAGCTGa                                      <  1:139813/251‑1 (MQ=255)                                                                                                                                                                                                       ggATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTCGACGATTTTCAGCAAACCACCCAGGAGCACGTCTGGCCGCGGCTGAATAAAgatgat                                     <  2:308823/251‑1 (MQ=255)                                                                                                                                                                                                                             cAGCGAACGGTTGCGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTCGACGATTTTCAGCAAACCACCCAGGAGCACGTCTGGCCGCGGCTGAATAAAGCTGATCTAATCGCCACCGCGCGTaaaa               <  2:119915/251‑1 (MQ=255)                                                                                                                                                                                                                                          cGAAAATCCGTCGGTGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTCACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTCGACGATTTTCAGCAAACCACCCAGGAGCACGTCTGGCCGCGGCTGAATAAAGCTGATCTAATCGCCACCGCGCGTAAAACCTGGGATGAAcg  >  2:251270/1‑251 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                             TTGTGATTACGCCGGCAGAAGTTCCGGCATTTTTACAAACCTTACCGCTGGCAAAAATCCCCGGCGTCGGCAAAGTCTCAGCGGCAAAACTGGAAGCGATGGGGCTGCGGACCTGCGGTGATGTACAAAAGTGTGATCTGGTGATGCTGCTTAAACGCTTTGGCAAATTTGGCCGCATTTTGTGGGAGCGTAGTCAGGGGATTGACGAACGCGATGTTAACAGCGAACGGTTGCGAAAATCCGTCGGCGTGGAACGCACGATGGCGGAAGATATTCATCACTGGTCTGAATGTGAAGCGATTATCGAGCGGCTGTATCCGGAACTTGAACGCCGTCTGGCAAAGGTAAAACCTGATTTACTGATTGCTCGCCAGGGGGTGAAATTAAAGTTCGACGATTTTCAGCAAACCACCCAGGAGCACGTCTGGCCGCGGCTGAATAAAGCTGATCTAATCGCCACCGCGCGTAAAACCTGGGATGAACG  >  NC_000913/251382‑251865

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