Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	964,019	G→A	7.1%	A210T (GCC→ACC)	phnC →	phosphonate ABC transporter ATP‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	964,019	0	G	A	7.1%	164.5 / 2.7	56	A210T (GCC→ACC)	phnC	phosphonate ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (27/25);  new base A (2/2);  total (29/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 = 9.34e-01

CTCGGGCGTTGTACCAGCGTCCGCAGGTGCTGCTGACCGACGAACCGGTGTCGGCCATGGACCC‑GGTGTTGGCCGACCATAGCCTGGCCTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTTGGCGCTGGCGCATTTTCCGCGGGTGATTGGCATTCGGGAAGGGCAGGTTGTGTTTGACTGCCCGGCCGAAGCGGTCACCGAGCAACTGCTGGACGCGCTATACGCCAACGAACAACTGGCTCCGCCGCCAACCCA  >  NC_002947/963876‑964164                                                                                                                                                 |                                                                                                                                                  cTCGGGCGAGGTACCGCCGGCCGCAGCACCTGCTGACCGCCCCACCGGTGCCGGCCCTGCTCCA‑GGTGTTGGCCGACCATAGCCTGGCCTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTa                                                                                                                                              <  5:245635/149‑1 (MQ=255)  tCGGGCGTTGTACCAGCGTCCGCAGGTGCTGCTGACCGACGAACCGGTGTCGGCCATGGACCC‑GGTGTTGGCCGACCATAGCCTGGCCTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACAcca                                                                                                                                                >  5:34706/1‑145 (MQ=255)  tCGGGCGTTGTACCAGCGTCCGCAGGTGCTGCTGACCGACGAACCGGTGTCGGCCATGGACCC‑GGTGTTGGCCGACCATAGCCTGGCCTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACAcca                                                                                                                                                >  5:211112/1‑145 (MQ=255)        gTTGTACCAGCGTCCGCAGGTGCTGCTGACCGACGAACCGGTGTCGGCCATGGACCC‑GGTGTTGGCCGACCATAGCCTGGCCTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTtgg                                                                                                                                       <  8:305242/149‑1 (MQ=255)          tGTACCAGCGTCCGCAGGTGCTGCTGACCGACGAACCGGTGTCGGCCATGGACCC‑GGTGTTGGCCGACCATAGCCTGGCCTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTtggcg                                                                                                                                     >  2:174198/1‑149 (MQ=255)             aCCAGCGTCCGCAGGTGCTGCTGACCGACGAACCGGTGTCGGCCATGGACCC‑GGTGTTGGCCGACCATAGCCTGGCCTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTtggcgctg                                                                                                                                  >  6:352969/1‑149 (MQ=255)                 gCGTCCGCAGGTGCTGCTGACCGACGAACCGGTGTCGGCCATGGACCC‑GGTGTTGGCCGACCATAGCCTGGCCTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGAGGTTGGTGGCGAGCCTGCACGCCGTAGAGTtggcgctggcg                                                                                                                               >  5:299723/1‑148 (MQ=255)                 gCGTCCGCAGGTGCTGCTGACCGACGAACCGGTGTCGGCCATGGACCC‑GGTGTTGGCCGACCATAGCCTGGCCTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGt                                                                                                                                          >  1:129176/1‑137 (MQ=255)                 gCGTCCGCAGGTGCTGCTGACCGACGAACCGGTGTCGGCCATGGACCC‑GGTGTTGGCCGACCATAGCCTGGCCTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGt                                                                                                                                          <  2:129176/137‑1 (MQ=255)                 gCGTCCGCAGGTGCTGCTGACCGACGAACCGGTGTCGGCCATGGACCC‑GGTGTTGGCCGACCATAGCCTGGCCTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTtggcgctggc                                                                                                                                <  1:290456/147‑1 (MQ=255)                 gCGTCCGCAGGTGCTGCTGACCGACGAACCGGTGTCGGCCATGGACCC‑GGTGTTGGCCGACCATAGCCTGGCCTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTtggcgctggc                                                                                                                                >  2:290456/1‑147 (MQ=255)                 gCGTCCGCAGGTGCTGCTGACCGACGAACCGGTGTCGGCCATGGACCC‑GGTGTTGGCCGACCATAGACTGGACTTGCTCAATCGCCATGCGCAGGCCACTGGCCTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTtggcgctggcg                                                                                                                               >  5:341157/1‑148 (MQ=255)                       gCAGGTGCTGCTGACCGACGAACCGGTGTCGGCCATGGACCCGGGTGTTGGCCGACCATAGCCTGGCCTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTTGGCGCTGGCGCAtttt                                                                                                                         >  1:71284/1‑149 (MQ=255)                          ggTGCTGCTGACCGACGAACCGGTGTCGGCCATGGACCC‑GGTGTTGGCCGACCATAGCCTGGCCTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTTGGCGCTGGCGCATTTTCcgc                                                                                                                     <  5:352969/149‑1 (MQ=255)                            tgctgcTGACCGACGAACCGGTGTCGGCCATGGACCC‑GGTGTTGGCCGACCATAGCCTGGCCTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTTGGCGCTGGCGCATTTTCCGCgg                                                                                                                   <  8:278368/149‑1 (MQ=255)                                    aCCGACGAACCGGTGTCGGCCATGGACCC‑GGTGTTGGCCGACCATAGCCTGGCCTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTTGGCGCTGGCGCATTTTCCGCGGGTGATTgg                                                                                                           >  7:205942/1‑149 (MQ=255)                                    aCCGACGAACCGGTGTCGGCCATGGACCC‑GGTGTTGGCCGACCATAGCCTGGCCTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACACCATAGAGTTGGCGCTGGCGCATTTTCCGCGGGTGATTg                                                                                                            <  6:34706/148‑1 (MQ=255)                                    aCCGACGAACCGGTGTCGGCCATGGACCC‑GGTGTTGGCCGACCATAGCCTGGCCTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACACCATAGAGTTGGCGCTGGCGCATTTTCCGCGGGTGATTg                                                                                                            <  6:211112/148‑1 (MQ=255)                                      cgacgaACCGGTGTCGGCCATGGACCC‑GGTGTTGGCCGACCATAGCCTGGCCTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGAAGAGTTGGCGCTGGCGCATTTTCCGCGGGTGATTGGCa                                                                                                         >  3:302683/1‑149 (MQ=255)                                        acgaACCGGTGTCGGCCATGGACCC‑GGTGTTGGCCGACCATAGCCTGGCCTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTTGGCGCTGGCGCATTTTCCGCGGGTGATTGGCAtt                                                                                                       >  8:113776/1‑149 (MQ=255)                                          gcaCCGGTGTCGTCCATGGCCCC‑GGTGTGGTCCGACCATAGCCTGTCCTGGCTCTATCGCCATGCGCAGGCCACTGGCGTGACGTTGGGGGCGAGCCGGCACGCCGTAGAGTTGGCGCGGGCGCTTTTTCCGCGGGGGAGTGGGATTCg                                                                                                     <  6:341157/147‑1 (MQ=255)                                              cGGTGTCGGCCATGGACCC‑GGTGTTGGCCGACCATAGCCTGGCCTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTa                                                                                                                                              <  5:70898/104‑1 (MQ=255)                                              cGGTGTCGGCCATGGACCC‑GGTGTTGGCCGACCATAGCCTGGCCTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTa                                                                                                                                              >  6:70898/1‑104 (MQ=255)                                              cGGTGTCGGCCATGGACCC‑GGTGTTGGCCGACCATAGCCTGGCCTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTTGGCGCTGGCGCATTTTCCGCGGGTGATTGGCATTCggg                                                                                                   >  2:48374/1‑147 (MQ=255)                                              cGGTGTCGGCCATGGACCC‑GGTGTTGGCCGACCATAGCCTGGCCTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTTGGCGCTGGCGCATTTTCCGCGGGTGATTGGCATTCGGGaa                                                                                                 <  6:151730/149‑1 (MQ=255)                                               ggTGTCGGCCATGGACCC‑GGTGTTGGCCGACCATAGCCTGGCCTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGTCGTTGGTGGCGAGCCTGCACGCCGTAGAGTTGGCGCTGGCGCATTTTCCGCGGGTGATTGGCATTCGGGAAg                                                                                                <  2:71284/149‑1 (MQ=255)                                                gtgtCGGCCATGGACCC‑GGTGTTGGCCGACCATAGCCTGGCCTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTTGGCGCTGGCGCATTTTCCGCGGGTGATTGGCATTCGGGAAgg                                                                                               >  2:292297/1‑149 (MQ=255)                                                         aTGGACCC‑GGTGTTGGCCGACCATAGCCTGGCCTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTTGGCGCTGGCGCATTTTCCGCGGGTGATTGGCATTCGGGAAGGGCAGGTtgt                                                                                      <  6:299723/149‑1 (MQ=255)                                                         aTGGACCC‑GGTGTTGGCCGACCATAGCCTGGCCTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTTGGCGCTGGCGCATTTTCCGCGGGTGATTGGCATTCGGGAAGGGCAGGTtgt                                                                                      <  1:261673/149‑1 (MQ=255)                                                             aCCC‑GGTGTTGGCCGACCATAGCCTGGCCTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTTGGCGCTGGCGCATTTTCCGCGGGTGATTGGCATTCGGGAAGGGCAGGTTGTGttt                                                                                  <  8:205942/149‑1 (MQ=255)                                                              ccc‑GGTGTTGGCCGACCATAGCCTGGCCTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTTGGCGCTGGCGCATTTTCCGCGGGTGATTGGCATTCGGGAAGGGCAGGtt                                                                                        <  6:35058/142‑1 (MQ=255)                                                              ccc‑GGTGTTGGCCGACCATAGCCTGGCCTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTTGGCGCTGGCGCATTTTCCGCGGGTGATTGGCATTCGGGAAGGGCAGGtt                                                                                        >  5:35058/1‑142 (MQ=255)                                                                         gCCGACCATAGCCTGGCCTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTTGGCGCTGGCGCATTTTCCGCGGGTGATTGGCATTCGGGAAGGGCAGGTTGTGTTTGACTGCCCGGc                                                                       >  1:87116/1‑149 (MQ=255)                                                                         gCCGACCATAGCCTGGCCTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTTGGCGCTGGCGCATTTTCCGCGGGTGATTGGCATTCGGGAAGGGCAGGTTGTGTTTGACTGCCCGGc                                                                       <  2:305672/149‑1 (MQ=255)                                                                                  aGCCTGGCCTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTTGGCGCTGGCGCATTTTCCGCGGGTGATTGGCATTCGGGAAGGGCAGGTTGTGTTTGACTGCCCGGCCGAAGCGGt                                                              >  5:353784/1‑149 (MQ=255)                                                                                   gCCTGGCCTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTAGGCGCTGGCGCATTTTCCGCGGGTGATTgg                                                                                                           <  1:28318/103‑1 (MQ=255)                                                                                   gCCTGGCCTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTAGGCGCTGGCGCATTTTCCGCGGGTGATTgg                                                                                                           >  2:28318/1‑103 (MQ=255)                                                                                         ccTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTTGGCGCTGGCGCATTTTCCGCGGGTGATTGGCATTCGGGAAGGGCAGGTTGTGTTTGACTGCCCGGCCGAAGCGGTCACCGAg                                                       >  1:325356/1‑149 (MQ=255)                                                                                          cTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTTGGCGCTGGCGCATTTTCCGCGGGTGATTGGCATTCGGGAAGGGCAGGTTGTGTTTGACTGCCCGGCCGAAGCGGTCACCGAg                                                       >  2:157086/1‑148 (MQ=255)                                                                                          cTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTTGGCGCTGGCGCATTTTCCGCGGGTGATTGGCATTCGGGAAGGGCAGGTTGTGTTTGACTGCCCGGCCGAAGCGGTCACCGAGc                                                      <  7:113776/149‑1 (MQ=255)                                                                                               tCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTTGGCGCTGGCGCATTTTCCGCGGGTGATTGGCATTCGGGAAGGGCAGGTTGTGTTTGACTGCCCGGCCGAAGCGGTCACCGAGCAActg                                                 >  5:340690/1‑149 (MQ=255)                                                                                                 aaTCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTTGGCGCTGGCGCATTTTCCGCGGGTGATTGGCATTCGGGAAGGGCAGGTTGTGTTTGACTGCCCGGCCGAAGCGGTCACCGAGCAActgct                                               <  1:157086/149‑1 (MQ=255)                                                                                                 aaTCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTTGGCGCTGGCGCATTTTCCGCGGGTGATTGGCATTCGGGAAGGGCAGGTTGTGTTTGACTGCCCGGCCGAAGCGGTCACCGAGCAActgc                                                <  6:353784/148‑1 (MQ=255)                                                                                                   tCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTTGGCGCTGGCGCATTTTCCGCGGGTGATTGGCATTCGGGAAGGGCAGGTTGTGTTTGACTGCCCGGCCGAAGCGGTCACCGAGCAACTGCTgg                                             <  1:292297/149‑1 (MQ=255)                                                                                                       cATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTTGGCGCTGGCGCATTTTCCGCGGGTGATTGGCATTCGGGAAGGGCAGGTTGTGTTTGACTGCCCGGCCGAAGCGGTCACCGAGCAACTGCTGGAcgc                                         >  3:52364/1‑149 (MQ=255)                                                                                                         tGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTTGGCGCTGGCGCATTTTCCGCGGGTGATTGGCATTCGGGAAGGGCAGGTTGTGTTTGACTGCCCGGCCGAAGCGGTCACCGAGCAACTGCTGGACgcgc                                       <  1:174198/149‑1 (MQ=255)                                                                                                          gcgcAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTTGGCGCTGGCGCATTTTCCGCGGGTGATTGGCATTCGGGAAGGGCAGGTTGTGTTTGACTGCCCGGCCGAAGCGGTCACCGAGCAACTGCTGGACgcgc                                       >  4:49876/1‑148 (MQ=255)                                                                                                              aGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTTGGCGCTGGCGCATTTTCCGCGGGTGATTGGCATTCGGGAAGGGCAGGTTGTGTTTGACTGCCCGGCCGAAGCGGTCACCGa                                                        <  7:237563/127‑1 (MQ=255)                                                                                                              aGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTTGGCGCTGGCGCATTTTCCGCGGGTGATTGGCATTCGGGAAGGGCAGGTTGTGTTTGACTGCCCGGCCGAAGCGGTCACCGa                                                        >  8:237563/1‑127 (MQ=255)                                                                                                                         gTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTTGGCGCTGGCGCATTTTCCGCGGGTGATTGGCATTCGGGAAGGGCAGGTTGTGTTTGACTGCCCGGCCGAAGCGGTCACCGAGCAACTGCTGGACGCGCTATACGCCAACGaac                        >  4:294616/1‑148 (MQ=255)                                                                                                                               ttGGTGGCGAGCCTGCACGCCGTAGAGTTGGCGCTGGCGCATTTTCCGCGGGTGATTGGCATTCGGGAAGGGCAGGTTGTGTTTGACTGCCCGGCCGAAGCGGTCACCGAGCAACTGCTGGACGCGCTATAc                                  >  7:43095/1‑132 (MQ=255)                                                                                                                                     gCGAGCCTGCACGCCGTAGAGTTGGCGCTGGCGCATTTTCCGCGGGTGATTGGCATTCGGGAAGGGCAGGTTGTGTTTGACTGCCCGGCCGAAGCGGTCAcc                                                          <  2:266264/102‑1 (MQ=255)                                                                                                                                     gCGAGCCTGCACGCCGTAGAGTTGGCGCTGGCGCATTTTCCGCGGGTGATTGGCATTCGGGAAGGGCAGGTTGTGTTTGACTGCCCGGCCGAAGCGGTCAcc                                                          >  1:266264/1‑102 (MQ=255)                                                                                                                                         gCCTGCACGCCGTAGAGTTGGCGCTGGCGCATTTTCCGCGGGTGATTGGCATTCGGGAAGGGCAGGtt                                                                                        >  4:22508/1‑68 (MQ=255)                                                                                                                                         gCCTGCACGCCGTAGAGTTGGCGCTGGCGCATTTTCCGCGGGTGATTGGCATTCGGGAAGGGCAGGtt                                                                                        <  3:22508/68‑1 (MQ=255)                                                                                                                                              cACGCCGTAGAGTTGGCGCTGGCGCATTTTCCGCGGGTGATTGGCATTCGGGAAGGGCAGGTTGTGTTTGACTGCCCGGCCGAAGCGGTCACCGAGCAACTGCTGGACGCGCTATACGCCAACGAACAACTGGCTCCGCCGCCAACCCa  <  1:48374/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                  CTCGGGCGTTGTACCAGCGTCCGCAGGTGCTGCTGACCGACGAACCGGTGTCGGCCATGGACCC‑GGTGTTGGCCGACCATAGCCTGGCCTTGCTCAATCGCCATGCGCAGGCCACTGGCGTGACGTTGGTGGCGAGCCTGCACGCCGTAGAGTTGGCGCTGGCGCATTTTCCGCGGGTGATTGGCATTCGGGAAGGGCAGGTTGTGTTTGACTGCCCGGCCGAAGCGGTCACCGAGCAACTGCTGGACGCGCTATACGCCAACGAACAACTGGCTCCGCCGCCAACCCA  >  NC_002947/963876‑964164

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