Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,361,789 G→A 7.2% intergenic (+24/‑81) 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,7890GA 7.2% 156.6 / 2.2 55intergenic (+24/‑81)PP_t45/yphAtRNA‑Leu/inner membrane protein
Reads supporting (aligned to +/- strand):  ref base G (26/25);  new base A (2/2);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.21e-01

AGCTGGTCGTTGCGACCCCTTCAAGCGTGCGAGTGTGGTGGAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGTTCGATGGTCAACGCGACATCATTATCCTGATCGCCCGCGTACTGTTGATG  >  NC_002947/4361653‑4361935
                                                                                                                                        |                                                                                                                                                  
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:178524/149‑1 (MQ=255)
                                                                                                              ccAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGTTCGATGGTCAACGCGACATCATTATc                          <  1:317841/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)
                                                                                                                                        gcggcgCCTTTTTGCTGCCTGACAGCTATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGTTCGATGGTCAACGCGACATCATTATCCTGATCGCCCGCGTACTGTTGATg  <  6:198861/147‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                                  
AGCTGGTCGTTGCGACCCCTTCAAGCGTGCGAGTGTGGTGGAATTGGTATACACAGCAGACTTAAAATCTGTCGGCGTGAGCCTTGCGGGTTCGAGTCCCGCCACTCGCACCAAACTTCGTTGAAAAGGCGCCCTCGCGGCGCCTTTTTGCTGCCTGACAGCAATGTCTGGCGCAGCTGCTAAAGTGGAGCTTGTCCCCGTCCACTGGAACGCCGCCATGCGCTACTCCCTGTTCGATGGTCAACGCGACATCATTATCCTGATCGCCCGCGTACTGTTGATG  >  NC_002947/4361653‑4361935

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