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

TTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGT‑GGTGGGTGGGCGGCCTGGGG‑AAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGG‑CCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAG  >  NC_002947/1544346‑1544635
                                                                                                                                                  |                                                                                                                                                  
ttggcggtggggggggggtggggggggggGTTGGCGCGGTGGCGGGTGTGGTTGGT‑GGTGGGGGGGCGGCCTGGGG‑AAGGGGCTGTAGCTCGATAAGCGGTATTGCCGGGGCCGG‑CCTCTTCGCGGGTAAACCCGCTCCTACACAggg                                                                                                                                                <  5:71554/141‑1 (MQ=37)
     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)
                                                                                                                                                  |                                                                                                                                                  
TTTGCCCTGGCGGCGGGGGTGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGT‑GGTGGGTGGGCGGCCTGGGG‑AAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGG‑CCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAG  >  NC_002947/1544346‑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: ‑