Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,601,098	C→A	S48S (TCG→TCT)	lsrR ←	DNA‑binding transcriptional repressor LsrR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,601,098	0	C	A	100.0%	50.9 / NA	18	S48S (TCG→TCT)	lsrR	DNA‑binding transcriptional repressor LsrR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (10/8);  total (10/8)

GGTTGCCTCGCCAAAACCAATCGCCAGCATCTGTTGTGGTTGAAGTAAACTCATCAACATATGCGCCGCGCCTATCCCCAGTCGCCCACCGACATCAGCATCCGCAAGCCCAGGGATCACCCGGACATGTTGCAGCGAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAGTAAAACCACGCGATCCGCGCGACCTGTTCTTCTTCACACATTCCCTGTTCTGAAATTGCCGAATCGTTGATTGTCATAATTCATTCTTCACTTTGAACATATTTAAATCTTTAATGCAATTGTTCAGTTCTTGCTCATTTATATCTGTGATGGCAACCACAGTTTGACTCTAC  >  NC_000913/1600858‑1601337                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                ggTTGCCTCGCCAAAACCAATCGCCAGCATCTGTTGTGGTTGAAGTAAACTCAGCAACATATGCGCCGCGCCTATCCCCAGTCGCCCACCGACATCAGCATCCGCAAGCCCAGGGATCACCCGGACATGTTGCAGCGAAAACTGACGACGTAATTGAGGTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGAGACATTTTc                                                                                                                                                                                                                                        <  1:247834/250‑1 (MQ=255)           ccAAAACCAATCGCCAGCATCTGTTGTGGTTGAAGTAAACTCATCAACATATGCGCCGCGCCTATCCCCAGTCGCCCACCGACATCAGCATCCGCAAGCCCAGGGATCACCCGGACATGTTGCAGCGAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGAGACATTTTCAAACgtgtn                                                                                                                                                                                                                               <  1:116528/249‑2 (MQ=255)                           gCATCTGTTGTGGTTGAAGTAAACTCATCAACATATGCGCCGCGCCTATCCCCAGTCGCCCACCGACATCAGCATCCGCAAGCCCAGGGATCACCCGGACATGTTGCAGCGAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGAGACATTTTCAAACGTGTCAGGCCGAGACGATCGCTg                                                                                                                                                                                                             >  1:19742/1‑251 (MQ=255)                                 gttgtGGTTGAAGTAAACTCATCAACATATGCGCCGCGCCTATCCCCAGTCGCCCACCGACATCAGCATCCGCAAGCCCAGGGATCACCCGGACATGTTGCAGCGAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGAGACATTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTNg                                                                                                                                                                                                       <  1:350527/251‑1 (MQ=255)                                   tgtgGTTGAAGTAAACTCATCAACATATGCGCCGCGCCTAACCCCAGTCGCACACCGACATAAGCATCAGCACGCCCAGGGATCACCCGGACATGTTGCAGCGAAAACTGACGACGTAATTGAGTTTAATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGAGACATTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGct                                                                                                                                                                                                     >  1:249049/1‑251 (MQ=255)                                                                                    ncccACCGACATCAGCATCCGCAAGCCCAGGGATCACCCGGACATGTTGCAGCGAAAACTGACGACGTCATTGAGTTTCATATTCCAGCCAGCCTTCAAAGCGAGAATTCATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGAGACATTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAGTAAAACCACGCGATCCGCGCGACCTGt                                                                                                                                                    >  1:53680/2‑251 (MQ=255)                                                                                          cGACATCAGCATCCGCAAGCCCAGGGATCACCCGGACATGTTGCAGCGAAAACTGACGACGAAATTGAGTTTACTATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGAGACATTTTCAAACGTGTCAGGcc                                                                                                                                                                                                                          >  1:451048/1‑175 (MQ=255)                                                                                                                                tGTTGCAGCGAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGAGACATTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAGTAAAAc                                                                                                                                                                         <  1:14519/186‑1 (MQ=255)                                                                                                                                     cAGCGAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGAGACATTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAGTAAAACCACGCGATCCGCGCGACCTGTTCTTCTTCACACATTCCCTGt                                                                                                                               >  1:75759/1‑223 (MQ=255)                                                                                                                                        cGAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGAGACATTTTCAAACGTGTCAGGCCGAGACGATCGCTGAACTCGCTCTGGGTCAGCCCGTCGTGATAGTAAAACCACGCGATCCGCGCGACCTGTTCTTCTTCACACATTCCCTGTTCTGAAATTGCCGAATCGTTGATTGTCATa                                                                                                 >  1:167897/1‑250 (MQ=255)                                                                                                                                                                         naGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGAGACATTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAGTAAAACCACGCGATCCGCGCGACCTGTTCTTCTTCACACATTCCCTGTTCTGAAATTGCCGAATCGTTGAt                                                                                                        >  1:262880/2‑210 (MQ=255)                                                                                                                                                                               gCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGAGACATTTTCAACCgtg                                                                                                                                                                                                                                 >  1:57805/1‑83 (MQ=255)                                                                                                                                                                                ccTTCAACGCGAGAATTAATCTGTACGCGAATCATGCCGGACTGATGCCCTTTCTCCAGCAATCGAGACATTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAGTAAAACCACGCGATCCGCGCGACCTGTTCTTCTTCACACATTCCCTGTTCTGAAATTGCCGAATCGtt                                                                                                           >  1:156228/1‑200 (MQ=255)                                                                                                                                                                                       aGCGAGAATTAATCTGTACGCGAATAATGGCGGACTGATGCGCTTTCTCCAGCAATCGAGACATTTTCAAACGTGTCAGGCCGAGACGATCGCtn                                                                                                                                                                                                             <  1:48583/95‑2 (MQ=255)                                                                                                                                                                                               ttAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGAGACATTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGc                                                                                                                                                                                           <  1:425082/105‑1 (MQ=255)                                                                                                                                                                                                                      ggACTGATGCCCTTTCTCCAGCAATCGAGACATTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAGTAAAACCACGCGATCCGCGCGACCTGTTCTTCTTCACACATTCCCTGTTCTGAAATTGCCGAATCGTTGATTGTCATAATTCATTCTTCACTTTGAACATATTTAAATCTTTAATGCAATTGTTCAGTTCTTGCTCAGTTATAGCTGTGATTGCAAc                  <  1:138623/251‑1 (MQ=255)                                                                                                                                                                                                                       gACTGATGCCCTTTCTCCAGCAATCGAGACATTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAGTAAAACCACGCGATCCGCGCGACCTGTTCTTCTTCACACATTCCCTGTTCTGAAATTGCCGAATCGTTGATTGTCATAATTCATTCTTCACTTTGAACATATTTAAATCTTTAATGCAATTGTTCAGTTCTTGCTCATTTATATCTGTGATGGCAAcc                 <  1:152394/251‑1 (MQ=255)                                                                                                                                                                                                                                         aGCAATCGAGACATTTTCAGAAGTGTCAGGCAGAGACGATCGATGATCGCGATGGGCGTCAGCCCGTGGAGATCGTAAAAGGACGCGATCCGCGCGCCCGGTTCTTCTTCACACATTCCGTGTTCTGAAATTGCCGAATCGTTGGTGGTCATAATGCATTCTTCACTTTGAACATATTTAAATCTTTAATGCAATTGTTCCGTTCGTGCTCATTTATATCTGTGATGGCAACCACGGTTTGCCTCTAc  >  1:232341/1‑248 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                GGTTGCCTCGCCAAAACCAATCGCCAGCATCTGTTGTGGTTGAAGTAAACTCATCAACATATGCGCCGCGCCTATCCCCAGTCGCCCACCGACATCAGCATCCGCAAGCCCAGGGATCACCCGGACATGTTGCAGCGAAAACTGACGACGTAATTGAGTTTCATATTCCAGACAGCCTTCAAAGCGAGAATTAATCTGTACGCGAATAATGCCGGACTGATGCCCTTTCTCCAGCAATCGCGACACTTTCAAACGTGTCAGGCCGAGACGATCGCTGATCTCGCTCTGGGTCAGCCCGTCGTGATAGTAAAACCACGCGATCCGCGCGACCTGTTCTTCTTCACACATTCCCTGTTCTGAAATTGCCGAATCGTTGATTGTCATAATTCATTCTTCACTTTGAACATATTTAAATCTTTAATGCAATTGTTCAGTTCTTGCTCATTTATATCTGTGATGGCAACCACAGTTTGACTCTAC  >  NC_000913/1600858‑1601337

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