Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	156,383	C→T	10.2%	intergenic (‑8/+35)	ygdQ ← / ← citN	transport protein/citrate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	156,383	0	C	T	10.2%	100.4 / 2.9	39	intergenic (‑8/+35)	ygdQ/citN	transport protein/citrate transporter
Reads supporting (aligned to +/- strand):  ref base C (19/16);  new base T (2/2);  total (21/18)
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.21e-01

CGGCTGACCAGAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGCCGAACTC  >  NC_002947/156239‑156528                                                                                                                                                 |                                                                                                                                                  aggCTGACCAGATTCGAGATGATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCTGTAAGGGTGTAAAAGGCAACCCAGATTTCCGGGGTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGt                                                                                                                                                <  8:248664/147‑1 (MQ=255)         cAGAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGcgcg                                                                                                                                       >  2:42144/1‑149 (MQ=255)          aGAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACACAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCTGACCGGGGcgcgt                                                                                                                                      <  4:172228/149‑1 (MQ=255)                        gatgatgTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTACCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATAc                                                                                                                        >  2:13802/1‑149 (MQ=255)                             tgttgtCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATg                                                                                                                   >  4:196375/1‑149 (MQ=255)                                  tCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGGCAGGGCGCGTCGCGTACGTGATACGGATGTATTa                                                                                                              <  3:96621/149‑1 (MQ=255)                                  tCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTCCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTa                                                                                                              <  1:13802/149‑1 (MQ=255)                                    gATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTaaa                                                                                                            <  7:249390/149‑1 (MQ=255)                                         cGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGacta                                                                                                       <  2:220500/149‑1 (MQ=255)                                         cGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGacta                                                                                                       >  1:20868/1‑149 (MQ=255)                                         cGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGa                                                                                                          >  8:192980/1‑146 (MQ=255)                                         cGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGa                                                                                                          <  7:192980/146‑1 (MQ=255)                                               cGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAAc                                                                                                 <  3:187620/149‑1 (MQ=255)                                                gATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGaa                                                                                                  >  7:221058/1‑147 (MQ=255)                                                     cAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCTGACCGGGGcgcgttg                                                                                                                                    <  8:45702/108‑3 (MQ=255)                                                     cAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCTGACCGGGGcgcgttg                                                                                                                                    >  7:45702/1‑106 (MQ=255)                                                     cAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCgg                                                                                            >  7:146920/1‑148 (MQ=255)                                                      aaGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGAt                                                                                                                          >  8:91729/1‑117 (MQ=255)                                                      aaGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGAt                                                                                                                          <  7:91729/117‑1 (MQ=255)                                                                  ggTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCAt                                                                                >  5:21295/1‑147 (MQ=255)                                                                   gTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATcag                                                                             >  3:11628/1‑149 (MQ=255)                                                                       aaaaGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGt                                                                         <  2:20868/149‑1 (MQ=255)                                                                              aaCCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATa                                                                                                                         <  7:93990/94‑1 (MQ=255)                                                                              aaCCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATa                                                                                                                         >  8:93990/1‑94 (MQ=255)                                                                                 ccAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCa                                                               >  3:30240/1‑149 (MQ=255)                                                                                            ggCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAgg                                                    >  5:9320/1‑149 (MQ=255)                                                                                                                  ttttCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGa                              <  1:42144/149‑1 (MQ=255)                                                                                                                      ccTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAg                                            <  1:265557/131‑1 (MQ=255)                                                                                                                      ccTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAg                                            >  2:265557/1‑131 (MQ=255)                                                                                                                          cGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGaa                       <  6:21295/148‑1 (MQ=255)                                                                                                                          cGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGaa                       <  4:11628/148‑1 (MQ=255)                                                                                                                                   aTGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCtgatga             >  6:275705/1‑149 (MQ=255)                                                                                                                                    tGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTgatgat            >  2:5359/1‑149 (MQ=255)                                                                                                                                     ggCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGcac                                                       >  3:228443/1‑105 (MQ=255)                                                                                                                                     ggCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGcac                                                       <  4:228443/105‑1 (MQ=255)                                                                                                                                       cgaCGCGCCCTGACCGGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGc         >  4:38005/3‑149 (MQ=255)                                                                                                                                       cAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTg                         >  2:228194/1‑133 (MQ=255)                                                                                                                                       cAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTg                         <  1:228194/133‑1 (MQ=255)                                                                                                                                       cAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGc         >  8:67979/1‑149 (MQ=255)                                                                                                                                              ccccGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGCCGAACTc  <  8:221058/149‑1 (MQ=255)                                                                                                                                              ccccGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGCCGAACTc  >  5:116104/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                  CGGCTGACCAGAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGCCGAACTC  >  NC_002947/156239‑156528

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