Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	157,484	G→T	7.1%	P81Q (CCG→CAG)	citN ←	citrate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	157,484	0	G	T	7.1%	175.8 / 4.1	56	P81Q (CCG→CAG)	citN	citrate transporter
Reads supporting (aligned to +/- strand):  ref base G (21/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 = 2.98e-01

CTGTACATCGGCAGCATGGCGGCGCAGCAGATCATGTAGGTGGTGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAGCATCATGGGGCCGATCTTGGCGGAAAAACCACCGAACAGGGCGAACAGGAT  >  NC_002947/157337‑157629                                                                                                                                                     |                                                                                                                                                  cTGTACATCGGCAGCATGGCGGCGCAGCAGATCATGTAGGTGGTGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCgg                                                                                                                                                  >  5:30110/1‑149 (MQ=255)        tCGGCAGCATGGCGGCGCAGCAGATCATGTAGGTGGTGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGaa                                                                                                                                            >  4:357672/1‑148 (MQ=255)        tCGGCAGCATGGCGGCGCAGCAGATCATGTAGGTGGTGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGaa                                                                                                                                            >  6:157751/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                                                                                                                                <  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                                                                                                            <  8:113276/149‑1 (MQ=255)                                       ggtggtGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATaga                                                                                                            <  6:76933/149‑1 (MQ=255)                                       ggtggtGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATaga                                                                                                            <  2:112090/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                                                                                                        <  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)                                                                                                                                                     |                                                                                                                                                  CTGTACATCGGCAGCATGGCGGCGCAGCAGATCATGTAGGTGGTGGCACCGTCACCGTCGAGCGAGACCACCAGGGCCA‑GTACAGCGGTGCCGACCGAGACTTTCAGCGGATCGCCCTTGACCAGCTTGAGGATCTTGCGAACGGCCGGGTCGAACAGGCCGGAATCGATCATCAGGGCGAAATAGAGGATGGCGAACATCAGCATCACGCCGGTAGGTGCGAGCTTGCTGATGCCCTCTAGCATCATGGGGCCGATCTTGGCGGAAAAACCACCGAACAGGGCGAACAGGAT  >  NC_002947/157337‑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: ‑