Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,479,111 C→G 7.3% F154L (TTC→TTG)  ybdR → Zn‑dependent oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,479,1110CG 7.3% 173.9 / 4.6 55F154L (TTC→TTG) ybdRZn‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (24/27);  new base G (2/2);  total (26/29)
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.75e-01

ACCACCAATACCGGGCGCGGATCGATCATCAACAAAAAGGGCATACCACCAGGCGCAGCGCTGTTCGGATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGT  >  NC_002947/4478965‑4479255
                                                                                                                                                  |                                                                                                                                                
accaccAATACCGGGCGCGGATCGATCATCAACAAAAAGGGCATACCACCAGGCGCAGCGCTGTTCGGATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCa                                                                                                                                                 >  4:170083/1‑148 (MQ=255)
  caccaATACCGGGCGCGGATCGATCATCAACAAAAAGGGCATACCACCAGGCGCAGCGCTGTTCGGATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAg                                                                                                                                              >  4:307758/1‑149 (MQ=255)
        tACCGGGCGCGGATCGATCATCAACAAAAAGGGCATACCACCAGGCGCAGCGCTGTTCGGATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCcga                                                                                                                                        <  7:335176/149‑1 (MQ=255)
        tACCGGGCGCGGATCGATCATCAACAAAAAGGGCATACCACCAGGCGCAGCGCTGTTCGGATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCcga                                                                                                                                        <  2:88879/149‑1 (MQ=255)
          ccGGGCGCGGATCGATCATCAACAAAAAGGGCATACCACCAGGCGCAGCGCTGTTCGGATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCcgacg                                                                                                                                      <  6:232607/149‑1 (MQ=255)
                  ggATCGATCATCAACAAAAAGGGCATACCACCAGGCGCAGCGCTGTTCGGATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTcgc                                                                                                                              >  4:290124/1‑149 (MQ=255)
                   gatcgatcATCAACAAAAAGGGCATACCACCAGGCGCAGCGCTGTTCGGATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTcgcg                                                                                                                             <  8:212835/149‑1 (MQ=255)
                   gatcgatcATCAACAAAAAGGGCATACCACCAGGCGCAGCGCTGTGCGGATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTcgcg                                                                                                                             <  8:157524/149‑1 (MQ=255)
                        atcatcAACAAAAAGGGCATACCACCAGGCGCAGCGCTGTTCGGATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGa                                                                                                                        >  2:146330/1‑149 (MQ=255)
                        atcatcAACAAAAAGGGCATACCACCAGGCGCAGCGCTGTTCGGATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGACTTTGAATGGCCcgacatt                                                                                                                                    <  6:168313/137‑4 (MQ=255)
                        atcatcAACAAAAAGGGCATACCACCAGGCGCAGCGCTGATCGGATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGACTTTGAATGGCCcgacatt                                                                                                                                    >  5:168313/1‑134 (MQ=255)
                          catcaACAAAAAGGGCATACCACCAGGCGCAGCGCTGTTCGGATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGAGCCATTCAAAGTCCCGACGACCCTCGCGGATGAc                                                                                                                       >  7:167686/1‑148 (MQ=255)
                           atcaACAAAAAGGGCATACCACCAGGCGCATCGCTGTTCGGATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGACTTTGAATGGCCcgac                                                                                                                                       <  4:109492/131‑1 (MQ=255)
                           atcaACAAAAAGGGCATACCACCAGGCGCAGCGCTGTTCGGATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGACTTTGAATGGCCcgac                                                                                                                                       >  3:109492/1‑131 (MQ=255)
                              aacaaAAAGGGCATACCACCAGGCGCAGCGCTGTTCGGATACCGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAgg                                                                                                                  <  6:40104/149‑2 (MQ=255)
                              aacaaAAAGGGCATACCACCAGGCGCAGCGCTGTTCGGATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAg                                                                                                                   <  8:317440/148‑1 (MQ=255)
                                caaAAAGGGCATACCACCAGGCGCAGCGCTGTTCGGATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGc                                                                                                                <  2:42994/149‑1 (MQ=255)
                                   aaaGGGCATACCACCAGGCGCAGCGCTGTTCGGATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGt                                                                                                             <  6:268951/149‑1 (MQ=255)
                                   aaaGGGCATACCACCAGGCGCAGCGCTGTTCGGATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGt                                                                                                             <  5:23117/149‑1 (MQ=255)
                                         cATACCACCAGGCGCAGCGCTGTTCGGATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTAt                                                                                                       >  3:151344/1‑149 (MQ=255)
                                            accaccAGGCGCAGCGCTGTTCGGATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATcc                                                                                                     >  4:197953/1‑148 (MQ=255)
                                             ccaccaGGCGCAGCGCTGTTCGGATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCg                                                                                                    <  1:146330/148‑1 (MQ=255)
                                                          cgcTGTTCGGATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCcgacg                                                                                                                                      >  5:187546/1‑101 (MQ=255)
                                                          cgcTGTTCGGATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCcgacg                                                                                                                                      <  6:187546/101‑1 (MQ=255)
                                                            cTGTTCGGATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGc                                                                                    <  4:151344/149‑1 (MQ=255)
                                                               ttCGGATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTg                                                                                 <  3:307758/149‑1 (MQ=255)
                                                                           cATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGtcat                                                                     <  3:197953/149‑1 (MQ=255)
                                                                                 tATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATg                                                                <  3:290124/148‑1 (MQ=255)
                                                                                   tGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGCCAAGCGGtgg                                                                      >  6:307900/1‑138 (MQ=255)
                                                                                           tCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGc                                                     >  5:207219/1‑149 (MQ=255)
                                                                                           tCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGc                                                     >  6:51709/1‑149 (MQ=255)
                                                                                           tCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGc                                                     <  6:207219/149‑1 (MQ=255)
                                                                                              cTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGc                                                                           >  4:139551/1‑124 (MQ=255)
                                                                                              cTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGc                                                                           <  3:139551/124‑1 (MQ=255)
                                                                                               tGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAg                                                 >  7:6570/1‑149 (MQ=255)
                                                                                                      cAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATc                                          >  1:141148/1‑149 (MQ=255)
                                                                                                       aaGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGCTTGGCGAAGGCTCATCg                                         >  5:110566/1‑149 (MQ=255)
                                                                                                       aaGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCg                                         >  5:222291/1‑149 (MQ=255)
                                                                                                             aTTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGc                                    >  6:157275/1‑148 (MQ=255)
                                                                                                                   ttAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTAt                             <  8:6570/149‑1 (MQ=255)
                                                                                                                      gggTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATgg                           <  3:170083/148‑1 (MQ=255)
                                                                                                                           cccAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCggg                     <  5:157275/149‑1 (MQ=255)
                                                                                                                            ccAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCTGATGACAAAGTGCTGTTACTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATTACGAGATTGGCGAAGGCTCATAGGTAGCGAACTATgtctcgtg                     >  6:52252/1‑141 (MQ=255)
                                                                                                                            ccAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCgggg                    <  5:52252/149‑1 (MQ=255)
                                                                                                                             cAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTa                                      <  7:327345/130‑1 (MQ=255)
                                                                                                                             cAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTa                                      >  8:327345/1‑130 (MQ=255)
                                                                                                                              aaaGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGt                                                                        <  2:185015/95‑1 (MQ=255)
                                                                                                                              aaaGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGt                                                                        >  1:185015/1‑95 (MQ=255)
                                                                                                                               aaGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGc                   >  7:160412/1‑147 (MQ=255)
                                                                                                                                aGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAg                >  6:123317/1‑149 (MQ=255)
                                                                                                                                aGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAg                <  6:222291/149‑1 (MQ=255)
                                                                                                                                    aaTGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCgg                      <  6:238905/139‑1 (MQ=255)
                                                                                                                                    aaTGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCgg                      >  5:238905/1‑139 (MQ=255)
                                                                                                                                            gCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGtt      >  4:72534/1‑147 (MQ=255)
                                                                                                                                              cATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGt  <  8:160412/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
ACCACCAATACCGGGCGCGGATCGATCATCAACAAAAAGGGCATACCACCAGGCGCAGCGCTGTTCGGATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGT  >  NC_002947/4478965‑4479255

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