Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,812,146	C→G	R138S (AGG→AGC)	yicR ←	RadC‑like JAB domain‑containing protein YicR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,812,146	0	C	G	96.8%	98.2 / ‑5.6	31	R138S (AGG→AGC)	yicR	RadC‑like JAB domain‑containing protein YicR
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base G (20/10);  total (20/11)
Fisher's exact test for biased strand distribution p-value = 3.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CATACTCTCCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGTTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTATGTGCAAGGATCAGCGCCGAGGCGTTTATTTTTATCGCTTCGCGGATAATTTCCCGAGGATGGACTTCAACATGGTTTAGCGTGCCGGAAAAAAGACGCCTATGCGTTATAACCCGGTGTTGGGAGTCGAGAAAGATCACCATAAAGATCTCCCGCTCCTCACCCGTGAGCTGGCTTTGTAAAAATTCCCTCGTCATCTCCGGGCTGAGTAAAGGGCTTTCTTCACGCATCCGCACGTTGTAGTAACGCCGCGCCAGTTCAGCAATCCCCTTTAACTGGGCAAATTTCGCCACGCCAATTCCATGAACGCCACTAAATTGTTC  >  NC_000913/3811920‑3812370                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                 cATACTCTCCACGCCCAATCACGCTATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGTTCAGTCATGAGTTTATCCGCTTTACTGGGTTCAGCACACCCCGAAGGGTGATTATGTGCAAGGATCAGCGCCGAGGCGTTTATTTTTATCGCTTCGCGGATAATTTCCCGAGGATGGACTTCAACATGGTTTAGCGTGCCGGAAAAAAGACGGCTATGCGTTATAACCCGGTGTTgt                                                                                                                                                                                                          >  2:373656/1‑250 (MQ=255)       ctcCACGCACAACCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGTTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTATGTGCAAGGCTCAGCGCCGAGGCGTTTATTTTTATCGCTTCGCGGATAATTTCCCGAGGATGGACTTCAACATGGTTTAGCGTGCCGGAAAAAAGACGGCTATGCGTTATAACCCGGTGTTGTGAGTCg                                                                                                                                                                                                    >  2:224290/1‑251 (MQ=255)        tcCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCAGGAACTGCCAACTCTTTATTATCCGTTCAGTAATGAGTTTCTCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTATGTGCAAGGATAAGCGCCGAGGCGTTTATTTTTATCGCTTCGCGGATAATTTCCCGAGGATGGACTTCAACATGGTTTAGCGTGCCGGAAAAAAGACGGCTATGCGTTATAACCCGGTGTTGTGAGTCgg                                                                                                                                                                                                   >  1:207024/1‑250 (MQ=255)         ccACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGTTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCTAAGGGTGATTATGTGCAAGGATCAGCGCCGAGGCGTTTATTTTTATCGCTTCGCGGATAATTTCCCGAGGAGGGACTTCAACATGGTTTAGCGTGCCGGAAAAAAGACGGCTATGCGTTATAACCCGGTGTTGTGAGTCgag                                                                                                                                                                                                  <  1:293267/251‑1 (MQ=255)         ccACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGTTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTATGTGCAAGGATCAGCGCCGAGGCGTTTATTTTTATCGCTTCGCGGATAATTTCCCGAGGATGGACTTCAACATGGTTTAGCGTGCCGGAAAAAAGACGGCTATGCGTTATAACCCGGTGTTGTGAGTc                                                                                                                                                                                                     <  1:273239/248‑1 (MQ=255)          cACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGTTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTATGTGCAAGGATCAGCGCCGAGGCGTTTATTTTTATCGCTTCGCGGATAATTTCCCGAGGATGGACTTCAACATGGTTTAGCGTGCCGGAAAAAAGACGGCTATGCGTTATAACCCGGTGTTGTGAGTCg                                                                                                                                                                                                    >  1:103838/1‑248 (MQ=255)                             gTCGAGCACGCGTAAATCCATGAACTTAAAACTCTTTAGTATCCTTTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCAGATCCCGAAGGGTGATTATGTGCGAGGATCAGCGCCGAGGTGTGTATTTTTAGCGCGTCGCGGATAATTTTCCGAGTATGGACTTCAACATCGTTTAGCGGGCCGGAAAAAAGACGCCGATGCCGTATAACCCGGTGTTGTGAGTCGAGAAAGATCAACATAAAGAtct                                                                                                                                                                              <  1:412057/251‑1 (MQ=255)                             gTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGTTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTATGTGCAAGGATCAGCGCCGAGGCGTTTATTTTTATCGCTTCGCGGATAATTTCCCGAGGATGGACTTCAACATGGTTTAGCGTGCCGGAAAAAAGACGGCTATGCGTTATAACCCGGTGTTGTGAGTCGAGAAAGATCACCATAAAGAtct                                                                                                                                                                              >  2:273905/1‑251 (MQ=255)                              tCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTCTCCGTTCCGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTATGTGCAAGGATCAGCGCCGAGGCGTTTATTTTTATCGCTTCGCGGATAATTTCCCGAGGATGGACTTCAACATGGTTTAGCGTGCCGGAAAAAAGACGGCTATGCGTTATAACCCGGTGTTGTGAGTCGAGAAAGATCACCATAAAGAtctc                                                                                                                                                                             >  2:32263/1‑251 (MQ=255)                                               cATGAACTGACAACTCTTTATTATCCGTTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTATGTGCAAGGATCAGCGCCGAGGCGTTTATTTTTATCGCTTCGCGGATAATTTCCCGAGGATGGACTTCAACATGGTTTAGCGTGCCGGAAAAAAGACGGCTATGCGTTATAACCCGGTGTTGTGAGTCGAGAAAGATCACCATAAAGATCTCCCGCTCCTCACCCGTg                                                                                                                                                             >  1:31970/1‑250 (MQ=255)                                                  gAACTGACAACTCTTTATTATCCGTTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTATGTGCAAGGATCAGCGCCGAGGCGTTTATTTTTATCGCTTCGCGGATAATTTCCCGAGGATGGACTTCAACATGGTTTAGCGTGCCGGAAAAAAGACGGCTATGCGTTATAACCCGGTGTTGTGAGTCGAGAAAGATCACCATAAAGATCTCCCGCTCCTCACCCGTGAGc                                                                                                                                                          <  2:118192/250‑1 (MQ=255)                                                                                    aGTTTATCCGCTTTACTGGGTTCAGCACATCCCGAAGGGTGATTATGTGCAAGGATCAGCGCCGAGGCGTTTATTTTTATCGCTTCGCGGATAATTTCCCGAGGATGGACTTCAACATGGTTTAGCGTGCCGGAAAAAAGACGGCTATGCGTTATAACCCGGTGTTGTGAGTCGAGAAAGATCACCATAAAGATCTCCCGCTCCTCACCCGTGAGCTGGCTTTGTAAAAATTCCCGCGTCATCTCCGGGCt                                                                                                                       <  1:410968/251‑1 (MQ=255)                                                                                                     tggTTCAGCACAACCCGAAGGTTGATTATGTGCAAGGATCAGCGCCGAGGCTTTTATTTTTATCGCTTCGCGGATAATTTCCCGAGGATGGACTTCAACATGGTTTAGCGTGCCGGAAAAAAGACGGCTATGCGTTATAACCCGGTGTTGTGAGTCGAGAAAGATCACCATAAAGATCTCCCGCTCCTCACCCGTGAGCTGGCTTTGTAAAAATTCCCGCGTCATCTCCGGGCTGAGTAAAGGGCTttc                                                                                                        <  2:237426/248‑1 (MQ=255)                                                                                                                     gAAGGGTGATTATTTGCCTGGATCAGCGCCGAGGCGTTTATTTTTATCGCTTCGCGGCTAATTTCCCGAGGTTGGCCTTCAACATGGTTTAGCGTGCCGGAAACAAGACGGCTATGCGTTATAACCCGGTGTTGTGAGTCGAGAAAGATCACCATAAAGATCTCCCGCTCCTCACCCGTGAGCTGGATTTGTAAAAATTCCCGCGTCATCTCCGGGCTGAGTAAAGGGCTTTCTTCACGCATCCGCACGtt                                                                                      >  2:181557/1‑251 (MQ=255)                                                                                                                         ggTGATTATGTGCAAGGATCAGCGCCGAGGCGTTTATTTTTATCGCTTCGCGGATAATTTCCCGAGGATGGACTTCAACATGGTTTAGCGTGCCGGAAAAAAGACGGCTATGCGTTATAACCCGGTGTTGTGAGTCGAGAAAGATCACCATAAAGATCTCCCGCTCCTCACCCGTGAGCTGGCTTTGTAAAAATTCCCGCGTCATCTCCGGGCTGAGTAAAGGGCTTTCTTCACGCATCCGCACGTTgtag                                                                                  <  1:181852/251‑1 (MQ=255)                                                                                                                              ttATGTGCAAGGATCAGCGCCGAGGCGTTTATTTTTATCGCTTCGCGGATAATTTCCCGAGGATGGACTTCAACATGGTTTAGCGTGCCGGAAAAAAGACGGCTATGCGTTATAACCCGGTGTTGTGAGTCGAGAAAGATCACCATAAAGATCTCCCGCTCCTCACCCGTGAGCTGGCTTTGTAAAAATTCCCGCGTCATCTCCGGGCTGAGTAAAGGGCTTTCTTCACGCATCCGCACGTTGTAGTAAcg                                                                             >  2:111213/1‑251 (MQ=255)                                                                                                                              ttATGTGCAAGGATCAGCGCCGAGGCGTTTATTTTTATCGCTTCGCGGATAATTTCCCGAGGATGGACTTCAACATGGTTTAGCGTGCCGGAAAAAAGACGGCTATGCGTTATAACCCGGTGTTGGGAGTCGAGAAAGATCACCATAAAGATCTCCCGCTCCTCACCCGTGAGCTGGCTTTGTAAAAATTCCCGCGTCATCTCCGGGCTGAGTAAAGGGCTTTCTTCACGCATCCGCACGTTGTAGTAAcg                                                                             >  1:309697/1‑251 (MQ=255)                                                                                                                                        ggATCAGCGCCGAGGCGTTTATTTTTATCGCTTCGCGGATAATTTCCCGAGGATGGACTTCAACATGGTTTAGCGTGCCGGAAAAAAGACGGCTATGCGTTATAACCCGGTGTTGTGAGTCGAGAAAGATCACCATAAAGATCTCCCGCGCCTCACCCGTGAGCTGGCTGTGTAAAAATTCCCGCGTCATCTCCGGGCTGAGTAAAGGGCTTTCTGCACGCATCCGCACGTTGTAGTAACGGCGCGCCAGt                                                                   <  1:145267/251‑1 (MQ=255)                                                                                                                                          aTCAGCGCCGAGGCGTTTATTTTTATCGCTTCGCGGATAATTTCCCGAGGATGGACTTCAACATGGTTTAGCGTGCCGGAAAAAAGACGGCTATGCGTTATAACCCGGTGTTGTGAGTCGAGAAAGATCACCATAAAGATCTCCCGCTCCTCACCCGTGAGCTGGCTTTGTAAAAATTCCCGCGTCATCTCCGGGCTGAGTAAAGGGCTTTCTTCACGCATCCGCACGTTGTAGTAACGTCGCGCCAGTTc                                                                 <  2:29045/251‑1 (MQ=255)                                                                                                                                             aGCGCCGCGGCGTTTATTTTTATCGCTTCGCGGATAATTTCCCGAGGATGGACTTCAACATGGTTTAGCGTGCCGGAAAAAAGACGGCTATGCGTTATAACCCGGTGTTGTGAGTCGAGAAAGATCACCATAAAGATCTCCCGCTCCTCACCCGTGAGCTGGCTTTGTAAAAATTCCCGCGTCATCTCCGGGCTGAGTAAAGGGCTTTCTTCACGCATCCGCACGTTGTAGTAACGTCGCGCCAGTTCAGc                                                              >  2:230655/1‑251 (MQ=255)                                                                                                                                                              ttttATCGCTTCGCGGATAATTTCCCGAGGATGGACTTCAACATGGTTTCGCGTGCCGGAAAAAAGACGGCTATGCGTTATAACCCGGTGTCGTGAGTCGAGAAAGATCACCATAAAGATCTCCCGCTCCTCACCCGTGAGCTGGCTTTGTAAAAATTCCCGCGTCATCTCCGGGCTGAGTAAAGGGCTTTCTTCACGCATCCGCACGTTGTAGTAACGTCGCGCCAGTTCAGCAATCCCCTTTAACTggg                                             >  1:146581/1‑251 (MQ=255)                                                                                                                                                                ttATCGCTTCGCGGATAATTTCCCGAGGATGGACTTCAACATGGTTTAGCGTGCCGGAAAAAAGACGGCTATGCGTTATAACCCGGTGTTGTGAGTCGAGAAAGATCACCATAAAGATCTCCCGCTCCTCACCCGTGAGCTGGCTTTGTAAAAATTCCCGCGTCATCTCCGGGCTGAGTAAAGGGCTTTCTTCACGCATCCGCACGTTGTAGTACCGTCGCGCCAGTTCAGCAATCCCCTTTAACTGGGc                                            >  2:337633/1‑250 (MQ=255)                                                                                                                                                                    cGCTTCGCGGATAATTTCCCGAGGATGGACTTCAACATGGTTTAGCGTGCCGGAAAAAAGACGGCTATGCGTTATAACCCGGTGTTGTGAGTCGAGAAAGATCACCATAAAGATCTCCCGCTCCTCACCCGTGAGCTGGCTTTGTAAAAATTACCGCGTCATCTCTGGGCCGAGTAAAGGGCTTTCTACACGCATCCGCACGTTGTAGTAACGTCGCGCCAGTTCAGCAATCCCCTTTAACTGGGCAAAtt                                       >  1:135618/1‑251 (MQ=255)                                                                                                                                                                      cTTCGCGGATAATTTCCCGAGGAAGGACTTCAACATGGGTTAGCGTGCCGGAAAAATGACGGCTATGCGTTATAACCCGGTGTTGTGAGTCGAGAAAGATCACCATAAAGATCTCCCGCTCCTCACCCGTGAGCTGGCTTTGTAAAAATTCCCGCGTCATCTCCGGGCTGAGTAAAGGGCTTTCTTCATGAATCCGCACGTTGTAGTAACGTCGCGACAGTTCAGCAATCCCCTTTAACTGGGCAAATTTc                                     >  2:97381/1‑251 (MQ=255)                                                                                                                                                                       ttCGCGGATCATTTTCCGAGGATGGAATTCAGCATGGTTTAGTTTGCAGGAAAAAAGACGGCTATGCGTTATAACCCAGTGTTGTGAGTCGAGAAAGATCAACATACAGATCTCCCGCTCCTCACCCGTGAGCTGGCTTTGTAATAATTCACGCGTCATCTCCGGGATGAGTACAGGGCTTTCTTCACGCCTCCGCAAGTTGTAGTAACGTCGCGCCAGTTCAGCAATACCCCTTAACTGGGCAAATTTcg                                    >  2:176346/1‑251 (MQ=255)                                                                                                                                                                                     cccGAGGATGGACTTCAACATGGTTTAGCGTGCCGGAAAAAAGACGGCTATGCGTTATAACCCGGTGTTGTGAGTCGAGAAAGATCACCATAAAGATCTCCCGCTCCTCACCCGTGAGCTGGCTTTGTAAAAATTCCCGCGTCATCTCCGGGCTGAGTAAAGGGCTTTCTTCACGCATCCGCACGTTGTAGTAACGTCGCGCCAGTTCAGCAATCCCCTTTAACTGGGCAAATTTCGCTACGCCAATTcc                       >  1:321576/1‑250 (MQ=255)                                                                                                                                                                                      ccGAGGATGGACTTCAACATGGGTGAGCGTGCCGGAAAAAAGACGGCTATGCGTTATAACCCGGTGTTGGGAGTCGAGAAAGATCACCATAAAGATCTCCCGCGCCTCACCCGTGAGCTGGCTTTGTAAAAATTCCCGCGTCATCTCCGGGCTGAGTAAAGGGCTGTCTTCACGCAGCCGCACGTTGTAGTAACGTCGCGCCAGTTCAGCAATCCCCTTTAACTGGGCAAATTTCGCTACGCCAATTCCAt                     <  1:181557/251‑1 (MQ=255)                                                                                                                                                                                                          gggttAGCGTGCCGGCAACAAACCGGCTATGCGTTATACCCCGGAGTTGTGCGTCGAGAAAGCTAAACATAACGATCTACCGCTCCTCACCCGTGAGCTCGCTTTGTAAAAATTCCCGCGTCATCTCCGGGCTGAGTAAAGGGCTTTCTTCACGAATCCGCACGTTGTCGTAACGTCGCGCCAGTTCAGCAATCCCCTTTAACTGGGCAAATTTCGCTACGCCAATTCCATGAACACCACTAAAATGTTc  >  2:268701/4‑250 (MQ=255)                                                                                                                                                                                                          ggTTTAGCGTGCCGGAAAAAAGACGGCTATGCGTTATAACCCGGTGTTGTGAGTCGAGAAAGATCACCATAAAGATCTCCCGCTCCTCACCCGTGAGCTGGCTTTGTAAAAATTCCCGCGTCATCTCCGGGCTGAGTAAAGGGCTTTCTTCACGCATCCGCACGTTGTAGTAACGTCGCGCCAGTTCAGCAATCCCCTTTAACTGGGCAAATTTCGCTACGCCAATTCCATGAACACCACTAAATTGTTc  >  1:143753/1‑250 (MQ=255)                                                                                                                                                                                                                      cGGAAAAAAGACGGCTATGCGTTATAACCCGGTGTTGTGAGTCGAGAAAGATCACCATAAAGATCTCCCGCTCCTCACCCGTGAGCTGGCTTTGTAAAAATTCCCGCGTCATCTCCGGGCTGAGTAAAGGGCTTTCTTCACGCATCCGCACGTTGTAGTAACGTCGCGCCAGTTCAGCAATCCCCTTTAACTGGGCAAATTTCGCTACGCCAATTc                        >  2:157397/1‑216 (MQ=255)                                                                                                                                                                                                                      cGGAAAAAAGACGGCTATGCGTTATAACCCGGTGTTGTGAGTCGAGAAAGATCACCATAAAGATCTCCCGCTCCTCACCCGTGAGCTGGCTTTGTAAAAATTCCCGCGTCATCTCCGGGCTGAGTAAAGGGCTTTCTTCACGCATCCGCACGTTGTAGTAACGTCGCGCCAGTTCAGCAATCCCCTTTAACTGGGCAAATTTCGCTACGCCAATTc                        <  1:157397/216‑1 (MQ=255)                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                 CATACTCTCCACGCCCAATCACGATATGGTCGAGCACGCGTAAATCCATGAACTGACAACTCTTTATTATCCGTTCAGTAATGAGTTTATCCGCTTTACTGGGTTCAGCACAACCCGAAGGGTGATTATGTGCAAGGATCAGCGCCGAGGCGTTTATTTTTATCGCTTCGCGGATAATTTCCCGAGGATGGACTTCAACATGGTTTAGCGTGCCGGAAAAAAGACGCCTATGCGTTATAACCCGGTGTTGGGAGTCGAGAAAGATCACCATAAAGATCTCCCGCTCCTCACCCGTGAGCTGGCTTTGTAAAAATTCCCTCGTCATCTCCGGGCTGAGTAAAGGGCTTTCTTCACGCATCCGCACGTTGTAGTAACGCCGCGCCAGTTCAGCAATCCCCTTTAACTGGGCAAATTTCGCCACGCCAATTCCATGAACGCCACTAAATTGTTC  >  NC_000913/3811920‑3812370

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