Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,544,499 C→G 17.6% intergenic (+58/+56) mdfA → / ← ampG multidrug MFS transporter/muropeptide permease AmpG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,544,4990CG17.6% 89.7 / 14.7 40intergenic (+58/+56)mdfA/ampGmultidrug MFS transporter/muropeptide permease AmpG
Reads supporting (aligned to +/- strand):  ref base C (13/20);  new base G (4/3);  total (17/23)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.62e-02

CGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGC‑‑TTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGATCAGCG  >  NC_002947/1544359‑1544643
                                                                                                                                              |                                                                                                                                                   
gggggggggCGGGGGGGTTGGCGCGGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCACTCCTACACAGGGCCGCGCCGCTTTGt                                                                                                                                             <  2:113784/148‑1 (MQ=255)
   gggTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAg                                                                                                                                           >  5:111745/1‑148 (MQ=255)
       ggggggggCGGGGGGGCTGTGGCCGCGGTGGGTTGGGGGGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCg                                                                                                                                      <  6:275365/146‑1 (MQ=17)
        ggggggggGTGGGCGCGGGGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCgg                                                                                                                                     <  4:296820/147‑1 (MQ=21)
         cGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCggg                                                                                                                                    <  6:111745/149‑1 (MQ=255)
         cGGGGGCGCTGGCGCTGTGGCCGCGGGGGGTTGGGGTGGGGGGGCGGCCGGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCggg                                                                                                                                    <  4:256655/149‑1 (MQ=17)
          gggggCGTTGGCGCTGTGGCCGCTGTGGTTAGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTTAAACCCCTCCGAAACAGGGCCGCGCCGCGTTTTAGGAgacggt                                                                                                                                   >  1:226627/1‑144 (MQ=255)
                   ggCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCAGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAAACCCCGCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCcgc                                                                                                                          >  2:82116/1‑148 (MQ=255)
                       cTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCCCTCCTAAACAGGGCCCCGCCGCATTTTAGGAGCGGGGTTACCcgcgca                                                                                                                       >  6:368757/1‑145 (MQ=39)
                            gCCGCTGGGGTTTGTGGGGGGGGGGCGGCCTGGGGAAGGGGCGGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCa                                                                                                                <  3:217007/149‑1 (MQ=255)
                                                 tgggCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGt                                                                                                                                             >  6:116603/1‑100 (MQ=255)
                                                 tgggCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGt                                                                                                                                             <  5:116603/100‑1 (MQ=255)
                                                   gatccccTGGGCAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGC‑‑TTCAg                                                                                       <  5:341972/144‑1 (MQ=255)
                                                     gggcCTGGGGATGGGGCTTTCGCTCTATATTCGTTTTTGCCGGGGCGGGCCTCTTCGCGGTTAACCCCGCTCCTACACCGGGCCGCGCCGTTTTTTCGTAGGGGGTTAACCCGCGAAGAGGCCAGCACAGGCAATACAAAGC‑‑TTCAggg                                                                                     <  8:53405/148‑1 (MQ=255)
                                                         cTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTATTAGCGGGTTAAACCCCGGCGAAACCGGGCCCCGCCGCGTTTTAGGAGCGGGGTTTACCGCGAAGAGGGCCGCCCAGGGAATAAAAATC‑‑GTGAgggga                                                                                   >  6:341972/1‑145 (MQ=255)
                                                            gggAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTTAAACCCCTCCGACACACGGCCCCGCCGCGTTGTAGGAGCGGGGTTAACCGCGCAGAGGCCCGCACAGGCAATACAAAGC‑‑TTCAGGGCCGcgcc                                                                              >  2:51575/1‑146 (MQ=255)
                                                            gggAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGc                                                                                                                  <  5:389123/115‑1 (MQ=255)
                                                                aGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGC‑‑TTCAGGGCCGCTCAGcc                                                                           <  6:383073/148‑1 (MQ=255)
                                                                   ggCGGTAGCTCGATACTCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACCCAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGC‑‑TTCAGGGCCGCTCAGCCACCg                                                                       <  5:343388/149‑1 (MQ=255)
                                                                           cTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGC‑‑TTCAGGGCCGCTCAGc                                                                            <  7:158895/136‑1 (MQ=255)
                                                                              gATAATCGGTATTGCCGGGGCCGGCGTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGC‑‑TTCAGGGCCGCTCAGCCACCGAGTGGTCGTcc                                                            <  2:226627/149‑1 (MQ=255)
                                                                              gATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGC‑‑TTCAGGGCCGCTCAGCCACCGAGTGGTCGTcc                                                            <  1:51575/149‑1 (MQ=255)
                                                                                              ggggCCGGCCTCTGCGCGGGAAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGGGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGC‑‑TTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTccccc                                            <  3:363946/149‑1 (MQ=255)
                                                                                                  ccggccTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAg                                                                                               >  6:356738/1‑96 (MQ=255)
                                                                                                  ccggccTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAg                                                                                               <  5:356738/96‑1 (MQ=255)
                                                                                                    ggccCCGTCCCGGGTAAACCCGCCCCAACCCCCGGCCCGCCCGGCTTTGTGGGCGGGGTTTACCCCGCGAAGAGGCCAGCACAGGCAATACAAAGC‑‑TTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCcccgcccg                                       <  3:2356/149‑1 (MQ=255)
                                                                                                        tctTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGC‑‑TTCAGGGCCGCTCAGCCTCCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTc                                  >  5:307322/1‑149 (MQ=21)
                                                                                                          ttCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGC‑‑TTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCt                                >  3:102581/1‑149 (MQ=255)
                                                                                                             gcgGGAAATCCCGCTCCTACAGTGGGCCGCGCCGCTTTGTTGGAGGGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGC‑‑TTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCGGCCCGGCTTCCTGAt                             <  5:328900/149‑1 (MQ=38)
                                                                                                             cggggTAAACCCGGCCCAACCCCGGGCCGCCGCGGGTTTTTAGCGGGGGTTTAACCCGCGAAGAGGCCAGCACAGGCAATACAAAGC‑‑TTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGa                              <  1:343529/147‑1 (MQ=255)
                                                                                                               gggTTAACCCTCTCCTTCCCAGGGCCGCGCGGCTTTTTAGGTGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTTTCTGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGAt                             <  1:298512/149‑1 (MQ=37)
                                                                                                                 gTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGC‑‑TTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCa                         >  2:290345/1‑149 (MQ=255)
                                                                                                                      cccGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCCGGCAATACGAAGC‑‑TTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTg                      >  2:166311/1‑147 (MQ=21)
                                                                                                                             cTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGC‑‑TTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTc                                  >  2:158676/1‑128 (MQ=25)
                                                                                                                             cTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGC‑‑TTCAGGGCCGCTCAGCCACCGAGTGGACGACCGCCACCGGCGTccccc                                            >  2:146497/1‑118 (MQ=21)
                                                                                                                                 cacaGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGC‑‑TTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatc            <  1:82116/147‑1 (MQ=255)
                                                                                                                                 cacaGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGC‑‑TTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatc            <  1:290345/147‑1 (MQ=255)
                                                                                                                                  acaGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGC‑‑TTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatcag          >  3:147636/1‑148 (MQ=255)
                                                                                                                                    agcggcgcggcccTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGC‑‑TTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatcaggat       >  1:154398/6‑149 (MQ=255)
                                                                                                                                    aGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGC‑‑TTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatcagga        >  6:52130/1‑148 (MQ=255)
                                                                                                                                      gtccgcgacgcTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGC‑‑TTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatcaggat       >  2:276791/3‑147 (MQ=255)
                                                                                                                                      ggCCGCGCCGCATTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGC‑‑TTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatcaggatca     >  6:136131/1‑149 (MQ=255)
                                                                                                                                         cgcgccgcTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGC‑‑TTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGATCAGCg  <  5:52130/149‑1 (MQ=255)
                                                                                                                                          gcgccgcTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGC‑‑TTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCg               >  1:275156/1‑135 (MQ=255)
                                                                                                                                          gcgccgcTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGC‑‑TTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCg               <  2:275156/135‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
CGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGC‑‑TTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGATCAGCG  >  NC_002947/1544359‑1544643

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: