Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,252,535 G→A 16.3% 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,5350GA16.3% 91.7 / 10.2 42intergenic (+22/+221)ureG/PP_2851urease accessory protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (21/14);  new base A (4/3);  total (25/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

AGCGCGACACCCAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCT  >  NC_002947/3252389‑3252683
                                                                                                                                                  |                                                                                                                                                    
aGCGCGACACCCAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCAt                                                                                                                                                   <  2:458448/150‑1 (MQ=255)
   gcgACACCCAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTc                                                                                                                                                >  1:1128569/1‑150 (MQ=255)
       cacCCAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGca                                                                                                                                            <  2:552795/150‑1 (MQ=255)
       cacCCAGCGCAAGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGca                                                                                                                                            >  2:50306/1‑150 (MQ=255)
          ccAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACag                                                                                                                                         >  1:326440/1‑150 (MQ=255)
           cAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAGCCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACaga                                                                                                                                        >  1:976871/1‑150 (MQ=255)
           cAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACaga                                                                                                                                        >  1:862562/1‑150 (MQ=255)
           cAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGACATGCCCGCACaga                                                                                                                                        <  1:1106624/150‑1 (MQ=255)
            aGCGCATTCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACagag                                                                                                                                       >  2:829454/1‑150 (MQ=255)
                 aTGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGTCTCTGTGCGGACATGCCCGCACagagtcctc                                                                                                                                  <  1:422809/150‑6 (MQ=255)
                         gCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCggg                                                                                                                                                      >  2:295398/1‑122 (MQ=255)
                         gCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCggg                                                                                                                                                      <  1:295398/122‑1 (MQ=255)
                             cGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGtggctgg                                                                                                                      >  2:238594/1‑150 (MQ=255)
                                 ccTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCAcc                                                                                                                  >  2:454380/1‑150 (MQ=255)
                                             gCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGACATGCCCGCACAGAGTcctcga                                                                                                                                >  2:863142/1‑120 (MQ=39)
                                             gCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATCCTATCGTCGAGCGCGGGATGCTGGGGGTGGGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCggg                                                                                                      >  2:204210/1‑149 (MQ=255)
                                                    gaagaaGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTcgcg                                                                                               <  1:829454/150‑1 (MQ=255)
                                                        aagGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGtt                                                                                           <  2:326440/150‑1 (MQ=255)
                                                               aGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCcgc                                                                                    <  1:50306/150‑1 (MQ=255)
                                                                 gggTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTTTGCCGGTGTTCGCGGGTTTACCCGCGa                                                                                  >  1:1028969/1‑150 (MQ=255)
                                                                  ggTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGt                                                                                            <  1:171749/139‑1 (MQ=255)
                                                                  ggTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGt                                                                                            >  2:171749/1‑139 (MQ=255)
                                                                          gCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCAc                                                                                                                   >  2:1011108/1‑108 (MQ=255)
                                                                          gCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCAc                                                                                                                   <  1:1011108/108‑1 (MQ=255)
                                                                           ccGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGc                                                                        <  2:862562/150‑1 (MQ=255)
                                                                            cggaaTCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCAAAGAGACCAACGCGGTCCATACCGAGATCGCTGACAGTGTTGTCGTTTTTACCCGCGAATACGGTAGCg                                                                       <  1:454380/146‑1 (MQ=255)
                                                                            cGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGACATGCCCGCACagag                                                                                                                                       >  2:131715/1‑86 (MQ=255)
                                                                            cGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGACATGCCCGCACagag                                                                                                                                       <  1:131715/86‑1 (MQ=255)
                                                                                   gATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGCGAGAAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGcaa                                                                >  1:509764/1‑150 (MQ=255)
                                                                                             tcgAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGcaca                                                                                                                                          >  2:207126/1‑66 (MQ=255)
                                                                                             tcgAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGcaca                                                                                                                                          <  1:207126/66‑1 (MQ=255)
                                                                                             tcgAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGAAAGAATGAAcgg                                                      >  1:654510/1‑150 (MQ=255)
                                                                                               gAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAAcggcg                                                    >  1:539987/1‑150 (MQ=255)
                                                                                                        aTGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAATGGCGGGTGGAAAt                                           <  1:238594/150‑1 (MQ=255)
                                                                                                            tGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGc                                       >  2:1031599/1‑150 (MQ=255)
                                                                                                                  tGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCcagcag                                 >  1:1019617/1‑150 (MQ=255)
                                                                                                                     gAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGggc                              <  1:1031599/150‑1 (MQ=255)
                                                                                                                            aCAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTAGCCAGCAGGGCCGGCctc                       >  1:857932/1‑150 (MQ=255)
                                                                                                                               gATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTc                    <  2:509764/150‑1 (MQ=255)
                                                                                                                                   gAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCgg                >  2:943129/1‑150 (MQ=255)
                                                                                                                                        ctctGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTaaa           >  1:713674/1‑150 (MQ=255)
                                                                                                                                         tctGTGCGGACATGCCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCTAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGACATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAAc          >  1:718146/1‑150 (MQ=255)
                                                                                                                                          ctGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAAcc         <  2:539987/150‑1 (MQ=255)
                                                                                                                                           tgtgCGGACATGCCCGCACAGAGGCCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAAccc        >  2:88035/1‑150 (MQ=255)
                                                                                                                                                 ggCATGTCCGCACAGAGACCAACGCGGTGGATGCGAGGTTCTCGGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCt  <  2:1188580/150‑1 (MQ=21)
                                                                                                                                                  |                                                                                                                                                    
AGCGCGACACCCAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTGCGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCT  >  NC_002947/3252389‑3252683

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