Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	157,486	G→T	7.1%	D80E (GAC→GAA)	citN ←	citrate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	157,486	0	G	T	7.1%	173.2 / 3.8	56	D80E (GAC→GAA)	citN	citrate transporter
Reads supporting (aligned to +/- strand):  ref base G (20/31);  new base T (2/2);  total (23/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

TCGGCAGCATGGCGGCGCAGCAGATCATGTAGGTGGTGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAGCATCATGGGGCCGATCTTGGCGGAAAAACCACCGAACAGGGCGAACAGGAT  >  NC_002947/157344‑157629                                                                                                                                                |                                                                                                                                                tCGGCAGCATGGCGGCGCAGCAGATCATGTAGGTGGTGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGaa                                                                                                                                            >  6:157751/1‑148 (MQ=255) tCGGCAGCATGGCGGCGCAGCAGATCATGTAGGTGGTGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGaa                                                                                                                                            >  4:357672/1‑148 (MQ=255)   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                                                                                                                                <  3:357672/149‑1 (MQ=255)            gcggcgCAGCAGATCATGTAGGTGGTGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATc                                                                                                                                <  5:157751/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                                                                                                            <  8:113276/149‑1 (MQ=255)                                ggtggtGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATaga                                                                                                            <  2:112090/149‑1 (MQ=255)                                ggtggtGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATaga                                                                                                            <  6:76933/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                                                                                                            >  2:14743/1‑130 (MQ=255)                                                   tcGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATaga                                                                                                            <  1:14743/130‑1 (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                                                                                                        <  2:322471/97‑1 (MQ=255)                                                                                         ccgaGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGAt                                                                                                        >  1:322471/1‑97 (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        >  3:242114/4‑148 (MQ=255)                                                                                                                                      ttccGGCCTGTTCGACCCGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAGCATCATGGGGCCGATCTTGGCGGAAAAACCACCGAACAGGGCGaa        >  6:14011/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)                                                                                                                                                |                                                                                                                                                TCGGCAGCATGGCGGCGCAGCAGATCATGTAGGTGGTGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAGCATCATGGGGCCGATCTTGGCGGAAAAACCACCGAACAGGGCGAACAGGAT  >  NC_002947/157344‑157629

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