Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,697,318 T→G 8.6% intergenic (+28/‑49) hslV → / → hslU ATP‑dependent HslVU protease peptidase subunit/protease HslVU ATPase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,697,3180TG 8.6% 162.8 / 7.7 58intergenic (+28/‑49)hslV/hslUATP‑dependent HslVU protease peptidase subunit/protease HslVU ATPase subunit
Reads supporting (aligned to +/- strand):  ref base T (31/22);  new base G (2/3);  total (33/25)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.38e-01

CAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGC  >  NC_002947/5697176‑5697457
                                                                                                                                               |                                                                                                                                           
cAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAggg                                                                                                                                       >  5:155229/1‑149 (MQ=255)
 aaGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCACCGGGAAgggg                                                                                                                                      >  1:307217/1‑149 (MQ=255)
   gaccgaccTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAgg                                                                                                                                        <  1:64613/145‑1 (MQ=255)
   gaccgaccTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAgg                                                                                                                                        >  2:64613/1‑145 (MQ=255)
   gaccgaccTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCACCGGGAAGGGGcg                                                                                                                                    <  5:341517/149‑2 (MQ=255)
   gaccgaccTCTCGCCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCt                                                                                                                                    <  4:237966/149‑1 (MQ=255)
    accgaccTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCtt                                                                                                                                   >  1:251877/1‑149 (MQ=255)
     ccgaccTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCACCGGGAAGGGGcgcc                                                                                                                                  <  6:241427/149‑4 (MQ=255)
      cgaccTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCACCGGGAAGGGGcgcca                                                                                                                                 <  2:307217/149‑5 (MQ=255)
          ctctCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCAc                                                                                                                             <  6:155229/149‑1 (MQ=255)
           tctcGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACg                                                                                                                            >  8:324893/1‑149 (MQ=255)
               ggCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGa                                                                                                                        <  3:111240/149‑1 (MQ=255)
               ggCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGa                                                                                                                        <  4:64512/149‑1 (MQ=255)
                            gCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGccc                                                                                                             >  5:362407/1‑147 (MQ=255)
                             ccGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAg                                                                                                          >  2:198964/1‑149 (MQ=255)
                                agaCCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGAc                                                                                                       <  4:88945/149‑1 (MQ=255)
                                     gCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTa                                                                                                   >  8:108103/1‑148 (MQ=255)
                                     gCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTa                                                                                                   <  7:108103/148‑1 (MQ=255)
                                         tGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAATGACCGTATTGAt                                                                                              <  2:251877/149‑1 (MQ=255)
                                               tCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTc                                                                                        >  7:265002/1‑149 (MQ=255)
                                                   cGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCt                                                                                                                          <  3:299064/111‑1 (MQ=255)
                                                   cGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCt                                                                                                                          >  4:299064/1‑111 (MQ=255)
                                                   cGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATg                                                                                    <  6:362407/149‑1 (MQ=255)
                                                                  gtTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCcgc                                                                           >  5:202406/1‑143 (MQ=255)
                                                                     cACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTc                                                                  >  3:335090/1‑149 (MQ=255)
                                                                     cACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTc                                                                  >  5:158146/1‑149 (MQ=255)
                                                                      accaaccaCAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAggg                                                                                                                                       <  5:62550/79‑1 (MQ=255)
                                                                      accaaccaCAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAggg                                                                                                                                       >  6:62550/1‑79 (MQ=255)
                                                                      accaaccaCAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTcc                                                                 >  6:224944/1‑149 (MQ=255)
                                                                        caaccaCACACTGACCATCGAGGATCAGGACCTGGCAGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCGCATGGACCGTATTTATCATGTCCTTGTCCCCCAGCGAGTTCGTCCAc                                                               <  8:66168/149‑1 (MQ=255)
                                                                         aaccaCAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCAc                                                               >  8:20613/1‑148 (MQ=255)
                                                                            cacaACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGaa                                                            >  8:223612/1‑148 (MQ=255)
                                                                            cacaACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCAAGACCCCCCGCGAGATCGTCCACGAAc                                                           <  1:198964/149‑1 (MQ=255)
                                                                                 ccTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGa                                                                                                                        >  6:350653/1‑83 (MQ=255)
                                                                                 ccTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGa                                                                                                                        <  5:350653/83‑1 (MQ=255)
                                                                                  cTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGTTCGTGGACGAtctcgcgg                                                     <  7:324893/149‑5 (MQ=255)
                                                                                     aCCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCTCGTTATGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGCCTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGc                                                   >  1:65549/1‑149 (MQ=255)
                                                                                      ccATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCca                                                 <  4:335090/149‑1 (MQ=255)
                                                                                        aTCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCcaca                                               <  7:223612/149‑1 (MQ=255)
                                                                                        aTCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGACATCGTCCACGAACTCACCCGCcac                                                >  4:7013/1‑148 (MQ=255)
                                                                                          cGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCGAcat                                              >  6:126035/1‑148 (MQ=255)
                                                                                            aggagCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCAcatcat                                           <  6:158146/149‑1 (MQ=255)
                                                                                            aggagCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCAcatcat                                           >  2:361662/1‑149 (MQ=255)
                                                                                            aggagCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCAcatcat                                           >  6:89427/1‑149 (MQ=255)
                                                                                            aggagCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCAcatcat                                           >  1:56860/1‑149 (MQ=255)
                                                                                                gCAGGACCTTGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTc                                                         <  8:208745/131‑1 (MQ=255)
                                                                                                gCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTc                                                         >  7:208745/1‑131 (MQ=255)
                                                                                                       cTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGgac                                 >  2:85045/1‑148 (MQ=255)
                                                                                                        tGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGCTCATGTCCATGACCCCCCGCGAGATCGTCCCCGAACTCAACCGCCACATCATCGGCCAGgacgc                               >  8:282972/1‑148 (MQ=255)
                                                                                                            cGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGc                            >  7:257372/1‑148 (MQ=255)
                                                                                                                 gATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAgcg                      >  1:2090/1‑149 (MQ=255)
                                                                                                                 gATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAgcg                      >  5:319365/1‑149 (MQ=255)
                                                                                                                  aTTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGACCTCAACCGCCACATCATCGGCCAGGACGACGCCAAgcgc                     >  4:297923/1‑149 (MQ=255)
                                                                                                                        cT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTg               <  1:361662/149‑1 (MQ=255)
                                                                                                                               tGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCAtt         >  7:262449/1‑149 (MQ=255)
                                                                                                                               aaaagcCCCTTCACCGGGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCAtt         >  1:351163/5‑149 (MQ=255)
                                                                                                                                  cgcCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGcc      <  7:20613/149‑1 (MQ=255)
                                                                                                                                      ccTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTgc  <  8:262449/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                           
CAAGACCGACCTCTCGGCCCGCGAAATCGCCGAGACCGCCCTGAACATCGCCGGCGACATCTGCGTGTTCACCAACCACAACCTGACCATCGAGGAGCAGGACCTGGCCGACTGATTCAGCT‑GTTCTGGCGCCCCTTCCCGGTGAAGGGGCTTTTCCACGCTGACTCACCTGCCCAAGGACCGTATTGATCATGTCCATGACCCCCCGCGAGATCGTCCACGAACTCAACCGCCACATCATCGGCCAGGACGACGCCAAGCGCGCCGTGGCCATTGCCCTGC  >  NC_002947/5697176‑5697457

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: ‑