Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,499,480 2 bp→AC intergenic (+57/‑106) trpS → / → zapE tryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,499,4800CA98.0% 147.4 / NA 50intergenic (+57/‑107)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (16/33);  total (16/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_0029471,499,4810AC98.0% 150.3 / NA 50intergenic (+58/‑106)trpS/zapEtryptophan‑‑tRNA ligase/nucleoside triphosphate hydrolase domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (15/34);  total (16/34)
Fisher's exact test for biased strand distribution p-value = 3.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTTTCGGGCACG‑‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTAT  >  NC_002947/1499339‑1499607
                                                                                                                                              ||                                                                                                                                   
cGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGCAACCGGCGCCGCCGTGCATGGCACCGGcgt                                                                                                                               >  8:5123/1‑147 (MQ=255)
  aCGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTg                                                                                                                             <  5:402378/149‑1 (MQ=255)
    gggCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGcc                                                                                                                           >  6:177790/1‑149 (MQ=255)
     ggCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCg                                                                                                                          <  5:177790/149‑1 (MQ=255)
     ggCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCg                                                                                                                          <  7:130337/149‑1 (MQ=255)
       ccGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGgt                                                                                                                        <  4:190120/149‑1 (MQ=255)
         gCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGgtgt                                                                                                                      <  3:115057/149‑1 (MQ=255)
          cAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGgtgt                                                                                                                      >  8:386743/1‑148 (MQ=255)
                gACAACGCTGAAGCCGGGGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTCAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCTAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgg                                                                                                                <  2:344956/148‑6 (MQ=37)
                gACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCg                                                                                                                          >  1:179962/1‑138 (MQ=255)
                gACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCg                                                                                                                          <  2:179962/138‑1 (MQ=255)
                gACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcggg                                                                                                               <  8:282482/149‑7 (MQ=255)
                gACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcggg                                                                                                               <  2:334515/149‑7 (MQ=255)
                gACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGttcgcgg                                                                                                                <  5:210617/148‑6 (MQ=255)
                 aCAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCCCGCAACACGCGCCGCCGGGCATGGCAACCGCGTTTGCGGTGttagcgggc                                                                                                              >  1:244194/1‑142 (MQ=255)
                     cGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCAcg‑‑‑c                                                                                                      <  6:143415/149‑1 (MQ=37)
                        tGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGCTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCg                                                                                                                          <  1:86309/131‑1 (MQ=255)
                        tGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCg                                                                                                                          >  2:86309/1‑130 (MQ=255)
                              cggTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCCGGACTGAGGCTGCCAAACGGGGGGCGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGC‑‑CCGCTCCca                                                                                              <  5:3116/147‑1 (MQ=21)
                               cGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGC‑‑CCGCTCCcac                                                                                             <  5:123918/149‑1 (MQ=25)
                                       ggCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGACCGCGAAACAGGCGCCGCCGTGCATGGCACCGGCGTTGCCGGTGTTAGCGGGCACGC‑‑CCGCTCCCACACACTGAt                                                                                     >  3:347383/1‑149 (MQ=255)
                                               agagCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCCTGCCCGCGCAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGC‑‑CCGCTCCCACAGACTGATTGCCAtt                                                                              >  5:179756/1‑148 (MQ=25)
                                                  gCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCgg                                                                                                                         <  8:199500/105‑1 (MQ=255)
                                                  gCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCgg                                                                                                                         >  7:199500/1‑105 (MQ=255)
                                                  gCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGC‑‑CCGCTCCCACAGACTGATTGCCATTAAGc                                                                          <  7:5123/149‑1 (MQ=21)
                                                    cTGAACGACGCCTTGCCACCACAGCAGGACTGAGGCTGCCAAACGGGGGGAGGGGGCGGGCCCGCGAAACAGGCGCCGCGGTGCAGGGCACGGGCTTTGCGGGTGTTCGCGGGCACGC‑‑CCGCTCCCACAGa                                                                                          <  1:369487/131‑1 (MQ=18)
                                                    cTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGGGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGC‑‑CCGCTCCCACAGACTGATTGCCATTAAGCgg                                                                        <  7:350340/149‑1 (MQ=21)
                                                     ttatttaaGCTTTGGCACCACAGCAGGACTGAGGCTGCCAAAGGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACGGGCTTTGCCGGTGTTCGCGGGCAGGC‑‑CCGCTCCCACAGACTGATTGCCATTAAGc                                                                          <  8:169742/140‑1 (MQ=18)
                                                     tGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCg                                                                                                                          >  2:420352/1‑101 (MQ=255)
                                                     tGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCg                                                                                                                          <  1:420352/101‑1 (MQ=255)
                                                          gAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGGGTTCGCGGGCACGC‑‑CCGCTCCCACAGACTGATTGCCATTAACCGGGCCTAt                                                                  <  3:337228/149‑1 (MQ=21)
                                                                      ccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACCGC‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGa                                                       <  1:315619/149‑1 (MQ=21)
                                                                      ccacAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGC‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGAc                                                      <  1:63180/149‑1 (MQ=21)
                                                                                      gcggCCAAACGGGGGGCGCGGGGGGGCCCGCGAAACTGGCGCCGCGGGGCAGGCCCCCGGCTTTGCGGGGGTGCGGGGGCCGCCCCCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgtt                                       <  1:334998/146‑1 (MQ=255)
                                                                                      gCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTAGCGGGCACGC‑‑CCGCTCCCACAGCCTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTttgttg                                      >  5:412635/1‑149 (MQ=21)
                                                                                          ccAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGC‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGc                                                               <  8:230629/120‑1 (MQ=17)
                                                                                          ccAAACGGGGGGAGCGGGCGTGCCCGCTAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGC‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGcc                                    <  2:244194/147‑1 (MQ=21)
                                                                                              aCGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGC‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTa                              <  1:330003/149‑1 (MQ=21)
                                                                                              aCGGTGGGAGCCGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGC‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCt                               >  3:302048/1‑148 (MQ=21)
                                                                                                       gCGGGCGTGCCCGCTAAACCGGCGCCGCGGTCAATGGCACGGGCTTTGCCGGTGTTCGCGGGAACGC‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCAt                     <  3:33342/149‑1 (MQ=25)
                                                                                                         gggCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGC‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGa                   >  6:6373/1‑149 (MQ=25)
                                                                                                                 cccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGC‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCtt            <  7:131089/149‑1 (MQ=37)
                                                                                                                  ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGC‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAAt                        <  3:133517/136‑1 (MQ=25)
                                                                                                                  ccGCGAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGC‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAAt                        >  4:133517/1‑136 (MQ=25)
                                                                                                                   cgcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGC‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCt             >  7:164368/1‑146 (MQ=37)
                                                                                                                    gcgAAACAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGC‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGATACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGa         <  1:109485/149‑1 (MQ=37)
                                                                                                                         aCAGGCCCCGCGGGGCATGGCACCGGCTTTGCCGGTGTTCGCGGCCACGC‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCt    <  6:241116/149‑1 (MQ=37)
                                                                                                                          cAGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGC‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTa   <  6:179756/149‑1 (MQ=37)
                                                                                                                           aGGCGCCGCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGC‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTAt  >  6:95521/1‑149 (MQ=37)
                                                                                                                           aGGCGACTCGGTGCATGGCACCGGCTTTGCCGGTGTTCGCGGGCACGC‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTAt  <  7:19295/149‑1 (MQ=37)
                                                                                                                                              ||                                                                                                                                   
CGACGGGCCGCAGTTCGACAACGCTGAAGCCCGTGATGCGGCCATTCAGAGCCTGAACGAAGCCTTGGCACCACAGCAGGACTGAGGCTGCCAAACGGTGGGAGCGGGCGTGCCCGCGAAACAGGCGCCGCGGTGCATGGCCACGGCTTTGCCGGTGTTTCGGGCACG‑‑‑CCGCTCCCACAGACTGATTGCCATTAAGCGGGCCTATCGCTACAGTGACGGCCCGTTTTTTGTTGCCTCGCTAACGAATCATGACTCCCTTAGAACGCTAT  >  NC_002947/1499339‑1499607

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: