Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	156,388	A→G	10.0%	intergenic (‑13/+30)	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,388	0	A	G	10.0%	105.9 / 3.3	40	intergenic (‑13/+30)	ygdQ/citN	transport protein/citrate transporter
Reads supporting (aligned to +/- strand):  ref base A (20/16);  new base G (2/2);  total (22/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 = 9.51e-01

CAGAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGCCGAACTC  >  NC_002947/156247‑156528                                                                                                                                              |                                                                                                                                             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                                                                                                       >  1:20868/1‑149 (MQ=255)                                 cGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGacta                                                                                                       <  2:220500/149‑1 (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                                                                                                                          <  7:91729/117‑1 (MQ=255)                                              aaGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGAt                                                                                                                          >  8:91729/1‑117 (MQ=255)                                                          ggTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCAt                                                                                >  5:21295/1‑147 (MQ=255)                                                           gTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATcag                                                                             >  3:11628/1‑149 (MQ=255)                                                               aaaaGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGt                                                                         <  2:20868/149‑1 (MQ=255)                                                                      aaCCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATa                                                                                                                         >  8:93990/1‑94 (MQ=255)                                                                      aaCCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATa                                                                                                                         <  7:93990/94‑1 (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  >  5:116104/1‑149 (MQ=255)                                                                                                                                      ccccGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGCCGAACTc  <  8:221058/149‑1 (MQ=255)                                                                                                                                              |                                                                                                                                             CAGAATCGAGATCATGATGATGTTGTCGATACCGAGGACGATCTCAAGCGCGGTAAGGGTGAAAAAGGCAACCCAGATTTCCGGGCTGGTCAGCCATTCCATGTATTTTCCTTCGAGCGGTAATGGCAACGCGCCCCGGTCAGGGCGCGTCGCGTACGTGATACGGATGTATTAAAGACTACTGAACAGCGGGAAAATCCCCATCAGCAGTGCGCCGAGCATTATGCACAGGCATACCAGCACTGCCCACTTGAGCGTGAAGCGCTGATGATCGCCGAACTC  >  NC_002947/156247‑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: ‑