Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,252,552 T→A 11.8% intergenic (+82/+161) 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,5520TA11.8% 86.9 / 4.1 34intergenic (+82/+161)ureG/PP_2851urease accessory protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (20/10);  new base A (2/2);  total (22/12)
Fisher's exact test for biased strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

GGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCAGGACAGTTGCCCCCCACA  >  NC_002947/3252404‑3252693
                                                                                                                                                    |                                                                                                                                             
gggCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGttt                                                                                                                                               <  6:123762/149‑1 (MQ=255)
    gAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTATTCGCGGGTAAAccc                                                                                                                                           >  7:19781/1‑149 (MQ=255)
       gggTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCcgc                                                                                                                                         <  4:68341/148‑1 (MQ=255)
         gTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCACCGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGaa                                                                                                                                      >  6:170841/1‑149 (MQ=255)
               ggCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTATTCGCGGGTAAACCCGCGAacaccgg                                                                                                                                <  8:19781/149‑4 (MQ=38)
                gCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTa                                                                                                                               <  6:48281/149‑1 (MQ=255)
                      aTCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTTACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTAGCGGGTTTACCCGCGAATACGGTAGCGGcc                                                                                                                         >  4:156576/1‑149 (MQ=255)
                                                gCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTAAACCCGCGaa                                                                                                                                      <  7:230177/110‑1 (MQ=255)
                                                gCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTAAACCCGCGaa                                                                                                                                      >  8:230177/1‑110 (MQ=255)
                                                  tGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTgg                                                                                             >  5:170308/1‑149 (MQ=255)
                                                      ggTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGt                                                                                                                                <  7:215469/110‑1 (MQ=255)
                                                      ggTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGt                                                                                                                                >  8:215469/1‑110 (MQ=255)
                                                        tGCGAGGCTGACAGATCGTGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACGGGCTCTGCCGGTGTGCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCcagca                                                                                       <  5:170841/149‑1 (MQ=255)
                                                                     gATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCtt                                                                          >  7:21760/1‑149 (MQ=255)
                                                                        cGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGcaa                                                                                                                     <  2:163886/103‑1 (MQ=255)
                                                                        cGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGcaa                                                                                                                     >  1:163886/1‑103 (MQ=255)
                                                                         gAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTg                                                                                              >  8:57556/1‑125 (MQ=255)
                                                                                        gCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTcgc                                                                       >  8:153909/1‑133 (MQ=255)
                                                                                        gCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTcgc                                                                       <  7:153909/133‑1 (MQ=255)
                                                                                                 cacaGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGa                                              <  8:21760/149‑1 (MQ=255)
                                                                                                 cacaGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGa                                              >  2:109746/1‑149 (MQ=255)
                                                                                                  acaGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGAc                                             >  5:163300/1‑149 (MQ=255)
                                                                                                          cAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGaa                                                                                                  >  1:141646/1‑88 (MQ=255)
                                                                                                          cAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGaa                                                                                                  <  2:141646/88‑1 (MQ=255)
                                                                                                          cAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAgg                                     >  5:240200/1‑149 (MQ=255)
                                                                                                                ggTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTc                                >  7:22620/1‑148 (MQ=255)
                                                                                                                 gTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGa                              <  1:109746/149‑1 (MQ=255)
                                                                                                                        acACCGGCAGAGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGa                       >  4:212716/2‑149 (MQ=255)
                                                                                                                         cACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAg                      <  3:156576/149‑1 (MQ=255)
                                                                                                                         cACCGGCAGAGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGACGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAg                      >  6:98629/1‑149 (MQ=255)
                                                                                                                                    ccGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCAGGACAGTTg           >  7:132978/1‑149 (MQ=255)
                                                                                                                                    ccGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCAGGACAGTTg           >  6:103019/1‑149 (MQ=255)
                                                                                                                                             gcgGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGtg                         >  8:86510/1‑125 (MQ=255)
                                                                                                                                             gcgGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCAGGACAGTTGCCCCCcaca  >  8:55126/1‑149 (MQ=255)
                                                                                                                                             gcgGGTTTACCAGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCAGGACAGTTGCCCCCcac   >  8:51443/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                             
GGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGGCGACCGCGCAGGCCGTCGAATGTTGAGCAGGACAGTTGCCCCCCACA  >  NC_002947/3252404‑3252693

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: