Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,385,015	C→G	7.7%	A43P (GCG→CCG)	nasA ←	nitrate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,385,015	0	C	G	7.7%	159.2 / 4.7	52	A43P (GCG→CCG)	nasA	nitrate transporter
Reads supporting (aligned to +/- strand):  ref base C (28/20);  new base G (2/2);  total (30/22)
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.23e-01

GGCCGATCAGGCCAGCCGTCTTGGGCGACAGGCGGTCGACCAGAACGCCCATGGCGAAGCGCAGGATTGCCCCGGCCA‑GGATCGGCGTGGCCACCATCAGGCCCCGTTGTTGTGCGCTCAGTTGCAGGTCGGCAGCAATCTGCACCGCCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTCGAAGTACAGAAACGCGGCGAACAGCGTGGGCACATGCCCGGATTTCCAGAAGCTGGCACTCATCGAACACCTCACTCGGTAGTACAACGGAACGAAAAAGACGCCGC  >  NC_002947/2384868‑2385161                                                                                                                                                     |                                                                                                                                                   ggCCGATCAGGCCAGCCGTCTTGGGCGACAGGCGGTCGACCAGAACGCCCATGGCGAAGCGCAGGATTGCCCCGGCCA‑GGATCGGCGTGGCCACCATCAGGCCCCGTTGTTGTGCGCTCAGTTGCAGGTCGGCAGCAATCTGCACCGcc                                                                                                                                                   >  6:316833/1‑149 (MQ=255) ggCCGATCAGGCCAGCCGTCTTGGGCGACAGGCGGTCGACCAGAACGCCCATGGCGAAGCGCAGGATTGCCCCGGCCA‑GGATCGGCGTGGCCACCATCAGGCCCCGTTGTTGTGCGCTCAGTTGCAGGTCGGCAGCAATCTGCACCGc                                                                                                                                                    >  7:155321/1‑148 (MQ=255)  gCCGATCAGGCCAGCCGTCTTGGGCGACAGGCGGTCGACCAGAACGCCCATGGCGAAGCGCAGGATTGCCCCGGCCA‑GGATCGGCGTGGCCACCATCAGGCCCCGTTGTTGTGCGCTCAGTTGCAGGTCGGCAGCAATCTGCACCGc                                                                                                                                                    <  5:172256/147‑1 (MQ=255)           gCCAGCCGTCTTGGGCGACAGGCGGTCGACCAGAACGCCCATGGCGAAGCGCAGGATTGCCCCGGCCA‑GGATCGGCGTGGCCACCATCAGGCCCCGTTGTTGTGCGCTCAGTTGCAGGTCGGCAGCAATCTGCACCGCCAACGGGCCca                                                                                                                                         <  1:8620/149‑1 (MQ=255)            ccAGCCGTCTTGGGCGACAGGCGGTCGACCAGAACGCCCATGGCGAAGCGCAGGATTGCCCCGGCCA‑GGATCGGCGTGGCCACCATCAGGCCCCGTTGTTGTGCGCTCAGTTGCAGGTCGGCAGCAATCTGCACCGCCAACGGGCCcag                                                                                                                                        >  4:62955/1‑149 (MQ=255)             cAGCCGTCTTGGGCGACAGGCGGTCGACCAGAACGCCCATGGCGAAGCGCAGGATTGCCCCGGCCAGGGATCGGCGTGGCCACCATCAGGCCCCGTTGTTGTGCGATCAGTTGCAGGTCGGCAGCAATCTGCACCGCCAACGGGCCcag                                                                                                                                        >  1:310121/1‑149 (MQ=255)                 cGTCTTGGGCGACAGGCGGTCGACCAGAACGCCCATGGCGAAGCGCAGGATTGCCCCGGCCA‑GGATCGGCGTGGCCACCATCAGGCCCCGTTGTTGTGCGCTCAGTTGCAGGTCGGCAGCAATCTGCACCGCCAACGGGCCCAGCAGGt                                                                                                                                   <  3:111492/149‑1 (MQ=255)                  gTCTTGGGCGACAGGCGGTCGACCAGAACGCCCATGGCGAAGCGCAGGATTGCCCCGGCCA‑GGATCGGCGTGGCCACCATCAGGCCCCGTTGTTGTGCGCTCAGTTGCAGGTCGGCAGCAATCTGCACCGCCAACGGGCCCAGCAGGTa                                                                                                                                  >  1:253735/1‑149 (MQ=255)                   tCTTGGGCGACAGGCGGTCGACCAGAACGCCCATGGCGAAGCGCAGGATTGCCCCGGCCA‑GGATCGGCGTGGCCACCATCAGGCCCCGTTGTTGTGCGCTCAGTTGCAGGTCGGCAGCAATCTGCACCGCCAACGGGCCCAGCAGGTAc                                                                                                                                 <  2:285920/149‑1 (MQ=255)                   tCTTGGGCGACAGGCGGTCGACCAGAACGCCCATGGCGAAGCGCAGGATTGCCCCGGCCA‑GGATCGGCGTGGCCACCATCAGGCCCCGTTGTTGTGCGCTCAGTTGCAGGTCGGCAGCAATCTGCACCGCCAACGGGCCCAGCAGGTAc                                                                                                                                 <  2:303604/149‑1 (MQ=255)                       gggCGACAGGCGGTCGACCAGAACGCCCATGGCGAAGCGCAGGATTGCCCCGGCCA‑GGATCGGCGTGGCCACCATCAGGCCCCGTTGTTGTGCGCTCAGTTGCAGGTCGGCAGCAATCTGCACCGCCAACGGGCCCAGCAGGTACCAGa                                                                                                                             >  7:145022/1‑149 (MQ=255)                          cGACAGGCGGTCGACCAGAACGCCCATGGCGAAGCGCAGGATTGCCCCGGCCA‑GGATCGGCGTGGCCACCATCAGGCCCCGTTGTTGTGCGCTCAGTTGCAGGTCGGCAGCAATCTGCACCGCCAACGGGCCCAGCAGGTACCAGACCa                                                                                                                          <  8:145022/149‑1 (MQ=255)                                 cGGTCGACCAGAACGCCCATGGCGAAGCGCAGGATTGCCCCGGCCA‑GGATCGGCGTGGCCACCATCAGGCCCCGTTGTTGTGCGCTCAGTTGCAGGTCGGCAGCAATCTGCACCGCCAACGGGCCCAGCAGGTACCAGACCATGAAGCt                                                                                                                   >  4:206168/1‑149 (MQ=255)                                      gACCAGAACGCCCATGGCGAAGCGCAGGATTGCCCCGGCCA‑GGATCGGCGTGGCCACCATCAGGCCCCGTTGTTGTGCGCTCAGTTGCAGGTCGGCAGCAATCTGCACCGCCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGt                                                                                                              >  6:143763/1‑149 (MQ=255)                                       aCCAGAACGCCCATGGCGAAGCGCAGGATTGCCCCGGCCA‑GGATCGGCGTGGCCACCATCAGGCCCCGTTGTTGTGCGCTCAGTTGCAGGTCGGCAGCAATCTGCACCGCCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTc                                                                                                             >  5:62874/1‑149 (MQ=255)                                                       cGAAGCGCAGGATTGCCCCGGCCA‑GGATCGGCGTGGCCACCATCAGGCCCCGTTGTTGTGCGCTCAGTTGCAGGTCGGCAGCAATCTGCACCGCCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTCGAAGTACAGAAAcgcg                                                                                             <  6:62874/149‑1 (MQ=255)                                                        gAAGCGCAGGATTGCCCCGGCCA‑GGATCGGCGTGGCCACCATCAGGCCCCGTTGTTGTGCGCTCAGTTTCAGGTCGGCAGCAATCTGCACCGCCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTCGAAGTACAGAAACgcgg                                                                                            <  1:93785/149‑1 (MQ=255)                                                             gcAGGATTGCCCCGGCCA‑GGATCGGCGTGGCCACCATCAGGCCCCGTTGTTGTGCGCTCAGTTGCAGGTCGGCAGCAATCTGCACCGCCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTCGAAGTACAGAAACGCGGCGaa                                                                                        >  3:301975/1‑148 (MQ=255)                                                             gcAGGATTGCCCCGGCCA‑GGATCGGCGTGGCCACCATCAGGCCCCGTTGTTGTGCGCTCAGTTGCAGGTCGGCAGCAATCTGCACCGCCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTCGAAGTACAGAAACGCGGCGaa                                                                                        >  3:307104/1‑148 (MQ=255)                                                             gcAGGATTGCCCCGGCCA‑GGATCGGCGTGGCCACCATCAGGCCCCGTTGTTGTGCGCTCAGTTGCAGGTCGGCAGCAATCTGCACCGCCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTCGAAGTACAGAAACGCGGCGAAc                                                                                       <  7:325777/149‑1 (MQ=255)                                                                   ttGCCCCGGCCA‑GGATCGGCGTGGCCACCATCAGGCCCCGTTGTTGTGCGCTCAGTTGCAGGTCGGCAGCAATCTGCACCGCCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTCGAAGTACAGAAACGCGGCGAACAGCGTg                                                                                 <  2:19795/149‑1 (MQ=255)                                                                      ccccGGCCA‑GGATCGGCGTGGCCACCATCAGGCCCCGTTGTTGTGCGCTCAGTTGCAGGTCGGCAGCAATCTGCACCGCCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTCGAACTACAGAAAGGCGGCGAACCGCGTGGGc                                                                              >  8:366270/1‑149 (MQ=255)                                                                        ccGGCCA‑GGATCGGCGTGGCCACCATCAGGCCCCGTTGTTGTGCGCTCAGTTGCAGGTCGGCAGCAATCTGCACCGCCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTCGAAGTACAGAAACGCGGCGAACAGCGTGGGca                                                                             >  4:118604/1‑148 (MQ=255)                                                                          ggCCA‑GGATCGGCGTGGCCACCATCAGGCCCCGTTGTTGTGCGCTCAGTTGCAGGTCGGCAGCAATCTGCACCGCCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTCGAAGTACAGAAACGCGGCGAACAGCGTGGGcac                                                                            >  2:300103/1‑147 (MQ=255)                                                                              a‑GGATCGGCGTGGCCACCATCAGGCCCCGTTGTTGTGCGCTCAGTTGCAGGTCGGCAGCAATCTGCACCGCCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTCGAAGTACAGAAACGCGGCGAACAGCGTGGGCACATGccc                                                                      <  2:253735/149‑1 (MQ=255)                                                                                   tCGGCGTGGCCACCATCAGGCCCCGTTGTTGTGCTCTCAGTTGCAGTTCGGCAGCAATCTGCACCGCCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTCGAAGTACAGAAACGCGGCGAACAGCGTGGGCACATGCCCGGAt                                                                  <  3:49840/149‑1 (MQ=255)                                                                                   tCGGCGTGGCCACCATCAGGCCCCGTTGTTGTGCGCTCAGTTGCAGGTCGGCAGCAATCTGCACCGCCAACGGGCCCAGTAGGTACCAGACCATGAAGCTGAGGTCGAAGTACAGAAACGCGGCGAACAGCGTGGGCACATGCCCGGAt                                                                  <  8:155321/149‑1 (MQ=255)                                                                                             caccaTCAGGCCCCGTTGTTGTGCGCTCAGTTGCAGGTCGGCAGCAATCTGCACCGCCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTCGAAGTACAGAAACGCGGCGAACAGCGTGGGCACATGCCCGGATTTCCAGAAGc                                                        >  2:231854/1‑149 (MQ=255)                                                                                                caTCAGGCCCCGTTGTTGTGCGCTCAGTTGCAGGTCCGCAGCAATCTGCACTGCCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTCGAAGTACAGAAACGCGGCGAACAGCGTGGGCACATGCCCGGATTTCCAGAAGCTgg                                                     <  5:363558/149‑1 (MQ=255)                                                                                                  tCAGGCCCCGTTGTTGTGCGCTCAGTTGCAGGTCGGCAGCAATCTGCACCGCCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTCGAAGTACAGAAACGCGGCGAACAGCGTGGGCACATGCCCGGATTTCCAGAAGCTGGCa                                                   >  3:96775/1‑149 (MQ=255)                                                                                                    aGGCCCCGTTGTTGTGCGCTCAGTTGCAGGTCGGCAGCAATCTGCACCGCCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTCGAAGTACAGAAACGCGGCGAACAGCGTGGGCACATGCCCGGATTTCCAGAAGCTGGCAc                                                  <  3:206168/148‑1 (MQ=255)                                                                                                    aGGCCCCGTTGTTGTGCGCTCAGTTGCAGGTCGGCAGCAATCTGCACCGCCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTCGAAGTACAGAAACGCGGCGAACAGCGTGGGCACATGCCCGGATTTCCAGAAGCTGGCACt                                                 >  6:196894/1‑149 (MQ=255)                                                                                                        cccGTTGTTGTGCGCTCAGTTGCAGGTCGGCAGCAATCTGCACCGCCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTCGAAGTACAGAAACGCGGCGAACAGCGTGGGCACATGCCCGGATTTCCAGAAGCTGGCACTCATc                                             <  8:8516/149‑1 (MQ=255)                                                                                                               ttgtGCGCTCAGTTGCAGGTCGGCAGCAATCTGCACCGCCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTCGAAGTACAGAAACGCGGCGAACAGCGTGGGCACATGCCCGGATTTCCAGAAGCTGGCACTCATCGAacac                                       <  5:196894/148‑1 (MQ=255)                                                                                                                       tCAGTTGCAGGTCGGCAGCAATCTGCACCGGCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTCGAAGTACAGAAACGCGGCGAACAGCGTGGGCACATGCCCGGATTTCCAGAAGCTGGCAc                                                  <  3:127940/129‑1 (MQ=255)                                                                                                                       tCAGTTGCAGGTCGGCAGCAATCTGCACCGGCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTCGAAGTACAGAAACGCGGCGAACAGCGTGGGCACATGCCCGGATTTCCAGAAGCTGGCAc                                                  <  3:133505/129‑1 (MQ=255)                                                                                                                       tCAGTTGCAGGTCGGCAGCAATCTGCACCGGCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTCGAAGTACAGAAACGCGGCGAACAGCGTGGGCACATGCCCGGATTTCCAGAAGCTGGCAc                                                  >  4:133505/1‑129 (MQ=255)                                                                                                                       tCAGTTGCAGGTCGGCAGCAATCTGCACCGGCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTCGAAGTACAGAAACGCGGCGAACAGCGTGGGCACATGCCCGGATTTCCAGAAGCTGGCAc                                                  >  4:127940/1‑129 (MQ=255)                                                                                                                       tCAGTTGCAGGTCGGCAGCAATCTGCACCGCCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTCGAAGTACAGAAACGCGGCGAACAGCGTGGGCACATGCCCGGATTTCCAGAAGCTGGCACTCATCGAACACCTCCCTCgg                              >  8:261010/1‑149 (MQ=255)                                                                                                                       tCAGTTGCAGGTCGGCAGCAATCTGCACCGCCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTCGAAGTACAGAAACGCGGCGAACAGCGTGGGCACATGCCCGGATTTCCAGAAGCTGGCACTCATCGAACACCTCACTCgg                              >  3:201543/1‑149 (MQ=255)                                                                                                                       tCAGTTGCAGGTCGGCAGCAATCTGCACCGCCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTCGAAGTACAGAAACGCGGCGAACAGCGTGGGCACATGCCCGGATTTCCAGAAGCTGGCACTCATCGAACACCTCACTCgg                              >  6:364885/1‑149 (MQ=255)                                                                                                                         aGTTGCAGGTCGGCAGCAATCTGCACCGCCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTCGAAGTACAGAAACGCGGCGAACAGCGTGGGCACATGCCCGGATTTCCAGAAGCTGGCACTCATCGAACACCTCACTCgg                              >  6:356709/1‑147 (MQ=255)                                                                                                                         aGTTGCAGGTCGGCAGCAATCTGCACCGCCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTCGAAGTACAGAAACGCGGCGAACAGCGTGGGCACATGCCCGGATTTCCAGAAGCTGGCACTCATCGAACACCTCACTCGgta                            <  2:310121/149‑1 (MQ=255)                                                                                                                             gCAGGTCGGCAGCAATCTGCACCGCCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTCGAAGTACAGAAACGCGGCGAACAGCGTGGGCACATGCCCGGATTTCCAGAAGCTGGCACTCATCGAACACCTCACTCGgtagta                         >  6:174655/1‑148 (MQ=255)                                                                                                                             gCAGGTCGGCAGCAATCTGCACCGCCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTCGAAGTACAGAAACGCGGCGAACAGCGTGGGCACATGCCCGGATTTCCAGAAGCTGGCACTCATCGAACACCTCACTCGgtagta                         >  2:77744/1‑148 (MQ=255)                                                                                                                               aGGTCGGCAGCAATCTGCACCGCCAACGGGCCCAGCAGGTACCAGACAATGAAGCTGAGGTAGAAGTAGAGAAACGCGGCGAACAGCGTGGGCACATGCCCGGATTTCCAGAAGCTGGCACTCATCGAACACCTCACTCGGTAGTACat                      >  5:311651/1‑148 (MQ=255)                                                                                                                                  tCGGCAGCAATCTGCACCGCCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTCGCAGTACAGAAACGCGGCGAACAGCGTGGGCACAAGCCCGGATTTCCAGAAGCTGGCACTCATCGAACACCACACTCGGTAGTAAAACg                    >  6:291468/1‑148 (MQ=255)                                                                                                                                    gtctgcaCTCTGCACCGCCAACGGGCCCAGCAGGTACCTGACCATGAAGCTGAGGTCGAAGTACAGAAACGCGGCGAACAGCGTGGGCACATGCCCGGATTTcc                                                              <  3:325182/100‑1 (MQ=255)                                                                                                                                    ggCAGCAATCTCCACCGCCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTCGAAGTACAGAAACGCGGCGACCAGCGTGGGCCCAGGCCCGGAtatac                                                              >  4:325182/1‑100 (MQ=255)                                                                                                                                        gcaATCTGCACCGCCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTCGAAGTACAGAAACGCGGCGAACAGCGTGGGCACATGCCCGGATTTCCAGAAGCTGGCACTCATCGAACACCTCACTCGGTAGTACAACGGAACGaa             <  3:62955/149‑1 (MQ=255)                                                                                                                                          aaTCTGCACCGCCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTCGAAGTACAGAAACGCGGCGAACAGCGTGGGCACATGCCCGGATTTCCAGAAGCTGGCACTCATCGAACACCTCACTCGGTAGTACAACGGAACGaaaa           <  5:316833/149‑1 (MQ=255)                                                                                                                                               gCACCGCCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTCGAAGTACAGAAACGCGGCGAACAGCGTGGGCACATGCCCGGATTTCCAGAAGCTGGCACTCATCGAACACCTCACTCGGTAGTACAACGGAACGAAAAAGAcg      >  5:133775/1‑149 (MQ=255)                                                                                                                                                cACCGCCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTCGAAGTACAGAAACGCGGCGAACAGCGTGGGCACATGCCCGGATTTCCAGAAGCTGGCACTCATCGAACACCTCACTCGGTAGTACAACGGAACGAAAAAGAcgc     >  8:42547/1‑149 (MQ=255)                                                                                                                                                  ccGCCAACGGGCCCAGCAGGTACCAGAACATGAAGCTGAGGTCGAAGTACAGAAACGCGGCGAACACCGTGGGCAAATGACCGGATTTCCAGAAGCTGGCACGCATCGAACACCTCACTCGGTAGTACAACTGAACGAAAAAGacgacg   >  2:97693/1‑149 (MQ=255)                                                                                                                                                   cGCCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTCGAAGTACAGAAACGCGGCGAACAGCGTGGGCACATGCCCGGATTTCCAGAAGCTGGCACTTATCGAACACCTCACTCGGTAGTACAACGGAACGAAAAAGAcgccgc  >  3:3046/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                   GGCCGATCAGGCCAGCCGTCTTGGGCGACAGGCGGTCGACCAGAACGCCCATGGCGAAGCGCAGGATTGCCCCGGCCA‑GGATCGGCGTGGCCACCATCAGGCCCCGTTGTTGTGCGCTCAGTTGCAGGTCGGCAGCAATCTGCACCGCCAACGGGCCCAGCAGGTACCAGACCATGAAGCTGAGGTCGAAGTACAGAAACGCGGCGAACAGCGTGGGCACATGCCCGGATTTCCAGAAGCTGGCACTCATCGAACACCTCACTCGGTAGTACAACGGAACGAAAAAGACGCCGC  >  NC_002947/2384868‑2385161

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑