Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,557,194 A→C 12.2% L282R (CTG→CGG)  ddpF ← putative D,D‑dipeptide ABC transporter ATP‑binding subunit DdpF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,557,1940AC12.2% 117.9 / 11.0 55L282R (CTG→CGG) ddpFputative D,D‑dipeptide ABC transporter ATP‑binding subunit DdpF
Reads supporting (aligned to +/- strand):  ref base A (25/21);  new base C (4/3);  total (31/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.57e-02

AAGCATTGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTAT  >  NC_000913/1557056‑1557328
                                                                                                                                          |                                                                                                                                      
aaGCATTGACGGAGGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGGGAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGCGGGCTGCCAGccg                                                                                                                                    >  2:1032103/1‑141 (MQ=255)
     ttGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCTGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCCGCTGGCAAc                                                                                                                               <  2:952556/143‑1 (MQ=255)
      tGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCCGCTGGCAAcc                                                                                                                              <  2:308600/143‑2 (MQ=255)
        ccGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTAATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGACAGCGGGCAACGtt                                                                                                                            >  2:99612/1‑143 (MQ=255)
          ggATGCCGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCCGCTGGCAAcccacg                                                                                                                          <  1:1050061/143‑6 (MQ=255)
            aTGTCGCGTCATGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTa                                                                                                                        <  2:291642/143‑1 (MQ=255)
                 gcgTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAGCGTGCGTAGaaaa                                                                                                                   >  2:286054/1‑143 (MQ=255)
                 gcgTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCGCACGTCCGTCCTCCCTTAGTGCTAATGATTGCCCGACGTCACATCCGTGGGTTGCCAGCGGGCACCGCTCGTAGcaac                                                                                                                   >  1:422972/1‑142 (MQ=255)
                  cgTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGGTGACGTCCTTCCTCCCTTATTGCTAATGATTGCCGGACTTCCCATCCGTGGGTGGCCGCCGGGCAACGTTCGTAGaaaca                                                                                                                  >  1:943359/1‑141 (MQ=255)
                   gTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCACCTCTCCGCCAGCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGaaaaaa                                                                                                                 <  1:336054/143‑1 (MQ=255)
                      aGGCGCGTGAATGCCTGATGCGTTGTTAGCAGCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAg                                                                                                              <  1:979484/143‑1 (MQ=255)
                                      gATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGc                                                                                              >  1:683182/1‑143 (MQ=255)
                                        tGCGGTGTTAGCACCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGGCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCgg                                                                                            <  1:446970/143‑1 (MQ=255)
                                           gttgttAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGATGCCCGCCGGCACCGTTCGTCGACCAAACAGCCATCCGGCAAAGTCCGCttg                                                                                         >  2:330726/1‑142 (MQ=255)
                                             tgttAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAACGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTcc                                                                                       <  1:529818/143‑1 (MQ=255)
                                                tAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCgg                                                                                    <  1:930214/143‑1 (MQ=255)
                                                       tctACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCACCGTGCGCATTCCCGGCAGATcc                                                                             >  1:964663/1‑143 (MQ=255)
                                                         tACAGCGCCCGCCAGCACTGGAGCTCACTTCCGTCCTCCCGTGTGCCTAAGGATCGCGGGCCTTCACAGCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGt                                                                           <  1:981004/143‑1 (MQ=255)
                                                           ccacgcccgccAGCCGCTGAGCTCACGCCCTTCCTCTCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGttt                                                                         <  2:105062/140‑1 (MQ=255)
                                                           cAGGGCGCGCAAGCGCCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGttt                                                                         <  2:162239/143‑1 (MQ=255)
                                                           cAGCGCCCGCCAGCACCGGCCCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACGTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGACGAAGCAGCCTTGCGACAACGTCCGGTCTCCCGGCAGATCCGttt                                                                         >  2:1036178/1‑143 (MQ=255)
                                                             gcgcCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGCACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTa                                                                       <  2:221164/143‑1 (MQ=255)
                                                               gcccgccAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCAGCAGATCCGTTTTACg                                                                     >  1:128161/1‑143 (MQ=255)
                                                               gcccgccAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGGGCGGTTTCCCGGCAGATCCGTTCTACg                                                                     >  2:168212/1‑143 (MQ=255)
                                                                 ccgccAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTa                                                                   >  1:866096/1‑143 (MQ=255)
                                                                 ccgccAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTa                                                                   >  1:919124/1‑143 (MQ=255)
                                                                   gccAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTact                                                                 >  1:918738/1‑141 (MQ=255)
                                                                     cAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCCGATCCGTTTTACGTAATgg                                                               >  1:861519/1‑142 (MQ=255)
                                                                           cGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTCCGTAATGCCCAttc                                                         >  2:994256/1‑143 (MQ=255)
                                                                             gAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCCCATAAGTGGCATGTAAGCGAGCAACATTCGTAGGCAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGCTTTACGTAATGCCCAttctt                                                       >  1:539479/1‑143 (MQ=255)
                                                                              aGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCAttcttc                                                      <  2:755812/143‑1 (MQ=255)
                                                                              aGCTCACGTCCGTCCGCCCTTATGGCTAATGAGTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCAttcttc                                                      <  1:71871/143‑1 (MQ=255)
                                                                                  cACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCCGATCCGCTTTCCGTCCTGCCCCTTCTtccccc                                                  >  1:6266/1‑143 (MQ=255)
                                                                                   aCGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCa                                                 >  1:625210/1‑143 (MQ=255)
                                                                                             tccCTTATTGCTAATGATTGCCGGACTTCACATCTATCGGTTCCCAGCGGGCGACGTACGTAGAAAAAGCCGCCTTGCGGCAACGTGCGGTTTCTCGGCAGGTCCGCGTTACGTAAGGCGTCTTCTTCCTACGACGGTGTAtt                                       >  2:774108/1‑142 (MQ=255)
                                                                                                   aTTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCAATTATTCCTTGAGCGGTTTATCAATGGc                                 >  1:378006/1‑143 (MQ=255)
                                                                                                     tGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCgg                               <  1:1013959/143‑1 (MQ=255)
                                                                                                     tGCTAATGATTGCCGGACTCCACATCCGTGGTTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCgg                               <  1:336544/143‑1 (MQ=255)
                                                                                                     ggCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCgg                               <  2:707023/142‑1 (MQ=255)
                                                                                                        tactGATTGCCGGATTTCACATCCGTGGGTTGCCAGCGGGCAGCGTTCGTAGAAAAAACAGCCTTGCGGCAATGTGCGGTTTCCCGGCAGACCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGggga                            <  2:694672/140‑1 (MQ=255)
                                                                                                        tAATGATTGCCGGACTTCACATCCGTGTTTTCCCTACGGGTAACGTTCGTAGAAAAATCAGCGCTGCGACAACGAGCGCTTTCCCGGCAGATCCGTTTTACATAATGCTCAATATTCCTCCAGCGGTTTATCAATGGAgggtg                            >  2:473988/1‑141 (MQ=255)
                                                                                                        tAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGggga                            >  1:248396/1‑143 (MQ=255)
                                                                                                        tAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTGCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGCTTGTCGAGGGCggggg                            >  2:1038518/1‑142 (MQ=255)
                                                                                                         aagtaccgCCGGACTTCACACCCGGGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGgggag                           <  1:610373/136‑1 (MQ=255)
                                                                                                          aTGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTATTGCCTATTCTTCCTCCAGCGGTTTATCAATGGCGgggagg                          <  1:544648/143‑1 (MQ=255)
                                                                                                              ttGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAAt                                                                 >  2:537244/1‑100 (MQ=255)
                                                                                                              ttGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAAt                                                                 <  1:537244/100‑1 (MQ=255)
                                                                                                               tGCCGGACTTAACATCCGTGGGTGGCCAGCGGGTAACGAAAGAAAAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATc                     <  2:554253/143‑1 (MQ=255)
                                                                                                                gCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCAGCGACGGGCGGCTTCCCGGCAGATCCGTGTTACGTAATGCCCGCGAGTCCTCCAGCGGTGTGTGAACGGCGaggcggggagcc                    >  2:265872/1‑138 (MQ=255)
                                                                                                                         cacaTCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCaataat           >  1:587365/1‑143 (MQ=255)
                                                                                                                            aTCCGTGGGTTGCCAGCGGGCAAGGTTCGTAGAAAAAACAGCCTTGCGGGAACGTGCGGCTCCACCGCACCCTTGGCTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAAGGGCGGGGAGGGAATCCAGCAATAATCgg        <  2:571914/143‑1 (MQ=255)
                                                                                                                                 aacgTTGCCCGCTGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGgtgt  >  1:304126/4‑142 (MQ=255)
                                                                                                                                 aacgTTGCCCGCTGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTa   >  1:452378/4‑143 (MQ=255)
                                                                                                                                  gggTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGAGCGGTGTCCCGGCAAATCCGCTGTACGGAACGCCCCGTCGGCCTCCAGCGCTTTCTCAATGGCGCGGAGCGAAGCCAGCAATACTCCGgcgtcg  >  2:133300/1‑141 (MQ=255)
                                                                                                                                  acgTTGCCCGCTGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGCGGGAATCCAGCAATAATCGGGTGTAt  >  1:568744/3‑143 (MQ=255)
                                                                                                                                          |                                                                                                                                      
AAGCATTGCCGGATGACGCGTCAGGCGCGTGAATGCCTGATGCGTTGTTAGCATCTCTACAGCGCCCGCCAGCACCGGAGCTCACGTCCGTCCTCCCTTATTGCTAATGATTGCCGGACTTCACATCCGTGGGTTGCCAGCGGGCAACGTTCGTAGAAAAAACAGCCTTGCGGCAACGTGCGGTTTCCCGGCAGATCCGTTTTACGTAATGCCCATTCTTCCTCCAGCGGTTTATCAATGGCGGGGAGGGAATCCAGCAATAATCGGGTGTAT  >  NC_000913/1557056‑1557328

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