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%	79.4 / NA	25	T103T (ACC→ACT)	amiB	N‑acetylmuramoyl‑L‑alanine amidase B
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (14/11);  total (14/11)

CGCGACGCTCTCTGATATTCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACCGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAG  >  NC_000913/4396130‑4396608                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            cgcgACGCGCTCTGATATTCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATAACAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACg                                                                                                                                                                                                                                      <  1:212218/251‑1 (MQ=255)            ctGATATTCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAAcc                                                                                                                                                                                                                              <  1:104265/248‑1 (MQ=255)            ctGATATTCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGa                                                                                                                                                                                                                            >  1:162546/1‑250 (MQ=255)                                                 gATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGATTACCCATCAAAGCAAACGCACCGTTGCGCTCGATATCAAACAAACTGTCGTGATTCAGGGACGGTCGTTGTGGTTCAGCTGCAATAAGCGGGTGAAGGCGATTCGCTCTGGAACTCCTAAAGATGCACAAACGCTACGGCTGGTGGTCTATCTTACTGAAAACTGTAAAACCGAAGCGGTAAAGCGGCAGAAGGGCAGCAATTACACTttcg                                                                                                                                                                                      <  2:131474/251‑1 (MQ=255)                                                                            tGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTAGCTCTGGAACGCCTAACGATGCACAAACGCTACGGCTGGAGGTCAATCTTACTGAAAACGGTAAAACAGAAGCGGTGAAGCGGCAGAATGGACGAAATTACACTGTCGTCTTTACGATTAACGCCGATGTgccgc                                                                                                                                                           >  1:3950/1‑251 (MQ=255)                                                                                 aTGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAAAGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCg                                                                                                                                                      >  1:135206/1‑251 (MQ=255)                                                                                                     aaaCGCACCGTGGCGCTCGATATCAAAAAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGt                                                                                                                                  >  2:13819/1‑251 (MQ=255)                                                                                                       aCGCACCGTGGCGCTCGATATCAAAAAAACGGGCGGGATTCAGGGACTGCCGTTGTTGGTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCGCCGCCGCCCGTGGTTg                                                                                                                                <  2:221411/251‑1 (MQ=255)                                                                                                              gTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGGGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGAAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAAcg                                                                                                                         >  2:117550/1‑251 (MQ=255)                                                                                                                                            aTTCAGGGACTGACGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACcgcg                                                                                           >  2:215213/1‑251 (MQ=255)                                                                                                                                             ttCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGAGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACcgcgc                                                                                          >  1:206224/1‑251 (MQ=255)                                                                                                                                                        ccGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAATATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAAcc                                                                               <  1:178237/251‑1 (MQ=255)                                                                                                                                                                    aGCGGCAATAATCTGGTGAAGGCGATGCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGGCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATcc                                                                   <  1:102797/251‑1 (MQ=255)                                                                                                                                                                            taatCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTaaa                                                            >  1:148613/1‑250 (MQ=255)                                                                                                                                                                                  ggTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGGTACGGCTGGTGTTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTTTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTTGTTGCGAAACGCGTTGAAACGTCTGCGGGTGTCGCACTGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACg                                                         <  1:38634/247‑1 (MQ=255)                                                                                                                                                                                  ggTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGa                                                       <  1:117550/249‑1 (MQ=255)                                                                                                                                                                                       aGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAAc                                                >  1:100389/1‑251 (MQ=255)                                                                                                                                                                                        ggCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGGGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGAGAAGCGGCAGAATGGCAGCCATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGCAACGCGTTGAAACGCCCGCGGTAGTCGCCCCGAGCGTCAGCGCACCGGGGCGCAATCCGTTTAAACCGGAAAGTAAcc                                               >  2:82112/1‑251 (MQ=255)                                                                                                                                                                                          cGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGc                                             >  1:44709/1‑251 (MQ=255)                                                                                                                                                                                                gCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGTTGAAGCGGCAGAATGGCAGCAATGACACTTTCGGCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCGTCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCTCAATCCGTTTAAAACGGAAAGTAACCGCACTAc                                        <  2:130598/250‑1 (MQ=255)                                                                                                                                                                                                gCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTAc                                        >  1:130598/1‑250 (MQ=255)                                                                                                                                                                                                   ctGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACgg                                      <  1:73407/249‑1 (MQ=255)                                                                                                                                                                                                                                   cTGGTGGTTGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGc    >  2:58529/1‑251 (MQ=255)                                                                                                                                                                                                                                    tggtggTCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGc    >  1:49010/1‑250 (MQ=255)                                                                                                                                                                                                                                     ggtggtCGATCTTACTGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAg  <  2:44901/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            CGCGACGCTCTCTGATATTCAGGTTTCTAACGGTAATCAACAGGCGCGGATAACGTTGAGTTTTATTGGCGATCCTGATTATGCGTTTAGCCATCAAAGCAAACGCACCGTGGCGCTCGATATCAAACAAACGGGCGTGATTCAGGGACTGCCGTTGTTGTTCAGCGGCAATAATCTGGTGAAGGCGATTCGCTCTGGAACGCCTAAAGATGCACAAACGCTACGGCTGGTGGTCGATCTTACCGAAAACGGTAAAACCGAAGCGGTGAAGCGGCAGAATGGCAGCAATTACACTGTCGTCTTTACGATTAACGCCGATGTGCCGCCACCGCCTCCTCCGCCGCCCGTGGTTGCGAAACGCGTTGAAACGCCTGCGGTTGTCGCACCGCGCGTCAGCGAACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAG  >  NC_000913/4396130‑4396608

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