Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,783,187 Δ1 bp 8.5% coding (236/927 nt) PP_2436 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,783,1870A. 8.5% 139.1 / 7.6 47coding (236/927 nt)PP_2436LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (22/21);  new base . (2/2);  total (24/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 = 3.63e-01

GACTGCCGGCGCCAGGCGGCGCAGGCGCTGGAACAGCGGCGGCAAGATGACCGCCTCGGTCAGGTCGGTCATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAGGGCATTGGACACTGCCGGCTGGGTAATGC  >  NC_002947/2783045‑2783331
                                                                                                                                              |                                                                                                                                                 
gACTGCCGGCGCCAGGCGGCGCAGGCGCTGGAACAGCGGCGGCAAGATGACCGCCTCGGTCAGGTCGGTCATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAAcagcgc                                                                                                                                               >  5:216252/1‑147 (MQ=255)
gACTGCCGGCGCCAGGCGGCGCAGGCGCTGGAACAGCGGCGGCAAGATGACCGCCTCGGTCAGGTCGGTCATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAAcagcgc                                                                                                                                               >  5:70417/1‑147 (MQ=255)
 aCTGCCGGCGCCAGGCGGCGCAGGCGCTGGAACAGCGGCGGCAAGATGACCGCCTCGGTCAGGTCGGTCATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAAcagcgccag                                                                                                                                            >  4:293787/1‑149 (MQ=255)
   tGCCGGCGCCAGGCGGCGCAGGCGCTGGAACAGCGGCGGCAAGATGACCGCCTCGGTCAGGTCGGTCATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAAcagcgccagcg                                                                                                                                          >  3:82203/1‑149 (MQ=255)
       ggCGCCAGGCGGCGCAGGCGCTGGAACAGCGGCGGCAAGATGACCGCCTCGGTCAGGTCGGTCATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGCGtt                                                                                                                                     <  3:293787/149‑1 (MQ=255)
          gcCAGGCGGCGCAGGCGCTGGAACAGCGGCGGCAAGATGACCGCCTCGGTCAGGTCGGTCATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGACAGCGCGTTgcg                                                                                                                                  >  8:133543/1‑149 (MQ=255)
          gcCAGGCGGCGCAGGCGCTGGAACAGCGGCGGCAAGATGACCGCCTCGGTCAGGTCGGTCATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACGCGCTGGCGCTGTTgcg                                                                                                                                  <  4:45612/149‑1 (MQ=255)
                cggcgCAGGCGCTGGAACAGCGGCGGCAAGATGACCGCCTCGGTCAGGTCGGTCATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTGCGCAACGCGCTGGCGCTGTTGCGgacatc                                                                                                                            <  4:33605/149‑3 (MQ=255)
                      aGGCGCTGGAACAGCGGCGGCAAGATGACCGCCTCGGTCAGGTCGGTCATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCgat                                                                                                                      <  2:217486/149‑1 (MQ=255)
                          gcTGGAACAGCGGCGGCAAGATGACCGCCTCGGTCAGGTCGGTCATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGatgatg                                                                                                                  >  1:159325/1‑149 (MQ=255)
                                            aGATGACCGCCTCGGTCAGGTCGGTCATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGtgtg                                                                                                >  1:39814/1‑149 (MQ=255)
                                             gATGACCGCCTCGGTCAGGTCGGTCATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTgg                                                                                               <  2:39814/149‑1 (MQ=255)
                                              aTGACCGCCTCGGTCAGGTCGGTCATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTggg                                                                                              <  4:297666/149‑1 (MQ=255)
                                                 aCCGCCTCGGTCAGGTCGGTCATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCAc                                                                                           <  2:159325/149‑1 (MQ=255)
                                                        cGGTCAGGTCGGTCATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCt                                                                                    >  1:285370/1‑149 (MQ=255)
                                                          accgaggcggtcATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCtt                                                                                   <  1:328104/141‑1 (MQ=255)
                                                            cAGGTCGGTCATACTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTgcg                                                                                >  5:77178/1‑149 (MQ=255)
                                                                tcggtcATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGt                                                                            <  5:64558/149‑1 (MQ=255)
                                                                 cggtcATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTa                                                                           >  1:13651/1‑149 (MQ=255)
                                                                  ggtcATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTa                                                                           <  6:216252/148‑1 (MQ=255)
                                                                    tcATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACg                                                                         >  2:168696/1‑148 (MQ=255)
                                                                           gATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGaa                                                                   >  8:102677/1‑147 (MQ=255)
                                                                               cGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCgg                                                             <  4:288648/149‑1 (MQ=255)
                                                                               cGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCgg                                                             <  7:173581/149‑1 (MQ=255)
                                                                                ggAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCggg                                                            >  4:164387/1‑149 (MQ=255)
                                                                                  aaGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTc                                                          <  6:77178/149‑1 (MQ=255)
                                                                                    ggTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCAt                                                        >  3:111436/1‑149 (MQ=255)
                                                                                        ttgttgGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTg                                                      >  4:133329/1‑147 (MQ=255)
                                                                                           ttgGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGt                                                 >  1:131876/1‑149 (MQ=255)
                                                                                                   ctgcGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCAc                                                                      >  6:204291/1‑120 (MQ=255)
                                                                                                   ctgcGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCAc                                                                      <  5:204291/120‑1 (MQ=255)
                                                                                                             aaCGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAg                               >  5:251881/1‑149 (MQ=255)
                                                                                                              aCGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAgg                              <  4:82203/149‑1 (MQ=255)
                                                                                                              aCGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAgg                              <  7:102677/149‑1 (MQ=255)
                                                                                                                gTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGgcgc                                    <  7:176773/141‑1 (MQ=255)
                                                                                                                gTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGgcgc                                    >  8:176773/1‑141 (MQ=255)
                                                                                                                 tGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAGGGCa                           >  3:266175/1‑149 (MQ=255)
                                                                                                                    ggCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAGGGCATTg                        <  2:285370/149‑1 (MQ=255)
                                                                                                                      cTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAGGGCATTGGa                      <  2:131876/149‑1 (MQ=255)
                                                                                                                         tCCTGCACCGATGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAGGGCATTGGac                     <  2:13651/147‑1 (MQ=255)
                                                                                                                           cTGCACCGATGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAGGGCATTGGACACTg                 >  4:177740/1‑149 (MQ=255)
                                                                                                                            tGCACCGATGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAGGGCATTGGACACTGc                >  6:222081/1‑149 (MQ=255)
                                                                                                                                 cGATGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAGGGCATTGGACACTGCCGGCt           >  6:143419/1‑149 (MQ=255)
                                                                                                                                  gATGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAGGGCATTGGACACTGCCGGCTg          <  4:266175/149‑1 (MQ=255)
                                                                                                                                   ctgtgCGCAACGCGCTGGCGCTGTTGCGGACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAGGGCATTGGACACTGCCGGCTgg         >  5:301037/2‑149 (MQ=255)
                                                                                                                                   aTGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAGGGCATTGGACACTGCCGGCTgg         >  8:181732/1‑149 (MQ=255)
                                                                                                                                    tgtgCGCAACGCGCTGGCGCTGTTGCGGACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAGGGCATTGGACACTGCCGGCTggg        >  4:248116/1‑149 (MQ=255)
                                                                                                                                      tCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAGGGCATTGGACACTGCCGGCTGGGt       >  6:200063/1‑148 (MQ=255)
                                                                                                                                        cGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAGGGCATTGGACACTGCCGGCTGGGTAAt    >  1:174441/1‑149 (MQ=255)
                                                                                                                                        cGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAGGGCATTGGACACTGCCGGCTGGGTAAt    <  5:267087/149‑1 (MQ=255)
                                                                                                                                          cAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAGGGCATTGGACACTGCCGGCTGGGTAATGc  <  1:168696/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                 
GACTGCCGGCGCCAGGCGGCGCAGGCGCTGGAACAGCGGCGGCAAGATGACCGCCTCGGTCAGGTCGGTCATGCTGATGCGGAAGGTCTTGTTGGCCTGCTGCGGGTTGAACGTGCGGCTTTCCTGCACCGATGTCCGCAACAGCGCCAGCGCGTTGCGCACAGGGCCGATGATGTTCTGCGCCATGGGTGTGGGCACCATGCCTTGCGCGGTACGCACGAACAGCGGGTCATTGAAGGTTTCGCGCAGGCGCGACAGGGCATTGGACACTGCCGGCTGGGTAATGC  >  NC_002947/2783045‑2783331

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: