Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,544,499 C→G 31.2% 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,4990CG31.2% 20.7 / 9.7 16intergenic (+58/+56)mdfA/ampGmultidrug MFS transporter/muropeptide permease AmpG
Reads supporting (aligned to +/- strand):  ref base C (6/5);  new base G (3/2);  total (9/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

GGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGATCAGCG  >  NC_002947/1544354‑1544643
                                                                                                                                                 |                                                                                                                                                
ggggggggggggggggggggCGGTGGCGCGGTGGCCGCTGGGGTTTGTGGTGGGGGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGgccgcgccgc                                                                                                                                               <  7:103568/143‑1 (MQ=38)
gggggggggggTGGCGGGGGCTTTGGGTGTTGGGCCGCTTTGTTTTGTGGTGGGGTGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGTTATTGCCGGGGCCGGCCTCTTCGCTGGTAAACCCGCTCCTACACAGGgccgcgccgc                                                                                                                                               <  1:382709/143‑1 (MQ=255)
       ggggggggggggGCGTTGGGGCTGTGGCCGCTGTGGTTTGTGGTGGTTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAg                                                                                                                                        <  8:214974/149‑1 (MQ=255)
                                             tgtgGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTa                                                                                                                             <  6:82648/122‑1 (MQ=255)
                                             tgtgGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTa                                                                                                                             >  5:82648/1‑122 (MQ=255)
                                               gtgggggtgggCGGCCTGGTGAAGGGCCTGTAGCTCGATAATCGGTATTGCCGTGGCCGGCCTCTTCGAAGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATAc                                                                                                 <  5:338989/144‑1 (MQ=37)
                                                                                                 ccGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGgccgcgccg                                                                                                                                                <  5:132334/51‑1 (MQ=255)
                                                                                                 ccGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGgccgcgccg                                                                                                                                                >  6:132334/1‑51 (MQ=255)
                                                                                                                          aCCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGc                     >  1:152025/1‑149 (MQ=255)
                                                                                                                           cccGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTg                      >  4:60465/1‑147 (MQ=255)
                                                                                                                           cccGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTg                      >  4:8911/1‑147 (MQ=25)
                                                                                                                            ccGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAg                   >  4:403730/1‑149 (MQ=255)
                                                                                                                            ccGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCg                                                     >  7:71265/1‑115 (MQ=21)
                                                                                                                            ccGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCg                                                     <  8:71265/115‑1 (MQ=21)
                                                                                                                             cGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAgg                  >  8:191567/1‑149 (MQ=255)
                                                                                                                                tCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAggg                 >  4:233391/1‑147 (MQ=37)
                                                                                                                                tCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAggg                 <  3:233391/147‑1 (MQ=37)
                                                                                                                                  cTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgat             >  8:104954/1‑149 (MQ=255)
                                                                                                                                       caGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatcagga        <  7:104954/149‑1 (MQ=255)
                                                                                                                                        aGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatcaggat       <  3:60465/149‑1 (MQ=255)
                                                                                                                                              gcgccgcTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCAACTGCAGGGCGAACAGCATCACCg  >  3:311915/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
GGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGATCAGCG  >  NC_002947/1544354‑1544643

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: