Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	634,746	2 bp→AT		[ybdL]–[ybdM]	[ybdL], [ybdM]

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	634,746	0	G	A	100.0%	65.1 / NA	21	*387* (TAG→TAA)	ybdL	methionine aminotransferase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (7/14);  total (7/14)
*	NC_000913	634,747	0	C	T	100.0%	64.1 / NA	21	*210* (TAG→TAA)	ybdM	Spo0J family protein, ParB‑like nuclease domain
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (7/14);  total (7/14)

AGCACGGCGTAGCGGCGATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAGCTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGCTGTCCATTCCCAGCTCTTTGCCAATTTTATTATCGT  >  NC_000913/634616‑634871                                                                                                                                   ||                                                                                                                             aGCACGGCGTAGCGGCGATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTaa                                                                                                                       <  1:252780/139‑1 (MQ=255)      ggCGTAGCGGCGATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTc                                                                                                                  <  1:711275/139‑1 (MQ=255)          tAGCGGCGATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAATCTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAgg                                                                                                              <  1:349683/139‑1 (MQ=255)             cggcgATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGgcgc                                                                                                           <  2:282405/139‑1 (MQ=255)                cgATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCggg                                                                                                        >  1:630679/1‑139 (MQ=255)                              gtgtTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCg                                                                                          <  1:501375/139‑1 (MQ=255)                                   ttGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCaaa                                                                                     <  2:209515/139‑1 (MQ=255)                                    tGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAaac                                                                                    <  2:437526/139‑1 (MQ=255)                                           aTCCCGTCCCACATAAACTGATTCGTTTCGGTTTTGCCATGAAGGAAACGACGTTTCTGGCATCAGCTGAACGACTGCGCCAGTTTTAATTATTTAACTGTCCATGAGCCGGATTATTGACGGTCGGCAAACAACTCtt                                                                             >  1:189004/1‑139 (MQ=255)                                            tCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTg                                                                            <  2:135587/139‑1 (MQ=255)                                              ccTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCa                                                                          <  2:660584/139‑1 (MQ=255)                                                            cTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGc                                                            >  1:119140/1‑139 (MQ=255)                                                                                 aagaagGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTtcgtcg                                       <  2:630679/139‑1 (MQ=255)                                                                                  agaagGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGc                                      <  1:487793/139‑1 (MQ=255)                                                                                  agaagGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGc                                      >  2:487793/1‑139 (MQ=255)                                                                                    aagGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACtt                                            <  1:43813/131‑1 (MQ=255)                                                                                    aagGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACtt                                            >  2:43813/1‑131 (MQ=255)                                                                                                        agcagcTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGgcg                                                                                                            <  1:152969/47‑1 (MQ=255)                                                                                                        agcagcTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGgcg                                                                                                            >  2:152969/1‑47 (MQ=255)                                                                                                        agcagcTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGCTGTCCATTCCCAGCTCTTTGcc                >  2:79829/1‑139 (MQ=255)                                                                                                                      tGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCTGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGCTGTCCATTCCCAGCTCTTTGCTAATTTTATTATCGt  <  2:485862/139‑1 (MQ=255)                                                                                                                                   ||                                                                                                                             AGCACGGCGTAGCGGCGATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAGCTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGCTGTCCATTCCCAGCTCTTTGCCAATTTTATTATCGT  >  NC_000913/634616‑634871

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