Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,277,599 G→C 9.8% S174R (AGC→AGG)  ampC ← beta‑lactamase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,277,5990GC 9.8% 118.4 / 4.6 41S174R (AGC→AGG) ampCbeta‑lactamase
Reads supporting (aligned to +/- strand):  ref base G (16/21);  new base C (2/2);  total (18/23)
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.95e-01

CCCTTGGGCATAAAGCCCCTGGGCGCTGGCCGGAACCTGAAGGTAGGTGTTTTTCAGGCCCATCTTCGTCAGCAGGCCCTCACTCATCAGCGTGGAAAACGGTTGGTGCTGCGCACGGGCGGCCAGGTCACCGAACAAACCCAGGCTGGGGTTGGAGTAGCAGCGCTGGGTGCCAGGTTTGGCACGAGGTTGCCAGTGGCGGAAGAAATCGACTACCTGTTGCGGTGTTTGTACGGTGTCGGGGAACTGCAGTGGCAGGCAGTCGGCGCTATAGGCGCCCAGTTCCAGC  >  NC_002947/3277455‑3277743
                                                                                                                                                |                                                                                                                                                
cccTTGGGCATAAAGCCCCTGGGCGCTGGCCGGAACCTGAAGGTAGGTGTTTTTCAGGCCCATCTTCGTCAGCAGGCCCTCACTCATCAGCGTGGAAAACGGTTGGTGCTGCGCACGGGCGGCCAGGTCACCGAACAAACCCAGGCTg                                                                                                                                               <  7:349343/148‑1 (MQ=255)
      ggCATAAAGCCCCTGGGCGCTGGCCGGAACCTGAAGGTAGGTGTTTTTCAGGCCCATCTTCGTCAGCAGGCCCTCACTCATCAGCGTGGAAAACGGTTGGTGCTGCGCACGGGCGGCCAGGTCACCGAACAAACCCAGGCTGGGGTTgg                                                                                                                                        <  3:108907/149‑1 (MQ=255)
      ggCATAAAGCCCCTGGGCGCTGGCCGGAACCTGAAGGTAGGTGTTTTTCAGGCCCATCTTCGTCAGCAGGCCCTCACTCATCAGCGTGGAAAACGGTTGGTGCTGCGCACGGGCGGCCAGGTCACCGAACAAACCCAGGCTGGGGTTg                                                                                                                                         <  1:225451/148‑1 (MQ=255)
       gCATAAAGCCCCTGGGCGCTGGCCGGAACCTGAAGGTAGGTGTTTTTCAGGCCCATCTTCGTCAGCAGGCCCTCACTCATCAGCGTGGAAAACGGTTGGTGCTGCGCACGGGCGGCCAGGTCACCGAACAAACCCAGGCTGGGGTTGGa                                                                                                                                       <  4:161372/149‑1 (MQ=255)
       gCATAAAGCCCCTGGGCGCTGGCCGGAACCTGAAGGTAGGTGTTTTTCAGGCCCATCTTCGTCAGCAGGCCCTCACTCATCAGCGTGGAAAACGGTTGGTGCTGCGCACGGGCGGCCAGGTCACCGAACAAACCCAGGCTGGGGTTGGa                                                                                                                                       <  1:188190/149‑1 (MQ=255)
       gCATAAAGCCCCTGGGCGCTGGCCGGAACCTGAAGGTAGGTGTTTTTCAGGCCCATCTTCGTCAGCAGGCCCTCACTCATCAGCGTGGAAAACGGGTGGTGCTGCGCACGGGCGGCCAGGTCACCGTACAAACCCAGGCTGGGGTTGGa                                                                                                                                       <  6:274336/149‑1 (MQ=255)
          tAAAGCCCCTGGGCGCTGGCCGGAACCTGAAGGTAGGTGTTTTTCAGGCCCATCTTCGTCAGCAGGCCCTCACTCATCAGCGTGGAAAACGGTTGGTGCTGCGCACGGGCGGCCAGGTCACCGAACAAACCCAGGCGGGGGTTGGAg                                                                                                                                      >  2:279081/1‑147 (MQ=255)
             aGCCCCTGGGCGCTGGCCGGAACCTGAAGGTAGGTGTTTTTCAGGCCCATCTTCGTCAGCAGGCCCTCACTCATCAGCGTGGAAAACGGTTGGTGCTGCGCACGGGCGGCCAGGTCACCGAACAAACCCAGGCTGGGGTTGGAGTag                                                                                                                                   >  2:251387/1‑147 (MQ=255)
               ccccTGGGCGCTGGCCGGAACCTGAAGGTAGGTGTTTTTCAGGCCCATCTTCGTCAGCAGGCCCTCACTCATCAGCGTGGAAAACGGTTGGTGCTGCGCACGGGCGGCCAGGTCACCGAACAAACCCAGGCTGGGGTTGGAGTagcagc                                                                                                                               >  7:269204/1‑149 (MQ=255)
                     ggCGCTGGCCGGAACCTGAAGGTAGGTGTTTTTCAGGCCCATCTTCGTCAGCAGGCCCTCACTCATCAGCGTGGAAAACGGTTGGTGCTGCGCACGGGCGGCCAGGTCACCGAACAAACCCAGGCTGGGGTTGGAGTAGCAGCGCTg                                                                                                                           <  1:251387/147‑1 (MQ=255)
                            gCCGGAACCTGAAGGTAGGTGTTTTTCAGGCCCATCTTCGTCAGCAGGCCCTCACTCATCAGCGTGGAAAACGGTTGGTGCTGCGCACGGGCGGCCAGGTCACCGAACAAACCCAGGCTGGGGTTGGAGTAGCAGCGCTGGGTGCCAgg                                                                                                                  >  6:235020/1‑149 (MQ=255)
                               ggAACCTGAAGGTAGGTGTTTTTCAGGCCCATCTTCGTCAGCAGGCCCTCACTCATCAGCGTGGAAAACGGTTGGTGCTGCGCACGGGCGGCCAGGTCACCGAACAAACCCAGGCTGGGGTTGGAGTAGCAGCGCTGGGTGCCAGGttt                                                                                                               <  2:178555/149‑1 (MQ=255)
                                  aCCTGAAGGTAGGTGTTTTTCAGGCCCATCTTCGTCAGCAGGCCCTCACTCATCAGCGTGGAAAACGGTTGGTGCTGCGCACGGGCGGCCAGGTCACCGAACAAACCCAGGCTGGGGTTGGAGTAGCAGCGCTGGGTGCCAGGTTTGGc                                                                                                            <  1:156414/149‑1 (MQ=255)
                                    cTGAAGGTAGGTGTTTTTCAGGCCCATCTTCGTCAGCAGGCCCTCACTCATCAGCGTGGAAAACGGTTGGTGCTGCGCACGGGCGGCCAGGTCACCGAACAAACCCAGGCTGGGGTTGGAGTAGCAGCGCTGGGTGCCAGGTTTGGCAc                                                                                                          <  5:235020/149‑1 (MQ=255)
                                           taggtGTTTTTCAGGCCCATCTTCGTCAGCAGGCCCTCACTCATCAGCGTGGAAAACGGTTGGTGCTGCGCACGGGCGGCCAGGTCACCGAACAAACCCAGGCTGGGGTTGGAGTAGCAGCGCTGGGTGCCAGGTTTGGCACGAGGTTg                                                                                                   <  1:111156/149‑1 (MQ=255)
                                             ggtGTTTTTCAGGCCCATCTTCGTCAGCAGGCCCTCACTCATCAGCGTGGAAAACGGTTGGTGCTGCGCACGGGCGGCCAGGTCACCGAACAAACCCAGGCTGGGGTTGGAGTAGCAGCGCTGGGTGCCAGGTTTGGCACGAGGTTGc                                                                                                  <  6:148707/148‑1 (MQ=255)
                                                         gCCGATCTTCGTCAGCAGGCCCTCTCTCATCAGCGTGGAAAACGGTTGGTGCTGCGCACGGGCGGCCAGGTCACCGAACAAACCCAGGCTGGGGTTg                                                                                                                                         <  5:308749/97‑1 (MQ=255)
                                                           ccATCTTCGTCAGCAGGCCCTCACTCATCAGCGTGGAAAACGGTTGGTGCTGCGCACGGGCGGCCAGGTCACCGAACAAACCCAGGCTGGGGTTGGAGTAGCAGCGCTGGGTGCCAGGTTTGGCACGAGGTTGCCAGTGGCGGAAGaaa                                                                                   >  3:127777/1‑149 (MQ=255)
                                                                          ggCCCTCACTCATCAGCGTGGAAAACGGTTGGTGCTGCGCACGGGCGGCCAGGTCACCGAACAAACCCAGGCTGGGGTTGGAGTAGCAGCGCTGGGTGCCAGGTTTGGCACGAGGTTGCCAGTGGCGGAAGAAATCGACTACCTGTTGc                                                                    <  1:191681/149‑1 (MQ=255)
                                                                                 actcaTCAGCGTGGAAAACGGTTGGTGCTGCGCACGGGCGGCCAGGTCACCGAACAAACCCAGGCTGGGGTTGGAGTAGCAGCGCTGGGTGCCAGGTTTGGCACGAGGTTGCCAGTGGCGGAAGAAATCGACTACCTGTTGCGgtgt                                                               <  1:279081/147‑1 (MQ=255)
                                                                                  ctcaTCAGCGTGGAAAACGGTTGGTGCTGCGCACGGGCGGCCAGGTCACCGAACAAACCCAGGCTGGGGTTGGAGTAGCAGCGCTGGGTGCCAGGTTTGGCACGAGGTTGCCAGTGGCGGAAGAAATCGACTACCTGTTGCGGTGTTTg                                                            >  3:2786/1‑149 (MQ=255)
                                                                                      tcaGCGTGGAAAACGGTTGGTGCTGCGCACGGGCGGCCAGGTCACCGAACAAACCCAGGCTGGGGTTGGAGTAGCAGCGCTGGGTGCCAGGTTTGGCACGAGGTTGCCAGTGGCGGAAGAAATCGACTACCTGTTGCGGTGTTTGTACg                                                        >  1:338538/1‑149 (MQ=255)
                                                                                      tcaGCGTGGAAAACGGTTGGTGCTGCGCACGGGCGGCCAGGTCACCGAACAAACCCAGGCTGGGGTTGGAGTAGCAGCGCTGGGTGCCAGGTTTGGCACGAGGTTGCCAGTGGCGGAAGAAATCGACTACCTGTTGCGGTGTTTGTACg                                                        <  2:338538/149‑1 (MQ=255)
                                                                                             ggAAAACGGTTGGTGCTGCGCACGGGCGGCCAGGTCACCGAACAAACCCAGGCTGGGGTTGGAGTAGCAGCGCTGGGTGCCAGGTTTGGCACGAGGTTGCCAGTGGCGGAAGAAATCGACTACCTGTTGCGGTGTTTGTACGGTGTCgg                                                 >  5:253021/1‑149 (MQ=255)
                                                                                              gAAAACGGTTGGTGCTGCGCACGGGCGGCCAGGTCACCGAACAAACCCAGGCTGGGGTTGGAGTAGCAGCGCTGGGTGCCAGGTTTGGCACGAGGTTGCCAGTGGCGGAAGAAATCGACTACCTGTTGCGGTGTTTGTACGGTGTCggg                                                <  8:56329/149‑1 (MQ=255)
                                                                                                 aaCGGTTGGTGCTGCGCACGGGCGGCCAGGTCACCGAACAAACCCAGGCTGGGGTTGGAGTAGCAGCGCTGGGTGCCAGGTTTGGCACGAGGTTGCCAGTGGCGGAAGAAATCGACTACCTGTTGCGGTGTTTGTACGGTGTCGGGGaa                                             >  1:170566/1‑149 (MQ=255)
                                                                                                   cGGTTGGTGCTGCGCACGGGCGGCCAGGTCACCGAACAAACCCAGGCTGGGGTTGGAGTAGCAGCGCTGGGTGCCAGGTTTGGCACGAGGt                                                                                                     >  4:235241/1‑91 (MQ=255)
                                                                                                   cGGTTGGTGCTGCGCACGGGCGGCCAGGTCACCGAACAAACCCAGGCTGGGGTTGGAGTAGCAGCGCTGGGTGCCAGGTTTGGCACGAGGt                                                                                                     <  3:235241/91‑1 (MQ=255)
                                                                                                           gctgcGCACGGGCGGCCAGGTCACCGAACAAACCCAGGCTGGGGTTGGAGTAGCAGCGCTGGGTGCCAGGTTTGGCACGAGGTTGCCAGTGGCGGAAGAAATCGACTACCTGTTGCGGTGTTTGTACGGTGTCGGGGAACTGCAGt                                      <  4:255766/146‑1 (MQ=255)
                                                                                                           gctgcGCACGGGCGGCCAGGTCACCGAACAAACCCAGGCTGGGGTTGGAGTAGCAGCGCTGGGTGCCAGGTTTGGCACGAGGTTGCCAGTGGCGGAAGAAATCGACTACCTGTTGCGGTGTTTGTACGGTGTCGGGGAACTGCAGt                                      >  3:255766/1‑146 (MQ=255)
                                                                                                            ctgcGCACGGGCGGCCAGGTCACCGAACAAACCCAGGCTGGGGTTGGAGTAGCAGCGCTGGGTGCCAGGTTTGGCACGAGGTTGCCAGTGGCGGAAGAAATCGACTACCTGTTGCGGTGTTTGTACGGTGTCGGGGAACTGCAGTggca                                  <  4:127777/149‑1 (MQ=255)
                                                                                                                    gggCGGCCAGGTCACCGAACAAACCCAGGCTGGGGTTGGAGTAGCAGCGCTGGGTGCCAGGTTTGGCACGAGGTTGCCAGTGGCGGAAGAAATCGACTACCTGTTGCGGTGTTTGTACGGTGTCGGGGAACTGCAGTGGCAGGCAGTCg                          >  4:189065/1‑149 (MQ=255)
                                                                                                                              gTCACCGAACAAACCCAGGCTGGGGTTGGAGTAGCAGCGCTGGGTGCCAGGTTTGGCACGAGGTTGCCAGTGGCGGAAGAAATCGACTACCTGTTGCGGTGTTTGTACGGTGTCGGGGAACTGCAGTGGCAGGCAGTCGGCGCTATAgg                >  5:381561/1‑149 (MQ=255)
                                                                                                                                 aCCGAACAAACCCAGGCTGGGGTTGGAGTAGCAGCGCTGGGTGCCAGGTTTGGCACGAGGTTGCCAGTGGCGGAAGAAATCGACTACCTGTTGCGGTGTTTGTACGGTGTCGGGGAACTGCAGTGGCAGGCAGTCGGCGCTATAGgcgc             >  1:375540/1‑149 (MQ=255)
                                                                                                                                  ccGAACAAACCCAGGCTGGGGTTGGAGTAGCAGCGCTGGGTGCCAGGTTTGGCACGAGGTTGCCAGTGGCGGAAGAAATCGACTACCTGTTGCGGTGTTTGTACGGTGTCGGGGAACTGCAGTGGCAGGCAGTCGGCGCTATAGGCGcc            >  5:18140/1‑149 (MQ=255)
                                                                                                                                   cGAACAAACCCAGGCTGGGGTTGGAGTAGCAGCGCTGGGTGCCAGGTTTGGCACGAGGTTGCCAGTGGCGGAAGAAATCGACTACCTGTTGCGGTGTTTGTACGGTGTCGGGGAACTGCAGTGGCAGGCAGTCGGCGCTATAGGCGcc            >  8:48857/1‑148 (MQ=255)
                                                                                                                                     tccAACCCCAGCCT‑GGGTTGGAGTAGCAGCGCTGGGTGCCAGGTTTGGCACGAGGTTGCCAGTGGCGGAAGAAATCGACTACCTGTTGCGGTGTTTGTACGGTGTCg                                                  <  2:106039/105‑1 (MQ=255)
                                                                                                                                     tccAACCCCAGCCT‑GGGTTGGAGTAGCAGCGCTGGGTGCCAGGTTTGGCACGAGGTTGCCAGTGGCGGAAGAAATCGACTACCTGTTGCGGTGTTTGTACGGTGTCg                                                  >  1:106039/3‑107 (MQ=255)
                                                                                                                                     tccAACCCCAGCCT‑GGGTTGGAGTAGCAGCGCTGGGTGCCAGGTTTGGCACGAGGTTGCCAGTGGCGGAAGAAATCGACTACCTGTTGCGGTGTTTGTACGGTGTCGGGGAACTGCAGTGGCAGGCAGTCGGCGCta                    <  7:308829/135‑1 (MQ=255)
                                                                                                                                     tccAACCCCAGCCT‑GGGTTGGAGTAGCAGCGCTGGGTGCCAGGTTTGGCACGAGGTTGCCAGTGGCGGAAGAAATCGACTACCTGTTGCGGTGTTTGTACGGTGTCGGGGAACTGCAGTGGCAGGCAGTCGGCGCta                    >  8:308829/3‑137 (MQ=255)
                                                                                                                                            ccAGGCTGGGGTTGGAGTAGCAGCGCTGGGTGCCAGGTTTGGCACGAGGTTGCCAGTGGCGGAAGAAATCGACTACCTGTTGCGGTGTTTGTACGGTGTCGGGGAACTGCAGTGGCAGGCAGTCGGCGCTATAGGCGCCCAGTTCCAGc  <  6:253021/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CCCTTGGGCATAAAGCCCCTGGGCGCTGGCCGGAACCTGAAGGTAGGTGTTTTTCAGGCCCATCTTCGTCAGCAGGCCCTCACTCATCAGCGTGGAAAACGGTTGGTGCTGCGCACGGGCGGCCAGGTCACCGAACAAACCCAGGCTGGGGTTGGAGTAGCAGCGCTGGGTGCCAGGTTTGGCACGAGGTTGCCAGTGGCGGAAGAAATCGACTACCTGTTGCGGTGTTTGTACGGTGTCGGGGAACTGCAGTGGCAGGCAGTCGGCGCTATAGGCGCCCAGTTCCAGC  >  NC_002947/3277455‑3277743

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: ‑