Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,544,494 C→G 20.1% 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,4940CG20.1% 55.8 / 11.7 30intergenic (+53/+61)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 = 2.59e-01

CCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGT‑GGTGGGTGGGCGGCCTGGGG‑AAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGG‑CCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGATCAGC  >  NC_002947/1544351‑1544642
                                                                                                                                                  |                                                                                                                                                    
ccTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGT‑GGTGGGTGGGCGGCCTGGGG‑AAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGG‑CCTCTTAGCGGGTAAACCCGCTCCTACACAGGgccggcc                                                                                                                                                 >  4:54330/1‑146 (MQ=255)
ccTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGT‑GGTGGGTGGGCGGCCTGGGG‑AAGGGGCTGGAGCTAGATAATCGTTATTGCCTGGGCCGG‑CCTATTATCGGGTAAACCCGCTCCTACAAAGGgccgcgc                                                                                                                                                 >  2:257657/1‑149 (MQ=255)
        gggggggggggggggCGGTGGCGGTGGGGGCGGTG‑GGTGTGTGGGGGGGTGGGCGGCCTGGGG‑AAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGG‑CCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTgtgt                                                                                                                                         <  4:208076/148‑1 (MQ=18)
         gggggTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGT‑GGTGGGTGGGCGGCCTGGGG‑AAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGG‑CCTATTCGCGGGTAAACCCCCTACGACAACGGGCCGCGCCGCGTTg                                                                                                                                          >  6:99623/1‑147 (MQ=255)
                ggtgggcCGTTGGCGCTGTGGCGGCTGTGGTTTGT‑GGTGGGTGGGCGGCCGGGGG‑AAGGGTCTGTAGCTCGATAATCGGTTTTGCCGGGGCCGG‑CCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCgg                                                                                                                                 <  5:133614/146‑1 (MQ=38)
                        gTTGGCGCTGTGGCCGCTGTGGTTTGT‑GGTGGGTGGGCGGCCTGGGG‑AAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGG‑CCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTAccc                                                                                                                         >  1:83120/1‑149 (MQ=255)
                         ttGGCGCTGTGGCCGCTGTGGTTTGT‑GGTGGGTGGGCGGCCTGGGG‑AAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGG‑CCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCcg                                                                                                                        <  2:83120/149‑1 (MQ=255)
                         ttGGCGCTGTGGCCGCTGTGGTTTGT‑GGTGGGTGGGCGGCCTGGGG‑AAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGG‑CCGCTTCGCGGGTAAAACCCCGCCGACACACGGCCGCGCCGCGTTGTTGGGGGGGGgtttgccg                                                                                                                        >  4:323491/1‑144 (MQ=255)
                                         tgtgGGTTGG‑GGGGGGGGGGCCGGCTGGGGCACGGGGCTGTCGCCTGAAAATGGGTA‑TGCCGGGGGCGGCCCTCTTCGCGGCTAACCCCGCCCAAACCCCGGGCCGCGCCGTTTTGTTGCCGGGGTTTAACCCGCGAAGAGGCCAGcac                                                                                                         <  4:202719/149‑1 (MQ=255)
                                             gTTTGT‑GGTGGGTGGGCGGCCTGGGG‑AAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGG‑CCTCTTCGCGGGTAAACCCGCACCGACACAGGGCCGCGCCGCGTTGTAGGAGCGGGGTTGCCCGCGCAGAGGGCAGccccggg                                                                                                     >  6:353596/1‑142 (MQ=37)
                                              tttGT‑GGTGGGTGGGCGGCCTGGGG‑AAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGG‑CCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAgg                                                                                                      <  2:267657/146‑1 (MQ=255)
                                              tttGT‑GGTGGGTGGGCGGCCTGGGG‑AAGGGGCTGTAGCTCGATAATCAGTATTGCCGGGGCCGG‑CCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAgg                                                                                                      >  1:267657/1‑146 (MQ=255)
                                                  t‑ggtggGTGGGCGGCCTGGGG‑AAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGG‑CCTCTTCGCGGGTAAAACCGCTGCGACACAGGGCCGCGCCGCGTTGTTGGGGCGGGGTTACCCGCGCAGAGGCCCGCACAGGCAATaa                                                                                                >  7:15195/1‑147 (MQ=38)
                                                             gcggcgTGGGT‑AAGGGGCTGTCGCTCGATAATCGGTATTGCCGGGGCCGG‑CCTCTTCGCGGGTAACCCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAgg                                                                                     <  3:54330/149‑1 (MQ=255)
                                                                                  tAGCTCGATAATCGGTATTGCCGGGGCCGG‑CCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTgg                                                                 <  5:353596/149‑1 (MQ=255)
                                                                                   aGCTCGATCATCGGTATTGCCGGGGCCGG‑CCTATTCGCGTGTAAACCCGCTCCTACACAGGGCCGAGCCGCTTTGTAGGAGCGGTTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTgg                                                                 <  1:257657/148‑1 (MQ=38)
                                                                                        gATAATCGGTATTCCCGGGGCCGG‑CCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGCGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTcc                                                           <  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                             >  4:182258/1‑136 (MQ=25)
                                                                                                                                    gCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGa                             <  3:182258/136‑1 (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)
                                                                                                                                                  |                                                                                                                                                    
CCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGT‑GGTGGGTGGGCGGCCTGGGG‑AAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGG‑CCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGATCAGC  >  NC_002947/1544351‑1544642

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: ‑