Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,931,805	C→T	R251W (CGG→TGG)	ygdG →	flap endonuclease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,931,805	0	C	T	100.0%	100.2 / NA	34	R251W (CGG→TGG)	ygdG	flap endonuclease
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (8/26);  total (8/26)

GGGGACTTGCGGGGATCAGCAGTTCAAAGGTACCGGGTGTTGCGGGAATCGGACCAAAAAGCGCCACGCAGCTGCTGGTCGAGTTTCAGAGTCTGGAAGGGATATATGAGAATCTGGATGCGGTT‑GCCGAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAA‑CCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCA‑AACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCGGAATGTCTTCCTTGCGTACACCAACATCACGCAAATGTGGCGGAATACCGACAT  >  NC_000913/2931560‑2932019                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                        ggggACTTTCGGGGATCAGAAGTTCAAAGGTACCGGGTGTTGCGGGAATCGGACCAAAAAGGTCCACGCAGATGCTGGTCGAGTTTCAGAGTCTGGAAGGGATATATGATAATCTGGATGCGGTT‑GCCGAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCGGTGCCGCGATATTGCCCGCTTACAAA‑CCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTtggtatggtac                                                                                                                                                                                                                    <  2:604769/251‑2 (MQ=255)      cTTGCGGGGATCAGCAGTTCAAAGGTACCGGGGGTTGCGGGAATCGGTCCAAAAAGCGCCACGCAGCTGCTGGTCGAGTTTCAGAGTCTGGCAGGGATATATGAGAATCTGGATGCGGTT‑GCCGAAAAGTGGCGCAAAAA‑TTAGAAACCCATAACGAGAAGGCGTTTCTGTGCCGCGATATTGCCCACTTACAAA‑CCGATTTGAATATCGACGGCAATTTACCGCAATCGCGGTTGGTATGGTAACGGCg                                                                                                                                                                                                               >  1:403376/1‑250 (MQ=255)        tGCGGGGATCAGCAGTTCAAAGGGACCGGGTGTTGCGGGAATCGGACCAAAAAGCGCCACGCAGCTGCTGGGCGAGTTTCAGAGTCTGGAAGGGATATAGGAGAATCTGGATGGGGTT‑GCCGAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGGGCCGCGATAGTGCCCGCTTACAAA‑CCGATTTGCATATCGACGGCAATTTACAGCAATTTCGGTTGGTATGGTAACGGCGa                                                                                                                                                                                                              <  1:513792/250‑1 (MQ=255)         gCGGGGATCAGCAGTTCAAAGGTACCGGGTGTTGCGGGAATCGGACCAAAAAGCGCCACGCAGCTGCTGGTCGAGTTTCAGAGTCTGGAAGGGATATATGAGAATCTGGATGCGGTT‑GCCGAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAA‑CCGATTTGCATATCGACGGCAAGTTACAGCAATTGCGGTTGGTATGGTAACGGCGAg                                                                                                                                                                                                             <  2:301648/250‑1 (MQ=255)               aTCAGCAGTTCAAAGGTACCGGGTGTTGCGGGAATCGGACCAAAAAGCGCCACGCAGCTGCTGGTCGAGTTTCAGAGTCTGGAAGGGATATATGAGAATCTGGATGCGGTT‑GCCGAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAA‑CCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAGCCGGATa                                                                                                                                                                                                      <  2:169032/251‑1 (MQ=255)                  agcagTTCAAAGGTACCGGGTGTTGCGGGAATCGGACCAAAAAGCGCCACGCAGCTGCTGGTCGAGTTTCAGAGTCTGGAAGGGATATATGAGAATCTGGATGCGGTT‑GCCGAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAA‑CCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAGCCGGATAcgc                                                                                                                                                                                                   <  1:565858/251‑1 (MQ=255)                          aaaGGTACCGGGTGTTGCGGGAATCGGACCAAAAAGCGCCACGCAGCTGCTGGTCGAGTTTCAGAGTCTGGAAGGGATATATGAGAATCTGGATGCGGTT‑GCCGAAAAGTGGCGCAAAAAATTAGAAACACATAAAGAGATGGCGTTTATGTGCCGCGATATTGCCCGCTTACAAA‑CCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAGCCGGaagcgcggccaacg                                                                                                                                                                                           >  1:610057/1‑238 (MQ=255)                                 ccGGGTGTTGCGGGAATCGGACCAAAAAGCGCCACGCAGCTGCTGGTCGAGTTTCAGAGTCTGGAAGGGATATATGAGAATCTGGATGCGGTT‑GCCGAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAA‑CCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAACCGGGTACGCCGCA‑CACGTCGCATc                                                                                                                                                                                   >  1:209230/1‑251 (MQ=255)                                    ggTGTTGCGGGAAGCGGACCAAAAAGCGCCACGCAGCTGCTGGTCGAGTTTCAGAGTCTGGAAGGGATATATGAGAATCTGGATGCGGTT‑GCCGAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAA‑CCGATTTGCATAGCGACGGCAATGTACAGCAATTGCGGTTGGTATGGTAACGGCGAGCCGGATACGCCGCA‑AACGTCGTATCCgg                                                                                                                                                                                <  2:295415/251‑1 (MQ=255)                                              gAATCGGACCAAAAAGCGCCACGCAGCTGCTGGTCGAGTTTCAGAGTCTGGAAGGGATATATGAGAATCTGGATGCGGTT‑GCCGAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAA‑CCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAGCCCGATACGCCGCa‑aaagccggctcccgccttattac                                                                                                                                                                       >  1:594204/1‑229 (MQ=255)                                                                 ccccgcggttgcggttggggttccgagTCTGGGAAGTATATATGGGAATCTGGATGGGGTTTGCC‑AAAAGTGGCGCCAAAAATTAGGAACCCATAAAGAGATTGCCGTTCCGGGGCGCGATATTGCCCGCTTACAAAACGGATTT‑CATAACGACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAGCCGGATACGCCGCA‑AACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCg                                                                                                                                                   <  2:335903/227‑1 (MQ=255)                                                                        ctgctgGTCGAGTTTCAGAGTCTGGAAGGGATATATGAGAATCTGGATGCGGTT‑GCCGAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAA‑CCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAGCGGGAGAAGACGCA‑AACAACGTAACCGGCGTTATCAAATCAccgccatttcccggcgggagtgt                                                                                                                                            >  1:17315/1‑232 (MQ=255)                                                                          gctgGTCGAGTTTCAGAGTCTGGAAGGGATATATGAGAATCTGGATGCGGTT‑GCCGAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAA‑CCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAGCCGGGTAACCCCCC‑AACCACGTATTCGGCATTATTACATCAGCGCGTTTTCCAAGCggcctggga                                                                                                                                           >  2:193564/1‑243 (MQ=255)                                                                          gctgGTCGAGTTTCAGAGTCTGGAAGGGATATATGAGAATCTGGATGCGGTT‑GCCGAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAA‑CCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAGCCGGATACGCCGCA‑AACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTa                                                                                                                                           <  1:8145/250‑1 (MQ=255)                                                                             ggTCGAGTTTCAGAGTCTGGAAGGGATATATGAGAATCTGGATGCGGTT‑GCCGAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAA‑CCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGTTATGGTAACGGCGAGCCGGATACGCCGCA‑AACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGc                                                                                                                                       <  1:287533/251‑1 (MQ=255)                                                                                      gcAGAGTCTGGAAGTGATATATGAGAATCTGGATGCTGTT‑GCCGAAAAGTGGCGCAAAAAATCAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAA‑CCGATTTGCATATCTACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAGCCGGATACGCCGCA‑AACGTCGTATCCGGCGTTTTCACACCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAat                                                                                                                               <  1:54600/249‑1 (MQ=255)                                                                                               ggAAGGGATATATGAGAATCTGGATGCGGTT‑GCCGAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAA‑CCGATGTGCATATCGACTGCAATTTACAGCAATTGCGGTTGGTATGGTAAGGGCGAGCCGGATACGCCGCA‑AACGTCGTATCTGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCa                                                                                                                        <  2:235431/248‑1 (MQ=255)                                                                                                     gATATATGAGAATCTGGATGCGGTT‑GCCGAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAA‑CCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAGCCGGATACGCCGCA‑AACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCt                                                                                                               <  1:619430/251‑1 (MQ=255)                                                                                                                aTCTGGATGCGGTT‑GCCGAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAA‑CCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAGCCGGATACGCCGCACACCGCCG‑ATCCGCCGTTATCACAACATCGCATTTTCCACACCGGGTGGTAAAGCTCTACAAAATCCTCAACCGATGCTTCACGCggggt                                                                                                    >  1:473662/1‑249 (MQ=255)                                                                                                                                ccGAAAAGTGGCGAAAAAAAGTAGAAACCCATAAAGAGAGTGCGTTTCTGTGCCGCGATATGGCCCGCTTACAAA‑CCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAGCCGGATACGCCGCA‑AACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTAGGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGTTTGCCACCGGTACAAAc                                                                                     <  2:408241/251‑1 (MQ=255)                                                                                                                                 cGAAAAGTGGAGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGGGCCGCGATATTGCCCGCTTACAAA‑CCGATGTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAGCCGGATACGCCGCA‑AACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAAca                                                                                    <  2:266421/251‑1 (MQ=255)                                                                                                                                  gAAAAGTGGCGCAAAAAAATAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAA‑CCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAGCCGGATACGCCGCA‑AACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTAc                                                                                                                                  <  2:432981/204‑1 (MQ=255)                                                                                                                                             cAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAA‑CCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAGCCGGATACGCCGCA‑AACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTgcc                                                                        <  1:8542/251‑1 (MQ=255)                                                                                                                                             cAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAA‑CCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAGCCGGATACGCCGCA‑AACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTgcc                                                                        <  1:424002/251‑1 (MQ=255)                                                                                                                                                 aaaTTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAA‑CCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAGCCGGGTAACCCCCC‑AACATCGGATCCGGCCTTATCACAACATCGCATTTACCAAGCGGTATGGGGAAACTCTACAATATCCTCAAACGTTGCTTTCCGCGGGGTGCCACCCGGACCAAAACCATCCAGTgcccccc                                                                    >  1:18626/1‑247 (MQ=255)                                                                                                                                                 aaaTTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATAGTGCCCGCTTACAAA‑CCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAGCCGGATACGCCGCA‑AACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGTGTTGCCACCGGTACAAACATCATCCAGTGCCGCCt                                                                    <  1:128713/251‑1 (MQ=255)                                                                                                                                                     tAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAA‑CCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAGCCGGATACGCCGCA‑AACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTgcg                                                                 <  1:523803/250‑1 (MQ=255)                                                                                                                                                            ccATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAA‑CCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAGCCGGATACGCCGCA‑AACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCAGCCAGTGCCGCCTGCGCCAGTGcc                                                         <  1:161723/251‑1 (MQ=255)                                                                                                                                                                 aaGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAA‑CCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAGCCGGATACGCCGCA‑AACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCGGAAt                                                    <  1:168415/251‑1 (MQ=255)                                                                                                                                                                        ggCGTTTCTGTGCCGCGATATTGCCCGCTTACAAA‑CCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAGCCGGATACGCCGCA‑AACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCGGAATGTCTTc                                              <  1:351724/250‑1 (MQ=255)                                                                                                                                                                                 gtgCCGCGATATTGCCCGCTTACAAA‑CCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAGCCGGATACGCCGCA‑AACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCGGAATGTCTTCCTTGCGTac                                     <  1:33426/250‑1 (MQ=255)                                                                                                                                                                                        gATATTGCCCGCTTACAAA‑CCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAGCCGGATACGCCGCA‑AACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCGGAATGTCTTCCTTGCGTACACCAACAt                             <  1:252921/251‑1 (MQ=255)                                                                                                                                                                                                  gCTTACAAA‑CCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAGCCGGATACGCCGCA‑AACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCGGAATGTCTTCCTTGCGTACACCAACATCACGCAAAtg                   <  1:169956/251‑1 (MQ=255)                                                                                                                                                                                                      aCAAA‑CCGATTGGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAGCCGGATACGCCGCA‑AACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCGGAATGTCTTCCGTGCGTACACCAACATCACGCAAATGTgg                <  2:583092/250‑1 (MQ=255)                                                                                                                                                                                                                    cATATCGACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAGCCGGATACGCCGCA‑AACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCGGAATGTCTTCCTTGCGTACACCAACATCACGCAAATGTGGCGTAATACCGACAt  <  2:522735/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                        GGGGACTTGCGGGGATCAGCAGTTCAAAGGTACCGGGTGTTGCGGGAATCGGACCAAAAAGCGCCACGCAGCTGCTGGTCGAGTTTCAGAGTCTGGAAGGGATATATGAGAATCTGGATGCGGTT‑GCCGAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAA‑CCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCA‑AACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCGGAATGTCTTCCTTGCGTACACCAACATCACGCAAATGTGGCGGAATACCGACAT  >  NC_000913/2931560‑2932019

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