Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,450,059 C→A 10.0% R3S (AGG→AGT) ‡ PP_1268 ← HtrA suppressor protein SohA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,450,0590CA10.0% 145.9 / 7.1 50R3S (AGG→AGT) ‡PP_1268HtrA suppressor protein SohA
Reads supporting (aligned to +/- strand):  ref base C (27/18);  new base A (2/3);  total (29/21)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.27e-01

GTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGCCCATGA  >  NC_002947/1449911‑1450203
                                                                                                                                                    |                                                                                                                                                 
gTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGattgcattgc                                                                                                                                                   <  2:286136/149‑1 (MQ=255)
             cAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAg                                                                                                                                      >  2:58661/1‑149 (MQ=255)
                ggCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTg                                                                                                                                   <  3:240447/149‑1 (MQ=255)
                 gcgcgcCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGc                                                                                                                                  <  2:212222/149‑1 (MQ=255)
                         cACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCAc                                                                                                                          >  5:388367/1‑149 (MQ=255)
                          aCCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACg                                                                                                                         <  1:129209/149‑1 (MQ=255)
                            cGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGt                                                                                                                       >  1:382373/1‑149 (MQ=255)
                             gTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGt                                                                                                                       >  5:359048/1‑148 (MQ=255)
                                tCGGGAACATTGGATTTATAGCGATCTGTCATAGTTGACGCGACGTAGATGGTTGAAGCTATGAGTTATGTCCTGTGTGTTAGAGCGATACAACTCAAGGTAATTCGTTTGCTTTGCATTGTCATTTGAGTTGCATTACCACgtttat                                                                                                                    <  4:125996/148‑6 (MQ=255)
                                   ggAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGc                                                                                                               <  4:156394/149‑1 (MQ=255)
                                        aGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGt                                                                                                            >  5:223434/1‑147 (MQ=255)
                                          tGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCg                                                                                                              <  3:292275/143‑1 (MQ=255)
                                          tGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCg                                                                                                              >  4:292275/1‑143 (MQ=255)
                                                      gATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCACCTGCATTGACAAGGCAATGCAATCGCATTAcctt                                                                                                                          <  2:15597/120‑3 (MQ=21)
                                                      gATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCACCTGCATTGACAAGGCAATGCAATCGCATTAcctt                                                                                                                          <  2:8379/120‑3 (MQ=21)
                                                       aTCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGc                                                                                                               >  8:224459/1‑129 (MQ=255)
                                                       aTCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCaa                                                                                           >  2:314314/1‑149 (MQ=255)
                                                         cGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAaca                                                                                         >  4:400426/1‑149 (MQ=255)
                                                           gTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGAAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCACACCAACACt                                                                                       >  7:13718/1‑149 (MQ=255)
                                                             cAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTgg                                                                                     >  3:23706/1‑149 (MQ=255)
                                                                cGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCAt                                                                                  <  7:78136/149‑1 (MQ=255)
                                                                  tGGACTCGACGTCGAAGGTTGCAGCCATGAGTAATGTGATGTGGGTGACAGCGATACAACTCATGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATT‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGt                                                                                <  4:23706/149‑1 (MQ=255)
                                                                   ggACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGt                                                                                >  3:232469/1‑148 (MQ=255)
                                                                    gACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTc                                                                               >  2:93639/1‑148 (MQ=255)
                                                                       tCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTAcc                                                                           >  2:165935/1‑149 (MQ=255)
                                                                         gACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTg                                                                                      >  8:102723/1‑136 (MQ=255)
                                                                         gACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTg                                                                                      <  7:102723/136‑1 (MQ=255)
                                                                               gAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATGCAGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCaa                                                                    <  1:314314/149‑1 (MQ=255)
                                                                                aaGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAAt                                                                   >  1:174865/1‑148 (MQ=255)
                                                                                      gCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTc                                                            <  1:93639/149‑1 (MQ=255)
                                                                                                 tAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGa                                                 >  5:220332/1‑149 (MQ=255)
                                                                                                   aTGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGAc                                                <  6:388367/148‑1 (MQ=255)
                                                                                                   aTGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGAc                                                <  6:223434/148‑1 (MQ=255)
                                                                                                   aTGTCCTGAGGGTGACAGCGATACAACTAAAGGAAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGAAG‑AGCGGTCAGGCAGAGCACACCAACACTGGCATGTCTACCCGGTAAATCACCTTAGCTCTGATTGACg                                               >  1:85687/1‑149 (MQ=255)
                                                                                                      tCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACggg                                             <  4:232469/148‑1 (MQ=255)
                                                                                                       ccTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTc                                           >  7:113977/1‑149 (MQ=255)
                                                                                                           tgGGTGTCAGCGATACAACTGATGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGAATTACCACGGTGATG‑AGCGTTCAGTCAGAGCCCACCAACACTTGGATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAg                                        <  2:85687/148‑1 (MQ=255)
                                                                                                           tgGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCa                                         <  4:174561/147‑1 (MQ=255)
                                                                                                             ggTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCaa                                     <  6:220332/149‑1 (MQ=255)
                                                                                                             ggTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGACAATCACCTTCGCGCTGATTGACGGGTCCAGCaa                                     <  1:58661/149‑1 (MQ=255)
                                                                                                                      gATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGccacct                            >  8:264952/1‑145 (MQ=255)
                                                                                                                        tACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAAcc                          >  5:123236/1‑149 (MQ=255)
                                                                                                                              ccgtggTAATGCACCTGCATTGACAAGGCAATGCAATCGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGa                      <  6:171265/143‑1 (MQ=25)
                                                                                                                              ccgtggTAATGCACCTGCATTGACAAGGCAATGCAATCGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGa                      >  5:171265/5‑147 (MQ=25)
                                                                                                                               cgtggTAATGCACCTGCATTGACAAGGCAATGCAATCGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAg                   >  5:363043/4‑149 (MQ=37)
                                                                                                                                aaGGCAATGCAATCGCATTACCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAgg                  >  5:310740/1‑149 (MQ=255)
                                                                                                                                        gCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAggtggcggtg          >  7:122817/1‑149 (MQ=255)
                                                                                                                                        gCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAggtggcggtg          >  3:281601/1‑149 (MQ=255)
                                                                                                                                        gCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAggtggcggtg          >  2:318044/1‑149 (MQ=255)
                                                                                                                                          gATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAggtggcggtggc        >  5:395624/1‑149 (MQ=255)
                                                                                                                                          gATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAggtggcggtggc        >  4:16380/1‑149 (MQ=255)
                                                                                                                                               cattgcCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGCCCat    >  3:158580/1‑148 (MQ=255)
                                                                                                                                                attgcCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGCCCatga  >  3:195330/1‑149 (MQ=255)
                                                                                                                                                attgcCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGCCCatga  <  1:318044/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GTTTTTTCAAGCCCAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGCCCATGA  >  NC_002947/1449911‑1450203

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