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,0850CT86.2% 539.1 / 72.0 189A133T (GCA→ACA) ddpAputative D,D‑dipeptide ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand):  ref base C (10/16);  new base T (80/83);  total (90/99)
Fisher's exact test for biased strand distribution p-value = 3.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.60e-01

Only 100 of 194 total aligned reads displayed.

CACCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAAAAGAAAGTTT  >  NC_000913/1562038‑1562133
                                               |                                                
cACCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGttt                                                >  1:2683932/1‑50 (MQ=255)
 aCCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGttt                                                <  2:3428400/49‑1 (MQ=255)
 aCCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTc                                               >  2:1482075/1‑50 (MQ=255)
 aCCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTc                                               >  2:287268/1‑50 (MQ=255)
  ccGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTct                                              >  2:3243375/1‑50 (MQ=255)
  ccGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTct                                              <  2:136360/50‑1 (MQ=255)
  ccGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTct                                              >  2:4927377/1‑50 (MQ=255)
  ccGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTct                                              <  1:4435797/50‑1 (MQ=255)
   cGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctg                                             <  1:5228820/50‑1 (MQ=255)
   cGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctg                                             <  2:2257050/50‑1 (MQ=255)
    gTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctgc                                            >  1:1051876/1‑50 (MQ=255)
    gTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTctgc                                            >  1:3829998/1‑50 (MQ=255)
     tATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctgct                                           <  2:2585998/50‑1 (MQ=255)
     tATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctgct                                           >  1:4161514/1‑50 (MQ=255)
     tATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctgct                                           <  2:4022194/50‑1 (MQ=255)
      aTGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctgctg                                          >  1:126079/1‑50 (MQ=255)
       tGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTctgctg                                          <  2:4613452/49‑1 (MQ=255)
       tGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTgg                                         <  1:3016152/50‑1 (MQ=255)
        gTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGc                                        >  1:4947993/1‑50 (MQ=255)
        gTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGc                                        >  1:1087931/1‑50 (MQ=255)
         ttCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGcc                                       <  2:2439158/50‑1 (MQ=255)
           cgtcgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGcccc                                     >  2:1282110/1‑50 (MQ=255)
            gtcgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCt                                    >  1:881214/1‑50 (MQ=255)
            gtcgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCt                                    >  2:4344835/1‑50 (MQ=255)
            gtcgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCt                                    <  1:5359785/50‑1 (MQ=255)
             tcgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTg                                   <  2:1142478/50‑1 (MQ=255)
             tcgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTg                                   <  2:1900045/50‑1 (MQ=255)
             tcgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTg                                   <  1:5034499/50‑1 (MQ=255)
              cgGGAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTgg                                  >  2:5608784/1‑50 (MQ=255)
               gggAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGc                                 <  2:2368248/50‑1 (MQ=255)
               gggAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGc                                 >  2:5333815/1‑50 (MQ=255)
                ggAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGcc                                <  2:2423127/50‑1 (MQ=255)
                ggAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGcc                                <  2:3074030/50‑1 (MQ=255)
                 gAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCg                               <  1:2615817/50‑1 (MQ=255)
                 gAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCg                               <  2:5155268/50‑1 (MQ=255)
                 gAGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCg                               <  2:671710/50‑1 (MQ=255)
                  aGCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGa                              <  2:3623927/50‑1 (MQ=255)
                   gCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAt                             >  1:1187588/1‑50 (MQ=255)
                   gCATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAt                             >  2:42659/1‑50 (MQ=255)
                    cATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAtt                            >  2:2709586/1‑50 (MQ=255)
                    cATCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAtt                            <  1:371644/50‑1 (MQ=255)
                    cATCAATCTTTAAATCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGAtt                            <  2:1626874/50‑1 (MQ=255)
                     aTCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAttt                           <  2:5233259/50‑1 (MQ=255)
                     aTCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAttt                           >  2:875931/1‑50 (MQ=255)
                      tCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAtttt                          >  2:57527/1‑50 (MQ=255)
                      tCAATCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAtttt                          <  2:4759742/50‑1 (MQ=255)
                        aaTCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAttttt                         <  2:1986667/49‑1 (MQ=255)
                        aaTCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGAtttt                          >  2:1594816/1‑48 (MQ=255)
                        aaTCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTta                        <  1:2078674/50‑1 (MQ=255)
                        aaTCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTta                        <  1:2759868/50‑1 (MQ=255)
                        aaTCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTta                        <  1:3464214/50‑1 (MQ=255)
                        aaTCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTta                        <  2:2738774/50‑1 (MQ=255)
                         aTCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtag                       >  2:1336073/1‑50 (MQ=255)
                          tCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtag                       >  1:3386220/1‑49 (MQ=255)
                          tCTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtag                       >  2:3229845/1‑49 (MQ=255)
                           cTTTAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtag                       >  1:3569649/1‑48 (MQ=255)
                           cTTTAAATCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTtag                       >  1:5521128/1‑48 (MQ=255)
                            tttAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtagtag                    >  1:1331159/1‑50 (MQ=255)
                            tttAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtagtag                    >  1:2282094/1‑50 (MQ=255)
                            tttAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtagtag                    >  1:3142987/1‑50 (MQ=255)
                              tAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtagtag                    <  2:439900/48‑1 (MQ=255)
                              tAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtagtag                    <  1:509745/48‑1 (MQ=255)
                              tAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtagta                     <  2:2803448/47‑1 (MQ=255)
                              tAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtagta                     <  2:3667295/47‑1 (MQ=255)
                              tAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTtagt                      <  1:619856/46‑1 (MQ=255)
                              tAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGcc                  >  1:2703732/1‑50 (MQ=255)
                              tAAATCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGc                   <  2:2833835/49‑1 (MQ=255)
                              tAAATCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTtagta                     <  1:3109145/47‑1 (MQ=255)
                               aaaTCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGcc                  <  2:4816498/49‑1 (MQ=255)
                               aaaTCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCg                 >  2:1822249/1‑50 (MQ=255)
                               aaaTCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCg                 >  1:5104379/1‑50 (MQ=255)
                                aaTCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCg                 <  1:4996694/49‑1 (MQ=255)
                                aaTCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGc                >  1:1534717/1‑50 (MQ=255)
                                 aTCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCt               >  2:3851613/1‑50 (MQ=255)
                                 aTCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCt               >  2:3304206/1‑50 (MQ=255)
                                 aTCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCt               <  2:2024256/50‑1 (MQ=255)
                                 aTCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCt               >  1:4670756/1‑50 (MQ=255)
                                  tCTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCt               >  2:5401772/1‑49 (MQ=255)
                                   cTTTGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTc              >  2:5098656/1‑49 (MQ=255)
                                     ttGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCaaa           <  1:796476/50‑1 (MQ=255)
                                     ttGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCaa            >  2:1552182/1‑49 (MQ=255)
                                      tGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTc              <  2:120347/46‑1 (MQ=255)
                                      tGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCaa            >  1:871229/1‑48 (MQ=255)
                                      tGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCa             <  2:4703174/47‑1 (MQ=255)
                                      tGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCa             <  1:1975421/47‑1 (MQ=255)
                                      tGGGAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCa             <  2:2135922/47‑1 (MQ=255)
                                       gggAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCaaaa          <  1:4053817/49‑1 (MQ=255)
                                       gggAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCaaaa          <  2:2535398/49‑1 (MQ=255)
                                       gggAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCaaaa          <  2:2164041/49‑1 (MQ=255)
                                       gggAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCaaaa          <  2:1465327/49‑1 (MQ=255)
                                       gggAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCaaaa          <  1:4835050/49‑1 (MQ=255)
                                       gggAAATGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCaa            >  1:3415599/1‑47 (MQ=255)
                                         gAAATGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaagaa       <  2:2946004/50‑1 (MQ=255)
                                         gAAATGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaaga        >  1:3285535/1‑49 (MQ=255)
                                          aaaTGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaagaaa      <  1:3190105/50‑1 (MQ=255)
                                           aaTGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaagaaag     >  1:3064307/1‑50 (MQ=255)
                                           aaTGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaagaaag     <  1:4540429/50‑1 (MQ=255)
                                           aaTGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAaaagaaag     >  2:1935058/1‑50 (MQ=255)
                                            aTGTTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAAAAGAAAGt    >  2:4581789/1‑50 (MQ=255)
                                               cTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAAAAGAAAGttt  <  2:3571153/49‑1 (MQ=255)
                                               |                                                
CACCGTATGTTCGTCGGGAGCATCAATCTTTAAATCTTTGGGAAATGCTTCTGCTGGCCCCTGGCCGATTTTTAGTAGCCGCTCAAAAGAAAGTTT  >  NC_000913/1562038‑1562133

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