Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	700,928	C→A	G223* (GGA→TGA)	nagC ←	N‑acetylglucosamine‑inducible nag divergent operon transcriptional repressor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	700,928	0	C	A	100.0%	211.0 / NA	54	G223* (GGA→TGA)	nagC	N‑acetylglucosamine‑inducible nag divergent operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (18/36);  total (18/36)

GCGTTGGCAGCGATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTT  >  NC_000913/700780‑701074                                                                                                                                                     |                                                                                                                                                   gcgcTGGCAGCGATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCagc                                                                                                                                                  >  2:351434/1‑151 (MQ=255) gCGTTGGCAGCGATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCagc                                                                                                                                                  <  1:467163/151‑1 (MQ=255)     tGGCAGCGATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGAt                                                                                                                                               <  1:873015/150‑1 (MQ=255)         aGCGATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGaca                                                                                                                                          <  2:930288/151‑1 (MQ=255)            gATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGc                                                                                                                                       >  1:512361/1‑151 (MQ=255)                gTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCAcc                                                                                                                                    >  1:821906/1‑150 (MQ=255)                 tttCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCAcc                                                                                                                                    >  2:763630/1‑149 (MQ=255)                    ccAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAAt                                                                                                                               <  2:108925/151‑1 (MQ=255)                          aGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTc                                                                                                                         >  1:173630/1‑151 (MQ=255)                             cAAAGTTTCCGCAGTGGCAGCGTTCCCCCAGGGGTTCGACCTGAATTTTGCCAAATTCCCCGACGTTTCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATTATCCCGGCCCCGGTTCAGCGATGGACCCCCACCAGAATTGAGTCtt                                                                                                                       <  2:973525/150‑1 (MQ=255)                                 gTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCa                                                                                                                   >  2:690029/1‑150 (MQ=255)                                    gCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATc                                                                                                                <  2:965915/150‑1 (MQ=255)                                            ggCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTg                                                                                                        <  1:125504/150‑1 (MQ=255)                                             gCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGCGTCTTCGCAATCCTGACTTGc                                                                                                       >  2:301723/1‑150 (MQ=255)                                              cAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCAc                                                                                                     <  2:651593/151‑1 (MQ=255)                                               aGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCAcc                                                                                                    >  1:347618/1‑151 (MQ=255)                                                 cGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGa                                                                                                  <  1:952822/151‑1 (MQ=255)                                                    tCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAagt                                                                                               >  2:890101/1‑151 (MQ=255)                                                      aCCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAagt                                                                                               >  1:318193/1‑149 (MQ=255)                                                       cccAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAagtagt                                                                                            <  2:837654/151‑1 (MQ=255)                                                         cAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGc                                                                                          <  2:72158/151‑1 (MQ=255)                                                            cGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTcc                                                                                       <  2:611031/151‑1 (MQ=255)                                                              gTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCcgc                                                                                     <  2:167330/151‑1 (MQ=255)                                                                 cGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCcgcca                                                                                   <  2:787948/150‑1 (MQ=255)                                                                  gACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCcgccagc                                                                                 >  1:555834/1‑151 (MQ=255)                                                                       gAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCcgccagcgccag                                                                            <  1:247035/151‑1 (MQ=255)                                                                                cAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGat                                                                   <  2:923557/151‑1 (MQ=255)                                                                                cAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGat                                                                   <  2:306576/151‑1 (MQ=255)                                                                                    ttCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATc                                                                <  1:821519/150‑1 (MQ=255)                                                                                      cACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGt                                                              <  1:667624/150‑1 (MQ=255)                                                                                      cACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGt                                                              <  2:422157/150‑1 (MQ=255)                                                                                       aCCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGt                                                              <  1:118760/149‑1 (MQ=255)                                                                                         cGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGAcc                                                          <  1:100085/151‑1 (MQ=255)                                                                                         cGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGAcc                                                          <  1:226054/151‑1 (MQ=255)                                                                                         cGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGAcc                                                          >  2:445454/1‑151 (MQ=255)                                                                                          gccgtggccgttgcGGCCAATAAAAATGCGCTCGTTAGAGATGATCCCGAACCCGTATCAGCGATGGACGCGCCCCAGAATGGAGTCTTCGCAATCCTATCGGGCACCGAAGTGGTGCTCCGCCAGCGCCAGACttggg                                                                     <  2:268732/137‑5 (MQ=255)                                                                                             gttgccgttgcGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACg                                                      >  1:765538/1‑151 (MQ=255)                                                                                             gttgccgttgcGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACg                                                      <  1:143164/151‑1 (MQ=255)                                                                                              ttgccgttgcGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACg                                                      >  2:673507/1‑150 (MQ=255)                                                                                               tgccgttgcGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGa                                                     >  2:892120/1‑150 (MQ=255)                                                                                                       cGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCa                                            <  2:330459/151‑1 (MQ=255)                                                                                                          ccAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACtt                                         <  2:450644/151‑1 (MQ=255)                                                                                                            aaTAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTa                                       <  2:364800/151‑1 (MQ=255)                                                                                                                aaaaTGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAAc                                   >  1:293292/1‑151 (MQ=255)                                                                                                                   aTGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGtt                                >  2:16322/1‑151 (MQ=255)                                                                                                                   aTGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGtt                                <  2:417571/151‑1 (MQ=255)                                                                                                                     gcgcCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGttc                               <  1:263495/150‑1 (MQ=255)                                                                                                                       gcCCGTTAGAGATAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGttcttc                            <  1:163077/151‑1 (MQ=255)                                                                                                                                agaTAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTc                   >  1:234933/1‑151 (MQ=255)                                                                                                                                  aTAATCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTa                 <  2:347612/151‑1 (MQ=255)                                                                                                                                     aTCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCa              <  2:173626/151‑1 (MQ=255)                                                                                                                                     aTCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCa              <  2:565590/151‑1 (MQ=255)                                                                                                                                     aTCCCGGCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCa              <  1:689508/151‑1 (MQ=255)                                                                                                                                          ggCCCCGGTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGcccc         >  1:472702/1‑151 (MQ=255)                                                                                                                                                ggTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGttt   <  1:150866/151‑1 (MQ=255)                                                                                                                                                 gTTCAGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGtttt  <  1:106156/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GCGTTGGCAGCGATAGTTTCCAGGCAGCCAAAGTTGCCGCAGTGGCAGCGTTCACCCAGCGGTTCGACCTGAATATGGCCAATTTCACCGACGTTGCCGTTGCGGCCAATAAAAATGCGCCCGTTAGAGATAATCCCGGCCCCGGTTCCGCGATGGACACGCACCAGAATGGAGTCTTCGCAATCCTGACTTGCACCGAAGTAGTGCTCCGCCAGCGCCAGACTACGGATATCGTGACCAACGAAACAGGTCACTTTAAAACGTTCTTCCAGAGCTTCTACCAGCCCCCAGTTTT  >  NC_000913/700780‑701074

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