Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,915,089	G→T	H38N (CAT→AAT)	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,915,089	0	G	T	100.0%	195.4 / NA	52	H38N (CAT→AAT)	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 T (21/31);  total (21/31)

ACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATGCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTT  >  NC_000913/3914958‑3915232                                                                                                                                    |                                                                                                                                                aCCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCg                                                                                                                              <  2:515851/151‑1 (MQ=255)  cccAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCg                                                                                                                              <  2:605123/150‑1 (MQ=255)   ccAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGc                                                                                                                            >  2:715080/1‑151 (MQ=255)    cAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGc                                                                                                                            >  2:239966/1‑150 (MQ=255)    cAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCa                                                                                                                           <  2:167399/151‑1 (MQ=255)    cAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCa                                                                                                                           <  2:44634/151‑1 (MQ=255)       ctgctCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAgg                                                                                                                        >  1:526162/1‑151 (MQ=255)         gctCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAgg                                                                                                                        <  2:437314/149‑1 (MQ=255)         gctCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGt                                                                                                                      <  1:595680/151‑1 (MQ=255)              gCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATgca                                                                                                                 <  2:196122/151‑1 (MQ=255)                      cACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTc                                                                                                         <  2:2483/151‑1 (MQ=255)                               aaGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCg                                                                                                <  1:30629/151‑1 (MQ=255)                                aGGTTGTCGTCTTTCGGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCgg                                                                                               <  2:713943/151‑1 (MQ=255)                                  gTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGa                                                                                                      >  2:413064/1‑142 (MQ=255)                                  gTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGaa                                                                                             <  2:780959/151‑1 (MQ=255)                                         tcTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATgc                                                                                      >  2:650205/1‑151 (MQ=255)                                           tttCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATgcgc                                                                                    >  2:425697/1‑151 (MQ=255)                                            ttCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATgcgc                                                                                    <  1:747187/150‑1 (MQ=255)                                              cAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTACGGAAGATCGGAATACATGCGCGTg                                                                                 <  1:602949/151‑1 (MQ=255)                                               aGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGc                                                                                >  2:132367/1‑151 (MQ=255)                                                  gcCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCtt                                                                             >  1:643910/1‑151 (MQ=255)                                                        tttAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCt                                                                        <  2:559030/150‑1 (MQ=255)                                                        tttAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCt                                                                        <  1:435042/150‑1 (MQ=255)                                                         ttAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTg                                                                       <  1:337099/150‑1 (MQ=255)                                                                aTCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAgg                                                               <  1:220493/151‑1 (MQ=255)                                                                aTCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAgg                                                               <  2:18594/151‑1 (MQ=255)                                                                      gccgTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACt                                                         <  1:15024/151‑1 (MQ=255)                                                                       ccgTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACt                                                         <  1:724808/150‑1 (MQ=255)                                                                        cgTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTAc                                                       <  1:593485/151‑1 (MQ=255)                                                                        cgTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTAc                                                       >  2:114446/1‑151 (MQ=255)                                                                         gTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACg                                                      >  1:594191/1‑151 (MQ=255)                                                                            aCCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTc                                                   >  1:605061/1‑151 (MQ=255)                                                                             ccGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCa                                                  <  1:627822/151‑1 (MQ=255)                                                                                     cAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCatt                                           >  1:59394/1‑150 (MQ=255)                                                                                          gAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCaa                                     <  2:353981/151‑1 (MQ=255)                                                                                            aCGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAAc                                    >  1:81585/1‑150 (MQ=255)                                                                                            aCGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAAc                                    >  1:487891/1‑150 (MQ=255)                                                                                            aCGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAAc                                    >  1:87378/1‑150 (MQ=255)                                                                                                              tcattcGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTa                  >  1:715809/1‑150 (MQ=255)                                                                                                               cattcGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTaa                 <  1:66876/150‑1 (MQ=255)                                                                                                               cattcGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAAt                <  1:669523/151‑1 (MQ=255)                                                                                                                 ttcGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCCCGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAAttt              >  2:149679/1‑151 (MQ=255)                                                                                                                  tcGCAGCATCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTg             >  1:119537/1‑151 (MQ=255)                                                                                                                    gcagcaTCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTgag           <  1:725841/151‑1 (MQ=255)                                                                                                                    gcagcaTCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTg             >  1:563533/1‑149 (MQ=255)                                                                                                                     cagcaTCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTgaga          <  1:222519/151‑1 (MQ=255)                                                                                                                     cagcaTCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTg             >  2:452640/1‑148 (MQ=255)                                                                                                                      agcaTCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTgaga          >  1:745549/1‑150 (MQ=255)                                                                                                                       gcaTCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAAc        <  2:278782/151‑1 (MQ=255)                                                                                                                         aTCAATGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGa      <  1:413109/151‑1 (MQ=255)                                                                                                                            aaTGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAAtt   <  1:157858/151‑1 (MQ=255)                                                                                                                             aTGACATTCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAAttt  <  1:423112/151‑1 (MQ=255)                                                                                                                                    |                                                                                                                                                ACCCAGCTGCTCTGCCTGAAGCACCCAGTTAAGGTTGTCGTCTTTCAGCGCCTGTTTTAGCAGATCGCCGCCGTGACCGTACACCAGGTGAACGTGCGCTGCGCCTAATTCATTCGCAGCATCAATGACATGCTGAACCATCGCTTTCCCGGCAAGGGTATGCAGCACTTTCGGAAGATCGGAATACATGCGCGTGCCTTTGCCTGCGGCAAGGATCACTACGCTCATAGCATTATTCAACATACGCGTCCTGACTGTAATTTGAGAACGAATTT  >  NC_000913/3914958‑3915232

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