Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,479,114 A→T 7.4% K155N (AAA→AAT ybdR → Zn‑dependent oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,479,1140AT 7.4% 157.8 / 3.3 54K155N (AAA→AATybdRZn‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (23/27);  new base T (2/2);  total (25/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.54e-01

CACCAATACCGGGCGCGGATCGATCATCAACAAAAAGGGCATACCACCAGGCGCAGCGCTGTTCGGATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCG  >  NC_002947/4478967‑4479259
                                                                                                                                                   |                                                                                                                                                 
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                                                                                                                                          <  6:187546/101‑1 (MQ=255)
                                                        cgcTGTTCGGATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCcgacg                                                                                                                                          >  5:187546/1‑101 (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:222291/149‑1 (MQ=255)
                                                                                                                              aGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAg                    >  6:123317/1‑149 (MQ=255)
                                                                                                                                  aaTGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCgg                          >  5:238905/1‑139 (MQ=255)
                                                                                                                                  aaTGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCgg                          <  6:238905/139‑1 (MQ=255)
                                                                                                                                          gCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGtt          >  4:72534/1‑147 (MQ=255)
                                                                                                                                            cATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGt      <  8:160412/149‑1 (MQ=255)
                                                                                                                                                 aaaGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTgccg  <  3:72534/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CACCAATACCGGGCGCGGATCGATCATCAACAAAAAGGGCATACCACCAGGCGCAGCGCTGTTCGGATACAGCCATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCG  >  NC_002947/4478967‑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: