Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	606,748	A→G	S212S (AGT→AGC)	ybdK ←	carboxylate‑amine ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	606,748	0	A	G	100.0%	136.0 / NA	38	S212S (AGT→AGC)	ybdK	carboxylate‑amine ligase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (18/20);  total (18/20)

AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCT  >  NC_000913/606617‑606869                                                                                                                                    |                                                                                                                          aaTCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAGCTCATACa                                                                                                                   <  2:190903/140‑1 (MQ=255)                   gCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAGCTCAGACAGCGAAACAGGGCTTCaaa                                                                                                 >  2:102713/1‑140 (MQ=255)                     gTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAGCTCAGACAGCGAAACAGGGCTTCAAAtt                                                                                               >  2:178614/1‑140 (MQ=255)                      tGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGTGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAGCTCAGACAGCGAAACAGGGCTTCAAAttg                                                                                              >  2:387007/1‑140 (MQ=255)                       ggCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAGCTCAGACAGCGAAACAGGGCTTCAAAttgt                                                                                             >  2:64484/1‑140 (MQ=255)                        gCTAAGGGTTAACGGGGTATCCATCACCC‑AACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAGCTCAGACAGCGAAACAGGGCTTCAAAttgttg                                                                                           >  1:439777/1‑140 (MQ=255)                           aaGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAGCTCAGACAGCGAAACAGGGCTTCAAATTGTTGcc                                                                                         <  1:64484/140‑1 (MQ=255)                                  aaCGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAGCTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACt                                                                                  >  1:108604/1‑140 (MQ=255)                                  aaCGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAGCTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACt                                                                                  >  1:16813/1‑140 (MQ=255)                                        gTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAGCTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCa                                                                            >  2:443517/1‑140 (MQ=255)                                        gTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAGCTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCa                                                                            <  1:1055/140‑1 (MQ=255)                                         tATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAGCTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCAc                                                                           <  2:518902/140‑1 (MQ=255)                                                cACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAGCTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCg                                                                    >  2:44347/1‑140 (MQ=255)                                                    cGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAGCTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGcc                                                                <  2:108604/140‑1 (MQ=255)                                                        ccTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAGCTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTa                                                            <  1:284257/140‑1 (MQ=255)                                                             aCCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAGCTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAgg                                                       >  2:493845/1‑140 (MQ=255)                                                                 tGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAGCTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAg                                                   >  1:384303/1‑140 (MQ=255)                                                                                  gggCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAGCTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCgg                                  >  2:441241/1‑140 (MQ=255)                                                                                   ggCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAGCTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGt                                 >  1:93699/1‑140 (MQ=255)                                                                                       aaTATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAGCTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTg                             >  1:133596/1‑140 (MQ=255)                                                                                       aaTATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAGCTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTg                             <  2:384303/140‑1 (MQ=255)                                                                                         tatCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAGCTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTgag                           >  2:486915/1‑140 (MQ=255)                                                                                               aGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAGCTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCa                     <  1:486915/140‑1 (MQ=255)                                                                                                    aGATCTTTAATGCTGTCGATCATCGTGGTGTAGCTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAAcg                <  1:199871/140‑1 (MQ=255)                                                                                                    aGATCTTTAATGCTGTCGATCATCGTGGTGTAGCTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAAcg                <  1:226890/140‑1 (MQ=255)                                                                                                    aGATCTTTAATGCTGTCGATCATCGTGGTGTAGCTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAAcg                >  2:199871/1‑140 (MQ=255)                                                                                                        cTTTAATGCTGTCGATCATCGTGGTGTAGCTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTa            >  2:409172/1‑140 (MQ=255)                                                                                                         tttAATGCTGTCGATCATCGTGGTGTAGCTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTAt           <  1:387007/140‑1 (MQ=255)                                                                                                            aaTGCTGTCGATCATCGTGGTGTAGCTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCg        <  1:249482/140‑1 (MQ=255)                                                                                                            aaTGCTGTCGATCATCGTGGTGTAGCTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCg        <  1:178614/140‑1 (MQ=255)                                                                                                              tGCTGTCGATCATCGTGGTGTAGCTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGtt      <  1:493845/140‑1 (MQ=255)                                                                                                              tGCTGTCGATCATCGTGGTGTAGCTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGtt      <  2:93699/140‑1 (MQ=255)                                                                                                               gCTGTCGATCATCGTGGTGTAGCTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTc     <  2:439777/140‑1 (MQ=255)                                                                                                                cTGTCGATCATCGTGGTGTAGCTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTcc    <  2:133596/140‑1 (MQ=255)                                                                                                                  gTCGATCATCGTGGTGTAGCTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCt  <  1:409172/140‑1 (MQ=255)                                                                                                                  gTCGATCATCGTGGTGTAGCTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCt  <  1:441241/140‑1 (MQ=255)                                                                                                                  gTCGATCATCGTGGTGTAGCTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCt  <  1:443517/140‑1 (MQ=255)                                                                                                                  gTCGATCATCGTGGTGTAGCTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTAGTGAGGAGGCAAAACGCGTATCCGTTCCCt  >  1:285589/1‑140 (MQ=255)                                                                                                                                    |                                                                                                                          AATCCCGCCATATTTACTGCGTGGCTAAGGGTTAACGGGGTATCCATCACCCGAACCTCCACCGTGCCAAAATGAGGACTGGGGCGAATATCCCAGTGCAGATCTTTAATGCTGTCGATCATCGTGGTGTAACTCAGACAGCGAAACAGGGCTTCAAATTGTTGCCAGTTACTGACCCACGGCATCGGGCCATTATCAGGAAAGGCGGAAAAAATATTCGGTCGTGAGGAGGCAAAACGCGTATCCGTTCCCT  >  NC_000913/606617‑606869

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