Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,182,072 A→C 11.6% intergenic (‑130/+24) PP_4560 ← / ← csbD ribonuclease BN/stress response protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,182,0720AC11.6% 127.5 / 9.4 51intergenic (‑130/+24)PP_4560/csbDribonuclease BN/stress response protein
Reads supporting (aligned to +/- strand):  ref base A (29/15);  new base C (4/2);  total (33/18)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.64e-01

GCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCG‑CCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAACACCCTGTTT  >  NC_002947/5181944‑5182218
                                                                                                                                |                                                                                                                                                   
gCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCtc                                                                                                                                 <  5:132347/149‑1 (MQ=255)
         tgcgcctcTTGATCTTAACACAGCAAATTCTATTACCCCGTGCCCGTGCTGTGCTGCTCCCGGCCGTGGTTTAATAATTTTTTTACTAATATGTGCCGGCCTCTTCGCGGGTTCCCCCTATATTTGGCCGGCACCGCCtcaagtcaag                                                                                                                         <  8:197568/146‑1 (MQ=255)
             ccgcTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGcttc                                                                                                                    >  1:137417/1‑149 (MQ=255)
             ccgcTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGcttc                                                                                                                    >  5:195652/1‑149 (MQ=255)
                cTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGcttctt                                                                                                                  >  5:36888/1‑148 (MQ=255)
                            acaGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCtt                                                                                                      >  4:277050/1‑148 (MQ=255)
                                        aaTAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTccc                                                                                            >  6:237766/1‑146 (MQ=255)
                                        aaTAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTccc                                                                                            <  5:237766/146‑1 (MQ=255)
                                           aGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAACAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGAcc                                                                                      >  2:124913/1‑149 (MQ=255)
                                               cGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGACAAGGCATCTTCACCGCATCCTTGACATTTc                                                                                              >  7:136905/1‑137 (MQ=255)
                                                gTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACctgctg                                                                                 <  2:3258/149‑1 (MQ=255)
                                                  gCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTgcg                                                                               >  3:87992/1‑149 (MQ=255)
                                                   cccGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTgcgc                                                                              >  4:101379/1‑149 (MQ=255)
                                                    ccGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGcgcc                                                                             >  3:84098/1‑149 (MQ=255)
                                                      gAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGcgcctc                                                                           >  6:316102/1‑149 (MQ=255)
                                                               tgctCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCG‑CCttt                                                                    >  3:71131/1‑146 (MQ=255)
                                                               tgctCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCG‑CCTTTGAg                                                                 >  5:275838/1‑149 (MQ=255)
                                                               tgctCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCG‑CCTTTGAg                                                                 >  6:155161/1‑149 (MQ=255)
                                                                      ggCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCG‑CCTTTGAGCTCCTgc                                                          >  2:200216/1‑149 (MQ=255)
                                                                         cGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCG‑CCTTTGAGCTCCTGCgcc                                                       >  7:7127/1‑149 (MQ=255)
                                                                            ggTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCACTTGACCTGCTGCGCCTCG‑CCTTTGAGCTCCTGCgccttg                                                    >  5:182929/1‑149 (MQ=255)
                                                                                 aagaaTTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCG‑CCTTTGAGCTCCTGCGCCTTGCCTTc                                               >  7:69257/1‑149 (MQ=255)
                                                                                 aagaaTTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCG‑CCTTTGAGCTCCTGCGCCTTGCCTTc                                               >  8:82742/1‑149 (MQ=255)
                                                                                  agaaTTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCG‑CCTTTGAGCTCCTGCGCCTTGCCTTCg                                              <  1:202730/149‑1 (MQ=255)
                                                                                   gaaTTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCG‑CCTTTGAGCTCCTGCGCCTTGCCTTCgg                                             <  6:275838/149‑1 (MQ=255)
                                                                                    aaTTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCG‑CCTTTGAGCTCCTGCGCCTTGCCTTCGGc                                            >  4:257464/1‑149 (MQ=255)
                                                                                     ttttGCTATCAAATATGTGCTGTCCTCTTGGCGGTTCAACCCGTTATTTGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCG‑CCTTTGAGCTCCTGCGCCTTGCCTTCGGcc                                           <  6:182929/148‑1 (MQ=255)
                                                                                     aTTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCG‑CCTTTGAGCTCCTGCGCCTTGCCTTCGGc                                            <  8:69257/148‑1 (MQ=255)
                                                                                      tttGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCtt                                                                                                      >  2:48702/1‑90 (MQ=255)
                                                                                      tttGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCtt                                                                                                      <  1:48702/90‑1 (MQ=255)
                                                                                          cTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCG‑CCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAg                                      <  5:155161/149‑1 (MQ=255)
                                                                                            tACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCG‑CCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCt                                    >  5:296906/1‑149 (MQ=255)
                                                                                              cGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCG‑CCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCttt                                  >  5:277277/1‑149 (MQ=255)
                                                                                               gAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGcgcctc                                                                           <  3:257585/108‑1 (MQ=255)
                                                                                               gAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGcgcctc                                                                           >  4:257585/1‑108 (MQ=255)
                                                                                                 aTAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCG‑CCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCg                               >  4:51657/1‑149 (MQ=255)
                                                                                                    aGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCG‑CCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTg                            >  2:70894/1‑149 (MQ=255)
                                                                                                        ccGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTc                          <  1:70894/148‑1 (MQ=255)
                                                                                                          ggCCTCTTCTCGGGTGAACCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCG‑CCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGt                       >  6:252127/1‑148 (MQ=255)
                                                                                                          ggCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCG‑CCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTc                      <  1:124913/149‑1 (MQ=255)
                                                                                                              tcttctCGGGTGAACCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCG‑CCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCttt                                  >  6:242804/1‑133 (MQ=255)
                                                                                                              tcttctCGGGTGAACCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCG‑CCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCttt                                  <  5:242804/133‑1 (MQ=255)
                                                                                                               ctTCTCGGGTGAACCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCG‑CCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCtt                 >  6:309033/1‑149 (MQ=255)
                                                                                                                 tctcgGGTGAACCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCG‑CCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGc               >  6:128190/4‑149 (MQ=255)
                                                                                                                  ctcgGGTGAACCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCG‑CCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCg                                         <  2:1691/120‑1 (MQ=255)
                                                                                                                   gcgGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCG‑CCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGcc              >  7:31181/1‑148 (MQ=255)
                                                                                                                        ttCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCG‑CCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAACAcc        <  4:84098/149‑1 (MQ=255)
                                                                                                                           aCCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCG‑CCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAACACCCTg     <  4:87992/149‑1 (MQ=255)
                                                                                                                             ccGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCG‑CCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTc                                >  3:168748/1‑120 (MQ=255)
                                                                                                                             ccGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCG‑CCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTc                                <  4:168748/120‑1 (MQ=255)
                                                                                                                             ccGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCG‑CCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAACAGCCTGtt   >  1:285348/1‑149 (MQ=255)
                                                                                                                               gagaAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCG‑CCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAACACCCTGttt  >  7:10865/1‑148 (MQ=255)
                                                                                                                                |                                                                                                                                                   
GCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCG‑CCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAACACCCTGTTT  >  NC_002947/5181944‑5182218

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