Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	414,163	G→A	R531* (CGA→TGA)	sbcC ←	ATP‑dependent structure‑specific DNA nuclease ‑ SbcC subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	414,163	0	G	A	100.0%	31.3 / NA	12	R531* (CGA→TGA)	sbcC	ATP‑dependent structure‑specific DNA nuclease ‑ SbcC subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (5/7);  total (5/7)

CGTTCGTGCTCATCTTGTGCATCCAGCCACGGTTGAATATCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGTCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCGTCGCACCTTCTTCACCGAGCTTTTTAACTTCGTTTTCCAGCGCCAGTAATCGAGACTGATTAACGCCAGGCTCCAGCGCCTGATACGCCTCGACCGCCGGGTGGCTGGTGGAACCACAAAGTGGGCAAGGCTGACCCGCCTGTAACTGTGCACGTTGAGCTTCCAGCGTTTTGATGCGCGCTTCCTGCTCGCAAATGGTTTTCACATCGGCAAGTTGCTGCGTCTTTTCTTTATAACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCG  >  NC_000913/413925‑414384                                                                                                                                                                                                                                               |                                                                                                                                                                                                                              cGTTCGTGCTCATCTTGTGCATCCAGCCACGGTTGAATATCGTCCAGTGGCTGCGAGGTGATTTTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGTCGGAGGCTTTGCACTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCGTCGCACCTTCTTCACCGAGCTTTTTAATTTCGTTTTCCAGCGCCAGTAATCAAGACTGATTAAc                                                                                                                                                                                                                   >  1:56239/1‑251 (MQ=255)                                gTTGAATATCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGTCGGCTTGCCATTGTGGAGTAAGTGCTTGCTCATCTTGTCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACTTAGCGTCGCACCTTCTTCACCGAGCTTTTTAACTTCGTTTTCCAGCGCCAGTAATCAATACTGATTAACGCCAGGCTCCAGCGCCTGATACGCCTCGAcc                                                                                                                                                                                    <  1:351808/251‑1 (MQ=255)                                  tGAATATCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGTCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCGTCGCACCTTCTTCACCGAGCTTTTTAACTTCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGCCAGGCTCCAGCGCCTGATACGCCTCGAccg                                                                                                                                                                                   <  1:22358/250‑1 (MQ=255)                                                      cAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTTCTTGCTCATCTTGTCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCGTCGCACCTTCTTCACCGAGCTTTTTAACTTCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGCCAGGCTCCAGCGCCTGATACGCCTCGACCGCCGGGTGGCTGGTGGAACcan                                                                                                                                                              <  1:40957/251‑2 (MQ=255)                                                            gATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGTCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCGTCGCACCTTCTTCACCGAGCTTTTTAACTTCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGCCAGGCTCCAGCGCCTGATACGCCTCGACCGCCGGGTGGCTGGTGGAACCACAAAgn                                                                                                                                                         <  1:92577/250‑2 (MQ=255)                                                                                                                                           tcagcacccTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCGTCGAACCTTCTTCACCGCGCTTTTTAACTTCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGCCAGGCTCCAGCGCCTGATACGCCTCgccc                                                                                                                                                                                    >  1:6846/9‑141 (MQ=255)                                                                                                                                                                                 cccACGTAGCGTCGCACCTTCTTCACCGAGCTTTTTAACTTCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGCCAGGCTCCa                                                                                                                                                                                                        >  1:8347/1‑86 (MQ=255)                                                                                                                                                                                 cccACGTAGCGTCGCACCTTCTTCACCGAGCTTTTTAACTTCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGCCAGGCTCCa                                                                                                                                                                                                        <  2:8347/86‑1 (MQ=255)                                                                                                                                                                                        aNNGTCGCANNTTCTTCACCGAGCTTTTTAACTTCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGCCAGGCTCCAGCGCCTGATACGCCTCGAcc                                                                                                                                                                                    <  2:6846/99‑1 (MQ=255)                                                                                                                                                                                        aGCGTCGCACCTTCTTCACCGAGCTTTTTAACTTCGTTTTCCAGCGCCAGTAATCAAGACTGATTAACGCCAGGCTCCAGCGCCTGAt                                                                                                                                                                                               >  1:351744/1‑88 (MQ=255)                                                                                                                                                                                                aCCTTCTTCACCGAGCTTTTTAACTTCGTTTTCCAGCGCCCGTAATCAAGACTGATTAACGCCAGGATCCAGCGCCTGATACGCCTCGACCGCCGGGTGGCTGGTGGAAACACAAAGTGGGCAAGGCTGACCCGCCTGTAACTGTGCACGTTGAGCTTCCAGCGTTTTGATGCGCGCTTCCTGCTCGCAAATGGTTTTCACATCGGCAAGTTGCTGCGTCTTTTCTTTATAACGCTGGCGCATTTCGTTa                     >  1:160334/1‑250 (MQ=255)                                                                                                                                                                                                                   ttaattgcttttcccAGCGCCATTAATCAAGACTGATTAACCCCGGCCTCCAGCGCCTGAGCCGCCCCGCCCGCCGGTTGGCTGGTGAACCCACAAAGTGGGCAAGGCTGACCCGCCTGTAACTGTGCACGTTGAGCTTCCAGCGTTTTGATGCGCGCTTCCTGCTCGCAAATGGTTTTCACATCGGCAAGTTGCTGCGTCTTTTCTTTATAACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCg  <  1:49982/237‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                              CGTTCGTGCTCATCTTGTGCATCCAGCCACGGTTGAATATCGTCCAGTGGCTGCAAGGTGATATTGAGGCTGGCCGTGACGGCTTGCCATTGTTGAGTAAGTGCTTGCTCATCTTGTCGGAGGCTTTGCGCTTCGTTTTCATCACGCTGAAGCTGCTTTGTTATGGCGTCCAGTTGCCCACGTAGCGTCGCACCTTCTTCACCGAGCTTTTTAACTTCGTTTTCCAGCGCCAGTAATCGAGACTGATTAACGCCAGGCTCCAGCGCCTGATACGCCTCGACCGCCGGGTGGCTGGTGGAACCACAAAGTGGGCAAGGCTGACCCGCCTGTAACTGTGCACGTTGAGCTTCCAGCGTTTTGATGCGCGCTTCCTGCTCGCAAATGGTTTTCACATCGGCAAGTTGCTGCGTCTTTTCTTTATAACGCTGGCGCATTTCGTTAAGTGCGGCGTTACGTTGCG  >  NC_000913/413925‑414384

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