Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,703,250 G→C 18.8% intergenic (+99/‑68) 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,2500GC18.8% 65.7 / 12.1 32intergenic (+99/‑68)lysS/PP_1497lysine‑‑tRNA ligase/TetR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (16/10);  new base C (3/3);  total (19/13)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

GTCGATCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAGGC  >  NC_002947/1703103‑1703392
                                                                                                                                                   |                                                                                                                                              
gTCGATCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGc                                                                                                                                               >  5:105913/1‑149 (MQ=255)
gTCGATCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGGGAGCTGAACGATCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTAGAAAGGc                                                                                                                                               <  6:206843/149‑1 (MQ=255)
   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                                                                                                                                   >  3:179488/1‑126 (MQ=255)
                                   gTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCg                                                                                                                                   <  4:179488/126‑1 (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)
                                                                                                                                                   |                                                                                                                                              
GTCGATCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAGGC  >  NC_002947/1703103‑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: