Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,562,085 C→T A133T (GCA→ACA)  ddpA ← putative D,D‑dipeptide ABC transporter periplasmic binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,562,0850CT96.8% 105.2 / ‑5.7 31A133T (GCA→ACA) ddpAputative D,D‑dipeptide ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (18/12);  total (19/12)
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

ACCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAAAAGAAAG  >  NC_000913/1562039‑1562130
                                              |                                             
aCCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTc                                            <  1:1760998/50‑1 (MQ=255)
 ccGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTct                                           <  1:1491556/50‑1 (MQ=255)
  cGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctg                                          <  1:1089679/50‑1 (MQ=255)
   gTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctgc                                         >  1:1068785/1‑50 (MQ=255)
     aTGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctgctg                                       >  1:2057789/1‑50 (MQ=255)
       gTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGc                                     >  1:2044425/1‑50 (MQ=255)
        ttCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGcc                                    <  1:800950/50‑1 (MQ=255)
        ttCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGcc                                    <  1:2162359/50‑1 (MQ=255)
          cgtcgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGcccc                                  >  1:1088606/1‑50 (MQ=255)
          cgtcgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGcccc                                  >  1:818086/1‑50 (MQ=255)
           gtcgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGcccc                                  <  1:444039/49‑1 (MQ=255)
           gtcgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCt                                 <  1:1789658/50‑1 (MQ=255)
             cgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTgg                               >  1:209991/1‑50 (MQ=255)
              gggAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGc                              >  1:2342906/1‑50 (MQ=255)
               ggAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGcc                             >  1:1078772/1‑50 (MQ=255)
                gAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGcc                             >  1:2128003/1‑49 (MQ=255)
                   cATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAtt                         <  1:1116154/50‑1 (MQ=255)
                      cAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAtttt                       >  1:2324148/1‑49 (MQ=255)
                         tCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtag                    >  1:835188/1‑49 (MQ=255)
                         tCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtag                    >  1:1834256/1‑49 (MQ=255)
                           tttAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtagtag                 >  1:346096/1‑50 (MQ=255)
                           tttAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtagtag                 >  1:657789/1‑50 (MQ=255)
                            ttAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGc                >  1:1451547/1‑50 (MQ=255)
                             tAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtagt                   <  1:2132679/46‑1 (MQ=255)
                              aaaTCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCg              >  1:1057997/1‑50 (MQ=255)
                                aTCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCt            >  1:338796/1‑50 (MQ=255)
                                  cTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTc           >  1:2303371/1‑49 (MQ=255)
                                     tGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCaaa        <  1:328301/49‑1 (MQ=255)
                                     tGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCaa         <  1:2205465/48‑1 (MQ=255)
                                       ggAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaaga     <  1:2103561/50‑1 (MQ=255)
                                          aaTGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaagaaag  >  1:1092007/1‑50 (MQ=255)
                                              |                                             
ACCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAAAAGAAAG  >  NC_000913/1562039‑1562130

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