Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,194,888	C→A	21.4%	intergenic (+30/+23)	PP_4574 → / ← PP_4575	hypothetical protein/urea amidolyase‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,194,888	0	C	A	21.4%	121.9 / 21.3	55	intergenic (+30/+23)	PP_4574/PP_4575	hypothetical protein/urea amidolyase‑like protein
Reads supporting (aligned to +/- strand):  ref base C (21/22);  new base A (7/5);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 7.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGCCTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGT‑‑ATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCGCCAGCGCCAGCGGCGTCAG  >  NC_002947/5194745‑5195033                                                                                                                                                  |                                                                                                                                                  cAGGCCTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTGGATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGcc                                                                                                                                                >  1:103247/1‑149 (MQ=255)    gCCTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGT‑‑ATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCGAGCCcgcg                                                                                                                                            >  4:348936/1‑148 (MQ=255)      cTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGT‑‑ATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGa                                                                                                                                           >  5:365597/1‑147 (MQ=255)      cTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGT‑‑ATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGa                                                                                                                                           >  6:290415/1‑147 (MQ=255)      cTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGT‑‑ATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCGAGCCCGCGaaa                                                                                                                                         <  3:348936/149‑1 (MQ=255)      cTGGGACTACACCGCCTCGACCCATCTGGTGCGCGACCTGAACGAGGACAAGATGT‑‑ATCCGCTGACCCTGCAAGGGTTGGAGGACGTGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGaaa                                                                                                                                         <  2:256093/149‑1 (MQ=255)         ggAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGT‑‑ATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGGTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAg                                                                                                                                        <  7:128751/147‑1 (MQ=255)           aGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGT‑‑ATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGcc                                                                                                                                    <  5:290415/149‑1 (MQ=255)             tACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGT‑‑ATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCcag                                                                                                                                  >  3:133295/1‑149 (MQ=255)              acacCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGT‑‑ATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCcag                                                                                                                                  >  5:129057/1‑148 (MQ=255)                acCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGT‑‑ATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCcagcag                                                                                                                               >  8:46893/1‑149 (MQ=255)                          cccATCTGGTGCGCGACCTGATCGAGGACAAGATGT‑‑ATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGAcc                                                                                                                     <  3:30445/149‑1 (MQ=255)                           ccATCTGGTGCGCGACCTGATCGAGGACAAGATGT‑‑ATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGAcc                                                                                                                     >  5:9159/1‑148 (MQ=255)                           ccATCTGGTGCGCGACCTGATCGAGGACAAGATGT‑‑ATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGAcc                                                                                                                     >  6:171148/1‑148 (MQ=255)                           ccATCTGGTGCGCGACCTGATCGAGGACAAGAAGA‑‑ATCCGCTCACCCTGCAAGGGATGGAGCACGGGATGCGCTACCTGTCCAAGTCAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCtaaagcac                                                                                                                     >  7:208545/1‑141 (MQ=255)                             aTCTGGTGCGCGACCTGATCGAGGACAAGATGT‑‑ATCCGCTGACCCTGCAAGGGATGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCATCAGTCTTACGACCacc                                                                                                                  >  5:360676/1‑147 (MQ=255)                                tGGTGCGCGACCTGATCGAGGACAAGATGT‑‑ATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCGAGCCCGCGaaa                                                                                                                                         >  3:245182/1‑123 (MQ=255)                                tGGTGCGCGACCTGATCGAGGACAAGATGT‑‑ATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCGAGCCCGCGaaa                                                                                                                                         <  4:245182/123‑1 (MQ=255)                                      gcgACCTGATCGAGGACAAGATGT‑‑ATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCGTGTACGGGCCCTTTCGCGGGCGAGCCCGGGaa                                                                                                                                          >  2:154235/1‑116 (MQ=255)                                      gcgACCTGATCGAGGACAAGATGT‑‑ATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGt                                                                                                         <  4:133295/149‑1 (MQ=255)                                         aCCTGATCGAGGACAAGATGT‑‑ATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTaaa                                                                                                      <  7:46893/149‑1 (MQ=255)                                              aTCGAGGACAAGATGT‑‑ATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAgctg                                                                                                  <  2:103247/148‑1 (MQ=255)                                                    gACAAGATGT‑‑ATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGc                                                                                                                                     <  7:69145/107‑1 (MQ=255)                                                    gACAAGATGT‑‑ATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGc                                                                                                                                     >  8:69145/1‑107 (MQ=255)                                                    gACAAGATGT‑‑ATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAAc                                                                                                                  <  4:68985/126‑1 (MQ=255)                                                    gACAAGATGT‑‑ATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAAc                                                                                                                  >  3:68985/1‑126 (MQ=255)                                                    gACAAGATGT‑‑ATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGctct                                                                                           <  4:37784/149‑1 (MQ=255)                                                    gACAAGATGT‑‑ATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGctct                                                                                           >  5:318129/1‑149 (MQ=255)                                                             t‑‑ATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGcc                                                                                  <  6:129057/149‑1 (MQ=255)                                                             t‑‑ATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGc                                                                                   <  6:318129/148‑1 (MQ=255)                                                                    gCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTc                                                                             <  6:365597/149‑1 (MQ=255)                                                                    gCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTACAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGACCGAGAAAGGGCCAGCAGTCGTACGACCAACGCCCAAGGTAAAGCTGCTGCTCACGCCAGGCCTCGTc                                                                             >  3:279273/1‑149 (MQ=255)                                                                        aCCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTg                                                                          <  5:171148/148‑1 (MQ=255)                                                                           cTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGcacc                                                                      <  6:9159/149‑1 (MQ=255)                                                                            tGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGACTCGCCCGCGAAAGGGCCGGCAGTCTTACGACCAACGCCCAAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGcacca                                                                     >  6:65218/1‑149 (MQ=255)                                                                                              ggATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGcaccac                                                                    >  5:331808/1‑132 (MQ=255)                                                                                              ggATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGcaccac                                                                    <  6:331808/132‑1 (MQ=255)                                                                                               gATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGcac                                                     >  1:221404/1‑146 (MQ=255)                                                                                               gATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGcac                                                     <  2:221404/146‑1 (MQ=255)                                                                                                                aGTTAGGTTGCCGTTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCTAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGGATGctct                                 <  8:208545/149‑1 (MQ=255)                                                                                                                      gTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGc                                                                         <  4:103231/103‑1 (MQ=255)                                                                                                                      gTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGc                                                                         >  3:103231/1‑103 (MQ=255)                                                                                                                          ccTGTACGGGCCCTTTCTCGGGCTCGCCCTCGAAAGGGGCATCAGTCTTACGACCAACTCGCCAGGTAAAGCTGCTGCTCTCGGCAGGCCTCGTCCTGCACCAACGCCCTGAACCTCACAACCGTCCCCGGCATGCTCTGCGCCAttt                        <  4:279273/148‑1 (MQ=255)                                                                                                                          ccTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCCGTCGTACGCCCCACGCCCCAGGTAAAGCTGCGGCTCTCGCCCGGCCGCGTCCTGCACCCCCGCCCTGCCCCGCCCCCCCGGCCCCGGCCTGCACTGCGGCAGGTg                       >  4:228624/1‑149 (MQ=255)                                                                                                                                   gCCCTTTCGCGGGCGAGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAGGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCt                                                              >  7:249287/1‑101 (MQ=255)                                                                                                                                   gCCCTTTCGCGGGCGAGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAACGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCt                                                              <  8:249287/101‑1 (MQ=255)                                                                                                                                     ccTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTgcgccagcgcca            >  6:105647/1‑149 (MQ=255)                                                                                                                                     ccTTTCGCGGACTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCAAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTgcgccagcgcca            <  5:65218/149‑1 (MQ=255)                                                                                                                                       tttCGCGGGCGAGCCCGCGAAAGGGCCGGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCgccagcgccagc          >  1:340510/1‑149 (MQ=255)                                                                                                                                       tttCGCGGGCGAGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTccc                                            <  6:52274/115‑1 (MQ=255)                                                                                                                                       tttCGCGGGCGAGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTccc                                            >  5:52274/1‑115 (MQ=255)                                                                                                                                       tatCGCGGGCGAGCCCGCGAAAGGGCCGGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGca                                                      >  1:277837/3‑105 (MQ=39)                                                                                                                                       tatCGCGGGCGAGCCCGCGAAAGGGCCGGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGca                                                      <  2:277837/103‑1 (MQ=39)                                                                                                                                         tCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCGCCAGCGCCAgcgg        >  2:137434/1‑149 (MQ=255)                                                                                                                                               gCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCGCCAGCGCCAGCGGCGTCAg  <  1:137434/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                  CAGGCCTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGT‑‑ATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGTACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCGCCAGCGCCAGCGGCGTCAG  >  NC_002947/5194745‑5195033

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