Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,428,205 T→C 7.4% K292E (AAG→GAG)  yedA ← phenylalanine/threonine exporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,428,2050TC 7.4% 152.3 / 2.0 54K292E (AAG→GAG) yedAphenylalanine/threonine exporter
Reads supporting (aligned to +/- strand):  ref base T (27/23);  new base C (2/2);  total (29/25)
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.97e-01

CGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCGATGATCACTGCCATCGCCACACACTCTTCAGCGCCAATCTGCTCACCGGCAAAGACGATCCCTAGCAGCACCGCCACCGCCGGGTTGACGTAGGCGTAGCTGGTGGCCG  >  NC_002947/5428062‑5428351
                                                                                                                                               |                                                                                                                                                  
cgcAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTgcgc                                                                                                                                                >  5:200961/1‑148 (MQ=255)
  cagcagCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCAc                                                                                                                                             >  4:127075/1‑149 (MQ=255)
  cagcagCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCAc                                                                                                                                             >  4:6582/1‑149 (MQ=255)
                  aGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGATc                                                                                                                              <  2:142072/148‑1 (MQ=255)
                    gCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGAt                                                                                                                               >  6:336610/1‑145 (MQ=255)
                    gCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGAt                                                                                                                               <  5:336610/145‑1 (MQ=255)
                    gCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGATCAGc                                                                                                                           >  4:182164/1‑149 (MQ=255)
                      cGAGTTTCGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGAGCCTTGCGCCACTGCGGCAGCACGATCAGcac                                                                                                                         >  6:281039/1‑149 (MQ=255)
                          tttGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGGTTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGAt                                                                                                                               >  3:64277/1‑139 (MQ=255)
                          tttGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGGTTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGAt                                                                                                                               <  4:64277/139‑1 (MQ=255)
                             ggCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGATCAGCACCACcgc                                                                                                                   <  5:281039/148‑1 (MQ=255)
                                   aTTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCgatg                                                                                                            <  4:118063/149‑1 (MQ=255)
                                       aTGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCGATGATc                                                                                                         <  6:58861/148‑1 (MQ=255)
                                                ccGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCGATGATCACTGCCATc                                                                                                <  6:200961/148‑1 (MQ=255)
                                                    cAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGAt                                                                                                                               <  2:302777/113‑1 (MQ=255)
                                                    cAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGAt                                                                                                                               >  1:302777/1‑113 (MQ=255)
                                                           ttCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCGATGATCACTGCCATCGCCACACActct                                                                                    >  2:268598/1‑149 (MQ=255)
                                                           ttCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCGATGATCACTGCCATCGCCACACActct                                                                                    >  6:289749/1‑149 (MQ=255)
                                                               gCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGATCAGCaccacc                                                                                                                     <  4:234745/112‑1 (MQ=255)
                                                               gCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGATCAGCaccacc                                                                                                                     >  3:234745/1‑112 (MQ=255)
                                                                  aaGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCGATGATCACTGCCATCGCCACACACTCTTCAgcgc                                                                             <  3:127075/149‑1 (MQ=255)
                                                                  aaGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCGATGATCACTGCCATCGCCACACACTCTTCAgcgc                                                                             <  3:6582/149‑1 (MQ=255)
                                                                  aaGGCGGCAAGTTTACCGCATCGCCCCTCCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCGATGATCACTGCCATCGCCACACACTCTTCAgcgc                                                                             >  1:193707/1‑149 (MQ=255)
                                                                       ggcAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCGATGATCACTGCCATCGCCACACACTCTTCAGCGCCAATc                                                                        <  1:268598/149‑1 (MQ=255)
                                                                       ggcAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCGATGATCACTGCCATCGCCACACACTCTTCAGCGCCAATc                                                                        >  6:64830/1‑149 (MQ=255)
                                                                          aaGTTTACCGCATCGCCCCTTCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGATCAGCACTACCGCGCCGATGATCACTGCCATCGCCACACActc                                                                                     <  1:52418/133‑1 (MQ=255)
                                                                          aaGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCGATGATCACTGCCATCGCCACACActc                                                                                     >  2:52418/1‑133 (MQ=255)
                                                                                 ccGCATCGCCCCTGCCTACTCGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCGATGATCACTGCCATCGCCACACACTCTTCAGCGCCAATCTGCTCACCgg                                                              >  1:128623/1‑149 (MQ=255)
                                                                                    cATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCGATGATCACTGCCa                                                                                                  <  3:318858/110‑1 (MQ=255)
                                                                                    cATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCGATGATCACTGCCa                                                                                                  >  4:318858/1‑110 (MQ=255)
                                                                                         ccccTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCGATGATCACTGCCATCGCCACACACTCTTCAGCGCCAATCTGCTCACCGGCAAAGAc                                                       >  1:97505/1‑148 (MQ=255)
                                                                                          cccTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCGATGATCACTGCCATCGCCACACACTCTTCAGCGCCAATCTGCTCACCGGCAAAGACg                                                      >  4:131329/1‑148 (MQ=255)
                                                                                          cccTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCGATGATCACTGCCATCGCCACACACTCTTCAGCGCCAATCTGCTCACCGGCAAAGACGa                                                     <  2:193707/149‑1 (MQ=255)
                                                                                               ccTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCGATGATCACTGCCATCGCCACACACTCTTCAGCGCCAATCTGCTCACCGGCAAAGACGATccc                                                 <  2:97505/148‑1 (MQ=255)
                                                                                                   cTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCGATGATCACTGCCATCGCCACACACTCTTCAGCGCCAATCTGCTCACCGGCAAAGACGATCCCTag                                              <  5:289749/147‑1 (MQ=255)
                                                                                                      gCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCGATGATCACTGCCATCGCCACACACTCTTCAGCGCCAATCTGCTCACCGGCAAAGACGATCCCTagcag                                           <  2:128623/147‑1 (MQ=255)
                                                                                                             tCACCTTTCAACGGTTGTACAGGTACAGGCGCCTTGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCGATGAACACTGCCATCGCCACACACTCTTCAGCGCCAATCTGATCACCGGCAAAGACGATCCCTAGCAGcaccgccac                                  >  6:323404/1‑149 (MQ=255)
                                                                                                                  tttCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCGATGATCACTGCCATCGCCACACACTCTTCAGCGCCAATCTGCTCACCGGCAAAGACGATCCCTAGCAGCACCGCCAccgccg                             >  2:5606/1‑149 (MQ=255)
                                                                                                                     cAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCGATGATCACTGCCATCGCCacacac                                                                                       >  7:324803/1‑88 (MQ=255)
                                                                                                                     cAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCGATGATCACTGCCATCGCCacacac                                                                                       <  8:324803/88‑1 (MQ=255)
                                                                                                                      aaCGGTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCGATGATCACTGCCATCGCCACACACTCTTCAGCGCCAATCTGCTCACCGGCAAAGACGATCCCTAGCAGCACCGCCACCGCCGGGtt                         >  7:194634/1‑149 (MQ=255)
                                                                                                                         ggTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCGATGATCACTGCCATCGCCACACACTCTTCAGCGCCAATCTGCTCACCGGCAAAGACGATCCCTAGCAGCACCGCCACCGCCGGGTTGAc                      <  3:182164/149‑1 (MQ=255)
                                                                                                                         ggTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCGATGATCACTGCCATCGCCACACACTCTTCAGCGCCAATCTGCTCACCGGCAAAGACGATCCCTAGCAGCACCGCCACCGCCGGGTTGAc                      <  8:194634/149‑1 (MQ=255)
                                                                                                                               aCAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCGATGATCACTGCCATCGCCACACACTCTTCAGCGCCAATCTGCTCACCGGCAAAGACGATCCCTAGCAGCACCGCCACCGCCGGGTTGACGTAGGc                >  3:246249/1‑149 (MQ=255)
                                                                                                                                   gTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCGATGATCACTGCCATCGCCACACACTCTTCAGCGCCAATCTGCTCACCGGCAAAGACGATCCCTAGCAGCACCGCCACCGCCGGGTTGACGTAGGCGTAg            >  2:252061/1‑149 (MQ=255)
                                                                                                                                   gTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCGATGATCACTGCCATCGCCACACACTCTTCAGCGCCAATCTGCTCACCGGCAAAGACGATCCCTAGCAGCACCGCCACCGCCGGGTTGACGTAGGCGTAg            >  4:291141/1‑149 (MQ=255)
                                                                                                                                     tCAGGCGCCTCGCGCCACTGCGGCAGGCCGATCAGCACCACCTCGCCGATGATCACTGCCATCGCCACACACTCTTCAGCGCCAATCTGCTCACCGGCAAAGACGATCCCTAGCAGCACCGCCACCGCCGGGTTGACGTAg                  <  3:269140/141‑1 (MQ=255)
                                                                                                                                     tCAGGCGCCTCGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCGATGATCACTGCCATCGCCACACACTCTTCAGCGCCAATCTGCTCACCGGCAAAGACGATCCCTAGCAGCACCGCCACCGCCGGTTTGACGTAg                  <  3:273799/141‑1 (MQ=255)
                                                                                                                                     tCAGGCGCCTCGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCGATGATCACTGCCATCGCCACACACTCTTCAGCGCCAATCTGCTCACCGGCAAAGACGATCCCTAGCAGCACCGCCACCGCCGGGTTGACGTAg                  >  4:273799/1‑141 (MQ=255)
                                                                                                                                     tCAGGCGCCTCGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCGATGATCACTGCCATCGCCACACACTCTTCAGCGCCAATCTGCTCACCGGCAAAGACGATCCCTAGCAGCACCGCCACCGCCGGGTTGACGTAg                  >  4:269140/1‑141 (MQ=255)
                                                                                                                                       aGGCGCCTTGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCGATGATCACTGCCATCGCCACACACTCTTCAGCGCCAATCTGCTCACCGGCAAAGACGATCCCTAGCAGCACCGCCACCGCCGGGTTGACGTAGGCGTAGCtgg        >  3:163348/1‑149 (MQ=255)
                                                                                                                                         gcgcCTTGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCGATGATCACTGCCATCGCCACACACTCTTCAGCGCCAATCTGCTCACCGGCAAAGACGATCCCTAGCAGCACCGCCACCGCCGGGTTGACGTAGGCGTAGCtggtg      >  4:149298/1‑149 (MQ=255)
                                                                                                                                           gcCTTGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCGATGATCACTGCCATCGCCACACACTCTTCAGCGCCAATCTGCTCACCGGCAAAGACGATCCCTAGCAGCACCGCCACCGCCGGGTTGACGTAGGCGTAGCTGGTGGc    <  3:291141/149‑1 (MQ=255)
                                                                                                                                             cTTGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCGATGATCACTGCCATCGCCACACACTCTTCAGCGCCAATCTGCTCACCGGCAAAGACGATCCCTAGCAGCACCGCCACCGCCGGGTTGACGTAGGCGTAGCTGGTGGCCg  <  4:163348/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
CGCAGCAGCGGTTCGATCAGGCCGAGTTTGGCGAAATTCATGGCAATACCGTCAGGGGGTTCAGCGAAGGCGGCAAGTTTACCGCATCGCCCCTGCCTACTTGCAGATTTCACCTTTCAACGGTTGTACAGGTTCAGGCGCCTTGCGCCACTGCGGCAGGCCGATCAGCACCACCGCGCCGATGATCACTGCCATCGCCACACACTCTTCAGCGCCAATCTGCTCACCGGCAAAGACGATCCCTAGCAGCACCGCCACCGCCGGGTTGACGTAGGCGTAGCTGGTGGCCG  >  NC_002947/5428062‑5428351

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