Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	5,610,048	C→A	S173I (AGC→ATC)	PP_4929 ←	LysR family transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,610,048	0	C	A	100.0%	55.7 / NA	21	S173I (AGC→ATC)	PP_4929	LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (9/12);  total (9/12)

GGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGCTC  >  NC_002947/5609929‑5610194                                                                                                                        |                                                                                                                                                   ggTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCg                                                                                                                      >  2:763562/1‑150 (MQ=255)       cGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGc                                                                                                                >  2:189969/1‑150 (MQ=255)           ttgttgTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGCGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTc                                                                                                            <  2:379862/150‑1 (MQ=255)              ttgtGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTg                                                                                                         >  2:773445/1‑150 (MQ=255)                gtgGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTgga                                                                                                       <  1:607044/150‑1 (MQ=255)                 tgGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTggag                                                                                                      >  2:816753/1‑150 (MQ=255)                      tGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTggag                                                                                                      >  1:375156/1‑145 (MQ=255)                      tGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTggag                                                                                                      <  2:375156/145‑1 (MQ=255)                         tCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACa                                                                                              <  2:718304/150‑1 (MQ=255)                              aaCTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGccc                                                                                         <  2:220168/150‑1 (MQ=255)                                      gTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGgc                                                                                 <  1:464545/150‑1 (MQ=255)                                               cAATTGTTGGCGCGTGACCCGGCCTTCGCCCGCCAGCGGTGGGCCGACCGCGCGGACCGTGACCCTCTACACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTc                                                                        <  2:166059/150‑1 (MQ=255)                                                aaTTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTcc                                                                       >  2:5027/1‑150 (MQ=255)                                                                           cccgccAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTgg                                            >  2:104376/1‑150 (MQ=255)                                                                                   cGGGTGGCCGACCGGGCACACCGTGACCATCTCCACGATGCCCAACTCAGTGCGCTCCAAGTAAGCCGAGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGTCGCGCCCCTCCACCAGCTTGCCCGCCACATCGCCCTGGGCTCCCCt                                    <  2:654588/150‑1 (MQ=255)                                                                                           cGACCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGtt                            >  1:691627/1‑150 (MQ=255)                                                                                             aCCGCGCACACCGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCACCCTGGGCTCCGCTGGCCAGTTGc                          <  1:763562/150‑1 (MQ=255)                                                                                                        cGTGACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCACCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGaa               >  1:135755/1‑150 (MQ=255)                                                                                                           gACCATCTCCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACgg            <  1:189969/150‑1 (MQ=255)                                                                                                                 ctcCACGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCg      <  1:38103/150‑1 (MQ=255)                                                                                                                     aCGATGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGCTc  <  1:816753/150‑1 (MQ=255)                                                                                                                        |                                                                                                                                                   GGTAACCGCTTTGTTGTGGCGTGATCAGCAACTGCCGGTGGCGGGCCAATTGTTGGCGGGTGACCCGGCCTTCGCCCGCCAGCGGGTGGCCGACCGCGCACACCGTGACCATCTCCACGCTGCCCAACGCACGGCGCTCCAAGGAAGCCGGGATGCTTTCGTGGAGGAAGAACAGGCCCAGATCGGCGCGCCGCTCCACCAGCTTGCGCGCCACATCGCCCTGGGCTCCGCTGGCCAGTTGCACCTCCAGGAACGGGTAGCGGCTC  >  NC_002947/5609929‑5610194

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