Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,252,497 T→C 19.6% intergenic (+27/+216) 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,4970TC19.6% 93.9 / 19.8 51intergenic (+27/+216)ureG/PP_2851urease accessory protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (19/22);  new base C (5/5);  total (24/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

GACACCCAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACC  >  NC_002947/3252351‑3252633
                                                                                                                                                  |                                                                                                                                        
gACACCCAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTcc                                                                                                                                        >  7:138471/1‑149 (MQ=255)
  cacCCAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAAAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGc                                                                                                                                      <  3:256282/149‑1 (MQ=255)
     ccAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGACATGCCCGcaca                                                                                                                                   <  3:62396/149‑1 (MQ=255)
               cGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAAc                                                                                                                          <  3:337279/148‑1 (MQ=255)
                  ccGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCgg                                                                                                                      <  8:138471/149‑1 (MQ=255)
                             cTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCAcc                                                                                                           >  3:276955/1‑149 (MQ=255)
                              tGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCg                                                                                                          >  6:83477/1‑149 (MQ=255)
                                   cTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGtgg                                                                                                                   <  2:263024/135‑1 (MQ=255)
                                   cTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGtgg                                                                                                                   >  1:263024/1‑135 (MQ=255)
                                    ttCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTg                                                                                                    >  5:231980/1‑149 (MQ=255)
                                        gCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCgg                                                                                                >  4:50421/1‑149 (MQ=255)
                                        gCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGACATGCCCGCACaga                                                                                                                                 <  3:87302/116‑1 (MQ=255)
                                        gCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGACATGCCCGCACaga                                                                                                                                 >  4:87302/1‑116 (MQ=255)
                                         cAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGgt                                                                                               <  3:115769/149‑1 (MQ=255)
                                            ccTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGtt                                                                                            >  5:362095/1‑149 (MQ=255)
                                              tGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGcgca                                                                                                                                   >  6:329215/1‑108 (MQ=255)
                                                  gaagGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCggg                                                                                      <  6:231980/149‑1 (MQ=255)
                                                                    ggCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACgg                                                                     >  4:105799/1‑148 (MQ=255)
                                                                     gCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATGCgg                                                                     >  1:240051/1‑147 (MQ=255)
                                                                       cGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGaa                                                                          <  3:67132/140‑1 (MQ=255)
                                                                       cGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGaa                                                                          >  4:67132/1‑140 (MQ=255)
                                                                         tGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCgg                                                               <  3:50421/149‑1 (MQ=255)
                                                                          gATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGtt                                                                                    >  6:4882/1‑127 (MQ=255)
                                                                          gATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGtt                                                                                    <  5:4882/127‑1 (MQ=255)
                                                                           aTCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGcc                                                             >  7:320745/1‑149 (MQ=255)
                                                                                   tATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGgt                                                                                               >  5:114393/1‑107 (MQ=255)
                                                                                   tATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGgt                                                                                               <  6:114393/107‑1 (MQ=255)
                                                                                   tATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGcaacaa                                                      <  4:114822/148‑1 (MQ=255)
                                                                                    aTCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAAt                                                     <  5:83477/148‑1 (MQ=255)
                                                                                        tcgAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAAcg                                                >  1:205677/1‑149 (MQ=255)
                                                                                                gggATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTg                                         >  8:215293/1‑148 (MQ=255)
                                                                                                          gggTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATTCGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGcga                              <  2:240051/149‑3 (MQ=255)
                                                                                                                   gCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGggccg                     <  7:215293/149‑1 (MQ=255)
                                                                                                                   gCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAAATGGCCAGCAGggccg                     <  4:276955/149‑1 (MQ=255)
                                                                                                                       aCAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAAcggc                                              >  2:221015/1‑120 (MQ=255)
                                                                                                                       aCAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAAcggc                                              <  1:221015/120‑1 (MQ=255)
                                                                                                                       aCAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGggccg                     <  4:150261/145‑1 (MQ=255)
                                                                                                                       aCAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGggccg                     <  4:203216/145‑1 (MQ=255)
                                                                                                                       aCAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGggccg                     >  3:203216/1‑145 (MQ=255)
                                                                                                                       aCAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGggccg                     >  3:150261/1‑145 (MQ=255)
                                                                                                                         aGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCctct               <  3:105799/149‑1 (MQ=255)
                                                                                                                          gATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCtt              <  6:362095/149‑1 (MQ=255)
                                                                                                                               aGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCgg         <  7:293415/149‑1 (MQ=255)
                                                                                                                                 gtctctGTGCGGACATGCCCGCACAGAGTCCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCggg        >  2:216201/3‑148 (MQ=255)
                                                                                                                                   ctctGTGCGGACATGCCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTa      >  8:100185/1‑148 (MQ=255)
                                                                                                                                    tctGTGCGGACATGCCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGggcc                      <  4:259000/131‑1 (MQ=255)
                                                                                                                                    tctGTGCGGACATGCCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGggcc                      >  3:259000/1‑131 (MQ=255)
                                                                                                                                      tgtgCGGACATGCCCGCACAGAGTCCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAAcc  >  3:156590/1‑149 (MQ=255)
                                                                                                                                      tgtgCGGACATGCCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAAt                                                     <  4:43313/98‑1 (MQ=255)
                                                                                                                                        tgCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAAc   <  8:320745/146‑1 (MQ=255)
                                                                                                                                        tgCGGACATGCCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAAc   <  7:100185/146‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                        
GACACCCAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACC  >  NC_002947/3252351‑3252633

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