Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,029,214	G→C	9.8%	L289F (TTG→TTC)	pcaY →	aromatic acid chemoreceptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,029,214	0	G	C	9.8%	150.6 / 6.7	51	L289F (TTG→TTC)	pcaY	aromatic acid chemoreceptor
Reads supporting (aligned to +/- strand):  ref base G (21/25);  new base C (2/3);  total (23/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 = 7.61e-01

GAGTGTGCCGGTGCAAGGGCAGGGTAACAACGAAATCGGCCAACTGTTCCATGAGCTCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAA‑CCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATC  >  NC_002947/3029073‑3029360                                                                                                                                               |                                                                                                                                                    gAGTGTGCCGGTGCAAGGGCAGGGTAACAACGAAATCGGCCAACTGTTCCATGAGCTCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGCCg                                                                                                                                              >  7:271247/1‑149 (MQ=255)  aGTGTGCCGGTGCAAGGGCAGGGTAACAACGAAATCGGCCAACTGTTCCATGAGCTCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGcc                                                                                                                                               <  5:49315/147‑1 (MQ=255)   gtgtGCCGGTGCAAGGGCAGGGTAACAACGAAATCGGCCAACTGTTCCATGAGCTCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTCAACGCCGtg                                                                                                                                            <  5:314764/149‑2 (MQ=255)      tgCCGGTGCAAGGGCAGGGTAACAACGAAATCGGCCAACTGTTCCATGAGCTCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGCCGTCAcc                                                                                                                                         >  5:324959/1‑149 (MQ=255)      tgCCGGTGCAAGGGCAGGGTAACAACGAAATCGGCCAACTGTTCCATGAGCTCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGCCGTCAc                                                                                                                                          <  8:232068/148‑1 (MQ=255)       gCCGGTGCAAGGGCAGGGTAACAACGAAATCGGCCAACTGTTCCATGAGCTCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTg                                                                                                                                                     >  1:59363/1‑136 (MQ=255)       gCCGGTGCAAGGGCAGGGTAACAACGAAATCGGCCAACTGTTCCATGAGCTCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTg                                                                                                                                                     <  2:59363/136‑1 (MQ=255)        ccGGTGCAAGGGCAGGGTAACAACGAAATCGGCCAACTGTTCCATGAGCTCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGCCggcccc                                                                                                                                         >  2:278443/1‑142 (MQ=255)        ccGGTGCAAGGGCAGGGTAACAACGAAATCGGCCAACTGTTCCATGAGCTCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTCAACGCCGtgg                                                                                                                                           <  1:128522/145‑3 (MQ=255)        ccGGTGCAAGGGCAGGGTAACAACGAAATCGGCCAACTGTTCCATGAGCTCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTCAACGCCGtgg                                                                                                                                           >  2:128522/1‑143 (MQ=255)         cGGTGCAAGGGCAGGGTAACAACGAAATCGGCCAACTGTTCCATGAGCTCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGCCGTCACCGAg                                                                                                                                      <  8:271247/149‑1 (MQ=255)           gTGCAAGGGCAGGGTAACAACGAAATCGGCCAACTGTTCCATGAGCTCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGCCGTCACCGAgg                                                                                                                                     >  4:144818/1‑148 (MQ=255)                  ggCAGGGTAACAACGAAATCGGCCAACTGTTCCATGAGCTCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGCCGTCACCGAGGAAAgcgc                                                                                                                              >  1:246661/1‑148 (MQ=255)                     aGGGTAACAACGAAATCGGCCAACTGTTCCATGAGCTCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCcaac                                                                                                                          <  7:37217/149‑1 (MQ=255)                        gTAACAACGAAATCGGCCAACTGTTCCATGAGCTCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCaacaa‑c                                                                                                                      <  2:28099/149‑1 (MQ=255)                        gTAACAACGAAATCGGCCAACTGTTCCATGAGCTCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTCAACGCCGtggcgg                                                                                                                                        <  6:4954/132‑6 (MQ=255)                        gTAACAACGAAATCGGCCAACTGTTCCATGAGCTCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTCAACGCCGtggcgg                                                                                                                                        >  5:4954/1‑127 (MQ=255)                           acaacGAAATCGGCCAACTGTTCCATGAGCTCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAA‑CCt                                                                                                                    <  3:144818/148‑1 (MQ=255)                             aacGAAATCGGCCAACTGTTCAATGAGCTCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAA‑CCTGCg                                                                                                                 >  5:87032/1‑149 (MQ=255)                                       gCCAACTGTTCCATGAGCTCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAA‑CCTGCGTCAGCAgggg                                                                                                       <  6:324959/149‑1 (MQ=255)                                             tgctgcATGCGCTGGAGCTCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAA‑CCTGc                                                                                                                  <  1:76455/127‑1 (MQ=255)                                             tgctgcATGCGCTGGAGCTCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAA‑CCTGc                                                                                                                  >  2:76455/6‑132 (MQ=255)                                               ttCCATGAGCTCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCGACCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAA‑CCTGCGTCAGCAGGGGCAGGa                                                                                                  <  1:278443/147‑1 (MQ=255)                                                 ccATGAGCTCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAA‑CCt                                                                                                                    <  1:185384/126‑1 (MQ=255)                                                 ccATGAGCTCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAA‑CCt                                                                                                                    >  2:185384/1‑126 (MQ=255)                                                  cATGAGCTCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAA‑CCTGCGTCAGCAGGGGCAGGAGCTg                                                                                              >  3:240529/1‑147 (MQ=255)                                                  cATGAGCTCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAA‑CCTGCGTCAGCAGGGGCAGGAGCTGGa                                                                                            <  2:209884/149‑1 (MQ=255)                                                          ccAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAA‑CCt                                                                                                                    >  6:298024/1‑117 (MQ=255)                                                          ccAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAA‑CCt                                                                                                                    <  5:298024/117‑1 (MQ=255)                                                           cAGCTCATGCAGCATATCCAGCTCGAAACTCTGGTGCAGATCAACAACTGTGCACGGCAGCTGTACGACG‑ACGACACGTCTTTGAACGCCTTCACCGAGGAAAGCGCCAACAA‑CCTGCGTCAGCAGGGGCAGGAGCTGTAGCAgta                                                                                      <  4:240529/146‑3 (MQ=255)                                                             gcgcATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAA‑CCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCa                                                                                 <  5:150957/149‑1 (MQ=255)                                                                 aTGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAA‑CCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGc                                                                             >  4:43975/1‑149 (MQ=255)                                                                    cagcagAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCACGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAA‑CCTGCGTCATCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGt                                                                          <  6:87032/149‑1 (MQ=255)                                                                     agcagAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTCGACGCCG‑CCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAA‑CCAGCGTCAGCAGGGGCAGCAGCTGGAGCAGGCCGCCACTGCCGt                                                                          >  2:274747/1‑148 (MQ=255)                                                                              aGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAA‑CCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAAt                                                                 >  8:98237/1‑148 (MQ=255)                                                                                   gACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAA‑CCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCAc                                                           >  4:312959/1‑149 (MQ=255)                                                                                   gACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAA‑CCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCAc                                                           <  2:246661/149‑1 (MQ=255)                                                                                    acacTCTGTGGCAGATCCACAACTGTGCCCTGCAGCTGGACGCCG‑CCGACACGGCGTTGAACGCCGTCACCGAGGTAAGCGCCAACAACCCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGTCCGCCACTGCCGTGACCGAAATGACCAc                                                           <  1:274747/149‑1 (MQ=255)                                                                                          tGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAA‑CCTGCGTcagca                                                                                                           <  4:305211/94‑1 (MQ=255)                                                                                          tGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAA‑CCTGCGTcagca                                                                                                           >  3:305211/1‑94 (MQ=255)                                                                                            gggCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAA‑CCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTgga                                                  <  3:43975/149‑1 (MQ=255)                                                                                             ggCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAA‑CCTGCGTCAGCAGGGGCAGGAGCTGGAGCa                                                                                         >  6:102214/1‑109 (MQ=255)                                                                                             ggCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAA‑CCTGCGTCAGCAGGGGCAGGAGCTGGAGCa                                                                                         <  5:102214/109‑1 (MQ=255)                                                                                                aGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAA‑CCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGaa                                              <  8:87149/149‑1 (MQ=255)                                                                                                             tgCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAA‑CCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGt                                                                          <  2:319378/108‑1 (MQ=255)                                                                                                             tgCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAA‑CCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGt                                                                          >  1:319378/1‑108 (MQ=255)                                                                                                                ccGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAA‑CCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGa                              <  3:251416/149‑1 (MQ=255)                                                                                                                ccGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAA‑CCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGa                              >  4:251416/1‑149 (MQ=255)                                                                                                                                       aCGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAA‑CCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCaa                                   >  2:75515/1‑122 (MQ=255)                                                                                                                                       aCGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAA‑CCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCaa                                   <  1:75515/122‑1 (MQ=255)                                                                                                                                            gTTGAACGCCGTCACCGAGGAAAGCGCCAACAA‑CCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAAt   >  5:359299/1‑149 (MQ=255)                                                                                                                                            gTTGAACGCCGTCACCGAGGAAAGCGCCAACAA‑CCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAAt   >  1:264224/1‑149 (MQ=255)                                                                                                                                             ttGAACGCCGTCACCGAGGAAAGCGCCAACAA‑CCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATc  <  7:98237/149‑1 (MQ=255)                                                                                                                                               gAACGCCGTCACCGAGGAAAGCGCCAACAA‑CCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAac              >  2:290178/1‑135 (MQ=255)                                                                                                                                               gAACGCCGTCACCGAGGAAAGCGCCAACAA‑CCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAac              <  1:290178/135‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                    GAGTGTGCCGGTGCAAGGGCAGGGTAACAACGAAATCGGCCAACTGTTCCATGAGCTCCAGCGCATGCAGCAGAGCCAGCGCGACACTCTGGGGCAGATCAACAACTGTGCCCGGCAGCTGGACGCCG‑CCGCCACGGCGTTGAACGCCGTCACCGAGGAAAGCGCCAACAA‑CCTGCGTCAGCAGGGGCAGGAGCTGGAGCAGGCCGCCACTGCCGTGACCGAAATGACCACGGCAGTGGAGGAAGTTGCGCGCAACGCGATCACCACCTCGCAAACCACCAGTGAATC  >  NC_002947/3029073‑3029360

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