Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,249,667 G→T 7.4% L429L (CTG→CTT ureC → urease subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,249,6670GT 7.4% 165.8 / 3.9 54L429L (CTG→CTTureCurease subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (22/27);  new base T (2/2);  total (25/29)
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.26e-01

GATGAAGCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATCGGCC  >  NC_002947/3249520‑3249804
                                                                                                                                                     |                                                                                                                                           
gATGAAGCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTgg                                                                                                                                            >  2:276439/1‑149 (MQ=255)
  tGAAGCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGc                                                                                                                                           <  3:79350/148‑1 (MQ=255)
     aGCTGCGTCGCGGGCCGCTGGCACCGGAGACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGAc                                                                                                                                       <  7:257263/149‑1 (MQ=255)
      gCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGgacgaca                                                                                                                                      >  6:235708/1‑148 (MQ=255)
       cTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCCAGCTTGCCGACCt                                                                                                                                     <  1:242231/149‑1 (MQ=255)
               cgGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCACACGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCtg                                                                                                                              <  1:276439/149‑1 (MQ=255)
                        ggCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAg                                                                                                                    >  5:255374/1‑149 (MQ=255)
                                    cccGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGc                                                                                                                       >  1:251041/1‑134 (MQ=255)
                                    cccGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGc                                                                                                                       <  2:251041/134‑1 (MQ=255)
                                    cccGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCgg                                                                                                        >  5:93205/1‑149 (MQ=255)
                                    cccGTAAAGCTACAACTTCCGGGTCAACCGCTACATCGCCAAGTACACCATCAACCCGCCAATCACCCATCGCATAGGTCACGAAGTCGGCTCCGAGCAGGTCGGCAACCTCGCCGACCTGCTCATCTGGTCGCCACCGATATTTGCg                                                                                                         >  8:193256/1‑148 (MQ=255)
                                     ccGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTc                                                                                                                         >  8:298600/1‑131 (MQ=255)
                                     ccGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTc                                                                                                                         <  7:298600/131‑1 (MQ=255)
                                       gTACAGCGACAACTTCCGGTTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCa                                                                                                     <  7:269961/149‑1 (MQ=255)
                                          cAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGc                                                                                                  <  1:240961/149‑1 (MQ=255)
                                              gACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGc                                                                                              <  5:62056/149‑1 (MQ=255)
                                               aCAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACATCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGc                                                                                              >  1:88857/1‑149 (MQ=255)
                                                  aCTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTg                                                                                          <  6:255374/149‑1 (MQ=255)
                                                    ttCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGt                                                                                        >  2:107697/1‑149 (MQ=255)
                                                    ttCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGt                                                                                        >  4:25873/1‑149 (MQ=255)
                                                              aGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGc                                                                            >  3:177167/1‑149 (MQ=255)
                                                               gcgcTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGc                                                                                      <  3:78935/140‑1 (MQ=255)
                                                               gcgcTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGc                                                                                      >  4:78935/1‑140 (MQ=255)
                                                               gcgcTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCg                                                                           <  5:308936/149‑1 (MQ=255)
                                                               gcgcTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCg                                                                           >  8:60125/1‑149 (MQ=255)
                                                                     cATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGgatga                                                                     >  5:155452/1‑149 (MQ=255)
                                                                     cATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGgatga                                                                     >  3:121181/1‑149 (MQ=255)
                                                                      tcccccAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGgatgat                                                                    <  5:130671/145‑1 (MQ=255)
                                                                       tCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGTGATGATc                                                                   <  2:54325/149‑1 (MQ=255)
                                                                       tCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATc                                                                   >  5:255378/1‑149 (MQ=255)
                                                                           cAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCa                                                               <  6:255378/149‑1 (MQ=255)
                                                                            aaGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCAc                                                              <  6:93205/149‑1 (MQ=255)
                                                                                  acCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGcc                                                        >  4:40806/1‑149 (MQ=255)
                                                                                        aaCCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATgg                                                   <  3:25873/148‑1 (MQ=255)
                                                                                        aaCCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATgg                                                   <  4:121181/148‑1 (MQ=255)
                                                                                         aCCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATGGGc                                                 <  2:88857/149‑1 (MQ=255)
                                                                                          cccGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGTTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATGGGCg                                                >  7:259882/1‑149 (MQ=255)
                                                                                           ccGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATGGGCGa                                               >  5:98993/1‑149 (MQ=255)
                                                                                            cGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATGGGCGAc                                              <  7:60125/149‑1 (MQ=255)
                                                                                             ggCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATGGGCGACa                                             >  5:349024/1‑149 (MQ=255)
                                                                                               cACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATGGGCGACATc                                           <  6:349024/149‑1 (MQ=255)
                                                                                               cACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGGGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATGGGCGACATc                                           <  4:310036/149‑1 (MQ=255)
                                                                                                     aCCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAAc                                        <  6:155452/147‑1 (MQ=255)
                                                                                                     aCCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACgg                                      >  1:236743/1‑149 (MQ=255)
                                                                                                                                cTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCgcgc                                                         >  1:127446/1‑104 (MQ=255)
                                                                                                                                cTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCgcgc                                                         <  2:127446/104‑1 (MQ=255)
                                                                                                                                cTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGc            <  6:98993/149‑1 (MQ=255)
                                                                                                                                cTCGGTGGAGGTAGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGGTCAACCTGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGt              <  1:107697/149‑1 (MQ=255)
                                                                                                                                    gtgtAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTa        >  3:24254/1‑149 (MQ=255)
                                                                                                                                    gTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTa        >  7:32927/1‑149 (MQ=255)
                                                                                                                                         ggTCGGCCAGCTTGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGc            <  1:6415/140‑1 (MQ=255)
                                                                                                                                         ggTCGGCCAGCTTGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGc            >  2:6415/1‑140 (MQ=255)
                                                                                                                                         ggTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATCGGc   <  7:193256/149‑1 (MQ=255)
                                                                                                                                          gTCGGCCAGCTTGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATCGGcc  >  7:189838/1‑149 (MQ=255)
                                                                                                                                          gTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATCGGcc  <  3:40806/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                           
GATGAAGCTGCGTCGCGGGCCGCTGGCACCGGACACCCCGTACAGCGACAACTTCCGGGTCAAGCGCTACATCGCCAAGTACACCATCAACCCGGCACTCACCCATGGCATCGGTCACGAAGTGGGCTCGGTGGAGGTCGGCAAGCTGGCCGACCTGGTGCTGTGGTCGCCAGCGTTCTTTGCGGTCAAGCCGGCCCTGGTGCTCAA‑‑GGGCGGGATGATCGTCACCGCGCCCATGGGCGACATCAACGGCTCCATACCCACCCCCCAGCCGGTGCACTATCGGCC  >  NC_002947/3249520‑3249804

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: