Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,330,357 +T coding (1395/1434 nt) dacB → D‑alanyl‑D‑alanine carboxypeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,330,3531.T96.5% 87.1 / ‑5.0 29R464L (CGT→CTT) dacBD‑alanyl‑D‑alanine carboxypeptidase
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base T (14/14);  minor base A (1/0);  total (15/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGGTTCGTTGCAGGGGGTATATAACCTGGCGGGATTCATTACCACAGCGAGCGGGCAACGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTAT  >  NC_000913/3330215‑3330484
                                                                                                                                           |                                                                                                                                   
ccGGTTCGTTGCAGGGGGTATATAACCTGGCGGGATTCATTACCACAGCGAGCGGGCAACGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCg                                                                                                                                      <  1:559300/139‑1 (MQ=255)
              ggggTATATAACCTGGCGGGATTCATTACCACAGCGAGCGGGCAACGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGTGTAATCGCCGTATTCCGTTAGTGCGTTTTTGAAAGCCGt                                                                                                                        <  1:35456/139‑1 (MQ=255)
              ggggTATATAACCTGGCGGGATTCATTACCACAGCGAGCGGGCAACGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTTGAAAGCCGt                                                                                                                        <  1:519874/139‑1 (MQ=255)
              ggggTATATAACCTGGCGGGATTCATTACCACAGCGAGCGGGCAACGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTTGAAAGCCGt                                                                                                                        <  1:310128/139‑1 (MQ=255)
                 gTATATAACCTGGCGGGATTCATTACCACAGCGAGCGGGCAACGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTTGAAAGCCGTTTg                                                                                                                     <  2:570845/139‑1 (MQ=255)
                  tatataACCTGGCGGGATTCATTACCACAGCGAGCGGGCAACGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGtttttgaaa                                                                                                                            >  1:212271/1‑125 (MQ=255)
                  tatataACCTGGCGGGATTCATTACCACAGCGAGCGGGCAACGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGtttttgaaa                                                                                                                            <  2:212271/130‑6 (MQ=255)
                        aCCTGGCGGGATTCATTACCACAGCGAGCGGGCAACGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTTGAAAGCCGTTTGTATAAAg                                                                                                              >  1:574644/1‑139 (MQ=255)
                          cTGGCGGGATTCATTACCACAGCGAGCGGGCAACGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTTGAAAGCCGTTTGTATAAAGat                                                                                                            >  1:269949/1‑139 (MQ=255)
                          cTGGCGGGATTCATTACCACAGCGAGCGGGCAACGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTTGAAAGCCGTTTGTATAAAGat                                                                                                            >  2:221814/1‑139 (MQ=255)
                            ggCGGGATTCATTACCACAGCGAGCGGGCAACGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTTGAAAGCCGt                                                                                                                        <  2:220910/125‑1 (MQ=255)
                            ggCGGGATTCATTACCACAGCGAGCGGGCAACGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTTGAAAGCCGt                                                                                                                        >  1:220910/1‑125 (MQ=255)
                                             caGCGAGCGGGCAACGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTaat                                                                                         >  2:738852/1‑137 (MQ=255)
                                                 gagcgGGCAACGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCaaa                                                                                     >  1:735485/1‑139 (MQ=255)
                                                     gggCAACGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGa                                                                                 <  1:221814/139‑1 (MQ=255)
                                                     gggCAACGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGa                                                                                 >  1:146480/1‑139 (MQ=255)
                                                                       gtgCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCg                                                               >  2:222135/1‑139 (MQ=255)
                                                                              atCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTc                                                        <  2:269949/139‑1 (MQ=255)
                                                                                  ttCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTTGAAATCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGat                                                    >  1:670513/1‑139 (MQ=255)
                                                                                       gCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGGGCGATTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTg                                               >  1:331251/1‑139 (MQ=255)
                                                                                                  gAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAATAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTgcgc                                    <  2:574644/139‑1 (MQ=255)
                                                                                                       tGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGc                                                                 <  1:468205/105‑1 (MQ=255)
                                                                                                       tGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGc                                                                 >  2:468205/1‑105 (MQ=255)
                                                                                                       tGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGt                               >  1:53607/1‑139 (MQ=255)
                                                                                                        gCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTTGAAAGCCGTTTGTATAAAGATATTTAt                                                                                                      >  1:495731/1‑67 (MQ=255)
                                                                                                        gCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTTGAAAGCCGTTTGTATAAAGATATTTAt                                                                                                      <  2:495731/67‑1 (MQ=255)
                                                                                                        gCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGcac                                                                      >  1:93512/1‑99 (MQ=255)
                                                                                                        gCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGcac                                                                      <  2:93512/99‑1 (MQ=255)
                                                                                                                          cGTATTCCGTTAGTGCGTTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTa            <  1:222135/139‑1 (MQ=255)
                                                                                                                                 cGTTAGTGCGTTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCg                                                               >  1:425037/1‑81 (MQ=255)
                                                                                                                                 cGTTAGTGCGTTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCg                                                               <  2:425037/81‑1 (MQ=255)
                                                                                                                                    tAGTGCGTTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAATCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTAt  <  1:738852/139‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                   
CCGGTTCGTTGCAGGGGGTATATAACCTGGCGGGATTCATTACCACAGCGAGCGGGCAACGAATGGCGTTTGTGCAATATCTTTCTGGCTATGCAGTAGAACCTGCGGATCAGCGTAATCGCCGTATTCCGTTAGTGCGTTTTGAAAGCCGTTTGTATAAAGATATTTATCAGAACAATTAGTCAAAAAGAAACCCCGGCACATGGCCGGGGCTTCAGATTATTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTAAACGCCTTAT  >  NC_000913/3330215‑3330484

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: