Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	313,851	T→C	6.9%	intergenic (+22/+29)	ygaU → / ← yrfG	murein cross‑linking regulator/purine nucleotidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	313,851	0	T	C	6.9%	169.2 / 2.5	58	intergenic (+22/+29)	ygaU/yrfG	murein cross‑linking regulator/purine nucleotidase
Reads supporting (aligned to +/- strand):  ref base T (29/25);  new base C (2/2);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

TTGCCATGGTCGTGTATGGCAATGCCAATCAGTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGT  >  NC_002947/313711‑314000                                                                                                                                             |                                                                                                                                                       ttGCCATGGTCGTGTATGGCAATGCCAATCAGTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTC‑TTcgc                                                                                                                                                 <  1:679777/147‑1 (MQ=255) ttGCCATGGTCGTGTATGGCAATGCCAATCAGTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTC‑TTcgc                                                                                                                                                 >  2:679777/1‑147 (MQ=255)      aTGGTCGTGTATGGCAATGCCAATCAGTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGgcgc                                                                                                                                           >  1:1039312/1‑148 (MQ=255)      aTGGTCGTGTATGGCAATGCCAATCAGTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGgcgc                                                                                                                                           <  2:1039312/148‑1 (MQ=255)                  ggCAATGCCAATCAGTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATgggg                                                                                                                             <  2:625385/150‑1 (MQ=255)                  ggCAATGCCAATCAGTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATgggg                                                                                                                             <  2:1205800/150‑1 (MQ=255)                    cAATGCCAATCAGTACAACAAGATTTTCGAGGCCACCAAGGAGATTGCCAAGCAGGAGGACAAGATCTATCCTGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTC‑TGCGCTGGCGCGCCCGCGAATGGGGcc                                                                                                                           >  1:209727/1‑150 (MQ=255)                         ccAATCAGTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACa                                                                                                                      >  1:708058/1‑150 (MQ=255)                           aaTCAGTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAgg                                                                                                                    <  1:567722/150‑1 (MQ=255)                               aGTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTa                                                                                                                <  1:637523/150‑1 (MQ=255)                                    aacaaGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGt                                                                                                           >  2:1127599/1‑150 (MQ=255)                                    aacaaGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGt                                                                                                           >  1:303336/1‑150 (MQ=255)                                     acaaGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTa                                                                                                          <  2:118039/150‑1 (MQ=255)                                       ctgttcatCGAGGCCAACAAGCCGATGCTCAAACACGCGGACAAGATCTATCCGGGGCAGGTGCTGAGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATAGGGCCTTACAGGCCTACCAGTAAc                                                                                                        >  2:199799/8‑150 (MQ=255)                                       aaGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTACTCCTTAGTTATCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAAc                                                                                                        >  1:749492/1‑150 (MQ=255)                                         gATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCATGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTc                                                                                                      >  2:1082099/1‑150 (MQ=255)                                             ttCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCATGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCgg                                                                                                  <  1:1082099/150‑1 (MQ=255)                                                   gCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGcc                                                                                                                           >  1:1098278/1‑119 (MQ=255)                                                   gCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGcc                                                                                                                           <  2:1098278/119‑1 (MQ=255)                                                   gCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCCCATTCGCGGGCGCGCCCGCGaaga                                                                                                                                >  2:891244/1‑113 (MQ=39)                                                   gCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCCCATTCGCGGGCGCGCCCGCGaaga                                                                                                                                <  1:891244/115‑3 (MQ=39)                                                     caacaaGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCgg                                                                                                  >  1:1074831/1‑142 (MQ=255)                                                     caacaaGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCgg                                                                                                  <  2:1074831/142‑1 (MQ=255)                                                      aacaaGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCg                                                                                         >  2:871692/1‑150 (MQ=255)                                                              gATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATGCCTTACTGATCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTGACAGGCCTACCAGTAACTCCCGGTAATCCCCGACcgccgc                                                                                 >  1:90705/1‑150 (MQ=255)                                                               aTGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCa                                                                                >  1:1120370/1‑150 (MQ=255)                                                                 gCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCaaa                                                                              <  2:1120370/150‑1 (MQ=255)                                                                      aaCAGCCGGACAAGATATTACGGAGTCAGGTGCGCCGTATGCCTGACTGATCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCCCGCCCGCGCATGGGGCCTTACAGGCCTACCGGTAACTCCCGGTAATCCTCGACGGCCGCAAACTCtt                                                                         <  2:90705/150‑1 (MQ=255)                                                                      aaCACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCtt                                                                         >  1:724915/1‑150 (MQ=255)                                                                         acCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATNGCTGAAGTCTGTACAGCCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAg                                                                      >  1:752133/1‑150 (MQ=255)                                                                           ccGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTa                                                                    >  2:989853/1‑150 (MQ=255)                                                                                   gATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCgg                                                            <  1:989853/150‑1 (MQ=255)                                                                                           cGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTC‑TTCGCAGGCGCGCCCCCGCCGATGACCTGACAGGCCGACCAGTATCTCCCGGTAATCCTCGACACCCGCAGACTCCGCGTGATCCCCCGGCCCTGCtt                                                    >  1:967627/1‑150 (MQ=255)                                                                                               gCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTa                                                                                                                        <  2:956954/78‑1 (MQ=255)                                                                                               gCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTa                                                                                                                        >  1:956954/1‑78 (MQ=255)                                                                                                      ctgcGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCgg                                                                                                  <  1:984228/93‑1 (MQ=255)                                                                                                      ctgcGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCgg                                                                                                  >  2:984228/1‑93 (MQ=255)                                                                                                         cGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTaa                                                                                               <  1:520396/93‑1 (MQ=255)                                                                                                         cGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTaa                                                                                               >  2:520396/1‑93 (MQ=255)                                                                                                              tCCTGACTGATCGCTGAAGTCTGTACAGGCCTC‑TTc                                                                                                                                                   <  1:39732/36‑1 (MQ=255)                                                                                                              tCCTGACTGATCGCTGAAGTCTGTACAGGCCTC‑TTc                                                                                                                                                   >  2:39732/1‑36 (MQ=255)                                                                                                                ctgactgaTCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCg                               >  2:421600/1‑150 (MQ=255)                                                                                                                   atctaTCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCGAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCa                            <  2:516074/146‑1 (MQ=255)                                                                                                                     tgaTCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCt                                        >  1:661224/1‑136 (MQ=255)                                                                                                                     tgaTCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCt                                        <  2:661224/136‑1 (MQ=255)                                                                                                                      gaTCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGcc                         <  1:1122575/150‑2 (MQ=255)                                                                                                                        tCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGAccg                                                                                     >  1:26838/1‑88 (MQ=255)                                                                                                                        tCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGAccg                                                                                     <  2:26934/88‑1 (MQ=255)                                                                                                                        tCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGAccg                                                                                     <  2:26838/88‑1 (MQ=255)                                                                                                                        tCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGAccg                                                                                     >  1:26934/1‑88 (MQ=255)                                                                                                                        tCGCTGAAGTCTGTACAGGCATC‑TTCGCGGTCTCGCCCGCGAATGGGGCCTTACAGGACTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCGCCGCCAgagat                       >  1:516074/1‑146 (MQ=255)                                                                                                                            tGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGcc                   <  2:303336/150‑1 (MQ=255)                                                                                                                                    gTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCa           >  1:37811/1‑150 (MQ=255)                                                                                                                                      aCAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCt                                        >  1:1029112/1‑119 (MQ=255)                                                                                                                                      aCAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCt                                        <  2:1029112/119‑1 (MQ=255)                                                                                                                                      aCAGGCCCC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCaaa         >  1:902502/1‑150 (MQ=255)                                                                                                                                       cAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAAc        <  2:752133/150‑1 (MQ=255)                                                                                                                                          gCCCC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAAcagc     <  2:902502/150‑1 (MQ=255)                                                                                                                                             tc‑tTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGt  >  2:1005175/1‑150 (MQ=255)                                                                                                                                             |                                                                                                                                                       TTGCCATGGTCGTGTATGGCAATGCCAATCAGTACAACAAGATTTTCGAGGCCAACAAGCCGATGCTCAAACACCCGGACAAGATCTATCCGGGGCAGGTGCTGCGTATTCCTGACTGATCGCTGAAGTCTGTACAGGCCTC‑TTCGCGGGCGCGCCCGCGAATGGGGCCTTACAGGCCTACCAGTAACTCCCGGTAATCCTCGACCGCCGCAAACTCTTCAGTATCCCGCGGCCCTGCTTTGCTATCCGGCTGGCGCACCGCCAGCAAGTGCCCCACACCAAACCGCCGT  >  NC_002947/313711‑314000

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