Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,450,062 T→A 9.4% T2T (ACA→ACT) ‡ 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,0620TA 9.4% 153.9 / 5.8 53T2T (ACA→ACT) ‡PP_1268HtrA suppressor protein SohA
Reads supporting (aligned to +/- strand):  ref base T (29/19);  new base A (2/3);  total (31/22)
Fisher's exact test for biased strand distribution p-value = 6.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.78e-01

CAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGCCCATGA  >  NC_002947/1449924‑1450203
                                                                                                                                          |                                                                                                                                              
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                                                                                                              >  4:292275/1‑143 (MQ=255)
                             tGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCg                                                                                                              <  3:292275/143‑1 (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)
                                                                                                                                          |                                                                                                                                              
CAAGGCGCGCCGCACCGTTTCGGGAACAGTGGTTTGATAGCGATCGGTCAGCGTGGACTCGACGTCGAAGGTGGCAGCCATGAGTAATGTCCTGTGGGTGACAGCGATACAACTCAAGGTAATGCGATTGCATTGCCTTGTCAATGCAGGTGCATTACCACGGTGATG‑AGCGGTCAGGCAGAGCCCACCAACACTGGCATGTCTACCCGGTCAATCACCTTCGCGCTGATTGACGGGTCCAGCAACCGCCCCAACCGCGACAGGTGGCGGTGGCCCATGA  >  NC_002947/1449924‑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: ‑