Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,703,279 G→A 11.1% intergenic (+128/‑39) lysS → / → PP_1497 lysine‑‑tRNA ligase/TetR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,703,2790GA11.1% 93.3 / 3.0 36intergenic (+128/‑39)lysS/PP_1497lysine‑‑tRNA ligase/TetR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (18/14);  new base A (2/2);  total (20/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.41e-01

GAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAGGCCGCAAGACTGCCCGCCAAGGCAGCGAACAGCGTC  >  NC_002947/1703159‑1703426
                                                                                                                         |                                                                                                                                                   
gAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCAt                                                                                                                           >  2:55479/1‑147 (MQ=255)
     aGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTTAAACCGCGCCTAAAAAGGCCGTGTTGGAGCGGGGTTTCCCGCGCAGGGGGCCGGACAGACGAAACCCCCCCCTGAGGtt                                                                                                                     >  5:23948/1‑147 (MQ=255)
      gccgccTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTTCAAAGGGCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATc                                                                                                                   >  7:97343/1‑149 (MQ=255)
            tttCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACGGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTAACCGCGAAGGGGCCAGCACTGACGAAACACCCCCAAGAGGTATACCCCGt                                                                                                             >  1:247252/1‑149 (MQ=255)
            tttCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCa                                                                                                                                             >  7:109781/1‑117 (MQ=255)
                        tttttCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCAt                                                                                                 >  2:85144/1‑149 (MQ=255)
                         ttttCTTGCATTATGTTTTTTTTTTTCTTTACCGCCCGGTGCGGGGGTAAACCGGCTCATACAATGGCTTTGTTGTTGGGGTTTAACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCAt                                                                                                 <  2:247252/148‑1 (MQ=255)
                                cccTGAAGATTTATGTGTTCATTCCCGCCCTCTGCGGGGGTAACCCCGCACCTACAACGGCCTTGTCGGTGCGGTTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGAc                                                                                          <  8:97343/147‑1 (MQ=255)
                                         tttATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGGCAGCACTGACGAAACACCCCCATGAGGGATCGCCCGTGACACCCGCCCTGCCTGCCCAAGGcgccg                                                                                >  3:115736/1‑149 (MQ=37)
                                           tATGTTGTCAGTACCGGCCTCTTCGCGGGTTAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAAAACCCCCAAGCGGTATCCCCCGTGACACCCGCCCTGCCGGGCCCAAGCgccgcc                                                                              >  4:297799/1‑149 (MQ=25)
                                             tgttgtCAGACCCGGCCTCTCCGCGGGAAACCCCGCCCCAACAACGCCCTTGTCGGAGCGGTTTAACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGccgccg                                                                             <  3:297799/148‑1 (MQ=255)
                                                                  tccgGGAAAACCCCCCCCCACCAAGGCCGTTGTCGTGGGGGGTTACCCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGcc                                                         <  6:23948/146‑1 (MQ=255)
                                                                  cgcgGGTAAACCCGCTCCTACAAGGCCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCAt                                                                         <  8:67317/131‑1 (MQ=25)
                                                                  cgcgGGTAAACCCGCTCCTACAAGGCCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCAt                                                                         >  7:67317/1‑131 (MQ=25)
                                                                             ccGCTCCTACAAGGCCTTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCa                                                                                        <  4:317121/105‑1 (MQ=255)
                                                                             ccGCTCCTACAAGGCCTTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCa                                                                                        >  3:317121/1‑105 (MQ=255)
                                                                             ccGCTCCTACAAGGCCTTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCaa                                             >  7:74236/1‑148 (MQ=255)
                                                                             ccGCTCCTACAAGGCCTTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCaa                                             >  3:61518/1‑148 (MQ=255)
                                                                             ccGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCaa                                             >  2:293834/1‑148 (MQ=255)
                                                                              cGCTCCTACAAGGCCTTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATa                                           <  8:74236/149‑1 (MQ=255)
                                                                               gCTCCTACAAGGCCTTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATAc                                          >  6:189657/1‑149 (MQ=255)
                                                                                cTCCTACAAGGCCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGACAAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCaa                                                                                       <  4:255698/103‑1 (MQ=18)
                                                                                cTCCTACAAGGCCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGACAAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCaa                                                                                       >  3:255698/1‑103 (MQ=18)
                                                                                   cTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTg                                                >  4:214607/1‑139 (MQ=255)
                                                                                   cTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTg                                                <  3:214607/139‑1 (MQ=255)
                                                                                   cTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAg                                      <  1:85144/149‑1 (MQ=255)
                                                                                                  tAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAGGCCGCAAGACTGc                         <  1:55479/148‑1 (MQ=255)
                                                                                                       gCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCaa                                       >  2:243913/1‑129 (MQ=255)
                                                                                                       gCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCaa                                       <  1:243913/129‑1 (MQ=255)
                                                                                                       gCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAGGCCGCAAGACTGCCCGcc                    >  3:299860/1‑148 (MQ=255)
                                                                                                        cGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAGGCCGCAAGACTGCCCGCCaa                  <  4:61518/149‑1 (MQ=255)
                                                                                                          ggTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAGGCCGCAAGACTGCCCGCCAAgg                >  2:28115/1‑149 (MQ=255)
                                                                                                                cccGCGATGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAGGCCGCAAGACTGCCCGCCAAGGCAGCg           <  1:293834/148‑1 (MQ=255)
                                                                                                                 ccGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAGGCCGCAAGACTGCCCGCCAAGGCAGCGaa         <  4:115736/149‑1 (MQ=255)
                                                                                                                   gcgAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAGGCCGCAAGACTGCCCGCCAAGGCAGCGaa         >  4:24490/1‑147 (MQ=255)
                                                                                                                   gcgAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAGGCCGCAAGACTGCCCGCCAAGGCAGCGAACa       <  5:189657/149‑1 (MQ=255)
                                                                                                                        ggggCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAGGCCGCAAGACTGCCCGCCAAGGCAGCGAACAGCGTc  <  1:28115/149‑1 (MQ=255)
                                                                                                                        ggggCCAGCACTGACGAAACACACCCATGAGGTATACCCCGTGACACCCGCCATGCCTGACCAACGCGCCGCCGCCATCGCCACCGCCGTCGCCGAAAGCGAGCAAAACCAAGGCCGCAAGACTCCCCGCCAAGCCAGCGAACAGCgac  >  4:339388/1‑147 (MQ=255)
                                                                                                                         |                                                                                                                                                   
GAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAGGCCGCAAGACTGCCCGCCAAGGCAGCGAACAGCGTC  >  NC_002947/1703159‑1703426

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: