Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	284,586	T→A	15.0%	intergenic (+33/+30)	betT‑II → / ← atsK	choline/carnitine/betaine transporter/alpha‑ketoglutarate‑dependent sulfonate dioxygenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	284,586	0	T	A	15.0%	114.1 / 8.0	51	intergenic (+33/+30)	betT‑II/atsK	choline/carnitine/betaine transporter/alpha‑ketoglutarate‑dependent sulfonate dioxygenase
Reads supporting (aligned to +/- strand):  ref base T (23/20);  new base A (4/4);  total (27/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGA‑GGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCATTC  >  NC_002947/284446‑284730                                                                                                                                             |                                                                                                                                                  gAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCcgcg                                                                                                                                           <  3:33452/149‑1 (MQ=255) gAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGT‑CTGGCCTCTTCGCGGGCAAGCCCGCGa                                                                                                                                          <  4:8815/149‑1 (MQ=255)     gCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAaga                                                                                                                                       <  2:213110/149‑1 (MQ=255)     gCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCAAGCCCGCGAaga                                                                                                                                       <  6:173129/149‑1 (MQ=255)      cAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAATGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAaga‑g                                                                                                                                     <  5:103419/149‑1 (MQ=255)       aGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGAACAAGCTCTACAGTCCAGGCCACTTAGCGGGGAAGCCCGCGAAGA‑gg                                                                                                                                    >  7:2040/1‑149 (MQ=255)         ccAGAACTATCACCTGGTGGGCTACAGCAAGGAGCATATCAACAACCACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGA‑GGcc                                                                                                                                  <  1:374020/149‑1 (MQ=255)                tATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGaa                                                                                                                                         >  1:127118/1‑136 (MQ=255)                tATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGaa                                                                                                                                         <  2:127118/136‑1 (MQ=255)                   gACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGGCCGCGAAGC‑GGCCAGCCAGGCaa                                                                                                                        >  8:154455/1‑149 (MQ=255)                    aCCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGA‑GGCCAGACAGGCAAc                                                                                                                       >  6:134708/1‑149 (MQ=255)                                 cAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGACGGCCAGACAGGCAACCGCTAGAAAGGt                                                                                                           <  1:16877/149‑1 (MQ=255)                                 cAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGA‑GGCCAGACAGGCAACCGCTAGAAAGGTg                                                                                                          >  1:291891/1‑149 (MQ=255)                                   gCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGA‑GGCCAGACAGGCAACCGCTa                                                                                                                  <  7:49653/139‑1 (MQ=255)                                   gCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGA‑GGCCAGACAGGCAACCGCTa                                                                                                                  >  8:49653/1‑139 (MQ=255)                                    cAAGGAGCAGATCGTCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGA‑GGCCAGACAGGCAACCGCTAGAAAGGTGCAt                                                                                                       >  7:29976/1‑149 (MQ=255)                                        gAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCAGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGA‑GGCCAGACAGGCAACCGCTAGAAAGGTGCATccc                                                                                                    >  1:112204/1‑148 (MQ=255)                                            aGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGA‑GGCCAGACAGGCAACCGCTAGAAAGGt                                                                                                           <  5:141355/137‑1 (MQ=255)                                            aGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGA‑GGCCAGACAGGCAACCGCTAGAAAGGt                                                                                                           >  6:141355/1‑137 (MQ=255)                                             gATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGA‑GGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGa                                                                                              >  3:70394/1‑149 (MQ=255)                                              atcatcAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGA‑GGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGAt                                                                                             >  5:106161/1‑149 (MQ=255)                                                    aaCGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGA‑GGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGc                                                                                       >  1:318677/1‑149 (MQ=255)                                                          aTCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGT‑CTGGCCTCTTCGCGGGCAAGCCCGCGAaga‑g                                                                                                                                     <  3:222468/96‑1 (MQ=37)                                                          aTCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGT‑CTGGCCTCTTCGCGGGCAAGCCCGCGAaga‑g                                                                                                                                     >  4:222468/1‑96 (MQ=37)                                                             cTCGACCAGTACGAGCGGCACATTCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGA‑GGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCa                                                                               <  2:291891/148‑1 (MQ=255)                                                             cTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGA‑GGCCAGACAGGCAACCGCTa                                                                                                                  <  5:73514/113‑1 (MQ=255)                                                             cTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGA‑GGCCAGACAGGCAACCGCTa                                                                                                                  >  6:73514/1‑113 (MQ=255)                                                              tCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGA‑GGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCAt                                                                              >  2:121654/1‑148 (MQ=255)                                                                  ccAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGA‑GGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCAccc                                                                         <  2:256973/149‑1 (MQ=255)                                                                      tACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGA‑GGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACccgccg                                                                     <  8:29976/149‑1 (MQ=255)                                                                      tACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGA‑GGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACccgccg                                                                     <  6:106161/149‑1 (MQ=255)                                                                            cGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTGCTGGCCTCTTCGCGGGCTTGCCCGCGAAGA‑GGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCACCACCCGCCGGTtcgg                                                               <  4:70394/149‑1 (MQ=255)                                                                              ccgttctcAGCACCTGCCCGTGGGCCGCTCCCCATGCTCTCCTGTCCGGGCCTCTGCGGGGGCTTCCCCGCTAAGA‑GGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTtcggt                                                              <  1:373990/141‑1 (MQ=255)                                                                                       gTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGA‑GGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAatcat                                                    >  6:295084/1‑149 (MQ=255)                                                                                          ccTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGA‑GGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGt                                                                                  <  5:323545/116‑1 (MQ=255)                                                                                          ccTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGA‑GGCCAGACAGGCAACCGCTAGAAAGGAGCATCCCCAAGGATGGTCGCCCGGt                                                                                  >  6:323545/1‑116 (MQ=255)                                                                                             gCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGA‑GGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCg                                              >  7:78759/1‑149 (MQ=255)                                                                                             gCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGA‑GGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCg                                              >  2:8634/1‑149 (MQ=255)                                                                                                       cGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGA‑GGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTg                                    <  8:78759/149‑1 (MQ=255)                                                                                                              cAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGA‑GGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCACCCACCGCAAAATGCTGGGTCAcg                             >  4:165123/1‑149 (MQ=255)                                                                                                               aaGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGA‑GGCCAGACAGTCAACCGCTAGAAAGGTGCATCCCCAAGGAAGGTCGCCCGGTGCATCACCCGCCGGATCGGTCGGAAATCATCCACCGCAAAATGCTGGGTCAcgc                            >  8:169890/1‑149 (MQ=255)                                                                                                               aaGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGA‑GGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCAcgc                            >  8:174856/1‑149 (MQ=255)                                                                                                                       ctgtcctgGCCTCTTCGCGGGCTTGCCCGCGAAGA‑GGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTcc                    >  3:33865/1‑149 (MQ=255)                                                                                                                       cggcccgggCCTCTGCCCGGGCTTCCCCGCGAAGA‑GGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTcc                    <  4:33865/142‑1 (MQ=255)                                                                                                                            ctgGCCTCTTCCCGGGCATGCCCGCCAAGA‑GGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCCCCCGCCGCTTCGGTCGGTAATAATCCACCCCAAAATGCTGGGTAACGCGGTTGTCCCAGaa               >  6:242850/1‑149 (MQ=255)                                                                                                                             tgGCCTCTTCGCGGGCTTGCCCGCGAAGA‑GGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTcc                    <  4:26547/143‑1 (MQ=255)                                                                                                                             tgGCCTCTTCGCGGGCTTGCCCGCGAAGA‑GGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTcc                    >  3:26547/1‑143 (MQ=255)                                                                                                                                 ctctTCGCGGGCTTGCCCGCGAAGA‑GGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCa          >  2:31301/1‑149 (MQ=255)                                                                                                                                  tctTCGCGGGCAAGCCCGCGAAGA‑GGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGc                                             <  2:314093/113‑1 (MQ=255)                                                                                                                                  tctTCGCGGGCAAGCCCGCGAAGA‑GGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGc                                             >  1:314093/1‑113 (MQ=255)                                                                                                                                         gggCTTGCCCGCGAAGA‑GGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCATtc  <  1:31301/149‑1 (MQ=255)                                                                                                                                             |                                                                                                                                                  GAAGGCAGCCAGAACTATGACCTGGTGGGCTACAGCAAGGAGCAGATCATCAACGACATCCTCGACCAGTACGAGCGGCACATGCAGTACCTGCACCTGGTCCGCTGACCAAGCTCTCCTGTCCTGGCCTCTTCGCGGGCTTGCCCGCGAAGA‑GGCCAGACAGGCAACCGCTAGAAAGGTGCATCCCCAAGGATGGTCGCCCGGTGCATCACCCGCCGGTTCGGTCGGTAATCATCCACCGCAAAATGCTGGGTCACGCGGTTGTCCCAGAACGCCACGTCATTC  >  NC_002947/284446‑284730

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