Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,182,113 T→G 14.2% intergenic (‑123/+31) 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,1130TG14.2% 107.9 / 9.1 49intergenic (‑123/+31)PP_4560/csbDribonuclease BN/stress response protein
Reads supporting (aligned to +/- strand):  ref base T (21/20);  new base G (4/3);  total (26/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCGCAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAACA  >  NC_002947/5181969‑5182258
                                                                                                                                                |                                                                                                                                                 
cgcgCAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCAccc                                                                                                                                              >  2:458621/1‑150 (MQ=255)
   gcAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCgag                                                                                                                                           >  1:904270/1‑150 (MQ=255)
     aGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCTCGGGTGAACCCGcgaa                                                                                                                                         <  2:308593/150‑1 (MQ=255)
     aGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGagaa                                                                                                                                         <  2:684970/150‑1 (MQ=255)
     aGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGagaa                                                                                                                                         <  2:418547/150‑1 (MQ=255)
       gTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGGTCACCCGagaaga                                                                                                                                       >  1:401100/1‑150 (MQ=255)
                     tAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAAATTTCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCAGCACCGCCt                                                                                                                         >  1:143174/1‑150 (MQ=255)
                               tttCGCCGCTTGATCGTTACACAGGCAAGTCTTATGGACCGTGCCCGCGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAgg                                                                                                               <  2:401100/150‑1 (MQ=255)
                               aggagacgcTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAgg                                                                                                               <  1:1005707/144‑1 (MQ=255)
                                                tACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGATAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCAc                                                                                                      <  1:977513/142‑1 (MQ=255)
                                                tACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCAc                                                                                                      >  2:977513/1‑142 (MQ=255)
                                                   acaGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGGGGCCGGCACCGcc                                                                                                                          <  2:723112/119‑1 (MQ=255)
                                                   acaGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGGGGCCGGCACCGcc                                                                                                                          >  1:723112/1‑119 (MQ=255)
                                                         aaGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTc                                                                                     >  1:1149945/1‑150 (MQ=255)
                                                           gTCTAATAGACCGTGCCCGCGCTGTGCTGCTCCCGGCTGTGGTTGAAGACTTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTccc                                                                                   <  1:973822/150‑1 (MQ=255)
                                                               aaTAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGGACACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGa                                                                               >  1:380862/1‑150 (MQ=255)
                                                                    aCCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACctgc                                                                          >  1:190284/1‑150 (MQ=255)
                                                                      catGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACctgctg                                                                        >  2:540441/3‑150 (MQ=255)
                                                                             gAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCACCGTGAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGcgcctcg                                                                 >  1:317656/1‑150 (MQ=255)
                                                                                 tgtgCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCtt                                                             <  2:1149945/150‑1 (MQ=255)
                                                                                   tgctgcTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCg                                                                                                                                  <  2:728888/79‑1 (MQ=255)
                                                                                   tgctgcTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCg                                                                                                                                  >  1:728888/1‑79 (MQ=255)
                                                                                    gctggtCCGGGGCGTGGTTGAAGAATTTGCTTACGCATCAGTGCCGGCCTCTTCGCGGGTTCCCCCGGGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCTCATTTCCCTTGACCTGCTGCGCCTCGCCTTTGa                                                          <  2:143174/150‑1 (MQ=255)
                                                                                      tgctCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGACGCCACCGACTCCAGTCCAGGCCTCTTCACCGCATCTTTTACTGGCCCCTTGACCTGCTGCGCCTCGCCTTTGAGc                                                        <  2:740698/150‑1 (MQ=255)
                                                                                        ctCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTc                                                      <  1:593311/150‑1 (MQ=255)
                                                                                         tCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTcc                                                     <  2:904270/150‑1 (MQ=255)
                                                                                           ccGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACAGTTCCCTTGACCTGCTGCGCCTCGTCTGTGAGCTCCTg                                                   >  1:549185/1‑150 (MQ=255)
                                                                                             ggCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTgcg                                                 >  1:792631/1‑150 (MQ=255)
                                                                                              gCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTgcgc                                                >  1:667611/1‑150 (MQ=255)
                                                                                                     ttGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGc                                                        <  2:115379/135‑1 (MQ=255)
                                                                                                     ttGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGc                                                        >  1:115379/1‑135 (MQ=255)
                                                                                                        aagaaTTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTcc                                                                                    >  2:726029/1‑104 (MQ=255)
                                                                                                        aagaaTTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTcc                                                                                    <  1:726029/104‑1 (MQ=255)
                                                                                                                      gAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTc                        <  1:540441/150‑1 (MQ=255)
                                                                                                                      gAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTc                        <  1:211620/150‑1 (MQ=255)
                                                                                                                      gAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTc                        >  2:526681/1‑150 (MQ=255)
                                                                                                                         tAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGTCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGtt                     <  1:458621/150‑1 (MQ=255)
                                                                                                                         tAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGtt                     >  2:348866/1‑150 (MQ=255)
                                                                                                                                cGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGt                      <  2:126760/142‑1 (MQ=255)
                                                                                                                                cGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGt                      >  1:126760/1‑142 (MQ=255)
                                                                                                                                 ggCCTCTTCTCGGGTGAACCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTg                                                           >  1:1101809/1‑104 (MQ=37)
                                                                                                                                 ggCCTCTTCTCGGGTGAACCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTg                                                           <  2:1101809/104‑1 (MQ=37)
                                                                                                                                 ggCCTCTTCTCGGGTGAACCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCa             >  2:1027551/1‑150 (MQ=39)
                                                                                                                                   ccTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCAcc           <  2:190284/150‑1 (MQ=255)
                                                                                                                                    ctctTCTCGGGTGAACCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGc                                    >  1:1207937/1‑124 (MQ=39)
                                                                                                                                    ctctTCTCGGGTGAACCCGCGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGc                                    <  2:1207937/124‑1 (MQ=255)
                                                                                                                                       ttCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCGGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGc       >  1:327172/1‑150 (MQ=255)
                                                                                                                                          tcgGGTGAACCCGCGAAGAGGCCGGCACCGCCTCACGTCCAGGCTTCTTCACCGCATCCTTCACACTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCAGCACCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTCGccca    >  2:313913/2‑148 (MQ=255)
                                                                                                                                            gggTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAaca  >  1:794511/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
CGCGCAGGTCGGGGAAAATCATAGCGTCTCCTTTCGCCGCTTGATCGTTACACAGGCAAGTCTAATAGACCGTGCCCGAGCTGTGCTGCTCCCGGCCGTGGTTGAAGAATTTGCTTACGAATAAGTGCCGGCCTCTTCGCGGGTTCACCCGAGAAGAGGCCGGCACCGCCTCAAGTCAAGGCTTCTTCACCGCATCCTTCACATTTCCCTTGACCTGCTGCGCCTCGCCTTTGAGCTCCTGCGCCTTGCCTTCGGCCCGCAGCTTTTCGTTGTCGGTCACCTTGCCAACA  >  NC_002947/5181969‑5182258

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