Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,703,252 C→T 20.0% intergenic (+101/‑66) 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,2520CT20.0% 56.8 / 8.8 30intergenic (+101/‑66)lysS/PP_1497lysine‑‑tRNA ligase/TetR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (15/9);  new base T (3/3);  total (18/12)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.59e-01

GATCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAGGC  >  NC_002947/1703106‑1703392
                                                                                                                                                  |                                                                                                                                            
gATCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCt                                                                                                                                             <  6:113826/148‑1 (MQ=255)
  tCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGCCtttg                                                                                                                                          >  5:165916/1‑146 (MQ=255)
  tCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCCCTTCGCGGGTAAACCCGCTCCTACAAGGCCtttgt                                                                                                                                         <  4:265969/149‑4 (MQ=255)
    cgcgACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAg                                                                                                                                        <  6:244860/149‑1 (MQ=255)
      cgACCTCATCCTGTCCCCGCCCATGCCTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTGCTGGCCCCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGa                                                                                                                                      <  8:304641/149‑1 (MQ=38)
        aCGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTCGTCAGTGCTGGCCCCTTCGCGGGTAAACCCGCTCCTACAAGGCCTTTGTAGGAGc                                                                                                                                    <  4:49545/149‑1 (MQ=37)
         cGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCg                                                                                                                                   >  4:79810/1‑149 (MQ=255)
          gTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGCCTTTGTAGGAGCgg                                                                                                                                  <  6:165916/149‑1 (MQ=255)
           tCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCggg                                                                                                                                 <  3:79810/149‑1 (MQ=255)
                cTGTTACCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTAc                                                                                                                            <  6:105913/149‑1 (MQ=255)
                                gTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCg                                                                                                                                   <  4:179488/126‑1 (MQ=255)
                                gTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCg                                                                                                                                   >  3:179488/1‑126 (MQ=255)
                                   cGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCgg                                                                                                                                  <  2:305169/124‑1 (MQ=255)
                                   cGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCgg                                                                                                                                  >  1:305169/1‑124 (MQ=255)
                                                                                                      gtCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCg                                     >  5:272338/1‑149 (MQ=255)
                                                                                                          atcgcGGCCTCTGCGCGGGAAACCCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTAACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCa                                  <  5:147947/145‑1 (MQ=255)
                                                                                                               ggCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGCGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCCAGAGGTATCACCCGTGACCCCCCCCCAGCCGGCCCAAGGCGCCCCCGGCCTCGccccacct                             >  6:147947/1‑143 (MQ=21)
                                                                                                                        gcgGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGa                     >  2:115725/1‑148 (MQ=255)
                                                                                                                             tAAACCCGCTCCTACAAGGCGTTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGt               >  3:234579/1‑149 (MQ=255)
                                                                                                                             tAAACCCGCTCCTACAAGGCCTTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGt               >  1:14653/1‑149 (MQ=255)
                                                                                                                             tAAACCAGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTAACCGCGAAGGGGCCAGCACTGAAGAAAAACCACCATGATGTATCCCCAGTGACACCCGCCATGCCTGACAAATGCGCCGCCGGCATCGCCACCGCCGGCGCCGACAGCgg               >  3:222843/1‑148 (MQ=38)
                                                                                                                                 cccGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCaa           >  4:286731/1‑149 (MQ=255)
                                                                                                                                  ccGCTCCTACAAGGCCTTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGccgccg                                           >  7:291724/1‑116 (MQ=255)
                                                                                                                                  ccGCTCCTACAAGGCCTTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGccgccg                                           <  8:291724/116‑1 (MQ=255)
                                                                                                                                  ccGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCaa           >  7:202500/1‑148 (MQ=255)
                                                                                                                                    gCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATAc        >  2:110815/1‑149 (MQ=255)
                                                                                                                                      tCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATAcc       >  7:28679/1‑148 (MQ=255)
                                                                                                                                       ccTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATAcc       >  8:33822/1‑147 (MQ=255)
                                                                                                                                        cTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAg    >  1:166388/1‑149 (MQ=255)
                                                                                                                                          aCAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAGGc  >  1:303182/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
GATCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAGGC  >  NC_002947/1703106‑1703392

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: