Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,951,184	G→A	R118H (CGC→CAC)	yecD →	isochorismatase family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,951,184	0	G	A	100.0%	201.7 / NA	52	R118H (CGC→CAC)	yecD	isochorismatase family protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (32/20);  total (32/20)

CGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCGCCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAAC  >  NC_000913/1951037‑1951329                                                                                                                                                    |                                                                                                                                                  cGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGccaccg                                                                                                                                                >  2:571239/1‑151 (MQ=255)   gTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAg                                                                                                                                              <  1:595201/151‑1 (MQ=255)     cGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGt                                                                                                                                            >  1:271362/1‑151 (MQ=255)     cGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGt                                                                                                                                            <  2:573226/151‑1 (MQ=255)     cGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGt                                                                                                                                            <  1:909331/151‑1 (MQ=255)          cgcCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGa                                                                                                                                       <  2:261515/151‑1 (MQ=255)            cccTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATa                                                                                                                                     <  1:559027/151‑1 (MQ=255)                   cGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAgt                                                                                                                              <  1:509109/151‑1 (MQ=255)                    ggCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAgtg                                                                                                                             >  1:1064968/1‑151 (MQ=255)                     gCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAgtg                                                                                                                             >  1:1038088/1‑150 (MQ=255)                       tAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGtt                                                                                                                           >  2:579977/1‑150 (MQ=255)                             gtTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTggg                                                                                                                    <  1:350150/151‑1 (MQ=255)                              ttGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTggg                                                                                                                    >  2:966191/1‑150 (MQ=255)                              ttGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTggg                                                                                                                    >  1:997129/1‑150 (MQ=255)                                gCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTGGGAtc                                                                                                                 <  2:1049720/151‑1 (MQ=255)                                gCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTGGGAtc                                                                                                                 <  2:883250/151‑1 (MQ=255)                                gCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTGGGAt                                                                                                                  <  1:825492/150‑1 (MQ=255)                                      aaaTTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGAcc                                                                                                           >  2:65789/1‑151 (MQ=255)                                       aaTTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGAcc                                                                                                           >  2:1044666/1‑150 (MQ=255)                                       aaTTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCa                                                                                                          <  2:654507/151‑1 (MQ=255)                                          tggtggCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAat                                                                                                        >  1:740797/1‑150 (MQ=255)                                             tggCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCg                                                                                                    >  2:687375/1‑151 (MQ=255)                                               gCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGgt                                                                                                  >  1:796025/1‑151 (MQ=255)                                                    aTCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTg                                                                                              >  2:904815/1‑150 (MQ=255)                                                          ctgcATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCAc                                                                                       >  1:492126/1‑151 (MQ=255)                                                                tAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCAccgcccg                                                                                 >  1:924986/1‑151 (MQ=255)                                                                tAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCAccgcccg                                                                                 >  1:598850/1‑151 (MQ=255)                                                                   gTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCaa                                                                              <  2:979014/151‑1 (MQ=255)                                                                   gTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCaa                                                                              <  2:861554/151‑1 (MQ=255)                                                                    tGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAAt                                                                             >  2:428857/1‑151 (MQ=255)                                                                        accgaccgCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCt                                                                         >  1:184402/1‑151 (MQ=255)                                                                        aCCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCt                                                                         >  2:254536/1‑151 (MQ=255)                                                                               gCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTGGCATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACt                                                                  <  1:737477/151‑1 (MQ=255)                                                                                cGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTCGAATCCACCGCCCGCAATGCCTGGGAACtt                                                                 >  1:25575/1‑151 (MQ=255)                                                                                 gATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTg                                                                >  2:810920/1‑151 (MQ=255)                                                                                      cGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCttt                                                           <  2:356690/151‑1 (MQ=255)                                                                                           tcatcaAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCt                                                      >  2:490771/1‑151 (MQ=255)                                                                                                  aCGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGTCCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGAtt                                               >  1:697219/1‑151 (MQ=255)                                                                                                  aCGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGAtt                                               >  1:797664/1‑151 (MQ=255)                                                                                                      cAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCg                                           >  2:17772/1‑151 (MQ=255)                                                                                                       aaTGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGa                                          <  2:174357/151‑1 (MQ=255)                                                                                                        aTGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGaa                                         <  2:628649/151‑1 (MQ=255)                                                                                                           gggTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGCGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGAc                                      >  1:656884/1‑151 (MQ=255)                                                                                                               gCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCt                                  >  1:364596/1‑151 (MQ=255)                                                                                                                  ttttACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGt                                >  2:659293/1‑150 (MQ=255)                                                                                                                   tttACGGTACGGATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGt                                >  2:694633/1‑149 (MQ=255)                                                                                                                        ggtacggATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTgccg                         <  2:187084/151‑1 (MQ=255)                                                                                                                         gtacggATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTgccgc                        <  1:449721/151‑1 (MQ=255)                                                                                                                           acggATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTgccgcc                       >  2:1042954/1‑150 (MQ=255)                                                                                                                              gATCTGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAgcg                   >  2:875783/1‑151 (MQ=255)                                                                                                                                  tGGAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCg                >  1:500614/1‑150 (MQ=255)                                                                                                                                   ggAGTTGCAATTACGCCACCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCg                >  2:681892/1‑149 (MQ=255)                                                                                                                                              tACGCCACCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATa    <  2:277919/150‑1 (MQ=255)                                                                                                                                                cgccacCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAAc  >  1:707464/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                  CGGTCGATGCGCCCTCACCGGCTAAAGTGTTGCCCGAAAATTGGTGGCAACATCCTGCTGCATTAGGTGCAACCGACAGCGATATCGAAATCATCAAACGTCAATGGGGTGCGTTTTACGGTACGGATCTGGAGTTGCAATTACGCCGCCGAGGTATCGATACAATAGTGTTATGTGGGATCTCGACCAATATCGGTGTTGAATCCACCGCCCGCAATGCCTGGGAACTTGGCTTTAATCTGGTGATTGCCGAAGACGCCTGTAGTGCCGCCAGCGCCGAGCAGCACAATAAC  >  NC_000913/1951037‑1951329

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