Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	157,493	A→C	7.2%	L78R (CTG→CGG)	citN ←	citrate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	157,493	0	A	C	7.2%	166.6 / 4.7	56	L78R (CTG→CGG)	citN	citrate transporter
Reads supporting (aligned to +/- strand):  ref base A (20/32);  new base C (2/2);  total (22/34)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.33e-01

GGCAGCATGGCGGCGCAGCAGATCATGTAGGTGGTGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAGCATCATGGGGCCGATCTTGGCGGAAAAACCACCGAACAGGGCGAACAGGATCGGTACCAG  >  NC_002947/157346‑157638                                                                                                                                                     |                                                                                                                                                  ggCAGCATGGCGGCGCAGCAGATCATGTAGGTGGTGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACa                                                                                                                                                   <  5:373057/148‑1 (MQ=255)   cagcaTGGCGGCGCAGCAGATCATGTAGGTGGTGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGc                                                                                                                                                <  2:193415/149‑1 (MQ=255)    agcaTGGCGGCGCAGCAGATCATGTAGGTGGTGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCTGTTCGACCCGGcc                                                                                                                                               <  2:375974/149‑1 (MQ=255)     gcaTGGCGGCGCAGCAGATCATGTAGGTGGTGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCg                                                                                                                                              <  6:30110/149‑1 (MQ=255)        tGGCGGCGCAGCAGATCATGTAGGTGGTCGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCTGTTCGACCCGGCCgttc                                                                                                                                           <  7:109193/149‑4 (MQ=255)         ggcggaGAAGCAGCTCATGTAGGTGGTGGCACCGTCACCTTCGAGCGAGACCACAAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGAGGATAGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAAt                                                                                                                                          >  4:329187/1‑149 (MQ=255)          gcggcgCAGCAGATCATGTAGGTGGTGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATc                                                                                                                                         <  5:157751/149‑1 (MQ=255)          gcggcgCAGCAGATCATGTAGGTGGTGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATc                                                                                                                                         <  3:357672/149‑1 (MQ=255)            ggcgCAGCAGATCATGTAGGTGGTGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCAGGTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATc                                                                                                                                         >  1:112090/1‑148 (MQ=255)                    agATCATGTAGGTGGTGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAgg                                                                                                                               >  5:231638/1‑149 (MQ=255)                            tAGGTGGTGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGaaa                                                                                                                         >  8:260406/1‑147 (MQ=255)                              ggtggtGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATaga                                                                                                                     <  2:112090/149‑1 (MQ=255)                              ggtggtGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATaga                                                                                                                     <  6:76933/149‑1 (MQ=255)                              ggtggtGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATaga                                                                                                                     <  8:113276/149‑1 (MQ=255)                                tggtggCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAgg                                                                                                                   <  1:78462/149‑1 (MQ=255)                                 ggtggCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGa                                                                                                                  >  5:108054/1‑149 (MQ=255)                                  gtggCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGag                                                                                                                 <  5:271827/149‑2 (MQ=255)                                    ggCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCCGAATCGATCATAAGGGCGAAATAGAGGATgc                                                                                                               >  7:333040/1‑148 (MQ=255)                                     gCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGc                                                                                                              <  2:195104/149‑1 (MQ=255)                                          gtcaccgtcGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACa                                                                                                         >  2:41930/1‑149 (MQ=255)                                           tcaccgtcGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACAt                                                                                                        <  2:382088/149‑1 (MQ=255)                                                 tcGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATaga                                                                                                                     <  1:14743/130‑1 (MQ=255)                                                 tcGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATaga                                                                                                                     >  2:14743/1‑130 (MQ=255)                                                   gagcgagACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCa                                                                                                >  8:322292/1‑149 (MQ=255)                                                        agaCCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAAc                                                                                                          >  6:237740/1‑134 (MQ=255)                                                        agaCCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAAc                                                                                                          <  5:237740/134‑1 (MQ=255)                                                         gaCCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCgg                                                                                          <  7:322292/149‑1 (MQ=255)                                                          accaccAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGAGAAATAGAGGATGGCGAACATCAGCATCACGCCgg                                                                                          >  5:265593/1‑148 (MQ=255)                                                               caGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTg                                                                                    >  8:67759/1‑149 (MQ=255)                                                                  ggCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGa                                                                                 <  6:265593/149‑1 (MQ=255)                                                                         tACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGc                                                                           <  6:108054/149‑1 (MQ=255)                                                                          actggtgTGCCGACCGAGACTTTAAGCGGATCTCCCTTGACCAGCTTTAGGATTTTGCGTACGGCCGGGTCGAAAAGGCCGGAATCGATAATCAGTTCGAAATTGAGGATGGCGAACATCAGCATCACGCCGGTAGTTGCGTGCTTGc                                                                           <  8:333040/142‑1 (MQ=255)                                                                                gTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGcc                                                                    >  4:70628/1‑149 (MQ=255)                                                                                 tGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGccc                                                                   <  1:41930/149‑1 (MQ=255)                                                                                   ccgaccgaGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCctc                                                                 <  7:260406/149‑1 (MQ=255)                                                                                     gaccgaGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTa                                                               <  3:372839/149‑1 (MQ=255)                                                                                     gaccgaGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGACGGCTCGAACAGGCGGGAATCGATCATCAGGGCGAAAAAGAGGATGCCGAACATCAGCATCACGCCGGTTGGTGCGAGCTTCCTGATGCCctc                                                                 >  8:280521/1‑147 (MQ=255)                                                                                      accgaGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTac                                                              >  6:223415/1‑148 (MQ=255)                                                                                      accgaGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAg                                                              <  3:237590/149‑1 (MQ=255)                                                                                       ccgaGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGAt                                                                                                                 >  1:322471/1‑97 (MQ=255)                                                                                       ccgaGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGAt                                                                                                                 <  2:322471/97‑1 (MQ=255)                                                                                         gagaCTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGCTCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAGcat                                                           <  2:256103/149‑1 (MQ=255)                                                                                          agaCTTTCAGCGGTTCGCACTTGACCAGCTTGAGGATCTTGCGTACGGCCGGGTCGAACAGGCCGGAATCGATCATCATGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAGcatc                                                          <  6:195105/149‑1 (MQ=255)                                                                                                          gCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAGCATCATGGGGCCGATCTTg                                           <  7:67759/148‑1 (MQ=255)                                                                                                           cccTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAGCATCATGGGGCCGATCTTGGc                                         <  4:345714/149‑1 (MQ=255)                                                                                                                        tGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAGCATCATGGGGCCGATCTTGGCGGAAAAaccacc                             >  8:306597/1‑148 (MQ=255)                                                                                                                         gAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAGCATCATGGGGCCGATCTTGGCGGAAAAACCACCg                            >  2:158580/1‑148 (MQ=255)                                                                                                                                  gCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTTCGAGCTTGCTGATGCCCTCTAGCATCATGGGGCCGATCTTGGCGGAAAAACCACCGAACAg                       <  3:329187/144‑1 (MQ=255)                                                                                                                                  gCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAGCATCATGGGGCCGATCTTGGCGGAAAAACCACCGAACAgg                      <  7:306597/145‑1 (MQ=255)                                                                                                                                    ttccGGCCTGTTCGACCCGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAGCATCATGGGGCCGATCTTGGCGGAAAAACCACCGAACAGGGCGaa                 >  6:14011/4‑148 (MQ=255)                                                                                                                                    ttccGGCCTGTTCGACCCGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAGCATCATGGGGCCGATCTTGGCGGAAAAACCACCGAACAGGGCGaa                 >  3:242114/4‑148 (MQ=255)                                                                                                                                     aaCGGCCGGGTCTAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAGCATCATGGGGCCGATCTTGGCGGAAAAACCACCGAACAGGGCGAACa               <  6:231638/149‑1 (MQ=255)                                                                                                                                      aCGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAGCATCATGGGGCCGATCTTGGCGGAAAAACCACCGAACAGGGCGAACAg              <  7:16762/149‑1 (MQ=255)                                                                                                                                         gCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAGCATCATGGGGCCGAGCTTGGCGGAAAAACCACCGAACAGGGCGAACAGGag           >  6:232061/1‑148 (MQ=255)                                                                                                                                                 gatcAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAGCATCATGGGGCCGATCTTGGCGGAAAAACCACCGAACAGGGCGAACAGGATCGGTACCa   <  5:223415/146‑1 (MQ=255)                                                                                                                                                 gAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAGCATCATGGGGCCGATCTTGGCGGAAAAACCACCGAACAGGGCGAACAGGATCGGTACCa   <  3:70628/149‑1 (MQ=255)                                                                                                                                                  aaCAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAGCATCATGGGGCCGATCTTGGCGGAAAAACCACCGAACAGGGCGAACAGGATCAGtccagg  >  3:112251/1‑144 (MQ=255)                                                                                                                                                     |                                                                                                                                                  GGCAGCATGGCGGCGCAGCAGATCATGTAGGTGGTGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAGCATCATGGGGCCGATCTTGGCGGAAAAACCACCGAACAGGGCGAACAGGATCGGTACCAG  >  NC_002947/157346‑157638

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