Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,820,554	T→C	T221A (ACC→GCC)	chbC ←	N,N'‑diacetylchitobiose‑specific PTS enzyme IIC component

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,820,554	0	T	C	100.0%	109.0 / NA	34	T221A (ACC→GCC)	chbC	N,N'‑diacetylchitobiose‑specific PTS enzyme IIC component
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (17/17);  total (17/17)

AAAGGAGTCCAGCATCGGCTTGGCCCAGATATGGAAGGTCTTACCGGCTGCCAGCGCCGCTTCGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAGGTACGCTGTCGGGCAGTTTAATGACCCAATTGCGGCGGACAATAAAGGTAAACATTTCTGCCACCACCAGGCCAATAATAATCCCGGAGATGATATTCGCCCCACCTAACCAGTTTGCA  >  NC_000913/1820312‑1820795                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  aaaGGAGTCCAGCATCGGCTTGGCCCAGATATGAAAGGTCTTATCGGCTGCCAGCGCCGCTTCGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGCAGCCATaa                                                                                                                                                                                                                                           <  2:240383/251‑1 (MQ=255)     gAGTCCAGCATCGGCTTGGCCCAGATATGGAAGGTCTTACCGGCTGCCAGCGCCGCTTCGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGCATCCATAATGAt                                                                                                                                                                                                                                       <  1:278317/251‑1 (MQ=255)     gAGTCCAGCATCGGCTTGGCCCAGATATGGAAGGTCTTACCGGCTGCCAGCGCCGCTTCGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGCATCCATAATGAt                                                                                                                                                                                                                                       <  1:284028/251‑1 (MQ=255)      aGTCCAGCATCGGCTTGGCCCAGATATGGAAGGTCTTACCGGCTGCCAGCGCCGCTTCGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGCATCCATAATGATc                                                                                                                                                                                                                                      <  2:6632/251‑1 (MQ=255)         ccatcaccGGCTTGGCCAAGATATGGAAAGTCTCACCGGCTGCCAGCGCCGCCCCGCCGGCACCATATTGCTACTAGGTCGCGATAGTTTACATTGCCCACGGCGTCATAATGCCGCTGTCCAGTGCGCTCAGCCCCAGCCAGCCATGCCTACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGACCCCACCCACCACGCTACCCCACGACGCCAGTGGGGTTGAGATCGCATCCATACTGATCTg                                                                                                                                                                                                                                    >  2:348235/8‑250 (MQ=255)                cGGCTTGGCCCAGATATGGAAGGTCTTACCGGCTGCCAGCGCCGCTTCGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTAGCCATCGATGCCAGTGGGGTTGAGATGGCATCCATAATGATCTGAGGGAAGt                                                                                                                                                                                                                            <  2:43313/251‑1 (MQ=255)                 ggCTTGGCCCAGATATGGAAGGTCTTACCGGCTGCCAGCGCCGCTTCGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGCATCCATAATGATCTGATGGAAGtt                                                                                                                                                                                                                           >  2:261305/1‑251 (MQ=255)                       gCCCAGATATGGAAGGTCTTACCGGCTGCCAGCGCCGCTTCGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGCATCCATAATGATCTGATGGAAGTTGGTGcc                                                                                                                                                                                                                     <  1:91929/251‑1 (MQ=255)                        cccAGATATGGAAGGTCTTACCGGCTGCCAGCGCCGCTTCGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGCATCCATAATGATCTGATGGAAGTTGGTGccc                                                                                                                                                                                                                    >  1:219858/1‑251 (MQ=255)                         ccAGATATGGAAGGTCTTACCGGCTGCCAGCGCCGCTTCGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGCATCCATAATGATCTGATGGAAGTTGGTGcccc                                                                                                                                                                                                                   >  1:305234/1‑251 (MQ=255)                               atGGAAGGTCTTACCGGCTGCCAGCGCCGCTTCGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGCATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTa                                                                                                                                                                                                              <  2:144463/250‑1 (MQ=255)                                  gAAGGTCTTACCGGCTGCCAGCGCCGCTTCGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGCATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCa                                                                                                                                                                                                          <  1:314917/251‑1 (MQ=255)                                                cTGCCAGCGCCGCTTCGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGCATCCATACTGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCaat                                                                                                                                                                                             >  2:84119/1‑250 (MQ=255)                                                      gcgcCGCTTCGACGGAACCATATTGCTGATAGGTAGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGCATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATcccc                                                                                                                                                                                      >  1:341697/1‑251 (MQ=255)                                                        gccgcTTCGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGCATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCAt                                                                                                                                                                                    >  2:81924/1‑251 (MQ=255)                                                           gcTTCGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGCATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCAc                                                                                                                                                                                 >  1:179012/1‑251 (MQ=255)                                                            cTTCGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGCATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACg                                                                                                                                                                                >  1:224899/1‑251 (MQ=255)                                                               cGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGCATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGa                                                                                                                                                                              <  1:350870/250‑1 (MQ=255)                                                                                      gTCGCGATCTTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGCATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTa                                                                                                                                                       >  1:118414/1‑250 (MQ=255)                                                                                                      gTGCCCACGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGCATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGACTAATAAAGCCGGGAATTACTGCCGAGAAGGAAcgc                                                                                                                                      >  2:217654/1‑251 (MQ=255)                                                                                                         cccACGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGCATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGAt                                                                                                                                   <  2:296020/251‑1 (MQ=255)                                                                                                                      aaTGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGCATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAGGt                                                                                                                       <  2:179012/250‑1 (MQ=255)                                                                                                                         gcccttctCCAGTGCGGTCAGCCCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCCGTGGGGTTGAGCTGGCATCCATCATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAAAGCCCAGGCAATAGTCCCCATCACGGAAAGCATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGCAGCAGGTACGc                                                                                                                   >  2:238568/8‑251 (MQ=255)                                                                                                                              tGTCCAGTGCGGTCAGCGCCAGCGCGCAATGAATACCGAAGAACCAGCGCAGTGGAACAAAGATCACATAGGCCCAGCCCACAACGCTACCCAACGATGCCAGTGGGGTTGAGATGGCATCAATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCCATAATCCACATCACGGAAAGAATAATAAAGACGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAGGGACGCTGTCg                                                                                                              >  1:297769/1‑251 (MQ=255)                                                                                                                                               gccagcgcGCCATGAATACCGAAGAACCAGAGCAGGGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGCATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAGGTACGCTGTCGGGCAGTTTAATGACCCa                                                                                             <  2:269586/251‑1 (MQ=255)                                                                                                                                                    cgcgcCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGCATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAGGTACGCTGTCGGGCAGTTTAATGACCCAATTgc                                                                                        <  1:149236/251‑1 (MQ=255)                                                                                                                                                            gAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGCATCCATAAt                                                                                                                                                                                                                                          >  2:98494/1‑97 (MQ=255)                                                                                                                                                            gAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGCATCCATAAt                                                                                                                                                                                                                                          <  1:98494/97‑1 (MQ=255)                                                                                                                                                                tccGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGCATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAGGTACGCTGTCGGGCAGTTTAATGACCCAATTGCGGCGGACAATa                                                                             >  2:100850/2‑250 (MQ=255)                                                                                                                                                                       aaCCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGCATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAGGTACGCTGTCGGGCAGTTTAATGACCCAATTGCGGCGGACAATCAAGGTaa                                                                      >  2:145154/1‑250 (MQ=255)                                                                                                                                                                                                                cacGCTACCCAACGATGCCAGTGGGGTTGAGATGGCATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAGGTACGCTGTCGGGCAGTTTAATGACCCAATTGCGGCGGACAATAAAGGTAAACATTTCTGCCACCACCAGGCCAATAATAATCCCTGAgatga                            <  1:324373/251‑1 (MQ=255)                                                                                                                                                                                                                              aTGCCAGTGGGGTTGAGATGGCATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAGGTACGCTGTCGGGCAGTTTAATGACCCAATTGCGGCGGACAATAAAGGTAAACATTTCTGCCACCACCAGGCCAATAATAATCCCGGAGATGATATTCGCCCCACCt              <  2:110784/251‑1 (MQ=255)                                                                                                                                                                                                                                   aGTGGGGTTGAGATGGCATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAGGTACGCTGTCGGGCAGTTTAATGACCCAATTGCGGCGGACAATAAAGGTAAACATTTCTGCCACCACCAGGCCAATAATAATCCCGGAGATGATATTCGCCCCACCTAAcc          >  2:245378/1‑250 (MQ=255)                                                                                                                                                                                                                                           tGAGATGGCATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAGGTACGCTGTCGGGCAGTTTAATGACCCAATTGCGGCGGACAATAAAGGTAAACATTTCTGCCACCACCAGGCCAATAATAATCCCGGAGATGATATTCGCCCCACCTAACCAGTTTGCa  <  1:245378/250‑1 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                  AAAGGAGTCCAGCATCGGCTTGGCCCAGATATGGAAGGTCTTACCGGCTGCCAGCGCCGCTTCGACGGAACCATATTGCTGATAGGTCGCGATATTTTCCAGTGCCCACGGCGTCATAATGCCGTTGTCCAGTGCGGTCAGCGCCAGCGCGCCATGAATACCGAAGAACCAGAGCAGTGGAACAAAGATCACATAGGCCCAGCCCACCACGCTACCCAACGATGCCAGTGGGGTTGAGATGGTATCCATAATGATCTGATGGAAGTTGGTGCCCCAGGTATTCAACGCCCAGGCAATAATCCCCATCACGGAAAGAATAATAAAGCCGGGAATTAATGCCGAGAAGGAACGCGATACTGAAGCAGGTACGCTGTCGGGCAGTTTAATGACCCAATTGCGGCGGACAATAAAGGTAAACATTTCTGCCACCACCAGGCCAATAATAATCCCGGAGATGATATTCGCCCCACCTAACCAGTTTGCA  >  NC_000913/1820312‑1820795

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