Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,544,489 C→A 23.8% 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,4890CA23.8% 39.2 / 9.2 21intergenic (+48/+66)mdfA/ampGmultidrug MFS transporter/muropeptide permease AmpG
Reads supporting (aligned to +/- strand):  ref base C (7/9);  new base A (3/2);  total (10/11)
Fisher's exact test for biased strand distribution p-value = 6.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.43e-01

TGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGA  >  NC_002947/1544365‑1544637
                                                                                                                              |                                                                                                                                                    
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)
                                                ggcCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGgccgcgccgc                                                                                                                                         <  3:232818/90‑1 (MQ=255)
                                                ggcCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCATCGCGGGTAAACACGCTCCTACACAGGgacgcgccgc                                                                                                                                         >  4:232818/1‑90 (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)
                                                                                                                              |                                                                                                                                                    
TGGCGGGGGCGTTGGCGCTGTGGCCGCTGTGGTTTGTGGTGGGTGGGCGGCCTGGGGAAGGGGCTGTAGCTCGATAATCGGTATTGCCGGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACACAGGGCCGCGCCGCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACAGGCAATACAAAGCTTCAGGGCCGCTCAGCCACCGAGTGGTCGTCCGCCACCGGCGTCCCCCGCCCGGCTTCCTGATGCCACTGCAGGGCGATCAGGA  >  NC_002947/1544365‑1544637

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: