Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,194,887	T→G	13.0%	intergenic (+29/+24)	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,887	0	T	G	13.0%	130.7 / 12.5	53	intergenic (+29/+24)	PP_4574/PP_4575	hypothetical protein/urea amidolyase‑like protein
Reads supporting (aligned to +/- strand):  ref base T (25/21);  new base G (5/2);  total (30/23)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.00e-01

CCAAGCAGGCCTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCGCCAGCGCCAGCGGCGTCA  >  NC_002947/5194740‑5195032                                                                                                                                                     |                                                                                                                                                  ccAAGCAGGCCTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTc                                                                                                                                                  >  2:33446/1‑149 (MQ=255)       aGGCCTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCGAGCCcgc                                                                                                                                            >  4:305242/1‑149 (MQ=255)           cTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGCGa                                                                                                                                          >  1:296670/1‑147 (MQ=255)           cTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGCGa                                                                                                                                          >  5:288414/1‑147 (MQ=255)           cTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCGAGCCCGCGaaa                                                                                                                                        <  3:305242/149‑1 (MQ=255)              ggAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTGCTACCGGCCCTTTCGCGGGCTCGCCCGCGaaa                                                                                                                                        <  1:33446/147‑1 (MQ=255)                     acCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCcagcag                                                                                                                              >  7:42467/1‑149 (MQ=255)                          cTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTa                                                                                                                         >  2:133869/1‑149 (MQ=255)                             gACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCCGCAGTCTTAcac                                                                                                                      >  7:272720/1‑147 (MQ=255)                               cccATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGAc                                                                                                                     >  6:182339/1‑148 (MQ=255)                                ccATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGAcc                                                                                                                    >  4:136659/1‑148 (MQ=255)                                  accTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAAc                                                                                                                 <  5:279648/147‑1 (MQ=255)                                  aTCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGTATGACGGGATGCGCTACCTGTCCAAGTTAGGTTGCCTG‑TACCGGCCCTATCGCGGGCTCGCCAGCGAAATGGCCATCATTCTTACGaacaac                                                                                                                 >  2:303836/1‑149 (MQ=255)                                       gTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGcccc                                                                                                            <  5:182339/149‑1 (MQ=255)                                          cgcgACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAgg                                                                                                         <  6:288414/149‑1 (MQ=255)                                               ccTTATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGATGAGGGGATGCTCTACCTGTCCAAGTGAGGTTGCATG‑TACCGGCCCTTTCGCGTGCTCGCCCTCGAAAGGTCCAGCAGTCTTACGACCAACGCCCCAGGTaaa                                                                                                     <  1:303836/148‑1 (MQ=255)                                               ccTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAg                                                                                                    >  5:13650/1‑149 (MQ=255)                                                cTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAgc                                                                                                   <  1:133869/149‑1 (MQ=255)                                                   aTCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAgctgc                                                                                                <  3:136659/149‑1 (MQ=255)                                                   aTCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAgct                                                                                                  <  8:42467/147‑1 (MQ=255)                                                    tCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAgctgct                                                                                               >  1:168097/1‑149 (MQ=255)                                                       aGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAgctgctgc                                                                                             >  8:299078/1‑148 (MQ=255)                                                          aCAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCtctc                                                                                         <  7:299078/149‑1 (MQ=255)                                                              gATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCa                                                                                     <  4:274177/149‑1 (MQ=255)                                                               aTGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAg                                                                                    >  5:83786/1‑149 (MQ=255)                                                                    tCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTAc                                                                                                                        <  4:56139/108‑1 (MQ=255)                                                                    tCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTAc                                                                                                                        >  3:56139/1‑108 (MQ=255)                                                                    tCCGCTGACCCCGCAAGGGTTGGAGGCCGGGATGCGCTACCAGTCCAACTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTc                                                                               <  8:272720/149‑1 (MQ=255)                                                                        cTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTcc                                                                           <  2:168097/149‑1 (MQ=255)                                                                             ccTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAgctgctg                                                                                              <  6:226645/125‑1 (MQ=255)                                                                             ccTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGAGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAgctgctg                                                                                              >  5:226645/1‑125 (MQ=255)                                                                                 cAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCaccacc                                                                  >  5:3604/1‑149 (MQ=255)                                                                                     ggTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGccc                                                              >  1:247293/1‑149 (MQ=255)                                                                                           aggaCGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGca                                                                       <  8:91133/134‑1 (MQ=255)                                                                                           aggaCGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGca                                                                       >  7:91133/1‑134 (MQ=255)                                                                                               cGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGcac                                                    <  6:83786/149‑1 (MQ=255)                                                                                                gggATGCGCTACCTGTCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCATCGTCCTGCACCACCGCCCTGAACCGcac                                                    >  5:64486/1‑148 (MQ=255)                                                                                                       gcTACCTGTCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGcacca                                                                    <  3:293521/125‑1 (MQ=255)                                                                                                       gcTACCTGTCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGcacca                                                                    >  4:293521/1‑125 (MQ=255)                                                                                                           ccTGTCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCACCgg                                        >  5:186361/1‑149 (MQ=255)                                                                                                               tCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATg                                    >  8:35191/1‑149 (MQ=255)                                                                                                                  aaGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCAc                                 <  2:247293/149‑1 (MQ=255)                                                                                                                    gTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGc                                                                        <  3:136/108‑1 (MQ=255)                                                                                                                         gTTGCCTG‑TACCGGCCCTTTCGCGGGCGCGCCCGCGAATGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCa                          <  2:296670/149‑1 (MQ=255)                                                                                                                                       gCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTgcgccagcgc             >  6:208150/1‑149 (MQ=255)                                                                                                                                         ccTTTCGCGGGCGAGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTgcgccagcgc             >  4:270280/1‑147 (MQ=255)                                                                                                                                          cTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGcgccagcgccag          >  8:93120/1‑149 (MQ=255)                                                                                                                                           tttCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCgccagcgccagc         >  2:11686/1‑149 (MQ=255)                                                                                                                                           tttCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCgccagcgccagc         <  6:13650/149‑1 (MQ=255)                                                                                                                                           tttCGCGGGCGAGCCCGCGAAAGGGCCGGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCg                                                       >  2:105574/1‑103 (MQ=255)                                                                                                                                           tttCGCGGGCGAGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCgccagcgccagc         >  5:285450/1‑149 (MQ=255)                                                                                                                                           tttCGCGGGCGAGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCgccagcgccagc         >  3:77335/1‑149 (MQ=255)                                                                                                                                                  ggCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCGCCAGCGCCAGCGGCGTCa  <  6:186361/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  CCAAGCAGGCCTGGGAGTACACCGCCTCGACCCATCTGGTGCGCGACCTGATCGAGGACAAGATGTATCCGCTGACCCTGCAAGGGTTGGAGGACGGGATGCGCTACCTGTCCAAGTGAGGTTGCCTG‑TACCGGCCCTTTCGCGGGCTCGCCCGCGAAAGGGCCAGCAGTCTTACGACCAACGCCCCAGGTAAAGCTGCTGCTCTCGCCAGGCCTCGTCCTGCACCACCGCCCTGAACCGCACCACCGTCCCCGGCATGCACTGCGCCAGTTGCGCCAGCGCCAGCGGCGTCA  >  NC_002947/5194740‑5195032

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