Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,131,717	G→A	G628E (GGG→GAG)	metL →	fused aspartate kinase/homoserine dehydrogenase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,131,717	0	G	A	100.0%	174.8 / NA	50	G628E (GGG→GAG)	metL	fused aspartate kinase/homoserine dehydrogenase 2
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (30/20);  total (30/20)

CGCCAGATCCACGACGCCTTCGAAAAAACCGGGCGTCACTGGCTGTACAATGCCACCGTCGGTGCGGGCTTGCCGATCAACCACACCGTGCGCGATCTGATCGACAGCGGCGATACTATTTTGTCGATCAGCGGGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACGGTAGCGTGCCGTTTACCGAGCTGGTGGATCAGGCGTGGCAGCAGGGCTTAACCGAACCTGACCCGCGTGACGATCTCTCTGGCAAAGACG  >  NC_000913/4131584‑4131851                                                                                                                                      |                                                                                                                                       cGCCAGATCCACGACGCCTTCGAAAAAACCGGGCGTCACTGGCTGTACAATGCCACCGTCGGTGCGGGCTTGCCGATCAACCACACCGTGCGCGATCTGATCGACAGCGGCGATACTATTTTGTCGATCAGCGAGAtctt                                                                                                                                  >  1:107076/1‑140 (MQ=255) cGCCAGATCCACGACGCCTTCGAAAAAACCGGGCGTCACTGGCTGTACAATGCCACCGTCGGTGCGGGCTTGCCGATCAACCACACCGTGCGCGATCTGATCGACAGCGGCGATACTATTTTGTAGATCAGCGAGAtctt                                                                                                                                  >  1:747141/1‑140 (MQ=255)           acgacgCCTTCGAAAAAACCGGGCGTCACTGGCTGTACAATGCCACCGTCGGTGCGGGCTTGCCGATCAACCACACCGTGCGCGATCTGATCGACAGCGGCGATACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACg                                                                                                                        <  2:491400/140‑1 (MQ=255)           acgacgCCTTCGAAAAAACCGGGCGTCACTGGCTGTACAATGCCACCGTCGGTGCGGGCTTGCCGATCAACCACACCGTGCGCGATCTGATCGACAGCGGCGATACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACg                                                                                                                        <  2:355359/140‑1 (MQ=255)           acgacgCCTTCGAAAAAACCGGGCGTCACTGGCTGTACAATGCCACCGTCGGTGCGGGCTTGCCGATCAACCACACCGTGCGCGATCTGATCGACAGCGGCGATACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACg                                                                                                                        <  1:520010/140‑1 (MQ=255)                 ccTTCGAAAAAACCGGGCGTCACTGGCTGTACAATGCCACCGTCGGTGCGGGCTTGCCGATCAACCACACCGTGCGCGATCTGATCGACAGCGGCGATACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACGCTCTcc                                                                                                                  >  2:517426/1‑140 (MQ=255)                 ccTTCGAAAAAACCGGGCGTCACTGGCTGTACAATGCCACCGTCGGTGCGGGCTTGCCGATCAACCACACCGTGCGCGATCTGATCGACAGCGGCGATACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACGCTCTcc                                                                                                                  >  2:415012/1‑140 (MQ=255)                     cGAAAAAACCGGGCGTCACTGGCTGTACAATGCCACCGTCGGTGCGGGCTTGCCGATCAACCACACCGTGCGCGATCTGATCGACAGCGGCGATACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGc                                                                                                              <  2:790867/140‑1 (MQ=255)                     cGAAAAAACCGGGCGTCACTGGCTGTACAATGCCACCGTCGGTGCGGGCTTGCCGATCAACCACACCGTGCGCGATCTGATCGACAGCGGCGATACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGc                                                                                                              >  2:36545/1‑140 (MQ=255)                        aaaaaCCGGGCGTCACTGGCTGTACAATGCCACCGTCGGTGCGGGCTTGCCGATCAACCACACCGTGCGCGATCTGATCGACAGCGGCGATACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGt                                                                                                           <  2:179396/140‑1 (MQ=255)                        aaaaaCCGGGCGTCACTGGCTGTACAATGCCACCGTCGGTGCGGGCTTGCCGATCAACCACACCGTGCGCGATCTGATCGACAGCGGCGATACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGt                                                                                                           <  1:36545/140‑1 (MQ=255)                               gggCGTCACTGGCTGTACAATGCCACCGTCGGTGCGGGCTTGCCGATCAACCACACCGTGCGCGATCTGATCGACAGCGGCGATACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCa                                                                                                    >  2:425781/1‑140 (MQ=255)                                        tGGCTGTACAATGCCACCGTCGGTGCGGGCTTGCCGATCAACCACACCGTGCGCGATCTGATCGACAGCGGCGATACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACgg                                                                                           >  1:15004/1‑140 (MQ=255)                                                  aTGCCACCGTCGGTGCGGGCTTGCCGATCAACCACACCGTGCGCGATCTGATCGACAGCGGCGATACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACGGTAGCGTGCCg                                                                                 >  1:142861/1‑140 (MQ=255)                                                      cACCGTCGGTGCGGGCTTGCCGATCAACCACACCGTGCGCGATCTGATCGACAGCGGCGATACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACGGTAGCGTGCCGTTTa                                                                             >  1:351099/1‑140 (MQ=255)                                                      cACCGTCGGTGCGGGCTTGCCGATCAACCACACCGTGCGCGATCTGATCGACAGCGGCGATACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACGGTAGCGTGCCGTTTa                                                                             >  2:570709/1‑140 (MQ=255)                                                        ccGTCGGTGCGGGCTTGCCGATCAACCACACCGTGCGCGATCTGATCGACAGCGGCGATACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACGGTAGCGTGCCGTTTAcc                                                                           <  2:142861/140‑1 (MQ=255)                                                        ccGTCGGTGCGGGCTTGCCGATCAACCACACCGTGCGCGATCTGATCGACAGCGGCGATACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACGGTAGCGTGCCGTTTAcc                                                                           >  2:659226/1‑140 (MQ=255)                                                          gTCGGTGCGGGCTTGCCGATCAACCACACCGTGCGCGATCTGATCGACAGCGGCGATACTATTTTGTAGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACGGTAGCGTGCCGTTTACCGa                                                                         <  2:747141/140‑1 (MQ=255)                                                             ggTGCGGGCTTGCCGATCAACCACACCGTGCGCGATCTGATCGACAGCGGCGATACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACGGTAGCGTGCCGTTTACCGAGCt                                                                      >  1:733603/1‑140 (MQ=255)                                                                             tCAACCACACCGTGCGCGATCTGATCGACAGCGGCGATACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACGGTAGCGTGCCGTTTACCGAGCTGGTGGATCAGGCGTgg                                                      >  2:464577/1‑140 (MQ=255)                                                                                aCCACACCGTGCGCGATCTGATCGACAGCGGCGATACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACGGTAGCGTGCCGTTTACCGAGCTGGTGGATCAGGCGTGgcag                                                   <  2:107076/140‑1 (MQ=255)                                                                                aCCACACCGTGCGCGATCTGATCGACAGCGGCGATACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACGGTAGCGTGCCGTTTACCGAGCTGGTGGATCAGGCGTGgcag                                                   <  1:464577/140‑1 (MQ=255)                                                                                  cacaCCGTGCGCGATCTGATCGACAGCGGCGATACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACGGTAGCGTGCCGTTTACCGAGCTGGTGGATCAGGCGTGgcagca                                                 <  1:517426/140‑1 (MQ=255)                                                                                       cGTGCGCGATCTGATCGACAGCGGCGATACTATTTTGTTGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACGGTAGCGTGCCGTTTACCGAGCTGGTGGATCAGGCGTGGCAGCAGGGCt                                            >  2:361280/1‑140 (MQ=255)                                                                                       cGTGCGCGATCTGATCGACAGCGGCGATACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACGGTAGCGTGCCGTTTACCGAGCTGGTGGATCAGGCGTGGCAGCAGGGCt                                            >  2:783315/1‑140 (MQ=255)                                                                                       cGTGCGCGATCTGATCGACAGCGGCGATACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACGGTAGCGTGCCGTTTACCGAGCTGGTGGATCAGGCGTGGCAGCAGGGCt                                            >  1:431693/1‑140 (MQ=255)                                                                                        gTGCGCGATCTGATCGACAGCGGCGATACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACGGTAGCGTGCCGTTTACCGAGCTGGTGGATCAGGCGTGGCAGCAGGGCtt                                           >  2:275242/1‑140 (MQ=255)                                                                                         tGCGCGATCTGATCGACAGCGGCGATACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACGGTAGCGTGCCGTTTACCGAGCTGGTGGATCAGGCGTGgcag                                                   <  1:446996/131‑1 (MQ=255)                                                                                         tGCGCGATCTGATCGACAGCGGCGATACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACGGTAGCGTGCCGTTTACCGAGCTGGTGGATCAGGCGTGgcag                                                   >  2:446996/1‑131 (MQ=255)                                                                                         tGCGCGATCTGATCGACAGCGGCGATACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACGGTAGCGTGCCGTTTACCGAGCTGGTGGATCAGGCGTGGCAGCAGGGCTTa                                          >  1:569322/1‑140 (MQ=255)                                                                                           cgcgATCTGATCGACAGCGGCGATACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACGGTAGCGTGCCGTTTACCGAGCTGGTGGATCAGGCGTGGCAGCAGGGc                                             >  1:554397/1‑135 (MQ=255)                                                                                           cgcgATCTGATCGACAGCGGCGATACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACGGTAGCGTGCCGTTTACCGAGCTGGTGGATCAGGCGTGGCAGCAGGGc                                             <  2:554397/135‑1 (MQ=255)                                                                                           cgcgATCTGATCGACAGCGGCGATACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACGGTAGCGTGCCGTTTACCGAGCTGGTGGATCAGGCGTGGCAGCAGGGCTTAAc                                        <  2:157762/140‑1 (MQ=255)                                                                                           cgcgATCTGATCGACAGCGGCGATACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACGGTAGCGTGCCGTTTACCGAGCTGGTGGATCAGGCGTGGCAGCAGGGCTTAAc                                        <  2:569322/140‑1 (MQ=255)                                                                                           cgcgATCTGATCGACAGCGGCGATACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACGGTAGCGTGCCGTTTACCGAGCTGGTGGATCAGGCGTGGCAGCAGGGCTTAAc                                        >  1:821261/1‑140 (MQ=255)                                                                                           cgcgATCTGATCGACAGCGGCGATACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACGGTAGCGTGCCGTTTACCGAGCTGGTGGATCAGGCGTGGCAGCAGGGCTTAAc                                        >  1:80531/1‑140 (MQ=255)                                                                                               aTCTGATCGACAGCGGCGATACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACGGTAGCGTGCCGTTTACCGAGCTGGTGGATCAGGCGTGGCAGCAGGGCTTAACCGaa                                    >  1:767670/1‑140 (MQ=255)                                                                                                          aGCGGCGATACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACGGTAGCGTGCCGTTTACCGAGCTGGTGGATCAGGCGTGGCAGCAGGGCTTAACCGAACCTGACCcgcg                         >  2:669801/1‑140 (MQ=255)                                                                                                           gcggcgATACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACGGTAGCGTGCCGTTTACCGAGCTGGTGGATCAGGCGTGGCAGCAGGGCTTAACCGAACCTGACCCGCGt                        <  1:409442/140‑1 (MQ=255)                                                                                                             ggcgATACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACGGTAGCGTGCCGTTTACCGAGCTGGTGGATCAGGCGTGGCAGCAGGGCTTAACCGAACCTGACCCGCGTGa                      <  1:10224/140‑1 (MQ=255)                                                                                                                 aTACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACGGTAGCGTGCCGTTTACCGAGCTGGTGGATCAGGCGTGGCAGCAGGGCTTAACCGAACCTGAcc                             <  2:521207/129‑1 (MQ=255)                                                                                                                 aTACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACGGTAGCGTGCCGTTTACCGAGCTGGTGGATCAGGCGTGGCAGCAGGGCTTAACCGAACCTGAcc                             >  1:521207/1‑129 (MQ=255)                                                                                                                  tACTATTTTGTCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACGGTAGCGTGCCGTTTACCGAGCTGGTGGATCAGGCGTGGCAGCAGGGCTTAACCGAACCTGACCCGCGTGACGAtc                 <  2:656082/140‑1 (MQ=255)                                                                                                                   aCTATTTTGTCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACGGTAGCGTGCCGTTTACCGAGCTGGTGGATCAGGCGTGGCAGCAGGGCTTAACCGAACCTGACCCGCGTGACGAtct                <  1:425781/140‑1 (MQ=255)                                                                                                                       ttttGTCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCCGTTCCTGCAATTCGACGGTAGCGTGCCGTTTACCGAGCTGGTGGATCAGGCGTGGCAGCAGGGCTTAACCGAACCTGACCCGCGTGACGATctctct            >  2:266224/1‑140 (MQ=255)                                                                                                                            tCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACGGTAGCGTGCCGTTTACCGCGCTGGTGGATCAGGCGTGGCAGCAGGGCTTAACCGAACCTGACCCGCGTGACGATCTCTCTGGCaa       >  2:378156/1‑140 (MQ=255)                                                                                                                            tCGATCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACGGTAGCGTGCCGTTTACCGAGCTGGTGGATCAGGCGTGGCAGCAGGGCTTAACCGAACCTGACCCGCGTGACGATCTCTCTGGCaa       >  1:827333/1‑140 (MQ=255)                                                                                                                               aTCAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACGGTAGCGTGCCGTTTACCGAGCTGGTGGATCAGGCGTGGCAGCAGGGCTTAACCGAACCTGACCCGCGTGACGATCTCTCTGGCAAAGa    >  1:151213/1‑140 (MQ=255)                                                                                                                                 cAGCGAGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACGGTAGCGTGCCGTTTACCGAGCTGGTGGATCAGGCGTGGCAGCAGGGCTTAACCGAACCTGACCCGCGTGACGATCTCTCTGGCAAAGACg  >  1:172693/1‑140 (MQ=255)                                                                                                                                      |                                                                                                                                       CGCCAGATCCACGACGCCTTCGAAAAAACCGGGCGTCACTGGCTGTACAATGCCACCGTCGGTGCGGGCTTGCCGATCAACCACACCGTGCGCGATCTGATCGACAGCGGCGATACTATTTTGTCGATCAGCGGGATCTTCTCCGGCACGCTCTCCTGGCTGTTCCTGCAATTCGACGGTAGCGTGCCGTTTACCGAGCTGGTGGATCAGGCGTGGCAGCAGGGCTTAACCGAACCTGACCCGCGTGACGATCTCTCTGGCAAAGACG  >  NC_000913/4131584‑4131851

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