Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	553,146	A→G	intergenic (‑46/+72)	purE ← / ← lpxH	N(5)‑carboxyaminoimidazole ribonucleotide mutase/UDP‑2,3‑diacylglucosamine diphosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	553,146	0	A	G	100.0%	113.1 / NA	33	intergenic (‑46/+72)	purE/lpxH	N(5)‑carboxyaminoimidazole ribonucleotide mutase/UDP‑2,3‑diacylglucosamine diphosphatase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (23/10);  total (23/10)

TTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTAT  >  NC_000913/553017‑553282                                                                                                                                  |                                                                                                                                         ttCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCGCAGAGAATAg                                                                                                                                >  1:889230/1‑140 (MQ=255)   cggcggCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCGCAGAGAATAGCa                                                                                                                              <  2:889230/140‑1 (MQ=255)         cgAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCGCAgaga                                                                                                                                    >  2:857784/1‑128 (MQ=255)         cgAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCGCAgaga                                                                                                                                    <  1:857784/128‑1 (MQ=255)           aaCTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGcgca                                                                                                                                        >  2:368068/1‑122 (MQ=255)                     tAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCGCAGAGAATAGCACGGAAAGAGAGCAAGGaa                                                                                                            >  1:1049332/1‑140 (MQ=255)                        cccAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCGCAGAGAATAGCACGGAAAGAGAGCAAGGAAAAc                                                                                                         >  2:1015639/1‑140 (MQ=255)                           aGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCGCAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGt                                                                                                      >  2:694686/1‑140 (MQ=255)                             tCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCGCAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTg                                                                                                    >  1:783066/1‑140 (MQ=255)                               gCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCGCAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCg                                                                                                  >  2:296338/1‑140 (MQ=255)                                  tttGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCGCAGAGAATAGCACGGAAAGAGAGCAAGGaaa                                                                                                           >  2:42789/1‑128 (MQ=255)                                  tttGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCGCAGAGAATAGCACGGAAAGAGAGCAAGGaaa                                                                                                           <  1:42789/128‑1 (MQ=255)                                         cccATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCGCAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGaa                                                                                        >  1:354534/1‑140 (MQ=255)                                             tCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCTTGTGGTGCACAACTCTCGGCTTTAGAGGGCGCAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTg                                                                                                <  2:944276/128‑1 (MQ=255)                                             tCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCGCAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTg                                                                                                >  1:944276/1‑128 (MQ=255)                                                                      attGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCGCAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACg                                                           <  1:694686/140‑1 (MQ=255)                                                                             gAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCGCAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATg                                                    >  1:196686/1‑140 (MQ=255)                                                                             gAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCGCAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATg                                                    >  1:572128/1‑140 (MQ=255)                                                                               aGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCGCAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGaa                                                  >  2:280605/1‑140 (MQ=255)                                                                                gACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCGCAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAAt                                                 <  2:783066/140‑1 (MQ=255)                                                                                         cTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCGCAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAAc                                        >  2:522176/1‑140 (MQ=255)                                                                                           tAAAACTCCTGTGGCGCACAACTCTCGGCTTTAGAGGGCGCAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGt                                      >  2:491680/1‑140 (MQ=255)                                                                                                   cTGTGGTGCACAACTCTCGGCTTTAGAGGGCGCAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATccgcc                              >  2:344262/1‑140 (MQ=255)                                                                                                   cTGTGGTGCACAACTCTCGGCTTTAGAGGGCGCAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATccgcc                              >  1:683162/1‑140 (MQ=255)                                                                                                   cTGTGGTGCACAACTCTCGGCTTTAGAGGGCGCAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATccgcc                              >  1:6376/1‑140 (MQ=255)                                                                                                      tggtgCACAACTCTCGGCTTTAGAGGGCGCAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTg                           >  2:699461/1‑140 (MQ=255)                                                                                                           cacaACTCTCGGCTTTAGAGGGCGCAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACttt                      <  1:280605/140‑1 (MQ=255)                                                                                                           cacaACTCTCGGCTTTAGAGGGCGCAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACttt                      >  2:924561/1‑140 (MQ=255)                                                                                                           cacaACTCTCGGCTTTAGAGGGCGCAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACttt                      <  1:344262/140‑1 (MQ=255)                                                                                                           cacaACTCTCGGCTTTAGAGGGCGCAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACttt                      >  1:348188/1‑140 (MQ=255)                                                                                                                     ggCTTTAGAGGGCGCAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGaa            >  2:611475/1‑140 (MQ=255)                                                                                                                     ggCTTTAGAGGGCGCAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGaa            >  2:274032/1‑140 (MQ=255)                                                                                                                       cTTTAGAGGGCGCAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAAcc          <  2:683162/140‑1 (MQ=255)                                                                                                                               ggCGCAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTAt  <  2:1049332/140‑1 (MQ=255)                                                                                                                                  |                                                                                                                                         TTCGGCGGCGAACTGCATGGTAGCCCAGTCGCTTTTGGACCCCATCACGATGGCGACACGCGCCGGATTATTGCGGGAAGACATGCGTCTTAAAACTCCTGTGGTGCACAACTCTCGGCTTTAGAGGGCACAGAGAATAGCACGGAAAGAGAGCAAGGAAAACGGTTGCGTGGCTGTGAAATCAGCAAAGTTGCGGGTTTTTTAAAACGGAAAATGAATCAGCTCAACGTCATCCGCCGTGACTTTCACCATTGAACCTTCCGTAT  >  NC_000913/553017‑553282

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