Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,479,177 G→C 9.1% W176C (TGG→TGC)  ybdR → Zn‑dependent oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,479,1770GC 9.1% 154.6 / 5.0 53W176C (TGG→TGC) ybdRZn‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (29/18);  new base C (3/2);  total (32/21)
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.17e-01

CATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGCACACCGTCTTCATGGTGGATGACAACGACTACCGG  >  NC_002947/4479040‑4479318
                                                                                                                                         |                                                                                                                                             
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                                                                                                                                          >  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                                                                                                                                       >  1:185015/1‑95 (MQ=255)
                                                   aaaGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGt                                                                                                                                       <  2:185015/95‑1 (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)
                                                                        aaaGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTgccg                                                             <  3:72534/148‑1 (MQ=255)
                                                                                      ccTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGc                                              >  5:291234/1‑149 (MQ=255)
                                                                                           cgGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGgc                                         <  5:123317/149‑1 (MQ=255)
                                                                                             gATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGgcgc                                       <  6:291234/149‑1 (MQ=255)
                                                                                                  cAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGgctc                                       >  2:158817/1‑142 (MQ=255)
                                                                                                  cAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGcaca                                  >  2:235302/1‑149 (MQ=255)
                                                                                                  cAAAGTGATGTTCCTTTGCGACTTCCTTCGCTCCGCTTGTCTAGCGGTCATCAATGCGGAGTTTGGCGAAGGCTCATCGGTAGGGATCTATGTCGCGGGGCCAGTGGTTCTGTTGAGTGCCGCCTGCGCAGGGATGCTGGgcgc                                       <  1:158817/144‑1 (MQ=255)
                                                                                                           gTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGCACACCGTCTTCa                         >  7:167312/1‑149 (MQ=255)
                                                                                                              ccTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGCACACCGTCTTCAtgg                      >  2:318885/1‑149 (MQ=255)
                                                                                                               cTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGCACACCGTCTTCAtggt                     <  8:167686/149‑1 (MQ=255)
                                                                                                               cTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGCACACCGTCTTCAtggt                     >  6:75344/1‑149 (MQ=255)
                                                                                                               cTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGCACACCGTCTTCAtgg                      <  6:314647/148‑1 (MQ=255)
                                                                                                               cTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGCACACCGTCTTCAtgg                      >  5:314647/1‑148 (MQ=255)
                                                                                                                 aTCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGCACACCGTCTTCAtggtg                    >  6:261247/1‑148 (MQ=255)
                                                                                                                    cGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGAGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGCACACCGTCTTCATGGTGGAt                 >  4:112191/1‑148 (MQ=255)
                                                                                                                     gACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGCACACCGTCTTCATGGTGGATGa               >  4:250965/1‑149 (MQ=255)
                                                                                                                           cTTCCCACCGCGTGGCCAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCGAGTGGGGCTGTTGAGTGCCGCCTGCGCTCGGATGCTGGGCGCGGACAGCGTCGTCAGGGTGGATGACacc           >  2:281055/1‑145 (MQ=255)
                                                                                                                            ttCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGCACACCGTCTTCATGGTGGATGACAACGAc        >  4:145468/1‑149 (MQ=255)
                                                                                                                              cccACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGCACACCGTCTTCATGGTGGATGACAACGACTa      >  2:312188/1‑149 (MQ=255)
                                                                                                                              atgaCCGCTTGCCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATg                                                                                           <  6:372144/61‑1 (MQ=255)
                                                                                                                              atgaCCGCTTGCCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATg                                                                                           >  5:372144/4‑64 (MQ=255)
                                                                                                                              atgaCCGCTTGCCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGcac                                   >  6:139289/4‑120 (MQ=255)
                                                                                                                              atgaCCGCTTGCCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGcac                                   <  5:139289/117‑1 (MQ=255)
                                                                                                                                  ccGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGCACACCGTCTTCATGGTGGATGACAACGACTACCgg  >  3:231044/1‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                             
CATCTGTATGGAGGGATCCCTGGAGGCCAAGCTGATTACGTTAGGGTGCCCAAAGGCAATGTCGGGCCATTCAAAGTCCCGACGACCCTCGCGGATGACAAAGTGCTGTTCCTATCCGACATCCTTCCCACCGCTTGGCAAGCGGTCATCAATGCCGAGATTGGCGAAGGCTCATCGGTAGCGATCTATGGCGCGGGGCCAGTGGGGCTGTTGAGTGCCGCCTGCGCACGGATGCTGGGCGCGCACACCGTCTTCATGGTGGATGACAACGACTACCGG  >  NC_002947/4479040‑4479318

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