Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,329,743 G→T D261Y (GAT→TAT)  dacB → peptidoglycan DD‑endopeptidase DacB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,329,7430GT95.5% 55.9 / ‑4.0 22D261Y (GAT→TAT) dacBpeptidoglycan DD‑endopeptidase DacB
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (8/13);  total (9/13)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

GCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCAGATCCTGCGCCAGC  >  NC_000913/3329503‑3329989
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      
gCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGCAACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTACGCCGGTGCAATTCTGAAATATGAGTTaa                                                                                                                                                                                                                                               >  1:188658/1‑250 (MQ=255)
                                cgTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGGTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTGTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTGTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACActgc                                                                                                                                                                                                              <  1:160466/251‑1 (MQ=255)
                                               ccccGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAgg                                                                                                                                                                                               <  1:197353/251‑1 (MQ=255)
                                                cccGTTACGATGTTCAGTCAGGTACGCACCCTCCCCAGCGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTCAAACAGGCGGGTATCACCCGGGACGGGACACTGCTGCGCCCGACTCAggg                                                                                                                                                                                              >  1:129591/1‑250 (MQ=255)
                                                                         gCACCCTCCCCCGTGTTTCTGCCGAAGCGCAATACTGCGAACGGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTAGGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCTGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGcc                                                                                                                                                                     <  1:249316/251‑1 (MQ=255)
                                                                            cccTCCCCCGGGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGTTGCCAGGCGACCTTAACCGCTTTACGCTGACTTGATTCCTGCCACAACGTTCTGAGCCGGTCCCGTTGGCTTGTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGTAATTCTGAAATATGAGTTAAAGCAGGCGGGGATCACCGGGAGCGGAACACTGCTGCGTCAGACTCAGGTTATCGATACTGGAACGGTAGTTGCCAGt                                                                                                                                                                  <  1:304068/251‑1 (MQ=255)
                                                                                     gTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATTAGTTAAAACAGGCGGGTATCGCCTGGAGCGGAACACGGCTTCGCCAGACTCATGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCg                                                                                                                                                         <  1:234378/251‑1 (MQ=255)
                                                                                                                                 ggCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGGANCACNNNNNGcnn                                                                                                                                                                                                                      >  2:208778/1‑144 (MQ=255)
                                                                                                                                  gCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGCTGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGGCGGCCCCCCCGCACGATCTGCTTTAGATTATGCTGAAAAAGTCGGa                                                                                                             >  1:261674/1‑250 (MQ=255)
                                                                                                                                                           gggATGCCTGCCACACCGTTCTGCGCCGCTCACGTTCCCTTTTGCCGTGCAGGCTGGAGACCGCTATGACTGTGCAATTCTGAAAGATGAGTTAAA‑CCGGCGGCAATTCCCTGGAGCGGAACACTGCTGCGCCGGACTCAGGTTAACGAAACTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGGCCGATCTGCTTAAGATTATGCTGGAAAAGTCGGACAACCTGAtaagcgggtggaaggagtg                                                                                  >  1:290761/1‑233 (MQ=255)
                                                                                                                                                                 cctgccACAACGTTCTGAGCCGCTCCCGTGGGCTTTTGCCGTGCAGGAGGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAacac                                                                                                                                                                                                                 <  1:160942/119‑1 (MQ=255)
                                                                                                                                                                                                       ccGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCan                                                                                                                                                                   <  1:85562/127‑2 (MQ=255)
                                                                                                                                                                                                            cAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGGTAAc                                                                                                                                                                                          >  1:29189/1‑99 (MQ=255)
                                                                                                                                                                                                              ngATGGANNCAGCTATGCCGGTGCAATTCTGCAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAAc                                                                                                                                                                                          <  2:29189/96‑1 (MQ=255)
                                                                                                                                                                                                                  gnNGCCAGCNNTGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCt                                                                                                                                                                                    <  2:4955/99‑1 (MQ=255)
                                                                                                                                                                                                                     gCCAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACATGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGGCGGACGc                         >  1:279830/1‑251 (MQ=255)
                                                                                                                                                                                                                      ccAGCTATGCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCCATGTGCCTGGAACATGGCGGGCCGGGTCGGGCGcc                        >  1:191796/1‑251 (MQ=255)
                                                                                                                                                                                                                         gctagtccGGTGCAATTCTGAAATATGAGTTAAAACTGGCGGGTATCCCCTGGAGCTGAACACTGCTGCGCCAGACTCAGGTTAACGATCCTGGAACGGTAGTTGCAAGTAAACAGTAGGCCCCGCTGCACGAGCTGCTTAAGATTATGCGGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATTTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTg                     <  1:19668/245‑1 (MQ=255)
                                                                                                                                                                                                                              gCCGGTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTc                 <  1:98611/250‑1 (MQ=255)
                                                                                                                                                                                                                                 ggTGCAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACGGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACAGGACCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGTTCGGATGCCGTGCGTCAGAt             <  1:123282/251‑1 (MQ=255)
                                                                                                                                                                                                                                     cAATTCTGAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGGACACTGCTGCGCCAGAATCAGGGTAACGAACCTGGGACGGGNNTTGCCANNa                                                                                                                                                                 >  2:5499/1‑99 (MQ=255)
                                                                                                                                                                                                                                            gAAATATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAATAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCAGATCCTGCGCCAGc  <  1:5499/251‑1 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      
GCCCCAAAGCCTGGTGATATGGCTTTTATACGCGTGGCATCTTATTACCCCGTTACGATGTTCAGCCAGGTACGCACCCTCCCCCGTGGTTCTGCCGAAGCGCAATACTGCGAACTGGATGTGGTGCCAGGCGACCTGAACCGCTTTACGCTGACGGGATGCCTGCCACAACGTTCTGAGCCGCTCCCGTTGGCTTTTGCCGTGCAGGATGGAGCCAGCTATGCCGGTGCAATTCTGAAAGATGAGTTAAAACAGGCGGGTATCACCTGGAGCGGAACACTGCTGCGCCAGACTCAGGTTAACGAACCTGGAACGGTAGTTGCCAGTAAACAGTCGGCCCCGCTGCACGATCTGCTTAAGATTATGCTGAAAAAGTCGGACAACATGATCGCCGATACGGTTTTCCGCATGATAGGCCATGCGCGCTTCAATGTGCCTGGAACATGGCGGGCCGGGTCGGACGCCGTGCGTCAGATCCTGCGCCAGC  >  NC_000913/3329503‑3329989

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑