Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,087,770	G→C	12.0%	R259G (CGC→GGC) ‡	metR‑II ←	homocysteine‑binding DNA‑binding transcriptional regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,087,770	0	G	C	12.0%	138.8 / 10.1	50	R259G (CGC→GGC) ‡	metR‑II	homocysteine‑binding DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (22/22);  new base C (3/3);  total (25/25)
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.52e-01

CCGGGCCTGACGCCACGCTAACACCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGC‑TTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGCTCGACGTCGGCGGGGTC  >  NC_002947/3087623‑3087917                                                                                                                                                    |                                                                                                                                                     ccggtacTGAAGCCACGCTAACACATGCCAGAGTGGATAACGATATGTCCTTGGCAGTGAGAAGAAAATAACGCATGTAAGTCACATCAAGAATGTCACTGCGAAGTGCAGCATAGAGGTTTGCCTGCAGGCAGTGATCGCCAAGgagc                                                                                                                                                     >  6:154611/7‑149 (MQ=255)   gggCCTGACGCCACGCTAACACCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGC‑tt                                                                                                                                                  <  3:54017/149‑1 (MQ=255)       cTGACGCCACGCTAACACCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGC‑TTGGCg                                                                                                                                              <  4:36073/149‑1 (MQ=255)       cTGACGCCACGCTAACACCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGC‑TTGGCg                                                                                                                                              <  1:272853/149‑1 (MQ=255)         gACGCCACGCTAACACCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGC‑TTGGCGGt                                                                                                                                            >  3:275547/1‑149 (MQ=255)           cGCCACGCTAACACCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGC‑TTGGCGGTTa                                                                                                                                          <  2:19177/149‑1 (MQ=255)                                     tgAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGC‑TTGGCGGTTACATAGCCGCGCGAACTGTATTCGTg                                                                                                                 >  3:107765/1‑148 (MQ=255)                                      gAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGC‑TTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCa                                                                                                               >  3:244498/1‑149 (MQ=255)                                         cGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGC‑TTGGCGGTTACATAGCCGCGCGAAc                                                                                                                           <  2:158036/134‑1 (MQ=255)                                         cGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGC‑TTGGCGGTTACATAGCCGCGCGAAc                                                                                                                           >  1:158036/1‑134 (MQ=255)                                            gaTGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGC‑TTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGcc                                                                                                          >  7:96645/1‑148 (MQ=255)                                                tCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGCGCC‑TTGGCga                                                                                                                                             <  6:5829/109‑2 (MQ=255)                                                tCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGCGCC‑TTGGCga                                                                                                                                             <  6:1175/109‑2 (MQ=255)                                                tCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGCGCC‑TTGGCga                                                                                                                                             >  5:5829/1‑108 (MQ=255)                                                tCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGCGCC‑TTGGCga                                                                                                                                             >  3:177359/1‑108 (MQ=255)                                                tCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGCGCC‑TTGGCga                                                                                                                                             >  5:1175/1‑108 (MQ=255)                                                tCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGCGCC‑TTGGCga                                                                                                                                             <  4:177359/109‑2 (MQ=255)                                                  cTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGC‑TTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTgg                                                                                                   >  7:184889/1‑149 (MQ=255)                                                          tGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGC‑TTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCa                                                                                                        <  4:25279/136‑1 (MQ=255)                                                          tGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGC‑TTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCa                                                                                                        <  4:20979/136‑1 (MQ=255)                                                          tGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGC‑TTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCa                                                                                                        >  3:25279/1‑136 (MQ=255)                                                          tGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGC‑TTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCa                                                                                                        >  3:20979/1‑136 (MQ=255)                                                                    aTCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGC‑TTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGAcgccgc                                                                                 >  7:124254/1‑149 (MQ=255)                                                                      cACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGC‑TTGGCGGTTACATAGCCGCGCGaa                                                                                                                            >  3:274248/1‑104 (MQ=255)                                                                      cACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGC‑TTGGCGGTTACATAGCCGCGCGaa                                                                                                                            <  4:274248/104‑1 (MQ=255)                                                                       aCGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGC‑TTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGcc                                                                              >  2:181819/1‑149 (MQ=255)                                                                            tGTACGGCACCTCCATCATTTGACTGCGCAGTTCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGC‑TTGGCGGTTACATATCCGCGCTAACTTTATTCGTGCACGGCCCAGTTGTGCAGGCAGCAGACGCCGCGCCCGGaa                                                                         <  4:201914/149‑1 (MQ=255)                                                                             gTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGC‑TTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAg                                                                        >  7:294375/1‑149 (MQ=255)                                                                                      cTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGC‑TTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGtt                                                               <  8:166137/149‑1 (MQ=255)                                                                                      cTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGC‑TTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGtt                                                               >  1:139181/1‑149 (MQ=255)                                                                                      cTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGC‑TTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGtt                                                               <  4:275547/149‑1 (MQ=255)                                                                                      cTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGC‑TTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGtt                                                               <  1:181819/149‑1 (MQ=255)                                                                                         cagcagGTCACTGCGCAGGGCAGCTTAGCGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGC‑TTGGGGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGca                                                            <  7:97948/149‑1 (MQ=255)                                                                                               gTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGC‑TTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGcatcatca                                                      >  8:233116/1‑149 (MQ=255)                                                                                                tgaCTGCGAAGTTCAGCATAGTGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGC‑TTGGCGGTTACATTGCAGCGCGAACTGTATTCGTGGACGGCCCAGTTGGTCAGGCAGCAGACGGCGCGCCCGGAAGCAACCTGtt                                                               <  6:98469/137‑1 (MQ=255)                                                                                                tCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGC‑TTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAAACAGtt                                                               >  5:98469/1‑139 (MQ=255)                                                                                                  aCTGCGCAGTGCGGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGCTTTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGcatcatcat                                                     >  1:86649/1‑148 (MQ=255)                                                                                                          tgcgagcaTAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGC‑TTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAg                                            <  6:174565/144‑1 (MQ=255)                                                                                                             cagcaTAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGC‑TTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTc                                         <  8:96645/148‑1 (MQ=255)                                                                                                               gcaTAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGC‑TTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGggcagg                                                                                              <  6:187127/93‑1 (MQ=255)                                                                                                               gcaTAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGC‑TTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGggcagg                                                                                              >  5:187127/1‑93 (MQ=255)                                                                                                               gcaTAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGC‑TTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGa                                      <  8:124254/149‑1 (MQ=255)                                                                                                                 aTAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGC‑TTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAg                                    <  4:107765/149‑1 (MQ=255)                                                                                                                       gTTTGCCTGCAGGCCGTGCTCGCCAAGGCGC‑TTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCAGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAgtacgta                              >  2:315825/1‑149 (MQ=255)                                                                                                                         ttGCCTGCAGGCCGTGCTCGCCAAGGCGC‑TTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTAcc                            >  8:209993/1‑149 (MQ=255)                                                                                                                               gCAGGCCGTGCTCGCCAAGGCGC‑TTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAAcc                                                                   >  3:182323/1‑104 (MQ=255)                                                                                                                               gCAGGCCGTGCTCGCCAAGGCGC‑TTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAAcc                                                                   <  4:182323/104‑1 (MQ=255)                                                                                                                                          tCGCCAAGGCGC‑TTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGCTCGACGTc           >  3:81983/1‑149 (MQ=255)                                                                                                                                               aaGGCGC‑TTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGCTCGACGTcggcgg      <  7:68838/149‑1 (MQ=255)                                                                                                                                               aaGGCGC‑TTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGCTCGACGTcggcgg      >  8:68838/1‑149 (MQ=255)                                                                                                                                                aGGCGC‑TTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGCTCGACGTCGGCggg     <  8:294375/149‑1 (MQ=255)                                                                                                                                                    gc‑TTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGCTCGACGTcggcgg      <  3:11882/144‑1 (MQ=255)                                                                                                                                                    gc‑TTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGCTCGACGTcggcgg      >  4:11882/1‑144 (MQ=255)                                                                                                                                                    gc‑TTGGCCGTTACATAGCCGCGCGAACTGTATTCGTGCACGCCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCACTAGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGCTCGACGTCGGCcggggc  >  1:224192/1‑146 (MQ=255)                                                                                                                                                    |                                                                                                                                                     CCGGGCCTGACGCCACGCTAACACCTGCCAGCGTGGTGAACGAGATGTCCTTGGCCGTGAGCAGAAAATCACGCATGTACGGCACCTCCAGCATGTCACTGCGCAGTGCAGCATAGAGGTTTGCCTGCAGGCCGTGCTCGCCAAGGCGC‑TTGGCGGTTACATAGCCGCGCGAACTGTATTCGTGCACGGCCCAGTTGGGCAGGCAGCAGACGCCGCGCCCGGAAGCAACCAGTTGCATCATCATCACGGTCAGTTCTGAAGTACGTACCTGTGCGGGCTCGACGTCGGCGGGGTC  >  NC_002947/3087623‑3087917

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