Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,252,535 G→A 8.9% intergenic (+22/+221) ureG → / ← PP_2851 urease accessory protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,252,5350GA 8.9% 152.0 / 5.8 56intergenic (+22/+221)ureG/PP_2851urease accessory protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (28/23);  new base A (2/3);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.55e-01

CGACACCCAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCT  >  NC_002947/3252393‑3252683
                                                                                                                                              |                                                                                                                                                    
cgACACCCAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGCGGACTCTGTGCGGGCATGTcc                                                                                                                                               <  2:732522/150‑1 (MQ=255)
 gACACCCAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACGTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCg                                                                                                                                              >  2:275501/1‑150 (MQ=255)
  acacCCAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGc                                                                                                                                             >  1:191591/1‑150 (MQ=255)
      ccAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACag                                                                                                                                         >  2:602570/1‑150 (MQ=255)
      ccAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACag                                                                                                                                         >  2:580589/1‑150 (MQ=255)
       cAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGTGCGAGTACTCTGTGCGGACATGCCCGCACaga                                                                                                                                        <  2:1230987/150‑1 (MQ=255)
        aGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACagag                                                                                                                                       >  1:920638/1‑150 (MQ=255)
          cgcATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGAc                                                                                                                                     <  1:1112846/150‑1 (MQ=255)
                      cAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGtggc                                                                                                                         <  1:522955/150‑1 (MQ=255)
                           gCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCa                                                                                                                    <  2:531959/150‑1 (MQ=255)
                               tGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCgg                                                                                                                <  1:78080/150‑1 (MQ=255)
                                 gACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGct                                                                                                              >  1:869121/1‑150 (MQ=255)
                                    cTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTg                                                                                                           >  1:480147/1‑150 (MQ=255)
                                     ttCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGc                                                                                                          >  2:90710/1‑150 (MQ=255)
                                          cAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGtggctggc                                                                                                                     <  2:1029701/134‑1 (MQ=255)
                                          cAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGtggctggc                                                                                                                     >  1:1029701/1‑134 (MQ=255)
                                             ccTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCAGAGCCGgtgcca                                                                                                  <  1:1144726/150‑4 (MQ=255)
                                              ctaaagaatGGCGAGTGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTcg                                                                                                 >  2:20603/4‑150 (MQ=255)
                                               tGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTcgc                                                                                                >  2:807095/1‑150 (MQ=255)
                                               tGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTcgc                                                                                                >  1:336152/1‑150 (MQ=255)
                                                gaagaaGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTcgcg                                                                                               >  1:1060250/1‑150 (MQ=255)
                                                      gggCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTa                                                                                         <  2:869121/150‑1 (MQ=255)
                                                           aGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCcgc                                                                                    <  2:613897/150‑1 (MQ=255)
                                                           aGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTAAACCcgc                                                                                    <  1:90710/150‑1 (MQ=255)
                                                            ggggTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCcgcg                                                                                   >  1:1295815/1‑150 (MQ=255)
                                                              ggTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTAAACCCGCGaa                                                                                 <  2:336152/150‑1 (MQ=255)
                                                               gTTGCAGGCCGTGATCGGTTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAAt                                                                                >  1:1172369/1‑150 (MQ=255)
                                                                      gCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTAAACCCGCGaa                                                                                 >  1:217344/1‑142 (MQ=255)
                                                                      gCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTAAACCCGCGaa                                                                                 <  2:217344/142‑1 (MQ=255)
                                                                            aTCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCg                                                                   >  2:372372/1‑150 (MQ=255)
                                                                               gATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTcgc                                                                                                >  1:21949/1‑118 (MQ=255)
                                                                               gATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTcgc                                                                                                <  2:21949/118‑1 (MQ=255)
                                                                                 ttttATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGACATGCCCGCACaga                                                                                                                                        >  2:882558/1‑76 (MQ=38)
                                                                                 ttttATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGACATGCCCGCACaga                                                                                                                                        <  1:882558/76‑1 (MQ=38)
                                                                                                     tGCTGGGGGTGCGAGGCTCACAGATCGAGGTCTCTGTGCGGGGATGTTCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCGGCGAATACGGTAGCGGCCCCAACAATGCACGGCGGGTGGAAAtt                                          <  2:27706/150‑1 (MQ=39)
                                                                                                        tGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGc                                       <  2:1060250/150‑1 (MQ=255)
                                                                                                         gggggTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGcc                                      >  1:1214699/1‑150 (MQ=255)
                                                                                                          ggggTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGCTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCca                                     <  2:920638/150‑1 (MQ=255)
                                                                                                              tGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCcagcac                                 <  1:275501/150‑2 (MQ=255)
                                                                                                                     cTGACAGATCGAGGACTCTGTGGGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGggccggc                          <  2:796126/150‑1 (MQ=255)
                                                                                                                      tGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGggccggcc                         >  2:1299619/1‑150 (MQ=255)
                                                                                                                            aTCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCgg                                                                                              >  2:148151/1‑75 (MQ=255)
                                                                                                                            aTCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCgg                                                                                              <  1:148151/75‑1 (MQ=255)
                                                                                                                             tCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTcgc                  >  1:873551/1‑150 (MQ=255)
                                                                                                                               gAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCgg                >  1:1339235/1‑150 (MQ=255)
                                                                                                                               gAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCgg                >  1:1376664/1‑150 (MQ=255)
                                                                                                                               gAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCgg                >  2:846163/1‑150 (MQ=255)
                                                                                                                                 ggtctctGTGCGGACATGCCCGCACAGAGTCCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTg                                         <  1:1174885/120‑1 (MQ=37)
                                                                                                                                 ggtctctGTGCGGACATGCCCGCACAGAGTCCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTg                                         >  2:1174885/4‑123 (MQ=37)
                                                                                                                                  gACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCgg                <  1:208777/147‑1 (MQ=255)
                                                                                                                                  gACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCgg                >  2:208777/1‑147 (MQ=255)
                                                                                                                                     tctGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAAcggc                                                     >  1:327631/1‑107 (MQ=255)
                                                                                                                                     tctGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAAcggc                                                     <  2:327631/107‑1 (MQ=255)
                                                                                                                                     tctGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAAc          >  2:1042417/1‑150 (MQ=255)
                                                                                                                                       tgtgCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAAccc        <  1:1299619/150‑1 (MQ=255)
                                                                                                                                         tgCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGc      <  2:1214699/150‑1 (MQ=255)
                                                                                                                                             ggCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCt  <  2:1172369/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
CGACACCCAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCT  >  NC_002947/3252393‑3252683

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