Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,726 A→G 7.0% L22P (CTC→CCC)  dppB ← dipeptide/heme ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,7260AG 7.0% 180.7 / 3.7 57L22P (CTC→CCC) dppBdipeptide/heme ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (24/29);  new base G (2/2);  total (26/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

CCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTATAGTGACCCGCAGCCTGGAACTGCGTGTCTGTATTGCTCACAAATCAATTCCC  >  NC_000913/3705578‑3705858
                                                                                                                                                    |                                                                                                                                    
ccAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTga                                                                                                                                      >  3:415546/1‑149 (MQ=255)
ccAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGACCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTga                                                                                                                                      >  1:137969/1‑149 (MQ=255)
       tAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGt                                                                                                                               <  4:219068/149‑1 (MQ=255)
          tGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAgg                                                                                                                                 >  1:259789/1‑144 (MQ=255)
          tGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAgg                                                                                                                                 <  2:259789/144‑1 (MQ=255)
          tGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAAt                                                                                                                            >  3:386488/1‑149 (MQ=255)
           ggAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATa                                                                                                                           >  1:12023/1‑149 (MQ=255)
               aTACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGa                                                                                                                       <  1:335065/149‑1 (MQ=255)
                    gctgcCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGGGAAGGGTAATAc                                                                                                                          >  1:220709/1‑141 (MQ=255)
                    gctgcCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGGGAAGGGTAATAc                                                                                                                          <  2:64695/141‑1 (MQ=255)
                    gctgcCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGGGAAGGGTAATAc                                                                                                                          <  2:220709/141‑1 (MQ=255)
                    gctgcCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGGGAAGGGTAATAc                                                                                                                          >  1:64695/1‑141 (MQ=255)
                       gcCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTg                                                                                                               <  4:386488/149‑1 (MQ=255)
                        ccACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTgg                                                                                                              <  4:461868/149‑1 (MQ=255)
                         ctcaTCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTggg                                                                                                             <  2:137969/147‑1 (MQ=255)
                          acaTCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGAt                                                                                                                      <  4:118121/139‑1 (MQ=255)
                          acaTCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGAt                                                                                                                      >  3:118121/1‑139 (MQ=255)
                          acaTCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTgggg                                                                                                            >  3:268269/1‑149 (MQ=255)
                           caTCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAggg                                                                                                                                <  1:162268/128‑1 (MQ=255)
                           caTCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAggg                                                                                                                                >  2:162268/1‑128 (MQ=255)
                            aTCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGAt                                                                                                          <  4:193352/149‑1 (MQ=255)
                            aTCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGAt                                                                                                          >  4:403040/1‑149 (MQ=255)
                                     tCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTcc                                                                                                 >  4:441768/1‑149 (MQ=255)
                                      cTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAgg                                                                                                                                 >  4:362009/1‑116 (MQ=255)
                                      cTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAgg                                                                                                                                 <  3:362009/116‑1 (MQ=255)
                                        anGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCaa                                                                                              >  3:110948/1‑149 (MQ=255)
                                            cGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGt                                                                                          <  4:4178/149‑1 (MQ=255)
                                             gAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTc                                                                                         <  1:366555/149‑1 (MQ=255)
                                              aGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCg                                                                                        <  4:268269/149‑1 (MQ=255)
                                               gTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCgg                                                                                       >  3:356932/1‑149 (MQ=255)
                                                      cagcagCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATaa                                                                                <  3:441768/149‑1 (MQ=255)
                                                        gcagcTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAAc                                                                              >  2:75100/1‑149 (MQ=255)
                                                           gcTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGga                                                                                      <  1:457820/138‑1 (MQ=255)
                                                           gcTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGga                                                                                      >  2:457820/1‑138 (MQ=255)
                                                               cgcgTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACa                                                                       <  1:60717/149‑1 (MQ=255)
                                                                     aCGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAGCTGCAACATAcccc                                                                 <  4:258803/149‑1 (MQ=255)
                                                                        ctctGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGa                                                              <  3:215155/149‑1 (MQ=255)
                                                                               gagaTCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAAAGCGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCta                                                       >  4:205291/1‑149 (MQ=255)
                                                                                   tCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTATAGt                                                   >  3:295715/1‑149 (MQ=255)
                                                                                    cccACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTATAGTg                                                  <  1:20626/149‑1 (MQ=255)
                                                                                     ccACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATaaa                                                                               >  3:339574/1‑119 (MQ=255)
                                                                                     ccACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATaaa                                                                               <  4:339574/119‑1 (MQ=255)
                                                                                          ttCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTATAGTGACCCGc                                            >  1:409383/1‑149 (MQ=255)
                                                                                               ccgccATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTATAGTGa                                                 >  1:469924/1‑139 (MQ=255)
                                                                                               ccgccATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTATAGTGa                                                 <  2:469924/139‑1 (MQ=255)
                                                                                               ccgccATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTATAGTGACCCGCAGCCt                                       <  2:409383/149‑1 (MQ=255)
                                                                                                       atcatcACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTATAGTGACCCGCAGCCTGGAACTGc                               <  4:415546/149‑1 (MQ=255)
                                                                                                         catcaCCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTATAGTGACCCGCAGCCTGGAACTGCgt                             <  2:125072/149‑1 (MQ=255)
                                                                                                                 gATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCtata                                                     >  4:197702/1‑117 (MQ=255)
                                                                                                                 gATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCtata                                                     <  3:197702/117‑1 (MQ=255)
                                                                                                                 gATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTATAGTGACCCGCAGCCTGGAACTGCGTGTCTGTAt                     <  4:356932/149‑1 (MQ=255)
                                                                                                                  aTCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTATAGTGACCCGCAGCCTGGAACTGCGTGTCTGTAtt                    >  4:376057/1‑149 (MQ=255)
                                                                                                                   tCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGAc                                                                                                       <  1:311654/65‑1 (MQ=255)
                                                                                                                   tCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGAc                                                                                                       >  2:311654/1‑65 (MQ=255)
                                                                                                                           gATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTATAGTGACCCGCAGCCTGGAACTGCGTGTCTGTATTGCTCACaaa           <  4:110948/149‑1 (MQ=255)
                                                                                                                                 gtgGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTATAGTGACCCGCAGCCTGGAACTGCGTGTCTGTATTGCTCACAAATCAAtt     >  3:197718/1‑149 (MQ=255)
                                                                                                                                    gACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTATAGTGACCCGCAGCCTGGAACTGCGTGTCTGTATTGCTCACAAATCAATTccc  >  4:171611/1‑149 (MQ=255)
                                                                                                                                      cAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTATAGTGACCCGCAGCCTGGAACTGCGTGTCTGTa                      <  4:164011/127‑1 (MQ=255)
                                                                                                                                      cAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTATAGTGACCCGCAGCCTGGAACTGCGTGTCTGTa                      >  3:164011/1‑127 (MQ=255)
                                                                                                                                                    |                                                                                                                                    
CCAAATGTAATGGAGATACTGCTGCCACATCGGTTTATCTAAGCCGAGTTCAGCCAGCAGCTGCGCGTGACGCTCTGGGGAGATCCCACGTTCGCCCGCCATGATCATCACCGGATCGCCCGGGATCATGTGGACAAAGGCAAATGTGAGAAGGGTAATACCGATAAACGTGGGGATGACGAGTCCCAAACGTCGGAGAATAAACTGCAACATAACCCGGATTCTCTATAGTGACCCGCAGCCTGGAACTGCGTGTCTGTATTGCTCACAAATCAATTCCC  >  NC_000913/3705578‑3705858

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑