Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,914,638	G→C	T188S (ACC→AGC)	glmU ←	fused N‑acetyl glucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,914,638	0	G	C	100.0%	159.7 / NA	43	T188S (ACC→AGC)	glmU	fused N‑acetyl glucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyl transferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (23/20);  total (23/20)

CCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCACGGGTGATACGTC  >  NC_000913/3914495‑3914782                                                                                                                                                |                                                                                                                                                 ccAGACGGGAGAGTTGCAGGCGGTTAGTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGCTCAGCtt                                                                                                                                           <  1:405269/151‑1 (MQ=255)    gACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGCTCAGCTTCgc                                                                                                                                        <  2:137622/151‑1 (MQ=255)           gagTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGGGATGTAGTATTTGCCCTGAGCATTATTGTTGCTCAGCTTCgccagcc                                                                                                                                   >  2:331443/1‑149 (MQ=255)                         tATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGCTCAGCTTCGCCAGCCAGCGTTTCATATCTg                                                                                                                    <  2:592412/150‑1 (MQ=255)                                gCCTTCTACTTCGCTTAAACGGTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGAAACGCCAGCGCAATAATGTAGGTGATGTAGTATTCCCCCTGAGCATTATTGTTGCTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGtt                                                                                                             >  1:487342/1‑150 (MQ=255)                                       aCTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGCTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCaa                                                                                                        <  2:506187/148‑1 (MQ=255)                                       aCTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGCTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAAt                                                                                                       <  1:647951/149‑1 (MQ=255)                                          tCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGCTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGa                                                                                                   >  1:266940/1‑150 (MQ=255)                                           cGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGCTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGaa                                                                                                  <  2:376484/150‑1 (MQ=255)                                           cGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGCTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAAt                                                                                                 >  1:357544/1‑151 (MQ=255)                                                       gCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGCTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGAt                                                                                     <  2:106988/151‑1 (MQ=255)                                                        cGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGCTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGAtc                                                                                    <  1:155375/151‑1 (MQ=255)                                                          gATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGCTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGAtctc                                                                                  >  2:544398/1‑151 (MQ=255)                                                          gATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTACTCGCCCTGAGCATTATTGTTGCTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGAtctc                                                                                  >  2:488258/1‑151 (MQ=255)                                                           aTGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGCTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTcc                                                                                 >  1:470787/1‑151 (MQ=255)                                                               aCGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGCTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCTGGTGTTGATCTCCTGaa                                                                             <  2:666283/151‑1 (MQ=255)                                                               aCGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGCTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGa                                                                              >  2:483917/1‑150 (MQ=255)                                                                  gCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGCTCAGCTTCGCCGGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCt                                                                          <  1:597672/151‑1 (MQ=255)                                                                          ttCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGCTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACgctgc                                                                  >  2:179059/1‑151 (MQ=255)                                                                            cACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGCTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCtc                                                                <  2:324292/151‑1 (MQ=255)                                                                              cGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGCTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCtcgt                                                              >  2:232872/1‑151 (MQ=255)                                                                                      cTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGCTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGc                                                       >  2:249535/1‑150 (MQ=255)                                                                                        gATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGCTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCa                                                      <  2:343032/149‑1 (MQ=255)                                                                                        gATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGCTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCa                                                      <  1:483915/149‑1 (MQ=255)                                                                                             gCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGCTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTtg                                                >  1:264468/1‑150 (MQ=255)                                                                                                  cgcAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGCTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTc                                           >  1:469911/1‑150 (MQ=255)                                                                                                   gcAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGCTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTcaa                                         <  2:236142/151‑1 (MQ=255)                                                                                                      ataatGTCGGTTATGTAGTATTCGCCCTGAGCATTATTGTTGCTCAGGTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTcaacaa                                      <  2:566389/151‑1 (MQ=255)                                                                                                      ataatGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGCTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTcaacaa                                      <  1:217257/151‑1 (MQ=255)                                                                                                       taatGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGCTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAAt                                     >  1:80899/1‑151 (MQ=255)                                                                                                       taatGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGCTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAAt                                     >  2:84359/1‑151 (MQ=255)                                                                                                        aatGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGCTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATg                                    <  2:251580/151‑1 (MQ=255)                                                                                                        aatGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGCTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATg                                    >  2:10248/1‑151 (MQ=255)                                                                                                        aatGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGCTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGGGGCATCTTTGTGCTCAACAATg                                    >  1:69333/1‑151 (MQ=255)                                                                                                             cGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGCTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGt                               >  2:339973/1‑151 (MQ=255)                                                                                                                  aTGTAGTATTCGCCCTGAGCATTATTGTTGCTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACtt                          >  2:200776/1‑151 (MQ=255)                                                                                                                   tGTAGTATTCGCCCTGAGCATTATTGTTGCTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACttt                         >  1:26996/1‑151 (MQ=255)                                                                                                                      agtaTTCGCCCTGAGCATTATTGTTGCTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGcc                      >  2:638593/1‑151 (MQ=255)                                                                                                                                cTGAGCATTATTGTTGCTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCACgg            >  1:19702/1‑151 (MQ=255)                                                                                                                                     cATTATTGTTGCTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCACGGGTGAt       <  2:574967/151‑1 (MQ=255)                                                                                                                                     cATTATTGTTGCTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCACGGGTGAt       >  1:504878/1‑151 (MQ=255)                                                                                                                                      attattGTTGCTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCACGGGTGATa      <  1:232872/151‑1 (MQ=255)                                                                                                                                           tgttgCTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCACGGGTGATACGTc  <  2:172611/150‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                 CCAGACGGGAGAGTTGCAGGCGGTTATTCACGCCTTCTACTTCGCTTAAACGTTGCGGATGAACGGCGACGATTTCACGCCCTTCCTGATACGCCAGCGCAATAATGTCGGTGATGTAGTATTCGCCCTGAGCATTATTGTTGGTCAGCTTCGCCAGCCAGCGTTTCATATCTGCGCCGTTGGCAATCAGAATGCCGGTGTTGATCTCCTGAATCTGACGCTGCTCGTCGGTGGCATCTTTGTGCTCAACAATGCCGGTAACTTTGCCGTTTTCACGGGTGATACGTC  >  NC_000913/3914495‑3914782

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