Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	157,141	Δ1 bp	9.8%	coding (585/1308 nt)	citN ←	citrate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	157,141	0	A	.	9.8%	153.9 / 11.0	51	coding (585/1308 nt)	citN	citrate transporter
Reads supporting (aligned to +/- strand):  ref base A (21/25);  new base . (2/3);  total (23/28)
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

CCATCAGGGCAGCGGTCAGGGCACCGTTGAAGTACAGCAGCTTGGGCCGGCGTGCGTCCGGGTACTGCGAGACGGTGATTTCGCTGTGGTCGATGTCGTCGGTAGGCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTACGAAGATGTCTGACGGGTCTACGTGCAGGGCGCTGGCAGCACGCGCGGTCGGCCCACCCCATGGGGTCATGTTCATCACGCCGCCGG  >  NC_002947/156996‑157289                                                                                                                                                  |                                                                                                                                                      ccATCAGGGCAGTGGTCAGGGCACCGTTGAAGTACACCAGCTTGGGCCGGCGTGCGTCCGGGTACTGCGAGACGGTGATTTCGCTGTGGTCGATGTCGTCGGTAGGCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCCATAc                                                                                                                                                    <  2:353770/149‑1 (MQ=255)   aTCAGGGCAGCGGTCAGGGCACCGTTGAAGTACAGCAGCTTGGGCCGGCGTGCGTCCGGGTACTGCGAGACGGTGATTTCGCTGTGGTCGATGTCGTCGGTAGGCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCCATACAt                                                                                                                                                  >  4:97397/1‑149 (MQ=255)            gCGGTCAGGGCACCGTTGAAGTACAGCAGCTTGGGCCGGCGTGCGTCCGGGTACTGCGAGACGGTGATTTCGCTGTGGTCGATGTCGTCGGTAGGCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAAt                                                                                                                                        <  1:327384/149‑1 (MQ=255)             cGGTCAGGGCACCGTTGAAGTACAGCAGCTTGGGCCGGCGTGCGTCCGGGTACTGCGAGACGGTGATTTCGCTGTGGTCGATGTCGTCGGTAGGCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCC‑TACATGTATGGCaagc                                                                                                                                       <  1:163777/149‑3 (MQ=255)                   gggCACCGTTGAAGTACAGCAGCTTGGGCCGGCGTGCGTCCGGGTACTGCGAGACGGTGATTTCGCTGTGGTCGATGTCGTCGGTAGGCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCC‑TACATGTATGGCaagcgtga                                                                                                                                   >  8:30840/1‑141 (MQ=255)                   gggCACCGTTGAAGTACAGCAGCTTGGGCCGGCGTGCGTCCGGGTACTGCGAGACGGTGATTTCGCTGTGGTCGATGTCGTCGGTAGGCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCC‑TACATGTATGGCaagcgtga                                                                                                                                   <  7:30840/147‑7 (MQ=255)                      cACCGTTGAAGTACAGCAGCTTGGGCCGGCGTGCGTCCGGGTACTGCGAGACGGTGATTTCGCTGTGGTCGATGTCGTCGGTAGGCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATc                                                                                                                              >  6:352008/1‑149 (MQ=255)                          ggtGAAGGAGAGCAGCTTGGGCCGGCGTGCGTCCGGGTACTGCGAGACGGTGATTTCGCTGTGGTCGATGTCGTCGGTAGGCAGATGCAGTTCGCCCAGGCGGGCACGCTCTCGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCa                                                                                                                          <  7:264265/147‑1 (MQ=255)                           ttGAAGTACAGCAGCTTGGGCCGGCGTGCGTCCGGGTACTGCGAGACGGTGATTTCGCTGTGGTCGATGTCGTCGGTAGGCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAg                                                                                                                         <  5:287548/149‑1 (MQ=255)                             gAAGTACAGCAGCTTGGGCCGGCGTGCGTCCGGGTACTGCGAGACGGTGATTTCGCTGTGGTCGATGTCGTCGGTAGGCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCC‑TACATGTATGGCaagc                                                                                                                                       <  4:225689/133‑3 (MQ=255)                             gAAGTACAGCAGCTTGGGCCGGCGTGCGTCCGGGTACTGCGAGACGGTGATTTCGCTGTGGTCGATGTCGTCGGTAGGCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCC‑TACATGTATGGCaagc                                                                                                                                       >  3:225689/1‑131 (MQ=255)                                  aCAGCAGCTTGGGCCGGCGTGCGTCCGGGTACTGCGAGACGGTGATTTCGCTGTGGTCGATGTCGTCGGTAGGCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACAcc                                                                                                                   >  2:277422/1‑148 (MQ=255)                                    agcagcTTGGGCCGGCGTGCGTCCGGGTACTGCGAGACGGTGATTTCGCTGTGGTCGATGTCGTCGGTAGGCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGaa                                                                                                                <  7:349201/149‑1 (MQ=255)                                    agcagcTTGGGCCGGCGTGCGTCCGGGTACTGCGAGACGGTGATTTCGCTGTGGTCGATGTCGTCGGTAGGCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGaa                                                                                                                >  3:227748/1‑149 (MQ=255)                                      cagcTTGGGCCGGCGTGCGTCCGGGTACTGCGAGACGGTGATTTCGCTGTGGTCGATGTCGTCGGTAGGCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGaac                                                                                                               >  5:5943/1‑148 (MQ=255)                                      cagcTTGGGCCGGCGTGCGTACGGGTACTGCGAGACGGTGATTTCGCTGTGGTCGATGTCGTCGGTAGGCAGATGCAGTTAGCCCAGGAGGGAACGATCACGCTTGCCATACATGTA‑GGAAATGGCCAGGATCGCCAGTACACAGaac                                                                                                               >  4:161079/1‑148 (MQ=255)                                        gcTTGCGCCGGCGTGCGTCCGGGTACTGCGAGACGGTGATTACGCTGTGGTCGATGTCGTCTGTAGGCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGATTGCCAAACATGTA‑GGCAATCGCCAGGATCGCCAGTACAAAGaacaa                                                                                                             >  2:377664/1‑148 (MQ=255)                                         cTTGGGCCGGCGTGCGTCCGGGTACTGCGAGACGGTGATTTCGCTGTGGTCGATGTCGTCGGTAGGCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAacaaca                                                                                                           >  5:358559/1‑149 (MQ=255)                                         cTTGGGCCGGCGTGCGTCCGGGTACTGCGAGACGGTGATTTCGCTGTGGTCGATGTCGTCGGTAGGCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAacaaca                                                                                                           <  2:217860/149‑1 (MQ=255)                                                        gtCCGGGTACCGCGAGACGGTGATTTCGCTGTGGTCGATGTCGTCGGTAGGCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCg                                                                                            >  6:169998/1‑149 (MQ=255)                                                           cGGGTACTGCGAGACGGTGATTTCGCTGTGGTCGATGTCGTCGGTAGGCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTa                                                                                         <  4:227748/149‑1 (MQ=255)                                                             ggTACTGCGAGTCGGTGAATTCGCTGTGGTCGTTGTCGTCGGTACGCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTACg                                                                                       <  1:336359/149‑1 (MQ=255)                                                               tACTGCGAGACGGTGATTTCGCTGTGGTCGATGTCGTCGGTAGGCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTACGaa                                                                                     <  3:97397/149‑1 (MQ=255)                                                                    cGAGACGGTGATTTCGCTGTGGTCGATGTCGTCGGTAGGCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGTTACGAAGATGt                                                                                >  2:51246/1‑149 (MQ=255)                                                                         cGGTTATTTCGGTGTGGTCGATGTCGTCGTTAGGCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCTAACAACATGTCCTGTATCATCGGTACGAAGATGTCTGAc                                                                           <  1:377664/149‑1 (MQ=255)                                                                          ggTGATTTCGCTGTGGTCGATGTCGTCGGTAGGCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTACGAAGATGTCTGACg                                                                          >  3:290289/1‑149 (MQ=255)                                                                          ggTGATTTCGCTGTGGTCGATGTCGTCGGTAGGCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTACGAAGATGTCTGACg                                                                          >  8:354645/1‑149 (MQ=255)                                                                             gATTTCGCTGTGGTCGATGTCGTCGGTAGGCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTACGAAGATGTCTGACGGGt                                                                       >  8:91389/1‑149 (MQ=255)                                                                                  cGCTGTGGTCGATGTCGTCGGTAGGCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTACGAAGATGTCTGACGGGTCTAc                                                                   <  5:352008/148‑1 (MQ=255)                                                                                   gCTGTGGTCGATGTCGTCGGTAGGCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTACGAAGATGTCTGACGGGTCTACGt                                                                 <  1:277422/149‑1 (MQ=255)                                                                                    cTGTGGTCGATGTCGTCGGTAGGCAGATGGAGTTCGGCCACGCTGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGTCCGGGATCATCGGTACGAAGATGTCTGACGGGTCTACGTg                                                                <  7:209935/149‑1 (MQ=255)                                                                                             aTGTCGTCGGTAGGCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTACGAAGATGTCTGACGGGTCTACGTGCAGGGCGCt                                                       >  1:132356/1‑149 (MQ=255)                                                                                              tGTCGTCGGTAGGCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTACGAAGATGTCTGACGGGTCTACGTGCAGGGCGCTg                                                      <  7:310199/149‑1 (MQ=255)                                                                                                tcgtcgGTAGGCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTACGAAGATGTCTGACGGGTCTACGTGCAGGGCGCTGgc                                                    <  5:378661/149‑1 (MQ=255)                                                                                                       tAGGCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTACGAAGATGTCTGACGGGTCTACGTGCAGGGCGCTGGCAGCAcgc                                             <  4:266305/149‑1 (MQ=255)                                                                                                         ggCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTACGAAGATGTCTGACGGGTCTACGTGCAGGGCGCTGGCAGCAcgcg                                            <  1:51246/148‑1 (MQ=255)                                                                                                              aTGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTACGAAGATGTCTGACGGGTCTACGTGCAGGGCGCTGGCAGCACGCGCGGTCg                                      >  5:145691/1‑149 (MQ=255)                                                                                                                gCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTACGAAGATGTCTGACGGGTCTACGTGCAGGGCGCTGGCAGCACGCGCGGTCGGc                                    >  7:189252/1‑149 (MQ=255)                                                                                                                 cAGTTCGCCCAGGCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTACGAAGATGTCTGACGGGTCTACGTGCATGGCGCTGGCAGCACGCGCGGTCGGc                                    <  2:132356/148‑1 (MQ=255)                                                                                                                    ttCGCCCAGGCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTACGAAGATGTCTGACGGGTCTACGTGCAGGGCGCTGGCAGCACGCGCGGTCGGCCCAc                                >  4:366623/1‑149 (MQ=255)                                                                                                                          cAGGCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTACGAAGATGTCTGACGGGTCTACGTGCAGGGCGCTGGCAGCACGCGCGGTCGGCCCACCCCATg                          <  4:290289/149‑1 (MQ=255)                                                                                                                           aGGCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTACGAAGATGTCTGACGGGTCTACGTGCAGGGCGCTGGCAGCACGCGCGGTCGGCCCACCCCATgg                         <  3:304888/149‑1 (MQ=255)                                                                                                                            ggCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTACGAAGATGTCTGACGGGTCTACGTGCAGGGCGCTGGCAGCACGCGCGGTCGGCCCACCCCATggg                        <  3:366623/149‑1 (MQ=255)                                                                                                                            ggCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTACGAAGATGTCTGACGGGTCTACGTGCAGGGCGCTGGCAGCACGCGCGGTCGGCCCACCCCATggg                        <  7:9400/149‑1 (MQ=255)                                                                                                                            ggCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTACGAAGATGTCTGACGGGTCTACGTGCAGGGCGCTGGCAGCACGCGCGGTCGGCCCACCCCATggg                        <  7:185345/149‑1 (MQ=255)                                                                                                                                   acgctcacgctTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTACGAAGATGTCTGACGGGTCTACGTGCAGGGCGCTGGCAGCACGCGCGGTCGGCCCACCCCATGGGGTCATGt                 >  1:334245/1‑149 (MQ=255)                                                                                                                                     gctcacgctTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTACGAAGATGTCTGACGGGTCTACGTGCAGGGCGCTGGCAGCACGCGCGGTCGGCCCACCCCATGGGGTCATGtt                >  7:112983/1‑148 (MQ=255)                                                                                                                                      ctcacgctTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTACGAAGATGTCTGACGGGTCTACGTGCAGGGCGCTGGCAGCACGCGCGGTCGGCCCACCCCATGGGGTCATGTtca              <  6:5943/149‑1 (MQ=255)                                                                                                                                       tcacgctTGCCATACATGTA‑GGCAATGGCCAGGATCGACAGTACACCGAACAACATGGCCGGGATCATCGGTACGAAGATGTCTGACGGGTCTACGTGCAGGGCGCTGGCAGCACGCGCGGTCGGCCCACCCCATGGGGTCATGTtca              >  1:68602/1‑148 (MQ=255)                                                                                                                                               gCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTACGAAGATGTCTGACGGGTCTACGTGCAGGGCGCTGGCAGCACGCGCGGTCGGCCCACCCCATGGGGTCATGTTCATCAcgcc       <  3:163285/147‑1 (MQ=255)                                                                                                                                                 cATACATGTA‑GGCAATGGCCAGGATCGACAGTACACCGAACAACATGGCCGGGATCATCGGTACGAAGATGTCTGACGGGTCTACGTGCAGGGCGCTGGCAGCACTCGCGGTCGGCCAACCCCATGGGGTCATGTTCATCACGccgccg   >  6:20632/1‑149 (MQ=255)                                                                                                                                                  aTACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTACGAAGATGTCTGACGGGTCTACGTGCAGGGCGCTGGCCGCACGCGCGGTCGGCCCACCCCATGGGGTCATGTTCATAACGCCGCCgg  >  4:308948/1‑149 (MQ=255)                                                                                                                                                  aTACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTACGAAGATGTCTGACGGGTCTACGTGCAGGGCGCTGGCAGCACGCGCGGTCGGCCCACCCCATGGGGTCATGTTCATCACGCCGCCgg  >  6:50273/1‑149 (MQ=255)                                                                                                                                                  aTACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTACGAAGATGTCTGACGGGTCTACGTGCAGGGCGCTGGCAGCACGCGCGGTCGGCCCACCCCATGGGGTCATGTTCATCACGCCGCCgg  >  5:279969/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                      CCATCAGGGCAGCGGTCAGGGCACCGTTGAAGTACAGCAGCTTGGGCCGGCGTGCGTCCGGGTACTGCGAGACGGTGATTTCGCTGTGGTCGATGTCGTCGGTAGGCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTACGAAGATGTCTGACGGGTCTACGTGCAGGGCGCTGGCAGCACGCGCGGTCGGCCCACCCCATGGGGTCATGTTCATCACGCCGCCGG  >  NC_002947/156996‑157289

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