Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	556,630	A→G	intergenic (‑77/+31)	ybcI ← / ← ybcJ	conserved inner membrane protein YbcI/putative RNA‑binding protein YbcJ

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	556,630	0	A	G	100.0%	109.9 / NA	31	intergenic (‑77/+31)	ybcI/ybcJ	conserved inner membrane protein YbcI/putative RNA‑binding protein YbcJ
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/16);  total (15/16)

ACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGT  >  NC_000913/556502‑556749                                                                                                                                 |                                                                                                                        acccTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAGCCgg                                                                                                                     <  2:512640/133‑1 (MQ=255)       aaCCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAGCCGGCCTGtt                                                                                                               <  2:634195/133‑1 (MQ=255)        aCCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAGCCGGCCTGTTg                                                                                                              >  2:587424/1‑133 (MQ=255)          cAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAGCCGGCCTGTTGAt                                                                                                            >  1:53837/1‑133 (MQ=255)           aaTGCATAAGGGAACGGCGGCGCGCGTAATGACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAGCCGGCCTGTTGATa                                                                                                           <  1:200383/133‑1 (MQ=255)            aTGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAGCCGGCCTGTTGATa                                                                                                           <  2:136657/132‑1 (MQ=255)            aTGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAGCCGGCCTGTTGATa                                                                                                           >  1:136657/1‑132 (MQ=255)            aTGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAGCCGGCCTGTTGATAc                                                                                                          >  2:557525/1‑133 (MQ=255)             tGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAGCCGGCCTGTTGATACa                                                                                                         >  1:148018/1‑133 (MQ=255)                aTAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAGCCGGCCTGTTGATACAGGc                                                                                                      <  1:24882/133‑1 (MQ=255)                  aGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAGCCGGCCTGTTGATACAGGCCg                                                                                                    <  2:36419/133‑1 (MQ=255)                  aGGGGAACGGCGGCGGGCGTAAAAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAGGAGAGCCGGCCTGTTGATACAGGCCg                                                                                                    <  1:250485/133‑1 (MQ=255)                   ggggAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAGCCGGCCTGTTGATa                                                                                                           >  1:555499/1‑125 (MQ=255)                   ggggAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAGCCGGCCTGTTGATa                                                                                                           <  2:555499/125‑1 (MQ=255)                    gggAACGGCGGCGAGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAGCCGGCCTGTTGATACAGGCCggg                                                                                                  >  2:593033/1‑133 (MQ=255)                       aaCGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAGCCGGCCTGTTGATACAGGCCGGGaaa                                                                                               <  1:446860/133‑1 (MQ=255)                             ggcgTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAGCCGGCCTGTTGATACAGGCCGGGAAAGGGATc                                                                                         <  2:631764/133‑1 (MQ=255)                             ggcgTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAGCCGGCCTGTTGATACAGGCCGGGAAAGGGATc                                                                                         <  2:638593/133‑1 (MQ=255)                              gcgtgcgtAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAGCCGGCCTGTTGATACAGGCCGGGAAAGGGATCa                                                                                        >  1:579057/1‑133 (MQ=255)                                       ataaCGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAGCCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCaacaac                                                                               <  2:148018/133‑1 (MQ=255)                                        taaCGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAGCCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAAcc                                                                              >  1:48006/1‑133 (MQ=255)                                                ggCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAGCCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCt                                                                      <  2:48006/133‑1 (MQ=255)                                                ggCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAGCCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCt                                                                      >  1:27326/1‑133 (MQ=255)                                                         gTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAGCCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGc                                                             <  1:587424/133‑1 (MQ=255)                                                                cGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAGCCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTc                                                      >  1:539820/1‑133 (MQ=255)                                                                            gATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAGCCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCg                                          >  2:334447/1‑133 (MQ=255)                                                                                  aTCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAGCCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGcgacg                                    <  2:53837/133‑1 (MQ=255)                                                                                        aTGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAGCCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTg                              >  2:609039/1‑133 (MQ=255)                                                                                             tGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAGCCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTATGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCg                         >  2:630030/1‑133 (MQ=255)                                                                                               aaCGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAGCCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGtt                       >  2:132815/1‑133 (MQ=255)                                                                                                                       tGAAAAGAGAGCCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGt  >  2:51641/1‑130 (MQ=255)                                                                                                                       tGAAAAGAGAGCCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGt  <  1:51641/130‑1 (MQ=255)                                                                                                                                 |                                                                                                                        ACCCTAAACCAATGCATAGGGGAACGGCGGCGTGCGTAATAACGGTTGGCATGATGGTCGCTTCGGCAAAATGTCGATGCTATCAGCATGGATGAACGGGGCGTAGAGGGCAAAAGTCTGAAAAGAGAACCGGCCTGTTGATACAGGCCGGGAAAGGGATCAGGCAACAACCTGTACGCTGTGACCTGCAAAACTCACTGTCTGACCGGCGACGATTTTGCAGCGTTTGCGCGTTTCAACCGCACCGT  >  NC_000913/556502‑556749

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: ‑