Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	157,149:1	+T	8.9%	coding (577/1308 nt)	citN ←	citrate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	157,149	1	.	T	8.9%	207.1 / 15.7	56	coding (577/1308 nt)	citN	citrate transporter
Reads supporting (aligned to +/- strand):  ref base . (26/25);  new base T (2/3);  total (28/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 = 9.16e-01

GCGGTCAGGGCACCGTTGAAGTACAGCAGCTTGGGCCGGCGTGCGTCCGGGTACTGCGAGACGGTGATTTCGCTGTGGTCGATGTCGTCGGTAGGCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTACGAAGATGTCTGACGGGTCTACGTGCAGGGCGCTGGCAGCACGCGCGGTCGGCCCACCCCATGGGGTCATGTTCATCACGCCGCCGGCAAGGAT  >  NC_002947/157007‑157296                                                                                                                                                |                                                                                                                                                    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                                                                                                                                              >  3:225689/1‑131 (MQ=255)                  gAAGTACAGCAGCTTGGGCCGGCGTGCGTCCGGGTACTGCGAGACGGTGATTTCGCTGTGGTCGATGTCGTCGGTAGGCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCC‑TACATGTATGGCaagc                                                                                                                                              <  4:225689/133‑3 (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                                                                                                                  <  2:217860/149‑1 (MQ=255)                              cTTGGGCCGGCGTGCGTCCGGGTACTGCGAGACGGTGATTTCGCTGTGGTCGATGTCGTCGGTAGGCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAacaaca                                                                                                                  >  5:358559/1‑149 (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                               <  7:9400/149‑1 (MQ=255)                                                                                                                 ggCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTACGAAGATGTCTGACGGGTCTACGTGCAGGGCGCTGGCAGCACGCGCGGTCGGCCCACCCCATggg                               <  7:185345/149‑1 (MQ=255)                                                                                                                 ggCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTACGAAGATGTCTGACGGGTCTACGTGCAGGGCGCTGGCAGCACGCGCGGTCGGCCCACCCCATggg                               <  3:366623/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)                                                                                                                                          cATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTACGAAGATGTCTGACGGGTCTACGTGCAGGGCGCTGGCAGCACGCGCGGTCGGCCCACCCCATGGGGTCATGTTCATCACGCCGCCGGCaa      <  7:354645/149‑1 (MQ=255)                                                                                                                                             gTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTACGAAGATGTCTGACGGGTCTACGTGCAGGGCGCTGGCAGCACGCGCGGTCGGCCCACCCCATGGGGTCATGTTCATCACGCCGCCGGCAAGga   >  7:111166/1‑149 (MQ=255)                                                                                                                                             gTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTACGAAGATGTCTGACGGGTCTACGTGCAGGGCGCTGGCAGCACGCGCGGTCGGCCCACCCCATGGGGTCATGTTCATCACGCCGCCGGCAAGga   >  8:13615/1‑149 (MQ=255)                                                                                                                                              tA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTACGAAGATGTCTGACGGGTCTACGTGCAGGGCGCTGGCAGCACGCGCGGTCGGCCCACCCCATGGGGTCATGTTCATCACGCCGCCGGCAAGgat  >  8:5097/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                    GCGGTCAGGGCACCGTTGAAGTACAGCAGCTTGGGCCGGCGTGCGTCCGGGTACTGCGAGACGGTGATTTCGCTGTGGTCGATGTCGTCGGTAGGCAGATGCAGTTCGCCCAGGCGGGCACGCTCACGCTTGCCATACATGTA‑GGCAATGGCCAGGATCGCCAGTACACCGAACAACATGGCCGGGATCATCGGTACGAAGATGTCTGACGGGTCTACGTGCAGGGCGCTGGCAGCACGCGCGGTCGGCCCACCCCATGGGGTCATGTTCATCACGCCGCCGGCAAGGAT  >  NC_002947/157007‑157296

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