Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,544,504 T→G 15.4% intergenic (+63/+51) 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,5040TG15.4% 60.8 / 6.4 26intergenic (+63/+51)mdfA/ampGmultidrug MFS transporter/muropeptide permease AmpG
Reads supporting (aligned to +/- strand):  ref base T (11/11);  new base G (2/2);  total (13/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.22e-01

TGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGATCAGCGTAGGCAC  >  NC_002947/1544365‑1544650
                                                                                                                                             |                                                                                                                                                  
gggcgggcgCGTTGGCGCGGGGGCCGCTGGGGTTTGTGGTGGGGGGGCGGCCTGGGGAAGGGGCTGCACCTCGATAAACGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGc                                                                                                                                           <  6:240938/148‑1 (MQ=21)
 ggggggggCGTTGGCGCTGTGGCCGGTGGGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCTGTGCTGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCg                                                                                                                                          <  3:65770/146‑1 (MQ=18)
      gggCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGttt                                                                                                                                     <  1:76121/149‑1 (MQ=255)
                    gggCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAACCCCGCGCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAATAGGc                                                                                                                       <  5:222221/148‑1 (MQ=255)
                                                 gcCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTAGCGGGTAAAACCGCTCCTAAAACGGGACGCGCCGCGTTGTATGAGCGGGTTTAACCGCGAATAGGCCAGAACCTGCAATAAAAAGCTTCAGTGcc                                                                                          >  1:203844/1‑149 (MQ=255)
                                                           ggggCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCGTTGTAGGAGCGGGTTTTCCCGCGATGAGGCCCGCACAGGCAATTCAAAGCTTCAGGGCCGGTCAGCCAc                                                                                >  5:313533/1‑149 (MQ=37)
                                                              gCTGTAGCTCGATAATCGGTATTCCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCaa                                                                                                           <  1:85548/119‑1 (MQ=255)
                                                                                     gCCGGGGCCGGCCTCTGCGCGGGTAAACCCGCTCCTCCCCAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTc                                                       <  2:203844/148‑1 (MQ=255)
                                                                                                ccTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGCGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCt                                           <  5:66042/149‑1 (MQ=255)
                                                                                                 ctctTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGcc                                                                >  1:11226/1‑127 (MQ=255)
                                                                                                 ctctTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGcc                                                                <  2:11226/127‑1 (MQ=255)
                                                                                                      cgcgGGTAAACCCCCACCGAAACACGGCCGCGCCGCTTTGTAGGCGCGGGTTTACCCGCGAAGAGGGCAGCCCACGCAATACAAACGGTAAGGGCCGCTCAGCCACCGAGTGGTCGTCCGGCACCCGCGTCCCCCGCCCCGCTTc                                         >  4:298540/1‑145 (MQ=255)
                                                                                                        cgGGTAAACCCGCTTCGCTACACAGGCCCGCGCCGCTTTGTAGGAGGGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGa                                     <  1:108680/149‑1 (MQ=38)
                                                                                                                cccGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTg                             >  4:169867/1‑147 (MQ=255)
                                                                                                                cccGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTg                             >  8:227047/1‑147 (MQ=255)
                                                                                                                cccGCTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTg                             >  8:20023/1‑147 (MQ=25)
                                                                                                                 ccGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATCGCCCTGCAg                          >  4:288301/1‑149 (MQ=255)
                                                                                                                    cTCCTACAAAGCGGCGCGGCCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCCGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGc                       >  5:202748/1‑149 (MQ=21)
                                                                                                                         cTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCg                      <  3:169867/147‑1 (MQ=255)
                                                                                                                            caccGGGCCGCGCCGCTTTGTAGGCGGGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatca                  <  7:218427/148‑1 (MQ=255)
                                                                                                                            cacaGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCCGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCcatcag                 >  8:179738/1‑149 (MQ=255)
                                                                                                                              caGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatcagga               >  8:274000/1‑149 (MQ=255)
                                                                                                                                gggCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatcaggat              >  5:71013/1‑148 (MQ=255)
                                                                                                                                   ccgcgccgcTTTGTAGGCGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGATCAGc          <  3:298540/149‑1 (MQ=255)
                                                                                                                                   ccgcgccgcTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATCGCCCTGCAGTGGcatcaggaagccg          <  3:288301/149‑6 (MQ=38)
                                                                                                                                           cTTTGTAGAAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGATCAGCGTAGGCAc  >  1:138836/1‑149 (MQ=255)
                                                                                                                                            tttGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatcaggatca            <  5:189901/138‑1 (MQ=255)
                                                                                                                                            tttGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCgatcaggatca            >  6:189901/1‑138 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
TGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGATCAGCGTAGGCAC  >  NC_002947/1544365‑1544650

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: