Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,128,696	G→T	A9S (GCA→TCA)	metB →	O‑succinylhomoserine(thiol)‑lyase/O‑ succinylhomoserine lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,128,696	0	G	T	100.0%	136.7 / NA	39	A9S (GCA→TCA)	metB	O‑succinylhomoserine(thiol)‑lyase/O‑ succinylhomoserine lyase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (21/18);  total (21/18)

TGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATA  >  NC_000913/4128619‑4128771                                                                              |                                                                            tGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCTCAg                                                                          >  1:769078/1‑81 (MQ=255)   tttACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAGCTTCATCACATGACGCGTAAACAGGCCACCATCTCAGt                                                                         <  1:769024/80‑1 (MQ=255)      aCTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCTCAGTGCg                                                                      <  1:769025/80‑1 (MQ=255)      aCTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCTCAGTGCg                                                                      <  2:769025/80‑1 (MQ=255)       ctctGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCTCAGTGCGt                                                                     <  2:769026/80‑1 (MQ=255)         ctGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCTCAGTGCGTAGc                                                                  >  2:769044/1‑81 (MQ=255)          tGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCTCAGTGCGTAg                                                                   <  1:769028/79‑1 (MQ=255)            gTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCTCAGTGCGTAGCggg                                                               <  1:769029/81‑1 (MQ=255)             tGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCTCAGTGCGTAGCGGGt                                                              >  2:769066/1‑81 (MQ=255)               ccTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCTCAGTGCGTAGCGGGtt                                                             <  1:769030/80‑1 (MQ=255)               ccTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCTCAGTGCGTAGCGGGtt                                                             >  2:769061/1‑80 (MQ=255)                      tttCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCTCAGTTCGTAGCGGGTTAAATgac                                                      >  1:769072/1‑80 (MQ=255)                      tttCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCTCAGTGCGTAGCGGGTTAAATgac                                                      >  1:769069/1‑80 (MQ=255)                        tCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCTCAGTGCGTAGCGGGTTAAATgacgac                                                   >  2:769081/1‑81 (MQ=255)                        tCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCTCAGTGCGTAGCGGGTTAAATgacgac                                                   >  2:769085/1‑81 (MQ=255)                        tCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCTCAGTGCGTAGCGGGTTAAATgacgac                                                   >  1:769065/1‑81 (MQ=255)                          aCCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCTCAGTGCGTAGCGGGTTAAATGacgacg                                                  >  2:769065/1‑80 (MQ=255)                              aCAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCTCAGTGCGTAGCGGGTTAAATGACGACGAACAg                                             >  2:769082/1‑81 (MQ=255)                               cAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCTCAGTGCGTAGCGGGTTAAATGACGACGAACAGt                                            <  1:769031/81‑1 (MQ=255)                                aaaGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCTCAGTGCGTAGCGGGTTAAATGACGACGAACAGTa                                           >  2:769076/1‑81 (MQ=255)                                  aGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCTCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATg                                         <  1:769032/81‑1 (MQ=255)                                        gggAACTTCATCACATGACGCGTAAACAGGCCACCATCTCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCg                                   <  1:769033/81‑1 (MQ=255)                                         ggAACTTCATCACATGACGCGTAAACAGGCCACCATCTCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGt                                  >  1:769076/1‑81 (MQ=255)                                          gAACTTCATCACATGACGCGTAAACAGGCCACCATCTCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGtt                                 <  2:769033/81‑1 (MQ=255)                                           aaCTTCATCACATGACGCGTAAACAGGCCACCATCTCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTg                                >  2:769072/1‑81 (MQ=255)                                            aCTTCATCACATGACGCGTAAACAGGCCACCATCTCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGt                               >  1:769080/1‑81 (MQ=255)                                             cTTCATCACATGACGCGTAAACAGGCCACCATCTCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTc                              <  1:769035/81‑1 (MQ=255)                                                catcaCATGACGCGTAAACAGGCCACCATCTCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTccc                            <  1:769036/80‑1 (MQ=255)                                                  tcaCATGACGCGTAAACAGGCCACCATCTCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCAcc                         <  2:769036/81‑1 (MQ=255)                                                     caTGACGCGTAAACAGGCCACCATCTCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGAt                      <  1:769038/81‑1 (MQ=255)                                                          cgcgTAAACAGGCCACCATCTCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATc                 >  2:769080/1‑81 (MQ=255)                                                            cgTAAACAGGCCACCATCTCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCtt               <  2:769038/81‑1 (MQ=255)                                                             gTAAACAGGCCACCATCTCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCttt              >  2:769091/1‑81 (MQ=255)                                                             gTAAACAGGCCACCATCTCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCttt              >  1:769095/1‑81 (MQ=255)                                                               aaaCAGGCCACCATCTCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTc             >  2:769073/1‑80 (MQ=255)                                                                   aGGCCACCATCTCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCa        >  1:769090/1‑81 (MQ=255)                                                                        accaTCTCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCtat   <  2:769039/81‑1 (MQ=255)                                                                        accaTCTCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCta    >  1:769100/1‑80 (MQ=255)                                                                         ccaTCTCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCtata  <  2:769040/81‑1 (MQ=255)                                                                              |                                                                            TGTTTACTCTGGTGCCTGACATTTCACCGACAAAGCCCAGGGAACTTCATCACATGACGCGTAAACAGGCCACCATCGCAGTGCGTAGCGGGTTAAATGACGACGAACAGTATGGTTGCGTTGTCCCACCGATCCATCTTTCCAGCACCTATA  >  NC_000913/4128619‑4128771

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