Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,361,786 C→G 7.4% intergenic (+21/‑84) PP_t45 → / → yphA tRNA‑Leu/inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,361,7860CG 7.4% 166.8 / 4.2 54intergenic (+21/‑84)PP_t45/yphAtRNA‑Leu/inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (26/24);  new base G (2/2);  total (28/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01

AGCTGGTCGTTGCGACCCCTTCAAGCGTGCGAGTGTGGTGGAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGTTCGATGGTCAACGCGACATCATTATCCTGATCGCCCGCGTACTGTTGA  >  NC_002947/4361653‑4361933
                                                                                                                                     |                                                                                                                                                   
aGCTGGTCGTTGCGACCCCTTCAAGCGTGCGAGTGTGGTGGAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCttttt                                                                                                                                      <  4:175264/149‑1 (MQ=255)
    ggTCGTTGCGACCCCTTCAAGCGTGCGAGTGTGGTGGAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTtgctg                                                                                                                                  >  6:72475/1‑149 (MQ=255)
     gTCGTTGCGACCCCTTCAAGCGTGCGAGTGTGGTGGAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTtgctgc                                                                                                                                 >  5:111683/1‑149 (MQ=255)
       cGTTGCGACCCCTTCAAGCGTGCGAGTGTGGTGGAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCt                                                                                                                               >  2:178524/1‑149 (MQ=255)
         ttGCGACCCCTTCAAGCGTGCGAGTGTGGTGGAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGa                                                                                                                             <  5:41069/149‑1 (MQ=255)
            cGACCCCTTCAAGCGTGCGAGTGTGGTGTAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAg                                                                                                                          <  6:228191/149‑1 (MQ=255)
              aCCCCTTCAAGCGTGCGAGTGTGGTGGAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCGCGAGGGCGCCttttcaa                                                                                                                                    <  6:321173/137‑4 (MQ=255)
              aCCCCTTCAAGCGTGCGAGTGTGGTGGAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCGCGAGGGCGCCttttcaa                                                                                                                                    >  5:321173/1‑134 (MQ=255)
                cccTTCAAGCGTGCGAGTGTGGTGGAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAAt                                                                                                                      >  2:92890/1‑149 (MQ=255)
                  cTTCAAGCGTGCGAGTGTGGTGGAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCGCGAGGGCGCCttttcaa                                                                                                                                    <  6:84427/133‑4 (MQ=255)
                  cTTCAAGCGTGCGAGTGTGGTGGAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCGCGAGGGCGCCttttcaa                                                                                                                                    >  5:84427/1‑130 (MQ=255)
                   ttCAAGCGTGCGAGTGTGGTGGAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTc                                                                                                                   >  7:106199/1‑149 (MQ=255)
                        gCGTGCGAGTGTGGTGGAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGgcg                                                                                                              >  4:4151/1‑149 (MQ=255)
                          gTGCGAGTGTGTTGGAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGgcgc                                                                                                             <  5:249084/148‑1 (MQ=255)
                              gAGTGTGGTGGAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAgct                                                                                                         >  1:131090/1‑148 (MQ=255)
                               aGTGTGGTGGAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAgctgc                                                                                                       <  5:276972/149‑1 (MQ=255)
                               aGTGTGGTGGAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAgctgc                                                                                                       <  4:234591/149‑1 (MQ=255)
                                    ggtggATTTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAg                                                                                                  >  4:15029/1‑149 (MQ=255)
                                     gtggAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATCTCTGCCGCAGCTCCTAAAGt                                                                                                 >  4:162431/1‑149 (MQ=255)
                                      tggATTTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTg                                                                                                <  3:15029/149‑1 (MQ=255)
                                        gAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGa                                                                                              <  5:72475/149‑1 (MQ=255)
                                        gAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGa                                                                                              <  8:106199/149‑1 (MQ=255)
                                          aTTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGc                                                                                            >  2:78720/1‑149 (MQ=255)
                                                ataCACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTATTGCTGTCAGGCAGCaaaaa                                                                                                                    >  6:190093/1‑116 (MQ=255)
                                                ataCACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTATTGCTGTCAGGCAGCaaaaa                                                                                                                    <  5:190093/119‑4 (MQ=255)
                                                  acacAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTc                                                                                       >  3:258211/1‑146 (MQ=255)
                                                  acacAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTc                                                                                       <  4:258211/146‑1 (MQ=255)
                                                   cacaGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCg                                                                                   >  2:317841/1‑149 (MQ=255)
                                                     cagcagACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAg                                                                                                                          >  5:101900/1‑108 (MQ=255)
                                                     cagcagACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAg                                                                                                                          <  6:101900/108‑1 (MQ=255)
                                                      agcagACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTc                                                                                 >  5:71537/1‑148 (MQ=255)
                                                        cagACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCAc                                                                              <  5:302973/149‑1 (MQ=255)
                                                          gACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTg                                                                            <  3:4151/149‑1 (MQ=255)
                                                           aCTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTg                                                                            >  7:245375/1‑148 (MQ=255)
                                                                          gCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAAcgc                                                                      >  1:141741/1‑139 (MQ=255)
                                                                          gCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAAcgc                                                                      <  2:141741/139‑1 (MQ=255)
                                                                                gCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTcc                                                      <  6:111683/149‑1 (MQ=255)
                                                                                     gCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGtt                                                 >  7:262341/1‑149 (MQ=255)
                                                                                      cGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGtt                                                 <  8:245375/148‑1 (MQ=255)
                                                                                             gAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGTTCGATGGt                                          <  6:71537/148‑1 (MQ=255)
                                                                                               gTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTa                                                          >  5:15881/1‑130 (MQ=255)
                                                                                               gTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTa                                                          <  6:15881/130‑1 (MQ=255)
                                                                                                  ccGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGTTCGATGGTCAAcgc                                    >  1:4508/1‑149 (MQ=255)
                                                                                                         tCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGTTCGATGGTCAACGCGAcatca                             >  5:198861/1‑149 (MQ=255)
                                                                                                           gCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGTTCGATGGTCAACGCGAcatca                             <  1:92890/147‑1 (MQ=255)
                                                                                                              ccAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGTTCGATGGTCAACGCGACATCATTATc                        >  3:350528/1‑149 (MQ=255)
                                                                                                              ccAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGTTCGATGGTCAACGCGACATCATTATc                        <  1:317841/149‑1 (MQ=255)
                                                                                                              ccAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGTTCGATGGTCAACGCGACATCATTATc                        <  1:178524/149‑1 (MQ=255)
                                                                                                               cAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGTTCGATGGTCAACGCGACATCATTATcc                       >  2:323701/1‑149 (MQ=255)
                                                                                                                      cGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGTTCGATGGTCAACGCGACATCATTATCCTGATc                  <  1:323701/147‑1 (MQ=255)
                                                                                                                            aaaGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCCTGCGGTACTCCCTGATCGATGGTCAACGCGACATCATTATCCTGATCGCCcgcg           >  3:257715/1‑148 (MQ=255)
                                                                                                                            aaaGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGTTCGATGGTCAACGCGACATCATTATCCTGATCGCCcgcg           >  3:124300/1‑148 (MQ=255)
                                                                                                                                gcgcCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGTTCGATGGTCAACGCGACATCATTATCCTGATCGCCCGCGTACTg      <  2:4508/149‑1 (MQ=255)
                                                                                                                                gcgcCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGTTCGATGGTCAACGCGACATCATTATCCTGATCGCCCGCGTACTg      >  4:172500/1‑149 (MQ=255)
                                                                                                                                    ccTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGTTCGATGGTCAACGCGACATCATTATCCTGATCGCCCGCGTACTGTTGa  <  8:262341/149‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                                   
AGCTGGTCGTTGCGACCCCTTCAAGCGTGCGAGTGTGGTGGAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGTTCGATGGTCAACGCGACATCATTATCCTGATCGCCCGCGTACTGTTGA  >  NC_002947/4361653‑4361933

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