Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,812,302 G→C 7.8% 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,3020GC 7.8% 156.3 / 4.8 51A39G (GCC→GGC) ‡PP_1615LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (26/21);  new base C (2/2);  total (28/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

CAGGTCGCCCATGGCGCGCAGCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCTGTGTGGCAATAATGTTT  >  NC_002947/1812157‑1812445
                                                                                                                                                   |                                                                                                                                                
cAGGTCGCCCATGGCGCGCAGCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGa                                                                                                                                               >  2:5430/1‑149 (MQ=255)
      gCCCATGGCGCGCAGCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGc                                                                                                                                          <  1:5430/148‑1 (MQ=255)
        ccATGGCGCGCAGCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCgg                                                                                                                                       <  8:301091/149‑1 (MQ=255)
                        ttCCGCACGGCCATCCTGCAAGCGTTGGCAATGGTGCAGGAAGGGCTGGCCGGCTTCGCTCAGGGGTACCCGGCGGGTGCGGCGAAACAGCAGGCGGGTTTGCAGGCGGTCCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAg                                                                                                                        <  1:238377/149‑1 (MQ=255)
                              aCGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAAc                                                                                                                   <  8:10167/147‑1 (MQ=255)
                                ggCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAAccccc                                                                                                               <  2:145107/149‑1 (MQ=255)
                                   cATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAgg                                                                                                            >  6:20576/1‑149 (MQ=255)
                                      ccTGCAAGCGTTGGCAATGCTGCAGTAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGAGTGTGCTGCGATACATCAGGCGCGTTTGCAGGCGGTAATCCATACGTGAGATCTGACGGATGATGTGAGATGATGAAACCCCCAGGctt                                                                                                         >  4:46191/1‑147 (MQ=255)
                                      ccTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGt                                                                                                         >  8:220996/1‑149 (MQ=255)
                                       cTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGGGATCTGCCGGCt                                                                                                                                     >  5:304248/1‑120 (MQ=255)
                                       cTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCt                                                                                                                                     >  1:69133/1‑120 (MQ=255)
                                       cTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCt                                                                                                                                     <  2:69133/120‑1 (MQ=255)
                                       cTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCt                                                                                                                                     <  6:304248/120‑1 (MQ=255)
                                        ctgcaGCATTGCCAACGCTTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGtt                                                                                                        >  1:171945/5‑149 (MQ=255)
                                           aaGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCACTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAg                                                                                                             >  5:256234/1‑140 (MQ=255)
                                           aaGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCACTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAg                                                                                                             <  6:256234/140‑1 (MQ=255)
                                                     aaTGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAACCCCCCAGGCGTTCGGCGGCGGCAGt                                                                                           >  4:184363/1‑148 (MQ=255)
                                                     aaTGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGt                                                                                           >  6:236936/1‑148 (MQ=255)
                                                     aaTGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGt                                                                                           >  6:320218/1‑148 (MQ=255)
                                                             catgAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGTCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGCTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGAGGCGGCAGTGAACTGGc                                                                                   <  3:46191/146‑1 (MQ=255)
                                                             cAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGAtgtg                                                                                                                               <  6:221542/105‑1 (MQ=255)
                                                             cAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGAtgtg                                                                                                                               >  5:221542/1‑105 (MQ=255)
                                                             cAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGc                                                                                   <  7:220996/149‑1 (MQ=255)
                                                                 aaGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGa                                                                               <  5:20576/149‑1 (MQ=255)
                                                                     tCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCg                                                                           >  8:59483/1‑149 (MQ=255)
                                                                            ggCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTc                                                                            >  8:280364/1‑141 (MQ=255)
                                                                            ggCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTc                                                                            <  7:280364/141‑1 (MQ=255)
                                                                                   cTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACaaa                                                            >  4:19302/1‑149 (MQ=255)
                                                                                       gggTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCTGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATtcg                                                        >  6:279683/1‑149 (MQ=255)
                                                                                              ccGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGc                                                 >  5:132500/1‑149 (MQ=255)
                                                                                              ccGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGc                                                 >  3:251179/1‑149 (MQ=255)
                                                                                                    ggTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTc                                             >  5:306627/1‑147 (MQ=255)
                                                                                                    ggTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTc                                             <  6:306627/147‑1 (MQ=255)
                                                                                                         tgcGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCgg                                      >  1:187355/1‑149 (MQ=255)
                                                                                                         tgcGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCgg                                      >  5:299481/1‑149 (MQ=255)
                                                                                                            gATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGgctg                                   <  4:251179/149‑1 (MQ=255)
                                                                                                             ttACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGgctgc                                  >  4:148840/2‑149 (MQ=255)
                                                                                                                 agcagGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGc                                                                        >  6:293861/1‑107 (MQ=255)
                                                                                                                     ggCGCGTTTGCAGGCGGTCCTCCATCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCg                            <  2:301375/147‑1 (MQ=255)
                                                                                                                     ggCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGa                           <  1:72815/149‑1 (MQ=255)
                                                                                                                     ggCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCg                            >  1:301375/1‑147 (MQ=255)
                                                                                                                         cgTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATc                      <  7:59483/149‑1 (MQ=255)
                                                                                                                          gTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATcc                     <  6:299481/149‑1 (MQ=255)
                                                                                                                              gCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCtgt                 >  8:272731/1‑149 (MQ=255)
                                                                                                                               cAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCtgt                 <  5:236936/148‑1 (MQ=255)
                                                                                                                               cAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCtgt                 <  5:320218/148‑1 (MQ=255)
                                                                                                                                          atcAGCCGGCAGATCGCCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCTGTGTGGCaataat      >  4:285722/3‑149 (MQ=255)
                                                                                                                                            cAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATc                             >  4:181192/1‑124 (MQ=255)
                                                                                                                                            cAGCCGGCAGATCGCCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCTGTGTGGCAATAATGt    >  8:193685/1‑149 (MQ=255)
                                                                                                                                             aGCCGGCAGATCGCCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCTGTGTGGCAATAATGtt   <  3:285722/149‑1 (MQ=255)
                                                                                                                                              gCCGGCAGATCGCCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCTGTGTGGCAATAATGttt  <  7:193685/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CAGGTCGCCCATGGCGCGCAGCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCTGTGTGGCAATAATGTTT  >  NC_002947/1812157‑1812445

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: