Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,968,037	A→G	K41R (AAA→AGA)	wecA →	UDP‑GlcNAc:undecaprenylphosphate GlcNAc‑1‑phosphate transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,968,037	0	A	G	100.0%	164.7 / NA	48	K41R (AAA→AGA)	wecA	UDP‑GlcNAc:undecaprenylphosphate GlcNAc‑1‑phosphate transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/26);  total (22/26)

TGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAAACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGTTTACGCAGGGATTTGCTTCACGTTCGGAATTGTCGATTACTATATTCCGCATGCATCTCTCTATCTCGCTTGTGCCGGTGTGCTTGT  >  NC_000913/3967902‑3968166                                                                                                                                        |                                                                                                                                  actgtcagtaaatTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACg                                                                                                                                >  1:625953/12‑139 (MQ=255)         tgcatTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACGTCAc                                                                                                                            >  2:507601/1‑135 (MQ=255)         tgcatTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACGTCAc                                                                                                                            <  1:507601/135‑1 (MQ=255)               gtgAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACGTCACCAGGGATTGa                                                                                                                  >  2:698438/1‑139 (MQ=255)               gtgAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACGTCACCAGGGATTGa                                                                                                                  <  1:333551/139‑1 (MQ=255)               gtgAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACGTCACCAGGGATTGa                                                                                                                  >  2:400999/1‑139 (MQ=255)                tgAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACGTCACCAGGGATTGAt                                                                                                                 >  2:107994/1‑139 (MQ=255)                               tGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACg                                                                                                                                <  1:221274/109‑1 (MQ=255)                               tGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACg                                                                                                                                >  2:221274/1‑109 (MQ=255)                                 aGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACGTCACCAGGGATTGATACCTCTCGTTGGGGGGa                                                                                                >  1:280886/1‑139 (MQ=255)                                  gTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACGTCACCa                                                                                                                          >  2:38927/1‑112 (MQ=255)                                  gTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACGTCACCa                                                                                                                          <  1:38927/112‑1 (MQ=255)                                       gATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGAc                                                                                                                                 <  2:571934/100‑1 (MQ=255)                                       gATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGAc                                                                                                                                 >  1:571934/1‑100 (MQ=255)                                       gATCTCATCAGTATTATTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCgg                                                                                          >  1:73982/1‑139 (MQ=255)                                        aTCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACg                                                                                                                                <  1:97908/100‑1 (MQ=255)                                        aTCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACg                                                                                                                                >  2:97908/1‑100 (MQ=255)                                         tctcATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGtt                                                                                        >  2:497163/1‑139 (MQ=255)                                             atcaGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGTTTACg                                                                                    <  1:437273/139‑1 (MQ=255)                                             atcaGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGTTTACg                                                                                    >  2:64226/1‑139 (MQ=255)                                             atcaGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGTTTACg                                                                                    >  2:437273/1‑139 (MQ=255)                                               caGTATTATTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGTTTACGCa                                                                                  <  2:73982/139‑1 (MQ=255)                                                   aTTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGTTTACGCAGGGa                                                                              <  2:186485/139‑1 (MQ=255)                                                    tttttttATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACGTCACCAGGGATTGATAc                                                                                                               <  1:606884/105‑1 (MQ=255)                                                    tttttttATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACGTCACCAGGGATTGATAc                                                                                                               >  2:606884/1‑105 (MQ=255)                                                       ttttATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGTTTACGCAGGGATTTg                                                                          >  2:324429/1‑139 (MQ=255)                                                          tATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGTTTACGCAGGGATTTGCtt                                                                       <  1:132740/139‑1 (MQ=255)                                                                     ctgtttctgtttTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGTTTACGCAGGGATTTGCTTCACGTTCGGaa                                                            <  1:698438/139‑1 (MQ=255)                                                                      tgtttctgtttTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGTTTACGCAGGGATTTGCTTCACGTTCGGAAt                                                           >  2:409187/1‑139 (MQ=255)                                                                        tttctgtttTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGTTTACGCAGGGATTTGCTTCACGTTCGGAATTg                                                         >  1:525780/1‑139 (MQ=255)                                                                        tttctgtttTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGTTTACGCAGGGATTTGCTTCACGTTCGGAATTg                                                         >  2:725499/1‑139 (MQ=255)                                                                        tttctgtttTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGTTTACGCAGGGATTTGCTTCACGTTCGGAATTg                                                         <  1:400999/139‑1 (MQ=255)                                                                          tctgtttTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGTTTACGCAGGGATTTGCTTCACGTTCGGAATTGTc                                                       <  1:409187/139‑1 (MQ=255)                                                                            tgtttTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGTTTACGCAGGGATTTGCTTCACGTTCGGAATTGTCGa                                                     <  1:366389/139‑1 (MQ=255)                                                                            tgtttTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGTTTACGCAGGGATTTGCTTCACGTTCGGAATTGTCGa                                                     <  2:625953/139‑1 (MQ=255)                                                                                 tttGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGTTTACGCAGGGATTTGCTTCACGTTCGGAATTGTCGATTACt                                                <  2:759630/139‑1 (MQ=255)                                                                                  ttGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACGTCACCAGGGATTGATACCTCTCGTTGGGGGGAtt                                                                                              >  1:372029/1‑92 (MQ=255)                                                                                  ttGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACGTCACCAGGGATTGATACCTCTCGTTGGGGGGAtt                                                                                              <  2:372029/92‑1 (MQ=255)                                                                                   tGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGTTTACGCAGGGATTTGCTTCACGTTCGGAATTGTCGATTACtat                                              <  1:324429/139‑1 (MQ=255)                                                                                          aaGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGTTTACGCAGGGATTTGCTTCACGTTCGGAATTGTCGATTACTATATTCCgc                                       <  2:2285/139‑1 (MQ=255)                                                                                           aGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGTTTACGCAGGGATTTGCTTCACGTTCGGAATTGTCGATTACTATATTCCgca                                      <  1:240330/139‑1 (MQ=255)                                                                                              tGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGTTTACGCAGGGATTTGCTTCACGTTCGGAATTGTCGATTACTATATTCCgcatgc                                   <  2:406363/139‑1 (MQ=255)                                                                                               ggCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGTTTACGCAGGGATTTGCTTCACGTTCGGAATTGTCGATTACTATATTCCgcatgca                                  <  2:174489/139‑1 (MQ=255)                                                                                               ggCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACGTCACCAGGGATTGATACCTCTCGTTGGGCGGATTTCGGTTTACGCAGGGATTTGCTTCACGTTCGGAATTGTCGATTACTATATTCCgcatgca                                  <  2:525780/139‑1 (MQ=255)                                                                                                   aaaaaaGTCGGTTTAGTGGATAAACCAAACTTCCGCAGACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGTTTACGCAGGGATTTGCTTCACGTTCGGAATTGTCGATTACTATATTCCGCATGCAtctc                              <  1:64226/139‑1 (MQ=255)                                                                                                                           ccAAACTTCCAAAGACGTCACCAGGGATTGATACCTCTCGGTGGAGGGATTTCCGTTTACGCAGGTATTTGCTTCACCGTCGGAATTGTCGATTACTATATTCCGCTTGCATCTCTCTGCCTCGCTTGTCCCGGCGTGc      >  2:561046/1‑139 (MQ=255)                                                                                                                               aCTTCCGCAGACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGTTTACGCAGGGATTTGCTTCACGTTCGGAATTGTCGATTACTATATTCCGCATGCATCTCTCTATCTCGCTTGTGCCGGTGTGCTTGt  >  1:375649/1‑139 (MQ=255)                                                                                                                                     gCAGACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGTTTACGCAGGGATTTGCTTCACGTTCGGAATTGTCGATTACTATATTCCgc                                       <  1:444950/96‑1 (MQ=255)                                                                                                                                        |                                                                                                                                  TGAGAGCATGCATTGTGAATTTACTGACAGTGAGTACTGATCTCATCAGTATTTTTTTATTCACGACACTGTTTCTGTTTTTTGCCCGTAAGGTGGCAAAAAAAGTCGGTTTAGTGGATAAACCAAACTTCCGCAAACGTCACCAGGGATTGATACCTCTCGTTGGGGGGATTTCGGTTTACGCAGGGATTTGCTTCACGTTCGGAATTGTCGATTACTATATTCCGCATGCATCTCTCTATCTCGCTTGTGCCGGTGTGCTTGT  >  NC_000913/3967902‑3968166

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