Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,396,373	C→T	T103T (ACC→ACT)	amiB →	N‑acetylmuramoyl‑L‑alanine amidase B

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,396,373	0	C	T	100.0%	96.4 / NA	29	T103T (ACC→ACT)	amiB	N‑acetylmuramoyl‑L‑alanine amidase B
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (20/9);  total (20/9)

CAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACCGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGT  >  NC_000913/4396149‑4396601                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                     cAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGa                                                                                                                                                                                                             <  2:191230/250‑1 (MQ=255)   ggTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAAAGGTAAAACCGAAGCGGGGAAgg                                                                                                                                                                                                          >  2:207440/1‑250 (MQ=255)          aaCGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACTGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCTGCAg                                                                                                                                                                                                     <  1:50293/249‑1 (MQ=255)           aCGGTAATCAACAGGCGAGGATAACGTTGAGTTTTATTGGCGATCCTGAATATGCGTTTAGCAATCACAGCAAACGCACCGTGGCGCTCGATATCAAACAAAAGGGCGTGATTCAGGGCATGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGACTAAAGATGCACAAACGCTACGGCTGGAGGTCGATCTTACTGAAAACGGTAAACCCGAAGAGGGGAAGCGGCaaaa                                                                                                                                                                                                   >  2:230415/1‑247 (MQ=255)                          cgcgGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACt                                                                                                                                                                                   >  1:25523/1‑251 (MQ=255)                               aTAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTgtcgt                                                                                                                                                                              >  1:202266/1‑251 (MQ=255)                                   cGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCttt                                                                                                                                                                          >  2:100587/1‑251 (MQ=255)                                        aGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGAt                                                                                                                                                                     <  1:27015/251‑1 (MQ=255)                                           cttATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTa                                                                                                                                                                   >  1:148864/2‑250 (MQ=255)                                              aTTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACg                                                                                                                                                                <  1:150807/250‑1 (MQ=255)                                                    gATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGAtgt                                                                                                                                                         <  1:175765/251‑1 (MQ=255)                                                               tGCGTTTAGCCATCAAAGCAACCGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTAGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGc                                                                                                                                              >  2:205281/1‑251 (MQ=255)                                                                    ttAGCCATCAACGCAAACGCACCGTGGCGCTCGCTATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGcctcct                                                                                                                                         >  2:34676/1‑251 (MQ=255)                                                                                 cAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTgg                                                                                                                            >  2:220105/1‑251 (MQ=255)                                                                                                                                    gCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAAc                                                                         >  1:56890/1‑251 (MQ=255)                                                                                                                                       gttgttgttCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCgg                                                                      >  1:220816/1‑251 (MQ=255)                                                                                                                                          gttgttCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTCGAACGCCTAAAGATGCACAAACGCTACGGATGGTGGTCGATCTTACTGAAAACGGTAAACCCGAAGCGGTGCAGCGGCCGAATGGCAGCAATTAAACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGgcgc                                                                   >  1:25045/1‑251 (MQ=255)                                                                                                                                                          aatCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAAc                                                   >  2:85927/1‑251 (MQ=255)                                                                                                                                                                 tGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGt                                            >  1:48706/1‑251 (MQ=255)                                                                                                                                                                           tCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTCCACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACt                                   <  2:148864/250‑1 (MQ=255)                                                                                                                                                                           tCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACt                                   <  2:219732/250‑1 (MQ=255)                                                                                                                                                                                  ggAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGgtgt                           <  1:34676/251‑1 (MQ=255)                                                                                                                                                                                    aaCGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGtt                          >  1:157605/1‑250 (MQ=255)                                                                                                                                                                                        ccTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCTACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATc                       >  2:172705/1‑249 (MQ=255)                                                                                                                                                                                         cTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATcagc                    >  1:1151/1‑251 (MQ=255)                                                                                                                                                                                          tAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGGAACCGCACTACGGGCGTTAAcagc                    >  2:50025/1‑250 (MQ=255)                                                                                                                                                                                                 gCACAAACGCTACGGCTGGTGGTAGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAACGGCAGAAAGGAAGCAATTACACTGTCGCCTTTACGATTAACGCCGACGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGAACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAACGTAACCGCACTACGGGTGTTATCAGCAGTAATAc            >  1:126730/1‑251 (MQ=255)                                                                                                                                                                                                       aCGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAAc      <  2:1151/251‑1 (MQ=255)                                                                                                                                                                                                           tACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGt  >  1:120985/1‑251 (MQ=255)                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                     CAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACCGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGT  >  NC_000913/4396149‑4396601

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