Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,913,871	T→C	T444A (ACT→GCT)	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,913,871	0	T	C	100.0%	176.1 / NA	49	T444A (ACT→GCT)	glmU	fused N‑acetyl glucosamine‑1‑phosphate uridyltransferase/glucosamine‑1‑phosphate acetyl transferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (24/25);  total (24/25)

GGGCGCCCCGAGCCTTGTAGAGAGTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGTCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTACGCGTCACAGTTGTACCCGCAGCAATGGTCG‑CGCCTTTGCCTACTGTTACCGGGGCCACCAGCTGAGTGTCGGAACCAACAAACACATCGTCGCCGATAATG  >  NC_000913/3913723‑3914016                                                                                                                                                     |                                                                                                                                                   gggCGCCCCGAGCCTTGTAGAGAGTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGCCt                                                                                                                                                  >  2:349121/1‑151 (MQ=255)        ccGAGCCTTGTAGAGAGTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGCCTGCGGcac                                                                                                                                           <  1:624857/151‑1 (MQ=255)         cGAGCCTTGTAGAGAGTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGCCTGCGGcaca                                                                                                                                          <  2:859896/151‑1 (MQ=255)                    gagagTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGCCTGCGGCACACGGCTGATAGc                                                                                                                               >  1:23534/1‑151 (MQ=255)                     agagTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGCCTGCGGCACACGGCTGATAGCt                                                                                                                              >  1:392894/1‑151 (MQ=255)                       agTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGCCTGCGGAACACGGCTGATAGCTa                                                                                                                             >  2:786215/1‑150 (MQ=255)                        gTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGCCTGCGGCACACGGCTGATAGCTAAt                                                                                                                           <  2:232937/151‑1 (MQ=255)                             gTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGCCTGCGGCACACGGCTGATAGCTAATGCAtt                                                                                                                      <  1:415322/151‑1 (MQ=255)                             gTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGCCTGCGGCACACGGCTGATAGCTAATGCAtt                                                                                                                      <  1:847152/151‑1 (MQ=255)                                aTTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTc                                                                                                                   >  2:752150/1‑151 (MQ=255)                                         gTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACAtt                                                                                                          >  2:510927/1‑151 (MQ=255)                                            aCTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTAcg                                                                                                       <  1:453436/151‑1 (MQ=255)                                                cTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTACGCGTc                                                                                                   >  1:187091/1‑151 (MQ=255)                                                cTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTACGCGTc                                                                                                   >  1:829811/1‑151 (MQ=255)                                                   gtgGGAGGGAGATTATGTTATCTCCTCATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTACGCGTcac                                                                                                 <  1:899219/150‑1 (MQ=255)                                                   gtgGGAGGGAGATTATGTTATCTCCTCATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTACGCGTcac                                                                                                 >  2:662290/1‑150 (MQ=255)                                                     gggagggaGATTATGTTATCTCCTCATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTACGCGTCACAg                                                                                               <  2:20238/150‑1 (MQ=255)                                                           gagaTTATGTTATCTCCTCATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTACGCGTCACAGTTGTAcc                                                                                        <  2:342379/151‑1 (MQ=255)                                                             gaTTATGTTATCTCCTCATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTACGCGTCACAGTTGTACCCg                                                                                      >  2:127298/1‑151 (MQ=255)                                                                 aTGTTATCTCCTCATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTACGCGTCACAGTTGTACCCgcagc                                                                                  <  1:895511/151‑1 (MQ=255)                                                                 aTGTTATCTCCTCATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTACGCGTCACAGTTGTACCCgcag                                                                                   <  1:746747/150‑1 (MQ=255)                                                                   gTTATCTCCTCATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTACGCGTCACAGTTGTACCCgcagca                                                                                 >  1:831602/1‑150 (MQ=255)                                                                       tctcCTCATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTACGCGTCACAGTTGTACCCGCAGCAAt                                                                               <  2:209522/148‑1 (MQ=255)                                                                          cctcATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTACGCGTCACAGTTGTACCCGCAGCAATGGTcg‑c                                                                        <  1:478648/151‑1 (MQ=255)                                                                             cATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTACGCGTCACAGTTGTACCCGCAGCAATGGTCG‑CGcc                                                                     >  1:441497/1‑151 (MQ=255)                                                                                   aTGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTACGCGTCACAGTTGTACCCGCAGCAATGGTCG‑CGCCTTTGcc                                                               >  1:431509/1‑151 (MQ=255)                                                                                     gtgACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTACGCGTCACAGTTGTACCCGCAGCAATGGTCG‑CGCCTTTGCCt                                                              <  2:548352/150‑1 (MQ=255)                                                                                               tAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTACGCGTCACAGTTGTACCCGCAGCAATGGTCG‑CGCCTTTGCCTACTGTTACCg                                                    <  1:637711/150‑1 (MQ=255)                                                                                                        aGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTACGCGTCACAGTTGTACCCGCAGCAATGGTCG‑CGCCTTTGCCTACTGTTACCGGGGccacc                                            >  1:507865/1‑149 (MQ=255)                                                                                                             cACTTTTTCTTTACCGGACGACGCCAACCTTCTTTCTGAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTACGCGTCACAGTTGTACCCGCAGCAATGGTCG‑CGCCTTTGCCTACTGTTACCGGGGCCACCAGCTGAg                                     <  1:325027/151‑1 (MQ=255)                                                                                                              aCTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTACGCGTCACAGTTGTACCCGCAGCAATGGTCG‑CGCCTTTGCCTACTGTTACCGGGGCCACCAGCTGAgt                                    >  1:406247/1‑151 (MQ=255)                                                                                                               cTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTACGCGTCACAGTTGTACCCGCAGCAATGGTCG‑CGCCTTTGCCTACTGTTACCGGGGCCACCAGCTGAgt                                    <  2:756355/150‑1 (MQ=255)                                                                                                                   ttCTTTACCGGACGACGCCAGCCTTCTTTCTGAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTACGCGTCACAGTTGTACCCGCAGCAATGGTCG‑CGCCTTTGCCTACTGTTACCGGGGCCACCAGCTGAGTGTCgg                               >  2:480277/1‑151 (MQ=255)                                                                                                                     cTTTACCGGACGACGCCAGCCTTCTTTCTGAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTCCGCGTCACAGTTGTACCCGCAGCAATGGTCG‑CGCCTTTGCCTACTGTTACCGGGGCCACCAGCTGAGTGTCg                                >  2:305145/1‑148 (MQ=255)                                                                                                                      tttCCCGGCCGACGCAAGCTTTCTTCCGGACCCGGCGCCACCCGCCGGTTGGCTAATGCATTTTCGCCGACATTCCGCGCCACAGTTGTCCCCGCAG‑AAAGGTCGCCGCTTTTGCCTCCTGTT‑CCCGGGCCACCAGCTGAGTGTCGGaa                             <  2:219632/149‑1 (MQ=255)                                                                                                                         aCCGGACGACGCCAGCCTTCTTTCTGAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTACGCGTCACAGTTGTACCCGCAGCAATGGTCG‑CGCCTTTGCCTACTGTTACCGGGGCCACCAGCTGAGTGTCGGAACca                          >  1:19769/1‑150 (MQ=255)                                                                                                                           cGGACGACGCCAGCCTTCTTTCTGAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTACGCGTCACAGTTGTACCCGCAGCAATGGTCG‑CGCCTTTGCCTACTGTTACCGGGGCCACCAGCTGAGTGTCGGAACcaaca                       <  1:516589/151‑1 (MQ=255)                                                                                                                            ggACGACGCCAGCCTTCTTTCTGAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTACGCGTCACAGTTGTACCCGCAGCAATGGTCG‑CGCCTTTGCCTACTGTTACCGGGGCCACCAGCTGAGTGTCGGAACcaaca                       >  2:858241/1‑150 (MQ=255)                                                                                                                             gacgacGCCAGCCTTCTTTCTGAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTACGCGTCACAGTTGTACCCGCAGCAATGGTCG‑CGCCTTTGCCTACTGTTACCGGGGCCACCAGCTGAGTGTCGGAACcaacaa                      >  2:721903/1‑150 (MQ=255)                                                                                                                             gacgacGCCAGCCTTCTTTCTGAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTACGCGTCACAGTTGTACCCGCAGCAATGGTCG‑CGCCTTTGCCTACTGTTACCGGGGCCACCAGCTGAGTGTCGGAACcaacaa                      >  2:96782/1‑150 (MQ=255)                                                                                                                                  cgCCAGCCTTCTTTCTGAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTACGCGTCACAGTTGTACCCGCAGCAATGGTCG‑CGCCTTTGCCTACTGTTACCGGGGCCACCAGCTGAGTGTCGGAACCAACAAACACAt                <  2:109558/151‑1 (MQ=255)                                                                                                                                   gCCAGCCTTCTTTCTGAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTACGCGTCACAGTTGTACCCGCAGCAATGGTCG‑CGCCTTTGCCTACTGTTACCGGGGCCACCAGCTGAGTGTCGGAACCAACAAACACAtc               <  2:572744/151‑1 (MQ=255)                                                                                                                                      aGCCTTCTTTCTGAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTACGCGTCACAGTTGTACCCGCAGCAATGGTCG‑CGCCTTTGCCTACTGTTACCGGGGCCACCAGCTGAGTGTCGGAACCAACAAACACAtcgtc            <  2:406221/151‑1 (MQ=255)                                                                                                                                         cttcttTCTGAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTACGCGTCACAGTTGTACCCGCAGCAATGGTCG‑CGCCTTTGCCTACTGTTACCGGGGCCACCAGCTGAGTGTCGGAACCAACAAACACATCGTCGcc         >  1:542316/1‑151 (MQ=255)                                                                                                                                           tctttctGAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTACGCGTCACAGTTGTACCCGCAGCAATGGTCG‑CGCCTTTGCCTACTGTTACCGGGGCCACCAGCTGAGTGTCGGAACCAACAAACACATCGTCGcc         <  1:527349/149‑1 (MQ=255)                                                                                                                                           tctttctGAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTACGCGTCACAGTTGTACCCGCAGCAATGGTCG‑CGCCTTTGCCTACTGTTACCGGGGCCACCAGCTGAGTGTCGGAACCAACAAACACATCGTCGCCg        >  1:437852/1‑150 (MQ=255)                                                                                                                                              ttctGAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTACGCGTCACAGTTGTACCCGCAGCAATGGTCG‑CGCCTTTGCCTACTGTTACCGGGGCCACCAGCTGAGTGTCGGAACCAACAAACACATCGTCGCCGAt      <  2:392868/149‑1 (MQ=255)                                                                                                                                               tctGAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTACGCGTCACAGTTGTACCCGCAGCAATGGTCG‑CGCCTTTGCCTACTGTTACCGGGGCCACCAGCTGAGTGTCGGAACCAACAAACACATCGTCGCCGATAAt   >  2:837966/1‑151 (MQ=255)                                                                                                                                                  gAGCCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTACGCGTCACAGTTGTACCCGCAGCAATGGTCG‑CGCCTTTGCCTACTGTTACCGGGGCCACCAGCTGAGTGTCGGAACCAACAAACACATCGTCGCCGATAATg  <  1:119320/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GGGCGCCCCGAGCCTTGTAGAGAGTGGGGTTATTTTTATAGTTACTGCTTGTGGGAGGGAGATTATGTTATCTCCTCATCCCATGTGACCGGGTTAGCCGGCCAGAATCACTTTTTCTTTACCGGACGACGCCAGCCTTCTTTCTGAGTCTGCGGCACACGGCTGATAGCTAATGCATTTTCGCCGACATTACGCGTCACAGTTGTACCCGCAGCAATGGTCG‑CGCCTTTGCCTACTGTTACCGGGGCCACCAGCTGAGTGTCGGAACCAACAAACACATCGTCGCCGATAATG  >  NC_000913/3913723‑3914016

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