Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,118	G→A	A106A (GCG→GCA)	thrB →	homoserine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,118	0	G	A	100.0%	152.2 / NA	41	A106A (GCG→GCA)	thrB	homoserine kinase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (29/12);  total (29/12)

CCAGTGGCGATGACCCTGGAAAAGAATATGCCGATCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGCGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGA  >  NC_000913/3041‑3192                                                                              |                                                                           ccAGTGGCGATGACCCTGGAAAAGAATATGCCGATCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCACt                                                                          >  2:696/1‑80 (MQ=255)  cAGTGGCGATGACCCTGGAAAAGAATATGCCGATCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCACTGa                                                                        >  2:684/1‑81 (MQ=255)   aGTGGCGATGACCCTGGAAAAGAATATGCCGATCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCACTGAt                                                                       <  1:633/81‑1 (MQ=255)   aGTGGCGATGACCCTGGAAAAGAATATGCCGATCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCACTGAt                                                                       >  2:683/1‑81 (MQ=255)    gTGGCGATGACCCTGGAAAAGAATATGCCGATCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCACTGATg                                                                      >  1:681/1‑81 (MQ=255)      ggCGATGACCCTGGAAAAGAATATGCCGATCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCACTGATGGc                                                                    >  2:677/1‑81 (MQ=255)       gCGATGACCCTGGAAAAGAATATGCCGATCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCACTGATGGCg                                                                   <  1:435570/81‑1 (MQ=255)         gATGACCCTGGAAAAGAATATGCCGATCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCACTGATGGCGat                                                                 <  1:634/81‑1 (MQ=255)             aCCCTGGAAAAGAATATGCCGATCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCACTGATGGCGatgaa                                                              <  1:636/80‑1 (MQ=255)             aCCCTGGAAAAGAATATGCCGATCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCACTGATGGCGatgaa                                                              <  1:635/80‑1 (MQ=255)               ccTGGAAAAGAATATGCCGATCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCACTGATGGCGatgaatg                                                            >  2:654937/1‑80 (MQ=255)                       aGAATATGCCGATCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCACTGATGGCGATGAATGAACACTGCg                                                   >  2:674/1‑81 (MQ=255)                        gAATATGCCGATCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCACTGATGGCGATGAATGAACACTGCgg                                                  >  1:666/1‑81 (MQ=255)                         aaTATGCCGATCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCACTGATGGCGATGAATGAACACTGCGGc                                                 <  1:637/81‑1 (MQ=255)                           tatGCCGATCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCACTGATGGCGATGAATGAACACTGCGGCaa                                               <  2:638/81‑1 (MQ=255)                             tGCCGATCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCACTGATGGCGATGAATGCACACTGCGGCAAGc                                             <  1:639/81‑1 (MQ=255)                              gCCGATCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCACTGATGGCGATGAATGAACACTGCGGCAAGcc                                            >  2:680/1‑81 (MQ=255)                               ccGATCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCACTGATGGCGATGAATGAACACTGCGGCAAGcc                                            <  2:640/80‑1 (MQ=255)                                 gATCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCACTGATGGCGATGAATGAACACTGCGGCAAGCCGCt                                         >  2:686/1‑81 (MQ=255)                                 gATCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCACTGATGGCGATGAATGAACACTGCGGCAAGCCGCt                                         >  2:562885/1‑81 (MQ=255)                                 gATCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCACTGATGGCGATGAATGAACACTGCGGCAAGCCGCt                                         >  2:689/1‑81 (MQ=255)                                   tCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCACTGATGGCGATGAATGAACACTGCGGCAAGCCGCtt                                        >  1:699/1‑80 (MQ=255)                                   tCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCACTGATGGCGATGAATGAACACTGCGGCAAGCCGCtt                                        >  2:700/1‑80 (MQ=255)                                        tCGGGCTTGGGCTCCAGTGCCTGTTCGGTGGTCGCGGCACTGATGGCGATGAATGAACACTGCGGCAAGCCGCTTAATGa                                   >  2:690/1‑80 (MQ=255)                                        tCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCACTGATGGCGATGAATGAACACTGCGGCAAGCCGCTTAATGa                                   >  1:676/1‑80 (MQ=255)                                          gggCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCACTGATGGCGATGAATGAACACTGCGGCAAGCCGCTTAATGaca                                 <  1:641/80‑1 (MQ=255)                                            gCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCACTGATGGCGATGAATGAACACTGCGGCAAGCCGCTTAATGACACTc                              <  1:642/81‑1 (MQ=255)                                             cTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCACTGATGGCGATGAATGAACACTGCGGCAAGCCGCTTAATGACACTCg                             <  1:643/81‑1 (MQ=255)                                              ttAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCACTGATGGCGATGAATGAACACTGCGGCAAGCCGCTTAATGACACTCGt                            >  1:691/1‑81 (MQ=255)                                               tAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCACTGATGGCGATGAATGAACACTGCGGCAAGCCGCTTAATGACACTCGtt                           >  1:675/1‑81 (MQ=255)                                                aGGCTCCAGTGCCTGTTCGGTGGTCGCGGCACTGATGGCGATGAATGAACACTGCGGCAAGCCGCTTAATGACACTCGttt                          >  1:685/1‑81 (MQ=255)                                                aGGCTCCAGTGCCTGTTCGGTGGTCGCGGCACTGATGGCGATGAATGAACACTGCGGCAAGCCGCTTAATGACACTCGttt                          >  2:673/1‑81 (MQ=255)                                                 ggCTCCAGTGCCTGTTCGGTGGTCGCGGCACTGATGGCGATGAATGAACACTGCGGCAAGCCGCTTAATGACACTCGTTTg                         >  2:687/1‑81 (MQ=255)                                                       aGTGCCTGTTCGGTGGTCGCGGCACTGATGGCGATGAATGAACACTGCGGCAAGCCGCTTAATGACACTCGTTTGCTGGCt                   >  1:707/1‑81 (MQ=255)                                                        gTGCCTGTTCGGTGGTCGCGGCACTGATGGCGATGAATGAACACTGCGGCAAGCCGCTTAATGACACTCGTTTGCTGGCtt                  >  1:697/1‑81 (MQ=255)                                                                   gtggtCGCGGCACTGATGGCGATGAATGAACACTGCGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAg       >  1:701/1‑81 (MQ=255)                                                                    tggtCGCGGCACTGATGGCGATGAATGAACACTGCGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGc      >  2:695/1‑81 (MQ=255)                                                                     ggtCGCGGCACTGATGGCGATGAATGAACACTGCGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCt     >  1:692/1‑81 (MQ=255)                                                                      gtCGCGGCACTGATGGCGATGAATGAACACTGCGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTg    >  2:678/1‑81 (MQ=255)                                                                        cgcgGCACTGATGGCGATGAATGAACACTGCGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTgg   >  2:694/1‑80 (MQ=255)                                                                        cgcgGCACTGATGGCGATGAATGAACACTGCGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTgg   >  2:705/1‑80 (MQ=255)                                                                         gcggcaCTGATGGCGATGAATGAACACTGCGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGa  <  1:644/80‑1 (MQ=255)                                                                              |                                                                           CCAGTGGCGATGACCCTGGAAAAGAATATGCCGATCGGTTCGGGCTTAGGCTCCAGTGCCTGTTCGGTGGTCGCGGCGCTGATGGCGATGAATGAACACTGCGGCAAGCCGCTTAATGACACTCGTTTGCTGGCTTTGATGGGCGAGCTGGA  >  NC_000913/3041‑3192

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