Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,479,116 T→G 7.4% V156G (GTC→GGC)  ybdR → Zn‑dependent oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,479,1160TG 7.4% 154.1 / 4.4 54V156G (GTC→GGC) ybdRZn‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (22/28);  new base G (2/2);  total (24/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

TACCGGGCGCGGATCGATCATCAACAAAAAGGGCATACCACCAGGCGCAGCGCTGTTCGGATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCG  >  NC_002947/4478973‑4479259
                                                                                                                                               |                                                                                                                                               
tACCGGGCGCGGATCGATCATCAACAAAAAGGGCATACCACCAGGCGCAGCGCTGTTCGGATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCcga                                                                                                                                            <  2:88879/149‑1 (MQ=255)
tACCGGGCGCGGATCGATCATCAACAAAAAGGGCATACCACCAGGCGCAGCGCTGTTCGGATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCcga                                                                                                                                            <  7:335176/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:207219/149‑1 (MQ=255)
                                                                                   tCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGc                                                         >  6:51709/1‑149 (MQ=255)
                                                                                      cTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGc                                                                               <  3:139551/124‑1 (MQ=255)
                                                                                      cTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGc                                                                               >  4:139551/1‑124 (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:222291/149‑1 (MQ=255)
                                                                                                                        aGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAg                    >  6:123317/1‑149 (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)
                                                                                                                                           aaaGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTgccg  <  3:72534/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                               
TACCGGGCGCGGATCGATCATCAACAAAAAGGGCATACCACCAGGCGCAGCGCTGTTCGGATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCG  >  NC_002947/4478973‑4479259

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