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%	132.4 / NA	42	*387* (TAG→TAA)	ybdL	methionine aminotransferase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (17/25);  total (17/25)
*	NC_000913	634,747	0	C	T	100.0%	133.3 / NA	42	*210* (TAG→TAA)	ybdM	Spo0J family protein, ParB‑like nuclease domain
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/25);  total (17/25)

GGAGCACGGCGTAGCGGCGATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAGCTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGCTGTCCATTCCCAGCTCTTTGCCAATTTTATTATCGTCCCAT  >  NC_000913/634614‑634876                                                                                                                                     ||                                                                                                                                  ggAGCACGGCGTAGCGGCGATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGTTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTAttt                                                                                                                              >  2:788883/1‑139 (MQ=255)  gAGCACGGCGTAGCGGCGATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTa                                                                                                                             >  2:528564/1‑139 (MQ=255)   aGCACGGCGTAGCGGCGATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTaa                                                                                                                            <  2:499989/139‑1 (MQ=255)   aGCACGGCGTAGCGGCGATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTaa                                                                                                                            <  1:1518792/139‑1 (MQ=255)              gcggcgATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGgcg                                                                                                                 <  2:536026/139‑1 (MQ=255)                ggcgATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGcgcg                                                                                                               <  2:606718/139‑1 (MQ=255)                   gATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGa                                                                                                            <  1:365837/139‑1 (MQ=255)                            gTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACg                                                                                                   <  1:799337/139‑1 (MQ=255)                                      tGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAaac                                                                                         <  2:251392/139‑1 (MQ=255)                                        cgcCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCgg                                                                                                              >  2:320398/1‑116 (MQ=255)                                        cgcCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCgg                                                                                                              <  1:320398/116‑1 (MQ=255)                                                       acaTAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGtt                                                                        >  2:1059912/1‑139 (MQ=255)                                                         aTAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTaa                                                                      >  2:1343384/1‑139 (MQ=255)                                                                 aTTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTaa                                                                                                                            <  1:852500/77‑1 (MQ=255)                                                                 aTTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTaa                                                                                                                            >  2:852500/1‑77 (MQ=255)                                                                   tCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAg                                                            >  1:545554/1‑139 (MQ=255)                                                                   tCGGCTCTGTTTTTCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAg                                                            >  2:1202611/1‑139 (MQ=255)                                                                     gTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGgc                                                          <  2:1208381/139‑1 (MQ=255)                                                                         ctGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGTCCGTTAATTTGCTTCAGGCgcag                                                      <  1:1493830/139‑1 (MQ=255)                                                                                    agaagGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGGCGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGc                                           <  1:851674/139‑1 (MQ=255)                                                                                    agaagGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGc                                           <  1:1059912/139‑1 (MQ=255)                                                                                    agaagGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGc                                           <  2:490094/139‑1 (MQ=255)                                                                                     gaagGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGCt                                          <  1:528564/139‑1 (MQ=255)                                                                                      aagGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGc                                                                            <  1:1466516/104‑1 (MQ=255)                                                                                      aagGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGc                                                                            >  2:1466516/1‑104 (MQ=255)                                                                                           aTCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGCTGTCCAt                                    >  2:1202731/1‑139 (MQ=255)                                                                                                cGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGATGGTCGGCAaacaa                                                                                       >  1:1638312/1‑83 (MQ=255)                                                                                                cGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAaacaa                                                                                       <  2:1638312/83‑1 (MQ=255)                                                                                                    gCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGCTGTcc                                      >  2:1003507/1‑128 (MQ=255)                                                                                                    gCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGCTGTcc                                      <  1:1003507/128‑1 (MQ=255)                                                                                                        gcagcaGCTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGCTGTCCATTCCCAGCTCTTTg                       <  1:1202731/139‑1 (MQ=255)                                                                                                        gcagcaGCTGAACGCCTGCGCCAGCTTTAATTATGTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGCTGTCCATTCCCAGCTCTTTg                       <  1:1343384/139‑1 (MQ=255)                                                                                                          agcagcTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGCTGTCCATTCCCAGCTCTTTGcc                     <  2:694278/139‑1 (MQ=255)                                                                                                            cagcTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGCTGTCCATTCCCAGCTCTTTGcc                     >  1:508243/1‑137 (MQ=255)                                                                                                            cagcTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGCTGTCCATTCCCAGCTCTTTGcc                     <  2:508243/137‑1 (MQ=255)                                                                                                             agcTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCa                                                                                            <  2:1165310/65‑1 (MQ=255)                                                                                                             agcTGAACGCCTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCa                                                                                            >  1:1165310/1‑65 (MQ=255)                                                                                                                      ccTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGCTGTCCATTCCCAGCTCTTTGCCAATTTTATTATc         >  2:233511/1‑139 (MQ=255)                                                                                                                       cTGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGCTGTCCATTCCCAGCTCTTTGCCAATTTTATTATCg        >  1:635752/1‑139 (MQ=255)                                                                                                                        tGCGCCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGCTGTCCATTCCCAGCTCTTTGCCAATTTTATTATCGt       <  2:221642/139‑1 (MQ=255)                                                                                                                         gcgcCAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGCTGTCCATTCCCAGCTCTTTGCCAATTTTATTATCGtc      >  2:714411/1‑139 (MQ=255)                                                                                                                             cAGCTTTAATTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGCTGTCCATTCCCAGCTCTTTGCCAATTTTATTATCGtcccat  <  2:635752/139‑1 (MQ=255)                                                                                                                                     ||                                                                                                                                  GGAGCACGGCGTAGCGGCGATTCCGCTGTCGGTGTTTTGCGCCGATCCCTTCCCACATAAACTGATTCGTCTCTGTTTTGCCAAGAAGGAATCGACGTTGCTGGCAGCAGCTGAACGCCTGCGCCAGCTTTAGCTATTTAACCGTCCAGGCGCGGGAATATTGACGGTCGGCAAACAACTCTTGCAGGCCGTTAATTTGCTTCAGGCGCAGCACTTCGTCGCTGTCCATTCCCAGCTCTTTGCCAATTTTATTATCGTCCCAT  >  NC_000913/634614‑634876

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