Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,544,494 C→G 17.3% intergenic (+53/+61) 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,4940CG17.3% 86.0 / 12.2 40intergenic (+53/+61)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 = 3.26e-01

CCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGC‑‑TTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGATCA  >  NC_002947/1544351‑1544640
                                                                                                                                                |                                                                                                                                                      
ccTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCggcgcgg                                                                                                                                                  >  3:296820/1‑148 (MQ=255)
        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:290345/147‑1 (MQ=255)
                                                                                                                                         cacaGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGC‑‑TTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatc         <  1:82116/147‑1 (MQ=255)
                                                                                                                                          acaGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGC‑‑TTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatcag       >  3:147636/1‑148 (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)
                                                                                                                                                |                                                                                                                                                      
CCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGC‑‑TTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGATCA  >  NC_002947/1544351‑1544640

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: