Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,544,499 C→G 20.1% 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,4990CG20.1% 53.6 / 12.1 30intergenic (+58/+56)mdfA/ampGmultidrug MFS transporter/muropeptide permease AmpG
Reads supporting (aligned to +/- strand):  ref base C (12/12);  new base G (3/3);  total (15/15)
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.42e-01

CCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGATCAGCGTAG  >  NC_002947/1544351‑1544646
                                                                                                                                                       |                                                                                                                                                   
ccTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGGAGCTAGATAATCGTTATTGCCTGGGCCGGCCTATTATCGGGTAAACCCGCTCCTACAAAGGgccgcgc                                                                                                                                                     >  2:257657/1‑149 (MQ=255)
        gggggggggggggggCGGTGGCGGTGGGGGCGGTGGGTGTGTGGGGGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTgtgt                                                                                                                                             <  4:208076/148‑1 (MQ=18)
         gggggTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTATTCGCGGGTAAACCCCCTACGACAACGGGCCGCGCCGCGTTg                                                                                                                                              >  6:99623/1‑147 (MQ=255)
                ggtgggcCGTTGGCGCTGTGGCGGCTGTGGTTTGTGGTGGGTGGGCGGCCGGGGGAAGGGTCTGTAGCTCGATAATCGGTTTTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCgg                                                                                                                                     <  5:133614/146‑1 (MQ=38)
                        gTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTAccc                                                                                                                             >  1:83120/1‑149 (MQ=255)
                         ttGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCcg                                                                                                                            <  2:83120/149‑1 (MQ=255)
                         ttGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCGCTTCGCGGGTAAAACCCCGCCGACACACGGCCGCGCCGCGTTGTTGGGGGGGGgtttgccg                                                                                                                            >  4:323491/1‑144 (MQ=255)
                                         tgtgGGTTGGGGGGGGGGGGCCGGCTGGGGCACGGGGCTGTCGCCTGAAAATGGGTATGCCGGGGGCGGCCCTCTTCGCGGCTAACCCCGCCCAAACCCCGGGCCGCGCCGTTTTGTTGCCGGGGTTTAACCCGCGAAGAGGCCAGcac                                                                                                             <  4:202719/149‑1 (MQ=255)
                                             gTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCACCGACACAGGGCCGCGCCGCGTTGTAGGAGCGGGGTTGCCCGCGCAGAGGGCAGccccggg                                                                                                         >  6:353596/1‑142 (MQ=37)
                                              tttGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAgg                                                                                                          <  2:267657/146‑1 (MQ=255)
                                              tttGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCAGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAgg                                                                                                          >  1:267657/1‑146 (MQ=255)
                                                  tggtggGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAAACCGCTGCGACACAGGGCCGCGCCGCGTTGTTGGGGCGGGGTTACCCGCGCAGAGGCCCGCACAGGCAATaa                                                                                                    >  7:15195/1‑147 (MQ=38)
                                                             gcggcgTGGGTAAGGGGCTGTCGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAACCCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAgg                                                                                         <  3:54330/149‑1 (MQ=255)
                                                                                  tAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTgg                                                                     <  5:353596/149‑1 (MQ=255)
                                                                                   aGCTCGATCATCGGTATTGCCGGGGCCGGCCTATTCGCGTGTAAACCCGCTCCTACACAGGGCCGAGCCGCTTTGTAGGAGCGGTTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTgg                                                                     <  1:257657/148‑1 (MQ=38)
                                                                                        gATAATCGGTATTCCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGCGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTcc                                                               <  3:323491/149‑1 (MQ=255)
                                                                                                                     ttCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTc                                     <  7:23240/147‑1 (MQ=255)
                                                                                                                      tCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTAAGCCCCCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTg                                  >  7:185596/1‑149 (MQ=255)
                                                                                                                              ttaCCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCAc                           >  3:95703/3‑148 (MQ=255)
                                                                                                                                    gCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAggg                    >  7:264732/1‑149 (MQ=255)
                                                                                                                                    gCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGa                                 <  3:182258/136‑1 (MQ=25)
                                                                                                                                    gCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGa                                 >  4:182258/1‑136 (MQ=25)
                                                                                                                                      tCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAg                      >  4:110018/1‑145 (MQ=255)
                                                                                                                                      tCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAg                      <  3:110018/145‑1 (MQ=255)
                                                                                                                                        cTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCga                 <  5:99623/148‑1 (MQ=255)
                                                                                                                                        cTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgat                >  2:98178/1‑149 (MQ=37)
                                                                                                                                         tACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGAt                                >  1:213147/1‑132 (MQ=25)
                                                                                                                                         tACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGAt                                <  2:213147/132‑1 (MQ=25)
                                                                                                                                         aaaacaGGGCCGCGCCGGTTTGTTGGTGCGGGTTTACCCGCGAATAGGCCATCACAGGCAATAAAAAGCTTCAGGGCCGCTCAGCCACCGAGTTGGCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatc               <  8:15195/146‑1 (MQ=255)
                                                                                                                                          acacAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgat                >  3:291698/1‑147 (MQ=255)
                                                                                                                                           cacaGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGAAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatcag             >  8:203265/1‑149 (MQ=255)
                                                                                                                                                  ccgcgccgcTTTGTAGGAGCGGTTTAACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGATCAGc      <  3:26705/149‑1 (MQ=255)
                                                                                                                                                  ccgcgccgcTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGATCAGc      <  8:264732/149‑1 (MQ=255)
                                                                                                                                                      gccgcTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGGAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGATCAGCGTAg  >  1:260033/1‑149 (MQ=255)
                                                                                                                                                      gccgcTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGATCAGCGTAg  >  5:298084/1‑149 (MQ=255)
                                                                                                                                                       |                                                                                                                                                   
CCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGATCAGCGTAG  >  NC_002947/1544351‑1544646

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