Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,812,302 G→C 11.3% A39G (GCC→GGC) ‡ PP_1615 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,812,3020GC11.3% 150.9 / 10.1 53A39G (GCC→GGC) ‡PP_1615LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (21/26);  new base C (3/3);  total (24/29)
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.67e-01

GCCAGGTCGCCCATGGCGCGCAGCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGA‑TTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCTGTGTGGCAATAATGTTTTGTCT  >  NC_002947/1812155‑1812450
                                                                                                                                                   |                                                                                                                                                     
gCCAGGTCGCCCATGGCGCGCAGCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGc                                                                                                                                                      <  6:303855/149‑1 (MQ=255)
 ccAGGTCGCCCATGGCGCTCAGCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCg                                                                                                                                                     <  1:124102/149‑1 (MQ=255)
      tCGCCCATGGCGCGCAGCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTg                                                                                                                                                >  2:274327/1‑149 (MQ=255)
        gCCCATGGCGCGCAGCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGCAGATcgccc                                                                                                                                              >  8:72709/1‑145 (MQ=255)
         cccATGGCGCGCAGCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCg                                                                                                                                             <  3:221288/149‑1 (MQ=255)
          ccATGGCGCGCAGCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCgg                                                                                                                                            >  3:130190/1‑149 (MQ=255)
             tGGCGCGCAGCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGCAGATCGCCCGGCTg                                                                                                                                         <  4:29592/149‑1 (MQ=255)
               gcgcgcAGCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGAt                                                                                                                                       >  1:1402/1‑149 (MQ=255)
                cgcgcAGCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGCAGATCGCCCGGCTggag                                                                                                                                      <  7:72709/149‑4 (MQ=255)
                   gcAGCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGc                                                                                                                                   <  7:42145/149‑1 (MQ=255)
                      gcgcTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGAt                                                                                                                                >  4:46847/1‑149 (MQ=255)
                         cTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAg                                                                                                                             <  6:222393/149‑1 (MQ=255)
                                aCGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAAccc                                                                                                                      >  4:177001/1‑149 (MQ=255)
                                aCGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAAc                                                                                                                        <  8:61199/147‑1 (MQ=255)
                                  ggCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAAccccc                                                                                                                    <  5:179591/149‑1 (MQ=255)
                                           gCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGCAGATCGCCCGGCTggaggac                                                                                                                                   >  8:188593/1‑119 (MQ=255)
                                           gCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTCCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGCAGATCGCCCGGCTggaggac                                                                                                                                   <  7:188593/125‑7 (MQ=255)
                                                  ttGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTcggcggcgg                                                                                                    <  7:67578/149‑1 (MQ=255)
                                                     gCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCCGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAg                                                                                                 <  4:130190/149‑1 (MQ=255)
                                                       aaTGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTg                                                                                               <  3:46847/149‑1 (MQ=255)
                                                           ctgcAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAAc                                                                                            <  2:292572/148‑1 (MQ=255)
                                                           ctgcAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACt                                                                                           >  1:219251/1‑149 (MQ=255)
                                                           ctgcAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACt                                                                                           >  6:221689/1‑149 (MQ=255)
                                                                                   gCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGA‑TTCGGCCACGGCGACa                                                                   >  3:264293/1‑148 (MQ=255)
                                                                                    cTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGA‑TTCGGCCACGGCGACaaa                                                                 >  7:309422/1‑149 (MQ=255)
                                                                                          gTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGA‑TTCGGCCACGGCGACAAATtcgtc                                                           <  7:323462/149‑1 (MQ=255)
                                                                                             aCCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGaaa                                                                                                                         >  6:130679/1‑85 (MQ=255)
                                                                                             aCCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGaaa                                                                                                                         <  5:130679/85‑1 (MQ=255)
                                                                                             aCCCGGCGGGTGCTGCGATACAGCAGGCGCGGTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCACGCGGTCGGCGGCGGCAGTGAACTGGCCCGATTTCGGCCACGGCGACAAATtcgtcg                                                          <  7:220435/148‑1 (MQ=255)
                                                                                              cccGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGA‑TTCGGCCACGGCGACAAATTCGTCGAt                                                        <  2:219251/148‑1 (MQ=255)
                                                                                                         ctgcGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGA‑TTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCg                                            <  2:1402/149‑1 (MQ=255)
                                                                                                          tgcGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTc                                                                                                            <  2:299748/85‑1 (MQ=255)
                                                                                                          tgcGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTc                                                                                                            >  1:299748/1‑85 (MQ=255)
                                                                                                          tgcGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGA‑TTCGGCCACGGCGACAAATTCGTCAATGCCTTCCCATCgg                                           >  7:135504/1‑149 (MQ=255)
                                                                                                            cGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCTCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGA‑tt                                                                                 >  7:313356/1‑109 (MQ=255)
                                                                                                            cGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCTCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGA‑tt                                                                                 <  8:313356/109‑1 (MQ=255)
                                                                                                              aTACAGCAGGCGCGTTTGCAGTCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGA‑TTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGgctgc                                       >  4:44792/1‑149 (MQ=255)
                                                                                                                aCAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGA‑TTCGGCCAc                                                                          <  4:285937/112‑1 (MQ=255)
                                                                                                                aCAGCAGGCGCGGTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGA‑TTCGGCCAc                                                                          >  3:285937/1‑112 (MQ=255)
                                                                                                                  agcagGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGA‑TTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCAt                                   <  8:309422/149‑1 (MQ=255)
                                                                                                                   gcagGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGA‑TTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATc                                  <  4:264293/149‑1 (MQ=255)
                                                                                                                     agGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGA‑TTCGGTCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGa                                <  3:177001/149‑1 (MQ=255)
                                                                                                                     agGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGA‑TTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGa                                >  7:126015/1‑149 (MQ=255)
                                                                                                                         gcgTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGA‑TTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTAt                            >  7:333488/1‑149 (MQ=255)
                                                                                                                            tttGCAGTCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGA‑TTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATccc                         <  3:44792/149‑1 (MQ=255)
                                                                                                                               gCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGA‑TTCGGCCACGGCGACAAATTCGTCGATGCCTTGCCATCGGCTGCTCAt                                   >  3:170231/1‑136 (MQ=255)
                                                                                                                               gCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGA‑TTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCAt                                   <  4:170231/136‑1 (MQ=255)
                                                                                                                                 aGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGA‑TTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCtgtgt                    >  3:249691/1‑149 (MQ=255)
                                                                                                                                  ggCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGA‑GTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCtgtgtg                   >  2:284387/1‑149 (MQ=255)
                                                                                                                                     ggTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGA‑TTCGGCCACGGCGACAAATTCGTCGATGCCTGCCCATCGGCTGCTCATCGATTATCCCTGTGTGGCa                <  1:274327/149‑1 (MQ=255)
                                                                                                                                          atcAGCCGGCAGATCGCCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGA‑TTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCTGTGTGGCaataat           >  4:171705/3‑149 (MQ=255)
                                                                                                                                           ccAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGA‑TTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCTGTGTGGCAATAATg          >  6:237826/1‑149 (MQ=255)
                                                                                                                                              gCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGA‑TTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCTGTGTGGCAATAATGttt       <  5:221689/149‑1 (MQ=255)
                                                                                                                                                 gggCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGA‑TTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCTGTGTGGCAATAATGTTTTGt    <  8:126015/149‑1 (MQ=255)
                                                                                                                                                 gggCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGA‑TTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCTGTGTGGCAATAATGTTTTGt    <  8:333488/149‑1 (MQ=255)
                                                                                                                                                   gCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGA‑TTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCTGTGTGGCAATAATGTTTTGTCt  >  7:26207/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCCAGGTCGCCCATGGCGCGCAGCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGA‑TTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCTGTGTGGCAATAATGTTTTGTCT  >  NC_002947/1812155‑1812450

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