Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,545,353:1	+G	8.6%	coding (497/1284 nt)	hyuC ←	bifunctional N‑carbamoyl‑beta‑alanine amidohydrolase/allantoine amidohydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,545,353	1	.	G	8.6%	220.9 / 15.4	58	coding (497/1284 nt)	hyuC	bifunctional N‑carbamoyl‑beta‑alanine amidohydrolase/allantoine amidohydrolase
Reads supporting (aligned to +/- strand):  ref base . (27/26);  new base G (2/3);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.26e-01

TTGGCCTGGTCCTCAAGGATCGGCCCCTGCTCGATATGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTGGTCCACACCACCACTTCGAGCGACGCTTCGGTTTCCACGCCCAGGTCGTTGAG  >  NC_002947/4545218‑4545498                                                                                                                                         |                                                                                                                                                  tgggCCTGATCCTCAAGGATCGGCCCCTGCTCGATATGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGGAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGAAT‑CACGCTTGGCCAgg                                                                                                                                     >  2:161543/3‑150 (MQ=255)  tGGCCTGGTCCTCAAGGATCGGCCCCTGCTCGATATGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAggg                                                                                                                                    <  2:1367579/150‑1 (MQ=255)        ggTCCTCAAGGATCGGCCCCTGCTCGATATGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCtt                                                                                                                              >  1:722121/1‑150 (MQ=255)         gTCCTCAAGGATCGGCCCCTGCTCGATATGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTc                                                                                                                             >  2:1003347/1‑150 (MQ=255)           ccTCAAGGATCGGCCCCTGCTCGATATGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCa                                                                                                                           <  1:1003347/150‑1 (MQ=255)             tCAAGGATCGGCCCCTGCTCGATATGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGTACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAgg                                                                                                                         >  2:525526/1‑150 (MQ=255)             tCAAGGATCGGCCCCTGCTCGATATGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAgg                                                                                                                         <  1:820601/150‑1 (MQ=255)             tCAAGGATCGGCCCCTGCTCGATATGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAgg                                                                                                                         <  1:601078/150‑1 (MQ=255)              cAAGGATCGGCCCCTGCTCGATATGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAggg                                                                                                                        >  1:985510/1‑150 (MQ=255)                 ggATCGGCCCCTGCTCGATATGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAGGGTGa                                                                                                                     >  1:480782/1‑150 (MQ=255)                       gCCCCTGCTCGATATGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGc                                                                                                               >  1:51099/1‑150 (MQ=255)                       gCCCCTGCTCGATATGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGc                                                                                                               >  1:498341/1‑150 (MQ=255)                            tGCTCGATATGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCgg                                                                                                            >  1:1264950/1‑148 (MQ=255)                            tGCTCGATATGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCgg                                                                                                            <  2:1264950/148‑1 (MQ=255)                               tcgaagtggtggtgtggaccaacgaagaaggttcgcgttttgcgccatgcatgatggggtctggggtattcgccggcaagttcaccctggaaGACACCCTGGGCAT‑CGCGCTTGGCCAGGgt                                                                                                                                   >  2:69006‑M1/93‑122 (MQ=255)                               tcgaagtggtggtgtggaccaacgaagaaggttcgcgttttgcgccatgcatgatggggtctggggtattcgccggcaagttcaccctggaaGACACCCTGGGCAT‑CGCGCTTGGCCAGGgt                                                                                                                                   <  1:69006‑M1/30‑1 (MQ=255)                                 gATATGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAgg                                                                                                                                     >  2:317159/1‑118 (MQ=255)                                 gATATGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAgg                                                                                                                                     <  1:317159/118‑1 (MQ=255)                                   tatGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATAccc                                                                                                   <  2:51099/150‑1 (MQ=255)                                     tGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCa                                                                                                 <  1:406405/150‑1 (MQ=255)                                     tGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCa                                                                                                 <  2:722121/150‑1 (MQ=255)                                          ttCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATAccc                                                                                                   <  1:1090749/143‑1 (MQ=255)                                          ttCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATAccc                                                                                                   >  2:1090749/1‑143 (MQ=255)                                          ttCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGAccc                                                                                            >  2:734286/1‑150 (MQ=255)                                                 tACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATcatg                                                                                     >  1:1248640/1‑150 (MQ=255)                                                    gCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATcatgcat                                                                                  >  1:1370800/1‑150 (MQ=255)                                                      ccccACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCTCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATcatgcatg                                                                                 >  1:754694/1‑149 (MQ=255)                                                      ccccACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCTCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATcatgcatg                                                                                 <  2:754694/149‑1 (MQ=255)                                                      ccccACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATgg                                                                                >  1:129580/1‑150 (MQ=255)                                                                    ccGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAgg                                                                                                                         >  2:231306/1‑95 (MQ=255)                                                                    ccGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAgg                                                                                                                         <  1:231306/95‑1 (MQ=255)                                                                     cGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCa                                                                                                                           >  2:451866/1‑92 (MQ=255)                                                                     cGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCa                                                                                                                           <  1:451866/92‑1 (MQ=255)                                                                              aCGAGGACCGGCGTAGCCAATGGAGTGCAACGCTTCGCCAACGCTGACACCCTGAGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGCCCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGtt                                                        <  1:655116/150‑1 (MQ=255)                                                                               cgcgGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTg                                                       <  2:480782/150‑1 (MQ=255)                                                                                gcgGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAggg                                                                                                                                    <  1:1147700/72‑1 (MQ=255)                                                                                gcgGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAggg                                                                                                                                    >  2:1147700/1‑72 (MQ=255)                                                                                      cGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTGGTc                                                    <  1:1032423/146‑1 (MQ=255)                                                                                      cGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTGGTc                                                    >  2:1032423/1‑146 (MQ=255)                                                                                        gCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTGGTCCAcacc                                              >  2:1264651/1‑150 (MQ=255)                                                                                         cGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTGGTCCAcacca                                             <  2:129580/150‑1 (MQ=255)                                                                                          gTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGTCGCAAAACGCGAACCTTCTTCGTTGGTCCAcaccac                                            <  2:985510/150‑1 (MQ=255)                                                                                             gCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGaa                                                                                                        >  2:1195735/1‑87 (MQ=255)                                                                                             gCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGaa                                                                                                        <  1:1195735/87‑1 (MQ=255)                                                                                                 atgtggTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCGTCATGCATGGCGCAAAACGCGAACCTTCTTCGTTGGTCCACACCACCACTTCg                                     >  2:1066355/6‑150 (MQ=255)                                                                                                  tGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTGGTCCACACCACCACTTCGa                                    <  1:91336/150‑1 (MQ=255)                                                                                                  tGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTGGTCCACACCACCACTTCGa                                    <  1:620236/150‑1 (MQ=255)                                                                                                              gttcaccctggaaGACACCCTGGCCAAGCGCGAT‑GCCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATcatgcat                                                                                  >  1:535795‑M1/14‑92 (MQ=255)                                                                                                              gttcaccctggaaGACACCCTGGCCAAGCGCGAT‑GCCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATcatgcat                                                                                  <  2:535795‑M1/79‑1 (MQ=255)                                                                                                               ttCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTGGTCCacac                                               <  2:1032840/126‑1 (MQ=255)                                                                                                               ttCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTGGTCCacac                                               >  1:1032840/1‑126 (MQ=255)                                                                                                                      tggaaGACACCCTGGCCAAGCGCGAT‑GCCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTg                                                       >  2:476964‑M1/6‑111 (MQ=255)                                                                                                                      tggaaGACACCCTGGCCAAGCGCGAT‑GCCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTg                                                       <  1:476964‑M1/106‑1 (MQ=255)                                                                                                                      tggaaGACACCCTGGCCAAGCGCGAT‑GCCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTGGTCCACACCACCACTTCGAGCGACGCTTCGGTTTCCACg                <  2:182926‑M1/145‑1 (MQ=255)                                                                                                                       cGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTGGTCCACACCACCACTTCGAGCGACGCTTCGGTTTCCACGc               >  1:1103206/1‑150 (MQ=255)                                                                                                                              ACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTGGTCCACACCACCACTTCGCCCGACGCTTCGGTTTCTACCCCCAGGTc        <  2:139348‑M1/150‑1 (MQ=255)                                                                                                                               CCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTGGTCCACACCACCACTTCGAGCTACGCTTCGGTTTCCACGCCCAGGTCg       <  2:1008308‑M1/150‑1 (MQ=255)                                                                                                                                    GGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTGGTCCACACCACCACTTCGAGCGACGCTTCGGTTTCCACGCCCAGGTCGTTGAg  >  1:734286‑M1/1‑150 (MQ=255)                                                                                                                                         |                                                                                                                                                  TTGGCCTGGTCCTCAAGGATCGGCCCCTGCTCGATATGCGCTTCGAAATACGCCCCCACCGGGTGCCCCGACACAGCACGCGGACCGGCGTAGCCAATGGAGTTCAACGCTTCGCCAACGCTGACACCCTGGGCAT‑CGCGCTTGGCCAGGGTGTCTTCCAGGGTGAACTTGCCGGCGAATACCCCAGACCCCATCATGCATGGCGCAAAACGCGAACCTTCTTCGTTGGTCCACACCACCACTTCGAGCGACGCTTCGGTTTCCACGCCCAGGTCGTTGAG  >  NC_002947/4545218‑4545498

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