Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,811,676 T→A 7.3% T262S (ACG→TCG)  PP_1615 ← LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,811,6760TA 7.3% 166.9 / 4.3 55T262S (ACG→TCG) PP_1615LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (25/26);  new base A (2/2);  total (27/28)
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.33e-01

CCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCAGCGCACCACTGTTCAGGTGTTCCAGCACGTAGTAATCCGGCAACTGGCACAACCCGATTCCTTGCAGCGCCGCATCCAGCACTGCCTGGCCACTGTTGCAGCGCCAGTT  >  NC_002947/1811530‑1811820
                                                                                                                                                   |                                                                                                                                                
ccGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGCACCGATCCTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTAtt                                                                                                                                               >  1:414972/1‑150 (MQ=255)
        agGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGtt                                                                                                                                       >  2:368812/1‑150 (MQ=255)
         ggCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTg                                                                                                                                      <  1:148410/150‑1 (MQ=255)
         ggCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTg                                                                                                                                      <  2:598586/150‑1 (MQ=255)
          gCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTgg                                                                                                                                     <  2:576281/150‑1 (MQ=255)
          gCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTgg                                                                                                                                     >  2:554931/1‑150 (MQ=255)
          gCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGAATACGGCGGtatg                                                                                                                                     <  1:392387/150‑4 (MQ=255)
          gCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGAATACGGCGGtatg                                                                                                                                     >  2:392387/1‑147 (MQ=255)
             aGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCACCAACCGCCGAATACGGCGGtatgggc                                                                                                                                  <  1:523238/150‑7 (MQ=255)
                                gTCAGCCGAAACGATACTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCa                                                                                                               >  1:460331/1‑150 (MQ=255)
                                 tCAGCCGAAACGATACTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCAg                                                                                                              >  2:217008/1‑150 (MQ=255)
                                     ccGAAACGATACTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCAGCGca                                                                                                          >  1:389217/1‑150 (MQ=255)
                                       gAAACGATACTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTAtt                                                                                                                                               <  2:178085/111‑1 (MQ=255)
                                       gAAACGATACTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTAtt                                                                                                                                               >  1:178085/1‑111 (MQ=255)
                                        aaaCGATACTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCAGCGcacca                                                                                                       >  2:573391/1‑150 (MQ=255)
                                         aaCGATACTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCAGCGcaccac                                                                                                      <  1:127718/150‑1 (MQ=255)
                                            gATACTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCAGCGCACCACTGt                                                                                                   >  2:165154/1‑150 (MQ=255)
                                                cTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTGGTgc                                                                                                                                  >  2:441374/1‑115 (MQ=255)
                                                cTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTGGTgc                                                                                                                                  <  1:441374/115‑1 (MQ=255)
                                                cTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTGGTGCGCCTCCAGCAGCGacac                                                                                                                 >  2:402275/1‑132 (MQ=255)
                                                cTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTGGTGCGCCTCCAGCAGCGacac                                                                                                                 <  1:402275/132‑1 (MQ=255)
                                                cTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCATACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCAGCGCACCACTGTTCAg                                                                                               <  2:155479/150‑1 (MQ=255)
                                                    ggCAACCCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCAGCGCACCACTGTTCAGgtgt                                                                                           <  1:443886/150‑1 (MQ=255)
                                                           ccGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTcggc                                                                                                                                           <  1:338763/95‑1 (MQ=255)
                                                           ccGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTcggc                                                                                                                                           >  2:338763/1‑95 (MQ=255)
                                                             gCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCAGCGCACCACTGTTCAGGTGTTCCAGCACg                                                                                  <  1:10568/150‑1 (MQ=255)
                                                                aaCCCCTCCTTCAGGTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCAGCGCACCACTGTTCAGGTGTTCCAGCACgtag                                                                               <  1:165154/150‑1 (MQ=255)
                                                                       cctTCAGGTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCAGCGCACCACTGTTCAGGTGTTCCAGCACGTAGTAATCCg                                                                        <  1:176568/150‑1 (MQ=255)
                                                                           catgTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCAGCGCACCACTGTTCAGGTGTTCCAGCACGTAGTAATCCGGCaa                                                                    <  2:463089/147‑1 (MQ=255)
                                                                            aGGTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTCGTGCGCCTCCAGCAGCGACACCATCGCACCACTGTTCAGGTGTTCCAGCACGGAGTAACCCGGCAAc                                                                   >  2:22853/1‑150 (MQ=255)
                                                                             ggTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCAGCGCACCACTGTTCAGGTGTTCCAGCACGTAGTAATCCGGCAACt                                                                  >  1:33931/1‑150 (MQ=255)
                                                                             ggTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCAGCGCACCACTGTTCAGGTGTTCCAGCACGTAGTAATCCGGCAACt                                                                  <  2:414972/150‑1 (MQ=255)
                                                                                   ccaccaGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCAGCGcac                                                                                                         >  1:179667/1‑105 (MQ=255)
                                                                                   ccaccaGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCAGCGcac                                                                                                         <  2:179667/105‑1 (MQ=255)
                                                                                    caccaGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCAGCGCACCACTGTTCAGGTGTTCCAGCACGTAGTAATCCGGCAACTGGCACaa                                                           >  2:534923/1‑150 (MQ=255)
                                                                                       caGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCAGCGCACCACTGTTCAGGTGTTCCAGCACGTAGTAATCCGGCAACTGGCACAAccc                                                        <  2:389217/150‑1 (MQ=255)
                                                                                                tACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCAGCGCACCACTGTTCAGGTGTTCCAGCACGTAGTACTCCGGCAACTGGCACAACCCGATTCCTTg                                               >  1:191865/1‑150 (MQ=255)
                                                                                                tACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCAGCGCACCACTGTTCAGGTGTTCCAGCACGTAGTAATCCGGCAACTGGCACAACCCGATTCCTTg                                               <  1:368812/150‑1 (MQ=255)
                                                                                                   cTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCAGCGCACCACTGTTCAg                                                                                               <  2:504133/99‑1 (MQ=255)
                                                                                                   cTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCAGCGCACCACTGTTCAg                                                                                               >  1:504133/1‑99 (MQ=255)
                                                                                                               gTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCAGCGCACCACTGTTCAGGTGTTCCAGCACGTAGTAATCCGGCAACTGGCACAACCCGa                                                      <  2:76104/128‑1 (MQ=255)
                                                                                                               gTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCAGCGCACCACTGTTCAGGTGTTCCAGCACGTAGTAATCCGGCAACTGGCACAACCCGa                                                      >  1:76104/1‑128 (MQ=255)
                                                                                                               gTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCAGCGCACCACTGTTCAGGTGTTCCAGCACGTAGTAATCCGGCAACTGGCACAACCCGATTCCtt                                                >  2:48327/1‑134 (MQ=255)
                                                                                                               gTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCAGCGCACCACTGTTCAGGTGTTCCAGCACGTAGTAATCCGGCAACTGGCACAACCCGATTCCtt                                                <  1:48327/134‑1 (MQ=255)
                                                                                                               gTGCCTCCGATGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTGGTGCGCCTCCAGCTGCGCCGCCAGCGCACAACTGTTCAGGTGTTCCCGCCACTCGTACTCCGGTAACCGGCACAGCCCGATTCCTTGCAGCGCCGCATCCAg                                <  1:351758/150‑1 (MQ=255)
                                                                                                                     ctgctgCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCAGCGCACCACTGTTCAGGTGTTCCAGCACGTAGTAATCCGGCAACTGGCACAACCCGATTCCTTGCAGCGCCGCATCCAGCACTGc                          <  1:523433/150‑1 (MQ=255)
                                                                                                                                aCAGCGCCCATACCGCCGTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCAg                                                                                                              <  1:229843/55‑1 (MQ=255)
                                                                                                                                aCAGCGCCCATACCGCCGTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCAg                                                                                                              >  2:229843/1‑55 (MQ=255)
                                                                                                                                      gcaccaaCCGCCGAATACGGCGGTTGGTGCGCCTCCAGCAGCGACACCAGCGCACCACTGTTCAGGTGTTCCAGCACGTAGTAATCCGGCAACTGGCACAACCCGATTCCTTGCAGCGCCGCATCCAGCACTGCCTGGCCACTGTTGCAg          >  2:363673/7‑150 (MQ=255)
                                                                                                                                      cccATACCGCCGTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCAGCGCACCACTGTTCAGGTGTTCCAGCACGTACTAATCCGGCAACTGGCACAAACCGATTCCTTGAAGCGCCGCATCCAGCACTGCCTGGCCACTGTTGCAgc         >  2:281662/1‑150 (MQ=255)
                                                                                                                                              cgccgTATTCGGCTGTTGGTGCGCCTCCAGCAGCGACACCAGCGCACCACTGTTTAGGTGTTCCAGCACGTAGTAATCCGGCAACTGGCACAACCCGATTCCTTGCAGATCGCCATCCAGCACTGCCTGGCCACTGTTGCAGCGCCAGtt  >  2:155721/1‑150 (MQ=255)
                                                                                                                                              cgccgTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCAGCGCACCACTGTTCAGGTGTTCCAGCACGTAGTAATCCGGCAACTGGCACAACCCGATTCCTTGCAGcgcc                                        >  2:385812/1‑112 (MQ=255)
                                                                                                                                              cgccgTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCAGCGCACCACTGTTCAGGTGTTCCAGCACGTAGTAATCCGGCAACTGGCACAACCCGATTCCTTGCAGcgcc                                        <  1:385812/112‑1 (MQ=255)
                                                                                                                                                 cgTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCAGCGCACCACTGTTCAGGTGTTCCAGCACGTAGTAATCCGGCAACTGGCACaa                                                           <  2:466979/90‑1 (MQ=255)
                                                                                                                                                 cgTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCAGCGCACCACTGTTCAGGTGTTCCAGCACGTAGTAATCCGGCAACTGGCACaa                                                           >  1:466979/1‑90 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CCGCGAAGAGGCCAGTACCAGCAGCACATCCCGTCAGCCGAAACGATACTCCGGCAACCCCGCCAACCCCTCCTTCAGGTAATCCACCAGCCGCCGTACCTTCGGCGACAGGTGCCTCTGCTGCGGATACAGCGCCCATACCGCCGTATTCGGCGGTTGGTGCGCCTCCAGCAGCGACACCAGCGCACCACTGTTCAGGTGTTCCAGCACGTAGTAATCCGGCAACTGGCACAACCCGATTCCTTGCAGCGCCGCATCCAGCACTGCCTGGCCACTGTTGCAGCGCCAGTT  >  NC_002947/1811530‑1811820

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: