Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	28,886	A→G	A171A (GCA→GCG)	dapB →	4‑hydroxy‑tetrahydrodipicolinate reductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	28,886	0	A	G	100.0%	99.6 / NA	28	A171A (GCA→GCG)	dapB	4‑hydroxy‑tetrahydrodipicolinate reductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (16/12);  total (16/12)

TGGCGTTAACGTCATGCTTAAGCTGCTGGAGAAAGCAGCCAAAGTGATGGGTGACTACACCGATATCGAAATTATTGAAGCACATCATAGACATAAAGTTGATGCGCCGTCAGGCACCGCACTGGCAATGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTCGTGAAGGCCACACCGGTGAACGTGTGCCTGGCACCATTGGTTTTGCCACCGTGCGTGCAGGTGACAT  >  NC_000913/28766‑29020                                                                                                                         |                                                                                                                                       tGGCGTTAACGTCATGCTTAAGCTGCTGGAGAAAGCAGCCAAAGTGATGGGTGACTACACCGATATCGAAATTATTGAAGCACATCATAGACATAAAGTTGATGCGCCGTCAGGCACCGCGCTGGCAATGGGAGAGGCGa                                                                                                                     <  2:278959/140‑1 (MQ=255) tGGCGTTAACGTCATGCTTAAGCTGCTGGAGAAAGCAGCCAAAGTGATGGGTGACTACACCGATATCGAAATTATTGAAGCACATCATAGACATAAAGTTGATGCGCCGTCAGGCACCGCGCTGGCAATGGGAGAGGCGa                                                                                                                     <  1:46073/140‑1 (MQ=255)  ggCGTTAACGTCATGCTTAAGCTGCTGGAGAAAGCAGCCAAAGTGATGGGTGACTACACCGATATCGAAATTATTGAAGCACATCATAGACATAAAGTTGATGCGCCGTCAGGCACCGCGCTGGCAATGGGAGAGGCGAt                                                                                                                    <  1:110524/140‑1 (MQ=255)        aaCGTCATGCTTAAGCTGCTGGAGAAAGCAGCCAAAGTGATGGGTGACTACACCGATATCGAAATTATTGAAGCACATCATAGACATAAAGTTGATGCGCCGTCAGGCACCGCGCTGGCAATGGGAGAGGCGATCGCCCa                                                                                                              >  2:240410/1‑140 (MQ=255)                                    cagcCAAAGTGATGGGTGACTACACCGATATCGAAATTATTGAAGCACATCATAGACATAAAGTTGATGCGCCGTCAGGCACCGCGCTGGCAATGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTgc                                                                                  >  1:341338/1‑140 (MQ=255)                                     agcCAAAGTGATGGGTGACTACACCGATATCGAAATTATTGAAGCACATCATAGACATAAAGTTGATGCGCCGTCAGGCACCGCGCTGGCAATGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTgcg                                                                                 <  2:341338/140‑1 (MQ=255)                                       ccAAAGTGATGGGTGACTACACCGATATCGAAATTATTGAAGCACATCATAGACATAAAGTTGATGCGCCGTCAGGCACCGCGCTGGCAATGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTgcg                                                                                 >  2:509992/1‑138 (MQ=255)                                       ccAAAGTGATGGGTGACTACACCGATATCGAAATTATTGAAGCACATCATAGACATAAAGTTGATGCGCCGTCAGGCACCGCGCTGGCAATGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTgcg                                                                                 <  1:509992/138‑1 (MQ=255)                                                 gggTGACTACACCGATATCGAAATTATTGAAGCACATCATAGACATAAAGTTGATGCGCCGTCAGGCACCGCGCTGGCAATGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTc                                                                     >  2:406519/1‑140 (MQ=255)                                                              gATATCGAAATTATTGAAGCACATCATAGACATAAAGTTGATGCGCCGTCAGGCACCGCGCTGGCAATGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTCGTGAAGGc                                                             <  1:157614/135‑1 (MQ=255)                                                              gATATCGAAATTATTGAAGCACATCATAGACATAAAGTTGATGCGCCGTCAGGCACCGCGCTGGCAATGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTCGTGAAGGc                                                             >  2:157614/1‑135 (MQ=255)                                                               atatCGAAATTATTGAAGCACATCATAGACATAAAGTTGATGCGCCGTCAGGCACCGCGCTGGCAATGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTCGTg                                                                  <  2:170473/129‑1 (MQ=255)                                                               atatCGAAATTATTGAAGCACATCATAGACATAAAGTTGATGCGCCGTCAGGCACCGCGCTGGCAATGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTCGTg                                                                  >  1:170473/1‑129 (MQ=255)                                                               atatCGAAATTATTGAAGCACATCATAGACATAAAGTTGATGCGCCGTCAGGCACCGCGCTGGCAATGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTCGTGAAGGCCACAcc                                                       >  1:105499/1‑140 (MQ=255)                                                                 atCGAAATTATTGAAGCACATCATAGACATAAAGTTGATGCGCCGTCAGGCACCGCGCTGGCAATGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTCGTGAAGGCCACACCgg                                                     >  2:192917/1‑140 (MQ=255)                                                                                gCACATCATAGACATAAAGTTGATGCGCCGTCAGGCACCGCGCTGGCAATGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTCGTGAAGGCCACACCGGTGAACGTGTGCCTgg                                      >  2:224740/1‑140 (MQ=255)                                                                                gCACATCATAGACATAAAGTTGATGCGCCGTCAGGCACCGCGCTGGCAATGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTCGTGAAGGCCACACCGGTGAACGTGTGCCTgg                                      <  1:406519/140‑1 (MQ=255)                                                                                  acaTCATAGACATAAAGTTGATGCGCCGTCAGGCACCGCGCTGGCAATGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTCGTGAAGGCCACACCGGTGAACGTGTGCCTGGCa                                    >  2:193142/1‑140 (MQ=255)                                                                                   catcatAGACATAAAGTTGATGCGCCGTCAGGCACCGCGCTGGCAATGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTCGTGAAGGCCACACCGGTGAACgt                                              >  1:504870/1‑129 (MQ=255)                                                                                   catcatAGACATAAAGTTGATGCGCCGTCAGGCACCGCGCTGGCAATGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTCGTGAAGGCCACACCGGTGAACgt                                              <  2:504870/129‑1 (MQ=255)                                                                                         aGACATAAAGTTGATGCGCCGTCAGGCACCGCGCTGGCAATGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTCGTGAAGGCCACACCGGTGAACGTGTGCCTGGCACCATTgg                             >  2:211691/1‑140 (MQ=255)                                                                                         aGACATAAAGTTGATGCGCCGTCAGGCACCGCGCTGGCAATGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTCGTGAAGGCCACACCAGTGAACgt                                              <  2:218279/123‑1 (MQ=255)                                                                                         aGACATAAAGTTGATGCGCCGTCAGGCACCGCGCTGGCAATGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTCGTGAAGGCCACACCAGTGAACgt                                              >  1:218279/1‑123 (MQ=255)                                                                                           aCATAAAGTTGATGCGCCGTCATGCACCGCGCTGGCAATGGGAGAGGCGATCGCCCACGCCCTTGATAACGATCTGAAAGATTGCGCGGTCTACAGTCGTGAAGGCCACACCGGTGAACGTGTGCCTGGCACCATTGGtt                           >  2:199842/1‑140 (MQ=255)                                                                                              tAAAGTTGATGCGCCGTCAGGCACCGCGCTGGCAATGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTCGTGAAGGCCACACCGGTGAACGTGTGCCTGGCACCATTGGTTTTg                        >  1:537454/1‑140 (MQ=255)                                                                                              tAAAGTTGATGCGCCGTCAGGCACCGCGCTGGCAATGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTCGTGAAGGCCACACCGGTGAACGTGTGCCTGGCACCATTGGTTTTg                        <  1:211691/140‑1 (MQ=255)                                                                                              tAAAGTTGATGCGCCGTCAGGCACCGCGCTGGCAATGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTCGTGAAGGCCACACCGGTGAACGTGTGCCTGGCACCATTGGTTTTg                        >  1:5864/1‑140 (MQ=255)                                                                                                                    aCCGCGCTGGCAATGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTCGTGAAGGCCACACCGGTGAACGTGTGCCTGGCACCATTGGTTTTGCCACCGTGCGTGCAGGTGACAt  <  2:537454/140‑1 (MQ=255)                                                                                                                      cgcgCTGGCAATGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTCGTGAAGGCCACACCGGTGAACGTGTGCCTGGCACCATTGGTTTTGCCACCGTGCGTGCAGGTGACAt  <  2:448524/138‑1 (MQ=255)                                                                                                                         |                                                                                                                                       TGGCGTTAACGTCATGCTTAAGCTGCTGGAGAAAGCAGCCAAAGTGATGGGTGACTACACCGATATCGAAATTATTGAAGCACATCATAGACATAAAGTTGATGCGCCGTCAGGCACCGCACTGGCAATGGGAGAGGCGATCGCCCACGCCCTTGATAAAGATCTGAAAGATTGCGCGGTCTACAGTCGTGAAGGCCACACCGGTGAACGTGTGCCTGGCACCATTGGTTTTGCCACCGTGCGTGCAGGTGACAT  >  NC_000913/28766‑29020

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