Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,277,593 A→C 9.8% G176G (GGT→GGG ampC ← beta‑lactamase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,277,5930AC 9.8% 106.6 / 5.1 41G176G (GGT→GGGampCbeta‑lactamase
Reads supporting (aligned to +/- strand):  ref base A (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 = 4.68e-01

GCGTCATACCCTTGGGCATAAAGCCCCTGGGCGCTGGCCGGAACCTGAAGGTAGGTGTTTTTCAGGCCCATCTTCGTCAGCAGGCCCTCACTCATCAGCGTGGAAAACGGTTGGTGCTGCGCACGGGCGGCCAGGTCACCGAACAAACCCAGGCTGGGGTTGGAGTAGCAGCGCTGGGTGCCAGGTTTGGCACGAGGTTGCCAGTGGCGGAAGAAATCGACTACCTGTTGCGGTGTTTGTACGGTGTCGGGGAACTGCAGTGGCAGGCAGTCGGCGCTATAGGCGCC  >  NC_002947/3277447‑3277733
                                                                                                                                                  |                                                                                                                                            
gCGTCATACCCTTGGGCATAAAGCCCCTGGGCGCTGGCCGGAACCTGAAGGTAGGTGTTTTTCAGGCCCATCTTCGTCAGCAGGCCCTCACTCATCAGCGTGGAAAACGGTTGGTGCTGCGCACGGGCGGCCAGGTCACCGAACAAAcc                                                                                                                                            <  5:362666/149‑1 (MQ=255)
        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)
                                                                                                                                             tccAACCCCAGCCTGGGTTGGAGTAGCAGCGCTGGGTGCCAGGTTTGGCACGAGGTTGCCAGTGGCGGAAGAAATCGACTACCTGTTGCGGTGTTTGTACGGTGTCg                                        <  2:106039/105‑1 (MQ=255)
                                                                                                                                             tccAACCCCAGCCTGGGTTGGAGTAGCAGCGCTGGGTGCCAGGTTTGGCACGAGGTTGCCAGTGGCGGAAGAAATCGACTACCTGTTGCGGTGTTTGTACGGTGTCg                                        >  1:106039/3‑107 (MQ=255)
                                                                                                                                             tccAACCCCAGCCTGGGTTGGAGTAGCAGCGCTGGGTGCCAGGTTTGGCACGAGGTTGCCAGTGGCGGAAGAAATCGACTACCTGTTGCGGTGTTTGTACGGTGTCGGGGAACTGCAGTGGCAGGCAGTCGGCGCta          <  7:308829/135‑1 (MQ=255)
                                                                                                                                             tccAACCCCAGCCTGGGTTGGAGTAGCAGCGCTGGGTGCCAGGTTTGGCACGAGGTTGCCAGTGGCGGAAGAAATCGACTACCTGTTGCGGTGTTTGTACGGTGTCGGGGAACTGCAGTGGCAGGCAGTCGGCGCta          >  8:308829/3‑137 (MQ=255)
                                                                                                                                                  |                                                                                                                                            
GCGTCATACCCTTGGGCATAAAGCCCCTGGGCGCTGGCCGGAACCTGAAGGTAGGTGTTTTTCAGGCCCATCTTCGTCAGCAGGCCCTCACTCATCAGCGTGGAAAACGGTTGGTGCTGCGCACGGGCGGCCAGGTCACCGAACAAACCCAGGCTGGGGTTGGAGTAGCAGCGCTGGGTGCCAGGTTTGGCACGAGGTTGCCAGTGGCGGAAGAAATCGACTACCTGTTGCGGTGTTTGTACGGTGTCGGGGAACTGCAGTGGCAGGCAGTCGGCGCTATAGGCGCC  >  NC_002947/3277447‑3277733

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: