Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,970,974	A→G	N281S (AAC→AGC)	wecB →	UDP‑N‑acetylglucosamine 2‑epimerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,970,974	0	A	G	100.0%	83.7 / NA	26	N281S (AAC→AGC)	wecB	UDP‑N‑acetylglucosamine 2‑epimerase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/8);  total (18/8)

TCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAACCTGTGCTGGTGATGCGCGATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGTACGGTGCGT  >  NC_000913/3970848‑3971101                                                                                                                               |                                                                                                                                tCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGTCTGGCTGATGAGCCACGc                                                                                                                           <  2:255183/133‑1 (MQ=255) tCCAGATCGTATATCCGGTGCATCGCGACCCGTACGTCAGATAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGTCTGGCTGATGAGCCACGc                                                                                                                           <  2:230224/133‑1 (MQ=255)  ccAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGTCTGGCTGATGAGCCACGcc                                                                                                                          >  1:21265/1‑133 (MQ=255)       ttGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACTGTTTGTCTGGCGGAAGAGCCAAGCCTGGCt                                                                                                                     >  1:159641/1‑133 (MQ=255)           cTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTGGTCTGGCCGATGAGCCACGCCTGGCTGAtt                                                                                                                 >  1:230224/1‑133 (MQ=255)            tATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGTCTGGCTGATGAGCCACGCCTGGCTGAttt                                                                                                                >  2:3833/1‑133 (MQ=255)                            aCCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGTCTGACTGATGAGGCACGCCTAGCTGATTTCGACCTACTCAggcgg                                                                                                >  2:195818/1‑133 (MQ=255)                                                gTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGTCTGGCTGATGAGCCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCtt                                                                            >  2:289501/1‑133 (MQ=255)                                                     tCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGTCTGGCTGATGAGCCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTg                                                                       >  1:312434/1‑133 (MQ=255)                                                       gCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGTCTGGCTGATGAGCCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTggg                                                                     <  1:3833/133‑1 (MQ=255)                                                       gCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGTCTGGCTGATGAGCCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTggg                                                                     <  1:55022/133‑1 (MQ=255)                                                       gCATTCTGAGGCATGGGACAAATGTCATTCAGATCGATCCCCAGGAGTATTTACCGTTTGTCTGGCTGATGAGCCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTggg                                                                     <  1:195818/133‑1 (MQ=255)                                                        cATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGTCTGGCTGATGAGCCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTgggg                                                                    >  2:89354/1‑133 (MQ=255)                                                             tGGGGCATGTTAAAAGTGTCATTCTGATTGATCCCCAGGAGCATTTACCGTTTGTCTGGCTGATGAGCCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAAcc                                                               <  2:159641/133‑1 (MQ=255)                                                               gggCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGTCTGGCTGATGAAACTCGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAATCCtg                                                             <  2:295067/133‑1 (MQ=255)                                                                            aTGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGTCTGGCTGATGAGCCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAACCTGTGCTGGTGATgcg                                                <  1:89354/133‑1 (MQ=255)                                                                                       atcgatcCCCAGGAGTATTTACCGTTTGTCTGGCTGATGAGCCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAACCTGTGCTGGTGATGCGCGATACCACTg                                     >  2:51999/1‑133 (MQ=255)                                                                                                   gAGTATTTACCGTTTGTCTGGCTGATGAGCCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAACCTGTGCTGGTGATGCGCGATACCACTGAGCGTCCGGAAg                         >  1:332759/1‑133 (MQ=255)                                                                                                   gAGTATTTACCGTTTGTCTGGCTGATGAGCCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAACCTGTGCTGGTGATGCCCGATACCACTGAGCGTCCGGAAg                         >  2:51387/1‑133 (MQ=255)                                                                                                      tATTTACCGTTTGTCTGGCTGATGAGCCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAACCTGTGCTGGTGATGCGCGATACCACTGAGCGTCCGGAAGCgg                      >  1:116907/1‑133 (MQ=255)                                                                                                       aTTTACCGTTTGTCTGGCTGATGAGCCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAACCTGTGCTGGTGATGCGCGATACCACTGAGCGTCCGGAAGCGGt                     <  2:312434/133‑1 (MQ=255)                                                                                                        tttACCGTTTGTCTGGCTGATGAGCCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAACCTGTGCTGGTGATGCGCGATATCACTGAACGTCCGGAAGCGGTg                    >  1:153794/1‑133 (MQ=255)                                                                                                             cGTTTGTCTGGCTGATGAGCCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAACCTGTGCTGGTGATGCGTGAGCCTATTGAGCGTCTGGAACCGGTGAcgac               >  2:237439/1‑131 (MQ=255)                                                                                                             cGTTTGTCTGGCTGATGAGCCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAACCTGTGCTGGTGATGCGCGATACCACTGAGCGTCCGGAAGCGGTGAcggc               >  1:500356/1‑133 (MQ=255)                                                                                                             cGTTTGTCTGGCTGATGAGCCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAACCTGTGCTGGTGATGCGCGATACCACTCAGCCTCAGGAAGCGGTGAcggc               >  1:189118/1‑133 (MQ=255)                                                                                                              gTTTGTCTGGCTGATGAGCCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAACCTGTGCTGGTGATGCGCGAAACCACTGAGCGACCAGAAGCGCTGAcggcg              >  1:245230/1‑133 (MQ=255)                                                                                                                          gatgaGCCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAACCTGTGCTGGTGATGCGCGATACCACTGAGAGTCCGAATCCGGTGACGGCGGCAACGGTGccg  >  2:242340/1‑131 (MQ=255)                                                                                                                               |                                                                                                                                TCCAGATTGTCTATCCGGTGCATCTCAACCCGAACGTCAGAGAACCGGTCAATCGCATTCTGGGGCATGTGAAAAATGTCATTCTGATCGATCCCCAGGAGTATTTACCGTTTGTCTGGCTGATGAACCACGCCTGGCTGATTTTGACCGACTCAGGCGGCATTCAGGAAGAAGCGCCTTCGCTGGGGAAACCTGTGCTGGTGATGCGCGATACCACTGAGCGTCCGGAAGCGGTGACGGCGGGTACGGTGCGT  >  NC_000913/3970848‑3971101

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