Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,526,622 T→A 13.5% I256N (ATT→AAT)  ddlB → D‑alanine‑‑D‑alanine ligase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,526,6220TA13.5% 95.7 / 8.6 37I256N (ATT→AAT) ddlBD‑alanine‑‑D‑alanine ligase B
Reads supporting (aligned to +/- strand):  ref base T (19/13);  new base A (2/3);  total (21/16)
Fisher's exact test for biased strand distribution p-value = 6.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.60e-01

GCCGCCGATCGCCCTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACT  >  NC_002947/1526473‑1526769
                                                                                                                                                      |                                                                                                                                                   
gccgccGATCGCCCTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCAt                                                                                                                                                     <  2:473854/150‑1 (MQ=255)
        tCGCCCTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCAGCCTTCGATGCCAATGGCATc                                                                                                                                              <  2:450347/150‑1 (MQ=255)
            ccTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCCCCCAGCCTTCGATGCCAATGGCATCGCAg                                                                                                                                          <  2:308647/150‑1 (MQ=255)
              tGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCt                                                                                                                                       >  2:247873/1‑150 (MQ=255)
                ggTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTgg                                                                                                                                     >  1:158841/1‑150 (MQ=255)
                     gccgcACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCg                                                                                                                                >  1:203199/1‑150 (MQ=255)
                            gtgtTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGAc                                                                                                                         <  1:247873/150‑1 (MQ=255)
                               ttCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCACGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTg                                                                                                                      >  2:56005/1‑150 (MQ=255)
                                   aCGACTATGACGCCAAGTACATCTTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGc                                                                                                                  >  2:346259/1‑150 (MQ=255)
                                       cTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGTGGTTGGACGTGATGCAGGa                                                                                                              <  2:472279/150‑1 (MQ=255)
                                        tATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGAc                                                                                                             <  2:292193/150‑1 (MQ=255)
                                                cAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAggg                                                                                                     >  2:493016/1‑150 (MQ=255)
                                                      cATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGtt                                                                                               >  2:236764/1‑150 (MQ=255)
                                                          gTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTgg                                                                                           >  2:579152/1‑150 (MQ=255)
                                                             aaCGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGCTCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGgctg                                                                                        >  2:426429/1‑150 (MQ=255)
                                                               cGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGAGGTGATGCAGGACGAGCAGGGCCGGTTCTGgctgct                                                                                      >  2:330449/1‑150 (MQ=255)
                                                                   aCCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGCCCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGaa                                                                                  <  1:485491/150‑1 (MQ=255)
                                                                     ccAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGt                                                                                >  2:410957/1‑150 (MQ=255)
                                                                          aCCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAaca                                                                           <  1:236764/150‑1 (MQ=255)
                                                                          aCCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAaca                                                                           <  1:579152/150‑1 (MQ=255)
                                                                            cGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACAcc                                                                         >  2:64927/1‑150 (MQ=255)
                                                                              cATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGc                                                                       >  2:555665/1‑150 (MQ=255)
                                                                                           cTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACCCCGCACCAGGCATGACa                                                          <  2:269522/150‑1 (MQ=255)
                                                                                               aTAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGTCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATc                                                      <  1:449250/150‑1 (MQ=255)
                                                                                                           aacaggaacTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGc                                          <  1:555665/150‑1 (MQ=255)
                                                                                                              aggaacTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCAt                                                              <  2:307156/127‑1 (MQ=255)
                                                                                                              aggaacTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCAt                                                              >  1:307156/1‑127 (MQ=255)
                                                                                                                         gACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGgc                                                                                          <  2:432320/88‑1 (MQ=255)
                                                                                                                         gACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGgc                                                                                          >  1:432320/1‑88 (MQ=255)
                                                                                                                         gACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGt                            >  1:585120/1‑150 (MQ=255)
                                                                                                                          aCCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGt                                                                                                                              >  1:564010/1‑51 (MQ=255)
                                                                                                                          aCCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGt                                                                                                                              <  2:564010/51‑1 (MQ=255)
                                                                                                                                      ccagcCTTCGATGCCAATGGCATCGCAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTg                >  2:599993/4‑150 (MQ=255)
                                                                                                                                                  gCCAATGGCATCGCAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCAt                                      >  2:176910/1‑116 (MQ=255)
                                                                                                                                                  gCCAATGGCATCGCAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCAt                                      <  1:176910/116‑1 (MQ=255)
                                                                                                                                                    cATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTgg                                                                                           <  2:16631/61‑1 (MQ=255)
                                                                                                                                                    cATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTgg                                                                                           >  1:16631/1‑61 (MQ=255)
                                                                                                                                                    cATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACt  >  2:403033/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
GCCGCCGATCGCCCTGGGTACGCCGCACGTGTTCTACGACTATGACGCCAAGTACATCGTCAACGATACCCAGTACCGCATCCCGTGCGGGCTGGATAGCGTCAAGGAACAGGAACTGATCGACCTGACCGCGCGCGCCTGCGATGCCATTGGCATCGAAGGCTGGGGGCGGTTGGACGTGATGCAGGACGAGCAGGGCCGGTTCTGGCTGCTCGAAGTCAACACCGCACCAGGCATGACAGATCATAGCCTGGTGCCCATGGCGGCCCGTGCGGCTGGCCTGGACTTCCAGCAACT  >  NC_002947/1526473‑1526769

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: