Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,489,016	C→T	E350K (GAA→AAA)	yjgR ←	DUF853 domain‑containing protein YjgR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,489,016	0	C	T	100.0%	111.6 / NA	35	E350K (GAA→AAA)	yjgR	DUF853 domain‑containing protein YjgR
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (17/18);  total (17/18)

GCCTGAAAGCCTTTTTGCAACATCTCATAGGCGGATTCCCGGTCCACCTCATCCTCATATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATCTTATCCAGCAGTACCTGCGGGGCATCGTTAAACAGC  >  NC_000913/4488784‑4489254                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                               gCCTGAAAGCCTTTTTGCAACATCCCATAGGCGGATTCCCGGTCCACCTCATCCTCATATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTGCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGAtt                                                                                                                                                                                                                              <  1:38925/251‑1 (MQ=255)  ccTGAAAGCCTTTTTGCAACATCTCATAGGCGGATTCCCGGTCCACCTCATCCTCATATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCCCAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTg                                                                                                                                                                                                                             >  1:242390/1‑251 (MQ=255)                      aTCTCATAGGCGGAGTCCCGGTCCACCTCGTCCTCATATTTGCCATACACCGGGGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCCCCGGCCCCATCCGCGAACAAGGCGCGATCCCCATCGCACGCTCCACCACAGCAGGGCTTCCTTTCGCATCCAGCAAAGAGATCAACGCCTCCCCGGTGCCCCGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTATgcgcagc                                                                                                                                                                                                         >  2:51438/1‑248 (MQ=255)                              ggCGGATTCCCGGTCCACCTCATCCTCATATTGGCCACACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCa                                                                                                                                                                                                 <  2:133032/251‑1 (MQ=255)                                  gatgcccGGTCCACCTCTGCATCATATTCGCCATACACCGGAGGGTTATTAAACAAGACATTACGCTAATCTTCCGTCAACGGCCCCATCCGCGAACAAGGCGCGCACACCATCGCACGCTCCACCACAGAAGGGCTTCCGTTCGCAGCCAGAAAATAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATAAAGCTCCGGATGGGCCTGCATGGTTTGCGCCGCTGCCTTCACTGc                                                                                                                                                                                             <  2:47302/247‑1 (MQ=255)                                                  catcctcatATTTTCCAGTCACCGGATAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGc                                                                                                                                                                             <  1:138555/251‑1 (MQ=255)                                                  catcctcatATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTg                                                                                                                                                                               <  1:198793/249‑1 (MQ=255)                                                     cctcatATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGt                                                                                                                                                                           <  1:121249/250‑1 (MQ=255)                                                         atatTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTaaa                                                                                                                                                                        >  1:17152/1‑249 (MQ=255)                                                               gCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGc                                                                                                                                                                 <  1:11232/250‑1 (MQ=255)                                                               gCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGc                                                                                                                                                                 >  2:128523/1‑250 (MQ=255)                                                                                                             gtccccGGCCCCATCCGCGAACAAGGCGCGATCACCAGCGCACGCTCCACCTCAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCTGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGGGCTGCCCGAGCACATTATCCgg                                                                                                                  <  1:189326/247‑1 (MQ=255)                                                                                                             gTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCgg                                                                                                                  >  1:88714/1‑251 (MQ=255)                                                                                                                  cGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAatat                                                                                                             <  2:98580/251‑1 (MQ=255)                                                                                                                                   aaGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGaa                                                                                            <  2:17152/251‑1 (MQ=255)                                                                                                                                   aaGGCGCGATAACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGCCGGGTTTTGCgc                                                                                             >  2:176754/1‑249 (MQ=255)                                                                                                                                                                    aaGGGCTTCCTTTCGCATCCAGACAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAAt                                                            >  2:41680/1‑250 (MQ=255)                                                                                                                                                                     aGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATaa                                                          <  1:64857/251‑1 (MQ=255)                                                                                                                                                                        gCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGcc                                                       >  2:92502/1‑251 (MQ=255)                                                                                                                                                                            ccTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCTATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTAt                                                   <  1:161250/251‑1 (MQ=255)                                                                                                                                                                                       cAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCaa                                        <  1:128523/251‑1 (MQ=255)                                                                                                                                                                                        aGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAAt                                       >  1:158967/1‑251 (MQ=255)                                                                                                                                                                                         gAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATc                                      <  2:88714/251‑1 (MQ=255)                                                                                                                                                                                             agagATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATCTTAt                                  >  1:198994/1‑251 (MQ=255)                                                                                                                                                                                                gaTCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATCTTATCca                               <  1:197116/251‑1 (MQ=255)                                                                                                                                                                                                 aTCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGTACATTATCCGTAATATCAGACGGGTTTTGCGaa                                                                                            <  1:80218/189‑1 (MQ=255)                                                                                                                                                                                                 aTCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATCTTATcc                                >  1:150178/1‑249 (MQ=255)                                                                                                                                                                                                 aTCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATCTTATCcag                              <  1:2387/251‑1 (MQ=255)                                                                                                                                                                                                 aTCAACGCCTCGACGGAGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGACCGCAAAGCATGTTGAACGCGCTTACCGAGCTGCCCGAGCACATTATCCGGACTATCAGACGGGTTTTGCGaa                                                                                            >  2:80218/1‑189 (MQ=255)                                                                                                                                                                                                      cGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATCTTATCCAGCAGt                          >  1:85401/1‑250 (MQ=255)                                                                                                                                                                                                           cGCCGGTGCCCAGTTCCTGAATCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATCTTATCCAGCAGTACCTGc                    >  2:138023/1‑251 (MQ=255)                                                                                                                                                                                                                           cTGAAACGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATCTTATCCAGCAGTACCTGCGGGGCATCGTTAAAc     >  2:89860/1‑250 (MQ=255)                                                                                                                                                                                                                                tCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGGGCGAATAAGCCTTATCACCTGCt                                           <  1:74220/207‑1 (MQ=255)                                                                                                                                                                                                                                tCGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCt                                           >  2:74220/1‑207 (MQ=255)                                                                                                                                                                                                                                 cGCCTTTTTGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATCTTATCCAGCAGTACCTGCGGGGCATCGTTAAAcagc  >  2:8257/1‑247 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                               GCCTGAAAGCCTTTTTGCAACATCTCATAGGCGGATTCCCGGTCCACCTCATCCTCATATTTGCCATACACCGGAGAGTGATTAATCAAGCCATTACGCTCATCTTCCGTCACCGGCCCCATCCGCGAACAAGGCGCGATCACCATCGCACGCTCCACCACAGAAGGGCTTCCTTTCGCATCCAGAAAAGAGATCAACGCCTCGCCGGTGCCCAGTTCCTGAATCGCCTTTTCGGTATCAAACGCCGGATTGGCCCGCATGGTTTGCGCCGCTGCCTTCACTGCTTTCTGATCTTTTGGCGTAAAAGCCCGCAAAGCATGTTGAACGCGATTACCGAGCTGCCCGAGCACATTATCCGGAATATCAGACGGGTTTTGCGAAACGAACCAGACGCCCACGCCTTTTGAGCGAATAAGCCTTATCACCTGCTCAATCTTATCCAGCAGTACCTGCGGGGCATCGTTAAACAGC  >  NC_000913/4488784‑4489254

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