Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,544,489 C→A 22.4% intergenic (+48/+66) 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,4890CA22.4% 61.4 / 16.5 36intergenic (+48/+66)mdfA/ampGmultidrug MFS transporter/muropeptide permease AmpG
Reads supporting (aligned to +/- strand):  ref base C (10/17);  new base A (5/3);  total (15/21)
Fisher's exact test for biased strand distribution p-value = 2.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.38e-01

CCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGC‑‑TTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAG  >  NC_002947/1544351‑1544635
                                                                                                                                           |                                                                                                                                                      
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                                                                                                                                     <  5:116603/100‑1 (MQ=255)
                                                         tgggCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGt                                                                                                                                     >  6:116603/1‑100 (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                                                                                       <  5:356738/96‑1 (MQ=255)
                                                                                                          ccggccTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAg                                                                                       >  6:356738/1‑96 (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)
                                                                                                                                           |                                                                                                                                                      
CCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGC‑‑TTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAG  >  NC_002947/1544351‑1544635

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: