Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,560,033 C→T A299A (GCG→GCA ddpB ← putative D,D‑dipeptide ABC transporter membrane subunit DdpB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,560,0330CT92.9% 143.4 / 1.9 42A299A (GCG→GCAddpBputative D,D‑dipeptide ABC transporter membrane subunit DdpB
Reads supporting (aligned to +/- strand):  ref base C (3/0);  new base T (13/26);  total (16/26)
Fisher's exact test for biased strand distribution p-value = 4.88e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

AGCAAATGAAACCACGACGGCAAAGCCCATCACTGCCGGGAAGTCGAGCGCCTGTATTGATGTTACTACCCATGCACCCATTCCGGGCCAGGCAA  >  NC_000913/1559985‑1560079
                                                |                                              
aGCAAATGAAACCACGACGGCAAAGCCCATCACTGCCGGGAAGTCGAgc                                                <  1:2253124/49‑1 (MQ=255)
  cAAATGAAACCACGACGGCAAAGCCCATCACTGCCGGGAAGTCGAGTGcc                                             >  1:1614195/1‑50 (MQ=255)
    aaTGAAACCACGACGGCAAAGCCCATCACTGCCGGGAAGTCGAGTGCCTg                                           <  1:3667075/50‑1 (MQ=255)
     aTGAAACCACGACGGCAAAGCCCATCACTGCCGGGAAGTCGAGTGCCTGt                                          <  1:3212765/50‑1 (MQ=255)
     aTGAAACCACGACGGCAAAGCCCATCACTGCCGGGAAGTCGAGTGCCTGt                                          <  1:1044069/50‑1 (MQ=255)
      tGAAACCACGACGGCAAAGCCCATCACTGCCGGGAAGTCGAGTGCCTGTa                                         >  1:3312627/1‑50 (MQ=255)
      tGAAACCACGACGGCAAAGCCCATCACTGCCGGGAAGTCGAGTGCCTGTa                                         >  1:657341/1‑50 (MQ=255)
       gAAACCACGACGGCAAAGCCCATCACTGCCGGGAAGTCGAGTGCCTGTAt                                        >  1:919103/1‑50 (MQ=255)
        aaaCCACGACGGCAAAGCCCATCACTGCCGGGAAGTCGAGTGCCTGTAtt                                       <  1:107774/50‑1 (MQ=255)
        aaaCCACGACGGCAAAGCCCATCACTGCCGGGAAGTCGAGTGCCTGTAt                                        <  1:1416552/49‑1 (MQ=255)
         aaCCACGACGGCAAAGCCCATCACTGCCGGGAAGTCGAGTGCCTGTAtt                                       >  1:1485015/1‑49 (MQ=255)
          aCCACGACGGCAAAGCCCATCACTGCCGGGAAGTCGAGTGCCTGTATTGa                                     <  1:1637870/50‑1 (MQ=255)
            cACGACGGCAAAGCCCATCACTGCCGGGAAGTCGAGTGCCTGTATTGATg                                   >  1:3056313/1‑50 (MQ=255)
            cACGACGGCAAAGCCCATCACTGCCGGGAAGTCGAGTGCCTGTATTGATg                                   >  1:1553618/1‑50 (MQ=255)
              cgacgGCAAAGCCCATCACTGCCGGGAAGTCGAGTGCCTGTATTGATGtt                                 <  1:3833121/50‑1 (MQ=255)
               gacgGCAAAGCCCATCACTGCCGGGAAGTCGAGTGCCTGTATTGATGTta                                <  1:836496/50‑1 (MQ=255)
                       aGCCCATCACTGCCGGGAAGTCGAGTGCCTGTATTGATGTTACTACCCAt                        <  1:2510981/50‑1 (MQ=255)
                        gCCCATCACTGCCGGGAAGTCGAGTGCCTGTATTGATGTTACTACCCAtc                       <  1:276657/50‑2 (MQ=255)
                         cccATCACTGCCGGGAAGTCGAGTGCCTGTATTGATGTTACTACCCATGc                      <  1:2614594/50‑1 (MQ=255)
                           cATCACTGCCGGGAAGTCGAGTGCCTGTATTGATGCTACTACCCATGCAc                    >  1:1365260/1‑50 (MQ=255)
                            aTCACTGCCGGGAAGTCGAGTGCCTGTATTGATGTTACTACCCATGCAcc                   <  1:2266823/50‑1 (MQ=255)
                             tCACTGCCGGGAAGTCGAGTGCCTGTATTGATGTTACTACCCATGCAccc                  <  1:2386607/50‑1 (MQ=255)
                               aCTGCCGGGAAGTCGAGTGCCTGTATTGATGTTACTACCCATGCACCCa                 <  1:2944804/49‑1 (MQ=255)
                                cTGCCGGGAAGTCGAGTGCCTGTATTGATGTTACTACCCATGCACCCAtt               <  1:888842/50‑1 (MQ=255)
                                  gCCGGGAAGTCGAGTGCCTGTATTGATGTTACTACCCATGCACCCATTcc             >  1:474944/1‑50 (MQ=255)
                                   ccGGGAAGTCGAGTGCCTGTATTGATGTTACTACCCATGCACCCATTcc             <  1:1780349/49‑1 (MQ=255)
                                   ccGGGAAGTCGAGTGCCTGTATTGATGTTACTACCCATGCACCCATTc              <  1:3062783/48‑1 (MQ=255)
                                   ccGGGAAGTCGAGTGCCTGTATTGATGTTACTACCCATGCACCCATTCCg            <  1:3703418/50‑1 (MQ=255)
                                   ccGGGAAGTCGAGTGCCTGTATTGATGTTACTACCCATGCACCCATTCCg            <  1:2962088/50‑1 (MQ=255)
                                   ccGGGAAGTCGAGCGCCTGTATTGATGTTACTACCCATGCACCCATTCCg            >  1:580345/1‑50 (MQ=255)
                                    cGGGAAGTCGAGTGCCTGTATTGATGTTACTACCCATGCACCCATTCCgg           >  1:3076105/1‑50 (MQ=255)
                                       gAAGTCGAGTGCCTGTATTGATGTTACTACCCATGCACCCATTCCGGGcc        <  1:486301/50‑1 (MQ=255)
                                       gAAGTCGAGTGCCTGTATTGATGTTACTACCCATGCACCCATTCCGGGc         >  1:2388092/1‑49 (MQ=255)
                                        aaGTCGAGTGCCTGTATTGATGTTACTACCCATGCACCCATTCCGGGCCa       <  1:3957805/50‑1 (MQ=255)
                                          gTCGAGTGCCTGTATTGATGTTACTACCCATGCACCCATTCCGGGCCAgg     <  1:360356/50‑1 (MQ=255)
                                          gTCGAGTGCCTGTATTGATGTTACTACCCATGCACCCATTCCGGGCCAgg     <  1:2484399/50‑1 (MQ=255)
                                           tCGAGTGCCTGTATTGATGTTACTACCCATGCACCCATTCCGGGCCAGGc    <  1:316457/50‑1 (MQ=255)
                                            cGAGTGCCTGTATTGATGTTACTACCCATGCACCCATTCCGGGCCAGGc    >  1:1718711/1‑49 (MQ=255)
                                            cGAGTGCCTGTATTGATGTTACTACCCATGCACCCATTCCGGGCCAGGc    >  1:1078899/1‑49 (MQ=255)
                                            cGAGCGCCTGTATTGATGTTACTACCCATGCACCCATTCCGGGCCAGGc    >  1:998700/1‑49 (MQ=255)
                                             gAGTGCCTGTATTGATGTTACTACCCATGCACCCATTCCGGGCCAGGCaa  <  1:833821/50‑1 (MQ=255)
                                             gAGTGCCTGTATTGATGTTACTACCCATGCACCCATTCCGGGCCAGGCa   <  1:2981289/49‑1 (MQ=255)
                                              aGCGCCTGTATTGATGTTACTACCCATGCACCCATTCCGGGCCAGGCaa  >  1:1230665/1‑49 (MQ=255)
                                                |                                              
AGCAAATGAAACCACGACGGCAAAGCCCATCACTGCCGGGAAGTCGAGCGCCTGTATTGATGTTACTACCCATGCACCCATTCCGGGCCAGGCAA  >  NC_000913/1559985‑1560079

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