Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,247,257 T→C 8.0% L119P (CTG→CCG)  ureD → urease accessory protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,247,2570TC 8.0% 136.9 / 2.6 50L119P (CTG→CCG) ureDurease accessory protein
Reads supporting (aligned to +/- strand):  ref base T (23/23);  new base C (2/2);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.06e-01

GCCTGGAACTGGACATTCACCTGGAGCCGGGTAGCCACGCGCTGCTGACCATGCCCGGCGCCAGCAAGTTCTACCGCAGCATTGGCCCGACCGCGCGGTTGGCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCGCAGGACAGCATCTTCTTCTCCGGCGCCCGCGCCAGCCTCGCCAGCCGTTTTACCCTTGAACCGGGCGCCCGCCTGCTGGCCTGGGAAACCCTGTGCCTGGGCCGTCCGGTCATGCACGAGCGTTTCGACCACGGCGCCCTCGA  >  NC_002947/3247113‑3247404
                                                                                                                                                 |                                                                                                                                                   
gCCTGGAACTGGACATTCACCTGGAGCCGGGTAGCCACGCGCTGCTGACCATGCCCGGCGCCAGCAAGTTCTACCGCAGCATTGGCCCGACCGCGCGGTTGGCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCg                                                                                                                                                 >  6:280209/1‑149 (MQ=255)
gCCTGGAACTGGACATTCACCTGGAGCCGGGTAGCCACGCGCTGCTGACCATGCCCGGCGCCAGCAAGTTCTACCGCAGCATTGGCCCGACCGCGCGGTTGGCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCg                                                                                                                                                 >  4:220548/1‑149 (MQ=255)
 ccTGGAACTGGACATTCACCTGGAGCCGGGTAGCCACGCGCTGCTGACCATGCCCGGCGCCAGCAAGTTCTACCGCAGCATTGGCCCGACCGCGCGGTTGGCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCg                                                                                                                                                 >  1:84527/1‑148 (MQ=255)
 ccTGGAACTGGACATTCACCTGGAGCCGGGTAGCCACGCGCTGCTGACCATGCCCGGCGCCAGCAAGTTCTACCGCAGCATTGGCCCGACCGCGCGGTTGGCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCg                                                                                                                                                 >  3:189724/1‑148 (MQ=255)
  cTGGAACTGGACATTCACCTGGAGCCGGGTAGCCACGCGCTGCTGACCATGCCCGGCGCCAGCAAGTTCTACCGCAGCATTGGCCCGACCGCGCGGTTGGCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCg                                                                                                                                                 >  7:223618/1‑147 (MQ=255)
  cTGGAACTGGACATTCACCTGGAGCCGGGTAGCCACGCGCTGCTGACCATGCCCGGCGCCAGCAAGTTCTACCGCAGCATTGGCCCGACCGCGCGGTTGGCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCg                                                                                                                                                 >  3:8418/1‑147 (MQ=255)
        cTGGACATTCACCTGGAGCCGGGTAGCCACGCGCTGCTGACCATGCCCGGCGCCAGCGAATTTTAACCCCAGCATTGGCCCGACCCGGCGGTTGGCCCACCGCTTCCATGTGGCCCCGGGCGCCACCCGGGAGTGGCTGCCGCAgggac                                                                                                                                          >  1:90648/1‑146 (MQ=255)
         tGGACATTCACCTGGAGCCGGGTAGCCACGCGCTGCTGACCATGCCCGGCGCCAGCAAGTTCTACCGCAGCATTGGCCCGCCCGCGCGGTTGGCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCCGCCGCAGGACAg                                                                                                                                         >  8:115646/1‑148 (MQ=255)
         tGGACATTCACCTGGAGCCGGGTAGCCACGCGCTGCTGACCATGCCCGGCGCCAGCAAGTTCTACCGCAGCATTGGCCCGACCGCGCGGTTGGCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCCGCCGCAGGACAg                                                                                                                                         >  8:111072/1‑148 (MQ=255)
            aCATTCACCTGGAGCCGGGTAGCCACGCGCTGCTGACCATGCCCGGCGCCAGCAAGTTCTACCGCAGCATTGGCCCGACCGCGCGGTTGGCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCGCAGGACAGCAtc                                                                                                                                     <  2:84527/149‑1 (MQ=255)
               ttCACCTGGAGCCGGGTAGCCACGCGCTGCTGACCATGCCCGGCGCCAGCAAGTTCTACCGCAGCATTGGCCCGACCGCGCGGTTGGCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCGCAGGACAGCAtcttc                                                                                                                                  <  8:252086/149‑1 (MQ=255)
               ttCACCTGGAGCCGGGTAGCCACGCGCTGCTGACCATGCCCGGCGCCAGCAAGTTCTACCGCAGCATTGGCCCGACCGCGCGGTTGGCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCGCAGGACAGCAtcttc                                                                                                                                  <  2:5545/149‑1 (MQ=255)
                 cACCTGGAGCCGGGTAGCCACGCGCTGCTGACCATGCCCGGCGCCAGCAAGTTCTACCGCAGCATTGGCCCGACCGCGCGGTTGGCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCGCAGGACAGCAtcttctt                                                                                                                                >  4:175797/1‑149 (MQ=255)
                 cACCTGGAGCCGGGTAGCCACGCGCTGCTGACCATGCCCGGCGCCAGCAAGTTCTACCGCAGCATTGGCCCGACCGCGCGGTTGGCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCGCAGGACAGCAtcttct                                                                                                                                 <  4:28765/148‑1 (MQ=255)
                  aCCTGGAGCCGGGTAGCCACGCGCTGCTGACCATGCCCGGCGCCAGCAAGTTCTACCGCAGCATTGGCCCGACCGCGCGGTTGGCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCGCAGGACAGCAtcttcttc                                                                                                                               <  4:210658/149‑1 (MQ=255)
                       gAGCCGGGTAGCCACGCGCTGCTGACCATGCCCGGCGCCAGCAAGTTCTACCGCAGCATTGGCCCGACCGCGCGGTTGGCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCGCAGGACAGCATCTTCTTCTCCgg                                                                                                                          <  8:286068/149‑1 (MQ=255)
                            gggTAGCCACGCGCTGCTGACCATGCCCGGCGCCAGCAAGTTCTACCGCAGCATTGGCCCGACCGCGCGGTTGGCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCGCAGGACAGCAtcttcttc                                                                                                                               >  6:257720/1‑139 (MQ=255)
                            gggTAGCCACGCGCTGCTGACCATGCCCGGCGCCAGCAAGTTCTACCGCAGCATTGGCCCGACCGCGCGGTTGGCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCGCAGGACAGCAtcttcttc                                                                                                                               <  5:257720/139‑1 (MQ=255)
                                    aCGCGCTGCTGACCATGCCCGGCGCCAGCAAGTTCTACCGCAGCATTGGCCCGACCGCGCGGTTGGCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCGCAGGACAGCATCTTCTTCTCCGGCGCCCGCGCCAGc                                                                                                             <  4:8418/149‑1 (MQ=255)
                                               aCCATGCCCGGCGCCAGCAAGTTCTACCGCAGCATTGGCCCGACCGCGCGGTTGGCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCGCAGGACAGCATCTACTTCTCCGGCGCCCGCGCCAGCCTCGCCAGCCg                                                                                                  >  7:164255/1‑149 (MQ=255)
                                                      ccGGCGCCAGCAAGTTCTACCGCAGCATTGGCCCGACCGCGCGGTTGGCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCGCAGGACAGCATCTTCTTCTCCGGCGCCCGCGCCAGCCTCGCCAGCCGTTTTAcc                                                                                           <  4:189724/149‑1 (MQ=255)
                                                       cGGCGCCAGCAAGTTCTACCGCAGCATTGGCCCGACCGCGCGGTTGGCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCGCAGGACAGCATCTTCTTCTCCGGCGCCCGCGCCAGCCTCGCCAGCCGTTTTAccc                                                                                          <  3:175797/149‑1 (MQ=255)
                                                        ggCGCCAGCAAGTTCTACCGCAGCATTGGCCCGACCGCGCGGTTGGCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCGCAGGACAGCATCTTCTTCTCCGGCGCCCGCGCCAGCCTCGCCAGCCGTTTTACCCt                                                                                         >  3:215373/1‑149 (MQ=255)
                                                             cAGCAAGTTCTACCGCAGCATTGGCCCGACCGCGCGGTTGGCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCGCAGGACAGCATCTTCTTCTCCGGCGCCCGCGCCAGCCTCGCCAGCCGTTTTACCCTTGAAc                                                                                    >  4:158991/1‑149 (MQ=255)
                                                                    gTTCTACCGCAGCATTGGCCCGACCGCGCGGTTGGCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCGCAGGACAGCATCTTCTTCTCCGGCGCCCGCGCCAGCCTCGCCAGCCGTTTTACCCTTGAACCGGgcg                                                                              >  1:312300/1‑149 (MQ=255)
                                                                    gTTCTACCGCAGCATTGGCCCGACCGCGCGGTTGGCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCGCAGGACAGCATCTTCTTCTCCGGCGCCCGCGCCAGCCTCGCCAGCCGTTTTACCCTTGAACCGGgcg                                                                              <  5:280209/149‑1 (MQ=255)
                                                                         aCCGCAGCATTGGCCCGACCGCGCGGTTGGCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCGCAGGACAGCATCTTCTTCTCCGGCGCCCGCGCCAGCCTCGCCAGCCGTTTTACCCTTGAACCGGGcgcccgc                                                                         >  5:57670/1‑149 (MQ=255)
                                                                           cGCAGCCTTGGCCCGACCGCGCGGTTGGCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCGCAGGACAGCATCTTCTTCTCCGGCGCCCGCGCCAGCCTCGCCAGCCGTTTTACCCTTGAACCGGGCGCCCGCct                                                                       <  4:23204/149‑1 (MQ=255)
                                                                           cGCAGCATTGGCCCGACCGCGCGGTTGGCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCGCAGGACAGCATCTTCTTCTCCGgcg                                                                                                                        >  5:212002/1‑100 (MQ=255)
                                                                           cGCAGCATTGGCCCGACCGCGCGGTTGGCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCGCAGGACAGCATCTTCTTCTCCGgcg                                                                                                                        <  6:212002/100‑1 (MQ=255)
                                                                           cGCAGCATTGGCCCGACCGCGCGGTTGGCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCGCAGGACAGCATCTTCTTCTCCGGCGCCCGCGCCAGCCTCGCCAGCCGTTTTACCCTTGAACCGGGCGCCCGCct                                                                       >  4:54686/1‑149 (MQ=255)
                                                                             cagcaTTGGCCCGACCGCGCGGTTGGCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCGCAGGACAGCATCTTCTTCTCCGGCGCCCGCGCCAGCCTCGCCAGCCGTTTTACCCTTGAACCGGGCGCCCGCctgc                                                                     >  6:81733/1‑149 (MQ=255)
                                                                              agcaTTGGCCCGACCGCGCGGTTGGCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCGCAGGACAGCATCTTCTTCTCCGGCGCCCGCGCCAGCCTCGCCAGCCGTTTTACCCTTGAACCGGGCGCCCGCctgct                                                                    <  4:215373/149‑1 (MQ=255)
                                                                                     gCCCGACCGCGCGGTTGGCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCGCAGGACAGCATCTTCTTCTCCGGCGCCCGCGCCAGCCTCGCCAGCCGTTTTACCCTTGAACCGGGCGCCCGCCTGCTGGCCTgg                                                             <  3:27199/149‑1 (MQ=255)
                                                                                     gCCCGACCGCGCGGTTGGCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCGCAGGACAGCATCTTCTTCTCCGGCGCCCGCGCCAGCCTCGCCAGCCGTTTTACCCTTGAACCGGGCGCCCGCCTGCTGGCCTg                                                              <  8:223618/148‑1 (MQ=255)
                                                                                       ccGACCGCGCGGTTGGCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCGCAGGACAGCATCTTCTTCTCCGGCGCCCGCGCCAGCCTCGCCAGCCGTTTTACCCTTGAACCGGGCGCCCGCCTGCTGGCCTggg                                                            <  4:165128/148‑1 (MQ=255)
                                                                                       ccGACCGCGCGGTTGGCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCGCAGGACAGCATCTTCTTCTCCGGCGCCCGCGCCAGCCTCGCCAGCCGTTTTACCCTTGAACCGGGCGCCCGCCTGCTGGCCTggg                                                            >  3:165128/1‑148 (MQ=255)
                                                                                                      gCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCGCAGGACAGCATCTTCTTCTCCGGCGCCCGCGCCAGCCTCGCCAGCCGTTTTACCCTTGAACCGGGCGCCCGCCTGCTGGCCTGGGAAACCCTGTGCCTgg                                             <  3:54686/148‑1 (MQ=255)
                                                                                                      gCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCCGCCGCAGGACAGCATCTTCTTCTCCAGCGCCCGCGCCAGCCTCGCCAGCCGTTTTACCCTTGAACCGGGCGCCCGCCTGCTGGCCTGGGAAACCCTGTGCCTgg                                             <  7:111072/148‑1 (MQ=255)
                                                                                                      gCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCCGCCGCAGGACAGCATCTTCTTCTCCAGCGCCCGCGCCAGCCTCGCCAGCCGTTTTACCCTTGAACCGGGCGCCCGCCTGCTGGCCTGGGAAACCCTGTGCCTgg                                             <  7:115646/148‑1 (MQ=255)
                                                                                                               ttCCATCTGGGCGCGGGCAGCTCCCTGGTGTGGCTGCCGGAGGACAGCATCTTCTTCTCCGGCGCCCGCGCCAGCCTCGCCAGCCGTTTTACCCTTGAACCGGGCGCCCGCCTGCTGGCCTg                                                              <  6:176167/122‑1 (MQ=255)
                                                                                                               ttCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCGCAGGACAGCATCTTCTTCTCCGGCGCCCGCGCCAGCCTCGCCAGCCGTTTTACCCTTGAACCGGGCGCCCGCCTGCTGGCCTg                                                              >  5:176167/1‑122 (MQ=255)
                                                                                                                  cATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCGCAGGACAGCATCTTCTTCTCCGGCGCCCGCGCCAGCCTCGCCAGCCGTTTTACCCTTGAACCGGGCGCCCGCCTGCTGGCCTGGGAAACCCTGTGCCTGGGCCGTCCGGTCAt                                <  3:110884/149‑1 (MQ=255)
                                                                                                                        gCCGCGGGCAGCACCCTGGAGTGGCTGCCGCAGGACAGCATCTTCTTCTCCGGCGCCCGCGCCAGCCTCGCCAGCCGTTTTACCCTTGAACCGGGCGCCCGCCTGCTGGCCTGGGAAACCCTGTGCCTGGGCCGTCCGGTCATGc                              <  3:220548/145‑1 (MQ=255)
                                                                                                                         ccGCGGGCAGCACCCTGGAGTGGCTGCCGCAGGACAGCATCTTCTTCTCCGGCGCCCGCGCCAGCCTCGCCAGCCGTTTTACCCTTGAACCGGGCGCCCGCCTGCTGGCCTGGGAAACCCTGTGCCTGGGCCGTCCGGTCATGCACGAg                         >  5:268373/1‑149 (MQ=255)
                                                                                                                         ccGCGGGCAGCACCCTGGAGTGGCTGCCGCAGGACAGCATCTTCTTCTCCGGCGCCCGCGCCAGCCTCGCCAGCCGTTTTACCCTTGAACAGGGCGCCCGCCTGCTGGCCTGGGAAACCCTGTGCCTGGGCCGTCCGGTCATGCACGAg                         >  7:25613/1‑149 (MQ=255)
                                                                                                                                   caCCCTGGAGTGGCTGCCGCAGGACAGCATCTTCTTCTCCGGCGCCCGCGCCAGCCTCGCCAGCCGTTTTACCCTTGAACCGGGCGCCCGCCTGCTGGCCTGGGAAACCCTGTGCCTGGGCCGTCCgg                                    >  8:293059/1‑128 (MQ=255)
                                                                                                                                   caCCCTGGAGTGGCTGCCGCAGGACAGCATCTTCTTCTCCGGCGCCCGCGCCAGCCTCGCCAGCCGTTTTACCCTTGAACCGGGCGCCCGCCTGCTGGCCTGGGAAACCCTGTGCCTGGGCCGTCCgg                                    <  7:293059/128‑1 (MQ=255)
                                                                                                                                     cccTGGAGTGGCTGCCGCAGGACAGCATCTTCTTCTCCGGCGCCCGCGCCAGCCTCGCCAGCCGTTTTACCCTTGAACCGGGCGCCCGCCTGCTGGCCTGGGAAACCCTGTGCCTGGGCCGTCCGGTCATGCACGAGCGTTTCGACCAc             <  7:40984/149‑1 (MQ=255)
                                                                                                                                            gTGGCTGCCGCAGGACAGCATCTTCTTCTCCGGCGCCCGCGCCAGCCTCGCCAGCCGTTTTACCCTTGAACCGGGCGCCCGCCTGCTGGCCTGGGAAACCCTGTGCCTGGGCCGTCCGGTCATGCACGAGCGTTTCGACCACGGCGccc      >  5:149323/1‑149 (MQ=255)
                                                                                                                                                 tGCCGCAGGACAGCATCTTCTTCTCCGGCGCCCGCGCCAGCCTCGCCAGCCGTTTTACCCTTGAACCGGGCGCCCGCCTGCTGGCCTGGGAAACCCTGTGCCTGGGCCGTCCGGTCATGCACGAGCGTTTCGACCACGGCGCCCTCGa  >  5:192363/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GCCTGGAACTGGACATTCACCTGGAGCCGGGTAGCCACGCGCTGCTGACCATGCCCGGCGCCAGCAAGTTCTACCGCAGCATTGGCCCGACCGCGCGGTTGGCCCAACGCTTCCATCTGGCCGCGGGCAGCACCCTGGAGTGGCTGCCGCAGGACAGCATCTTCTTCTCCGGCGCCCGCGCCAGCCTCGCCAGCCGTTTTACCCTTGAACCGGGCGCCCGCCTGCTGGCCTGGGAAACCCTGTGCCTGGGCCGTCCGGTCATGCACGAGCGTTTCGACCACGGCGCCCTCGA  >  NC_002947/3247113‑3247404

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: