Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,212,507	A→G	7.3%	intergenic (+71/+23)	PP_4593 → / ← PP_4594	hypothetical protein/cystathionine gamma‑synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,212,507	0	A	G	7.3%	156.3 / 2.8	55	intergenic (+71/+23)	PP_4593/PP_4594	hypothetical protein/cystathionine gamma‑synthase
Reads supporting (aligned to +/- strand):  ref base A (23/28);  new base G (2/2);  total (25/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

CCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCG‑ACAACCGCAGCAGGTTCTCGCTGATCCCCATCACCCCTTTCTGCTCTGGG  >  NC_002947/5212360‑5212653                                                                                                                                                    |                                                                                                                                                    ccGTCAACCGACTGCGCAAGGAAAGGTTTCTGATCAAGGACAAGATCGAACGGATGATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTAc                                                                                                                                                    <  7:322120/149‑1 (MQ=255)  cGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCAAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCCTTTCGCGGGTAAACCCGCTCCTACGACGGTGTAGGAGCGGGCTTGcc                                                                                                                                                   >  2:229355/1‑149 (MQ=25)       aCCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGa                                                                                                                                              <  6:257245/149‑1 (MQ=255)        ccGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGa                                                                                                                                              >  8:14799/1‑148 (MQ=255)          gCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAg                                                                                                                                           >  7:277232/1‑149 (MQ=255)           ccTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAgg                                                                                                                                          <  2:111031/149‑1 (MQ=255)               cgcAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCTCCCAAATCCTGACAGGCATATATTGAATCCACCCGACTCTTCGCGGGCAAGCCCGCTCATACACCGTCGTAGGAGCCTGATTACCCGCGAAACGGccg                                                                                                                                      >  6:330861/1‑148 (MQ=255)                gcAAGGCAACGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCCTCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTg                                                                                                                                     <  6:4371/149‑1 (MQ=255)                gcAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTg                                                                                                                                     <  4:173690/149‑1 (MQ=255)                 cAAGGAAAGGTTGCTGATCAAGGACAAGATCAAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCCTTTCGCGGGTAAACCCGCTCCTACGACGGTGTAGGAGCGGGCTTGCCCGCGAAGAGGccggt                                                                                                                                    <  1:229355/149‑4 (MQ=18)                       aaGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCt                                                                                                                              <  8:159018/149‑1 (MQ=255)                              cTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGGCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGCAGGCCCtttcgcg                                                                                                                        <  5:330861/148‑7 (MQ=38)                              cTGATCAACGCCCAGCTCGCACGGATCATTCCTCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCGCAGGCCCtttcgcg                                                                                                                        <  6:247633/148‑7 (MQ=37)                                                 aaCGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGc                                                                                                    <  8:277232/149‑1 (MQ=255)                                                     gATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAg                                                                                                                       <  6:75203/126‑1 (MQ=255)                                                     gATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAg                                                                                                                       >  5:75203/1‑126 (MQ=255)                                                      aTCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTcacac                                                                                                                         <  6:42222/123‑1 (MQ=255)                                                      aTCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTcacac                                                                                                                         >  5:42222/1‑123 (MQ=255)                                                      aTCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAg                                                                                                                <  7:94800/132‑1 (MQ=255)                                                      aTCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAg                                                                                                                >  8:94800/1‑132 (MQ=255)                                                                 gCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTc                                                                                     >  7:122490/1‑148 (MQ=255)                                                                   ccAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGAc                                                                                  >  6:267823/1‑149 (MQ=255)                                                                     aaaTCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGAc                                                                                  >  1:300114/1‑147 (MQ=255)                                                                               aGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACg                                                                         <  4:184720/146‑1 (MQ=255)                                                                               aGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACg                                                                         >  3:184720/1‑146 (MQ=255)                                                                                 gCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAAt                                                                    >  1:208421/1‑149 (MQ=255)                                                                                   aTCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGTGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACCGTAATc                                                                   >  6:91403/1‑148 (MQ=255)                                                                                   aTCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGc                                                                                                    <  4:230656/115‑1 (MQ=255)                                                                                   aTCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGc                                                                                                    >  3:230656/1‑115 (MQ=255)                                                                                   aTCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACgg                                                                        <  3:310983/143‑1 (MQ=255)                                                                                   aTCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACgg                                                                        >  4:310983/1‑143 (MQ=255)                                                                                        ttGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATc                                                                   >  1:99308/1‑143 (MQ=255)                                                                                        ttGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATc                                                                   <  2:99308/143‑1 (MQ=255)                                                                                        ttGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCt                                                                 <  7:14799/145‑1 (MQ=255)                                                                                        ttGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGa                                                              <  2:208421/148‑1 (MQ=255)                                                                                        ttGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGAACAGCGCATTCAACGCCCCATCCAGGTACTCGACCAGGGCACGGTAAtagatag                                                              >  8:155265/1‑142 (MQ=38)                                                                                        ttGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTAGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTaa                                                                     <  7:155265/141‑1 (MQ=255)                                                                                            aTCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCt                                                         <  8:122490/149‑1 (MQ=255)                                                                                                     ccTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCG‑ACAAc                                               <  2:300114/149‑1 (MQ=255)                                                                                                               ggCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCG‑ACAACCGCAGCAGGt                                     <  5:267823/149‑1 (MQ=255)                                                                                                                     cccGATCCTACACCGTCGTAGGAGCGGTTTTACCCGCGAAAGGGTCTGCGGCCCTAACACAGCGACCAGCGCTTTCAAAGCCCCATCCAGGTCCTCTACCAGGTCACGTTAATCCTCGATCCCTACCG‑ACAACCGCAGCAGGTTCTCGc                               <  5:91403/149‑1 (MQ=255)                                                                                                                       cGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCG‑ACAACCGCAGCAGGTTCTCGCt                              >  2:26772/1‑148 (MQ=255)                                                                                                                        gCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCG‑ACAACCGCAGCAGGTTCTCGCTGa                            >  7:242095/1‑149 (MQ=255)                                                                                                                               cacCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCG‑ACAACCGCAGCAGGTTCTCGCTGATCCCCAt                     <  6:179891/149‑1 (MQ=255)                                                                                                                                acCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCG‑ACAACCGCAGCAGGTTCTCGCTGATCCCCAt                     >  2:39539/1‑148 (MQ=255)                                                                                                                                acCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCG‑ACAACCGCAGCAGGTTCTCGCTGATCCCCAt                     >  1:107996/1‑148 (MQ=255)                                                                                                                                  cgtcgtAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCG‑ACAACCGCAGCAGGTTCTCGCTGATCCCCATCAc                  >  6:98694/1‑149 (MQ=255)                                                                                                                                   ggtgtAGGAGCGGGCTTGCCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCG‑ACAACCGCAGCAGGTTCTCGCTGATCCCCATCAcc                 >  4:121791/4‑149 (MQ=255)                                                                                                                                   ggtgtAGGAGCGGGCTTGCCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCG‑ACAACCGCAGCAGGTTCTCGCTGATCCCCATCAcc                 <  3:121791/146‑1 (MQ=255)                                                                                                                                        tttcgCGGGTAAACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGCCCAGGTCACGGTAATCCTCGATCCCTACCG‑ACAACCGCAGCAGGTTCTCGCTGATCCCCATCACCCCtt             >  2:82783/5‑148 (MQ=255)                                                                                                                                               ggTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCG‑ACAACCGCAGCAGGTTCTCGc                               >  6:142061/1‑123 (MQ=255)                                                                                                                                               ggTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCG‑ACAACCGCAGCAGGTTCTCGc                               <  5:142061/123‑1 (MQ=255)                                                                                                                                                gTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCG‑ACAACCGCAGCAGGTTCTCGCTGATCCCCATCACCCCTTTCTGCTCTg    <  1:39539/149‑1 (MQ=255)                                                                                                                                                 tttACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCG‑ACAACCGCAGCAGGTTCTCGCTGATCCCCATCACCCCTTTCTGCTCTgg   <  5:98694/149‑1 (MQ=255)                                                                                                                                                   tACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTCCCG‑ACAACCGCAGCAGGTTCTCGCTGATCCCCATCACCCCTTTCTGCTCTggg  <  1:26772/148‑1 (MQ=255)                                                                                                                                                   aaCCCGCGAAAGGGCCTGCGGCCCTCACACATCGACCATCGCTTTCAACGCCCCATGCAGGTCCTCGACCAGGTCACGGGAATCCTCGATCCCTACCGCCCAACCGCAGCAGGTTCTCGCTGATCCCCATCACCCCTTTCTGCTCTggg  <  1:82783/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    CCGTCAACCGCCTGCGCAAGGAAAGGTTGCTGATCAAGGACAAGATCGAACGGATCATTCATCCGCCCAAATCCTGACAGGCATCTATTGAATCCACCGGCCTCTTCGCGGGCAAGCCCGCTCCTACACCGTCGTAGGAGCGGGTTTACCCGCGAAAGGGCCTGCGGCCCTCACACAGCGACCAGCGCTTTCAACGCCCCATCCAGGTCCTCGACCAGGTCACGGTAATCCTCGATCCCTACCG‑ACAACCGCAGCAGGTTCTCGCTGATCCCCATCACCCCTTTCTGCTCTGGG  >  NC_002947/5212360‑5212653

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