Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,812,308 T→G 8.0% Q37P (CAG→CCG) ‡ 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,3080TG 8.0% 149.2 / 4.8 50Q37P (CAG→CCG) ‡PP_1615LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (25/21);  new base G (2/2);  total (27/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 = 8.82e-01

GCCCATGGCGCGCAGCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCTGTGTGGCAATAATGTTT  >  NC_002947/1812163‑1812445
                                                                                                                                                   |                                                                                                                                          
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                                                                                                                                     <  6:304248/120‑1 (MQ=255)
                                 cTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCt                                                                                                                                     >  1:69133/1‑120 (MQ=255)
                                 cTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCt                                                                                                                                     <  2:69133/120‑1 (MQ=255)
                                  ctgcaGCATTGCCAACGCTTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGtt                                                                                                        >  1:171945/5‑149 (MQ=255)
                                     aaGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCACTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAg                                                                                                             <  6:256234/140‑1 (MQ=255)
                                     aaGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCACTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAg                                                                                                             >  5:256234/1‑140 (MQ=255)
                                               aaTGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAACCCCCCAGGCGTTCGGCGGCGGCAGt                                                                                           >  4:184363/1‑148 (MQ=255)
                                               aaTGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGt                                                                                           >  6:320218/1‑148 (MQ=255)
                                               aaTGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGt                                                                                           >  6:236936/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                                                                            <  7:280364/141‑1 (MQ=255)
                                                                      ggCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTc                                                                            >  8:280364/1‑141 (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                                             <  6:306627/147‑1 (MQ=255)
                                                                                              ggTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTc                                             >  5:306627/1‑147 (MQ=255)
                                                                                                   tgcGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCgg                                      >  5:299481/1‑149 (MQ=255)
                                                                                                   tgcGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCgg                                      >  1:187355/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)
                                                                                                                                                   |                                                                                                                                          
GCCCATGGCGCGCAGCGCTTCCTCACGGCCATCCTGCAAGCGTTGGCAATGCTGCAGGAAGGTCTGCCCGGCTTCGCTCAGGGTTACCCGGCGGGTGCTGCGATACAGCAGGCGCGTTTGCAGGCGGTCCTCCAGCCGGGCGATCTGCCGGCTGATGTGCGATGAGGAAACCCCCAGGCGTTCGGCGGCGGCAGTGAACTGGCCCGATTCGGCCACGGCGACAAATTCGTCGATGCCTTCCCATCGGCTGCTCATCGATTATCCCTGTGTGGCAATAATGTTT  >  NC_002947/1812163‑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: