Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,814,174 T→A 11.7% intergenic (+20/+44) emrB → / ← luxS multidrug efflux pump membrane subunit EmrB/S‑ribosylhomocysteine lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,814,1740TA11.7% 130.3 / 8.8 52intergenic (+20/+44)emrB/luxSmultidrug efflux pump membrane subunit EmrB/S‑ribosylhomocysteine lyase
Reads supporting (aligned to +/- strand):  ref base T (19/26);  new base A (2/4);  total (22/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

GATTATTTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATC  >  NC_000913/2814034‑2814313
                                                                                                                                            |                                                                                                                                           
gctgagttCGGCCAATGAGAGCCTCTGGATGTCAGCCGGGATAGTCCTCGCCCGGCTGGGGCTGGTGTGGGTTGCTAAACCGCCATTTGGCGCAGGTGGCGGTGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAttt                                                                                                                                           <  1:321694/137‑1 (MQ=255)
gaggtcttCCGCCCAGGGTATCTTCTGGATGTCAGCCGGGAGATGCCTCGTCCTGCGGGGGCTGGTGTGGGGTGCTACACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAttt                                                                                                                                           <  2:563965/137‑1 (MQ=255)
  ttattTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGa                                                                                                                                         <  1:130263/143‑1 (MQ=255)
   tattTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGaa                                                                                                                                        <  2:753841/143‑1 (MQ=255)
       tCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAAATGAACTgg                                                                                                                                    <  1:208320/143‑1 (MQ=255)
        ccGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAAATGAACTGGc                                                                                                                                   <  2:172237/143‑1 (MQ=255)
           ccAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAAATGAACTGGCtta                                                                                                                                <  1:61809/143‑2 (MQ=255)
            cAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCtttt                                                                                                                               >  1:766287/1‑143 (MQ=255)
                gagaTCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCa                                                                                                                           <  1:776872/143‑1 (MQ=255)
                gagaTCTTCTGGATGGCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCa                                                                                                                           <  1:27880/143‑1 (MQ=255)
                  gaTCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAAGTGGCTTTTTTCaat                                                                                                                         >  2:697053/1‑143 (MQ=255)
                      ttctGGCTGTCAGCCGGGGTATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTTTTCTGGCTTTTTTCaattaat                                                                                                                     <  2:497339/143‑1 (MQ=255)
                      ttctGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTACGTGCAACTAAGCCAGTTCAGTTGAACTGGCCTTTTTCaattcat                                                                                                                     >  2:783842/1‑143 (MQ=255)
                      ttctGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCaattaat                                                                                                                     <  1:621847/143‑1 (MQ=255)
                      ttctGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAAATGAACTGGCtt                                                                                                                                 >  1:33174/1‑131 (MQ=255)
                      ttctGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCAAATGAACTGGCtt                                                                                                                                 <  2:33174/131‑1 (MQ=255)
                        ctGGATGTCAGCCGGGATATTCCTCGCCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATtg                                                                                                                   <  1:396498/143‑1 (MQ=255)
                          gtctGTCTGCCGCGATACTCTTCGTCCTGCTCTGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGCGGCGGCGTATTCGGTGCGCACTAATTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATtgtg                                                                                                                 <  2:595280/140‑1 (MQ=255)
                          ggATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGGGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATtgtg                                                                                                                 <  1:323152/143‑1 (MQ=255)
                                   ccGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGtt                                                                                                        <  1:752121/143‑1 (MQ=255)
                                       gATATCCCTCGGCCTGCTGGGGCTGGTGTGGGTTGCTGAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACt                                                                                                    <  2:817191/143‑1 (MQ=255)
                                          atTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAATAGCCATGTGGCCGATGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGAc                                                                                                 <  2:575273/143‑1 (MQ=255)
                                           ttCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACGACTAAGGCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACt                                                                                                <  1:207058/143‑1 (MQ=255)
                                             ccTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGCTCATTTGACCTGGCTTTCTTCCATCAATTGTGAATATCGTTGAATGACTac                                                                                              >  2:273737/1‑142 (MQ=255)
                                                 gTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGAtgt                                                                                          >  1:759452/1‑143 (MQ=255)
                                                    ctgctgGGGTTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCa                                                                                       <  1:525821/143‑1 (MQ=255)
                                                      gctgGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGt                                                                                     >  1:689790/1‑143 (MQ=255)
                                                               ggTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCTGCGGAGGCTGTTCGCACTAAGTACAACTAAGCTAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAAc                                                                            <  1:250445/143‑1 (MQ=255)
                                                                   tgGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTtct                                                                        <  2:671845/143‑1 (MQ=255)
                                                                   ggattggCTAAACGGCCACTTGGGGAAGGTGGCGGCGGAGGCGGTGCGCACGAAGTACAATGAAGCCAGTTCATTTGAACGGGCTTTGTTCAATTAAGTGTGAAGCTAGGGTACTGACTAGACGTGCAGTTCCTGCAACTtct                                                                        <  2:22170/137‑1 (MQ=255)
                                                                      tttGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCtt                                                                     <  1:168646/143‑1 (MQ=255)
                                                                           tAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCCATTAATTGTGCAGCTAGTTTACTGACTAGATGTGCAGTTGCTGCGGCTTCTCTTTCGGc                                                                >  1:226674/1‑143 (MQ=255)
                                                                           tAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGc                                                                >  1:239131/1‑143 (MQ=255)
                                                                           tAAACCGCCATTTGGCGCAGGAGGCGGCGGAGGCGGTGCGCACTAAATACAACTAAGCCAAATCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGc                                                                >  1:665178/1‑143 (MQ=255)
                                                                              aCCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATGGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGt                                                             <  2:517881/143‑1 (MQ=255)
                                                                                cGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGc                                                           >  1:774416/1‑143 (MQ=255)
                                                                                  ccATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCa                                                         >  2:211637/1‑143 (MQ=255)
                                                                                         gcgcAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGttcttc                                                  <  1:779001/143‑1 (MQ=255)
                                                                                             cggTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTg                                              <  2:487156/142‑1 (MQ=255)
                                                                                                  gcggcgGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGGGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGtt                                         >  1:750967/1‑143 (MQ=255)
                                                                                                      cggAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATg                                     <  2:523456/143‑1 (MQ=255)
                                                                                                             ggTGCGCACTAAGTACAACTAAGCCAGATCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCGCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGGGTACGt                              >  2:734786/1‑143 (MQ=255)
                                                                                                                 cgcACTAAGTACAACTAAGCCAGGTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTCGATGTGCAGTACCTGCAACTTCTCTTTCGGCCGTGCCAGTGCGTCGTTGCTGTTGATACGTACGTCACg                          >  2:721883/1‑143 (MQ=255)
                                                                                                                     cTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTcc                      >  1:735663/1‑143 (MQ=255)
                                                                                                                      tAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATCGTGAAGAAAGATTACTGACTAGCTGTGCAGTTCCTGCAACTTCTCTTTCGGCAGCGCCAGTGCTTCGTTGCGGTTGATGCGTCCGTCACGCTCCa                     >  2:585876/1‑143 (MQ=255)
                                                                                                                         gTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGaa                  <  2:185656/143‑1 (MQ=255)
                                                                                                                             aaCTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCt              <  2:257990/143‑1 (MQ=255)
                                                                                                                              aCTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTCCGTTGCTGTTGATGCGTACGTCACGTTCCACGATGCtt             >  2:642984/1‑142 (MQ=255)
                                                                                                                                tAAGCGAGGTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTAcg           <  2:472942/143‑1 (MQ=255)
                                                                                                                                 aaGCCTGTTCATTTTAACTGGCTTTTTTCATTTAATTGTGAAGATAGTTTACTGACTAGTTGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTAcgc          <  2:215302/143‑1 (MQ=255)
                                                                                                                                    ccAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATCCGTACGTCCCGTTCCAGAAGCCCACgcgcc       >  2:104889/1‑142 (MQ=255)
                                                                                                                                    ccAGTTCAAATGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATCCTACgcgcc       >  2:29180/1‑142 (MQ=255)
                                                                                                                                        ttCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAatat   >  1:105393/1‑143 (MQ=255)
                                                                                                                                         tCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTAATGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTCCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATc  >  1:99601/1‑143 (MQ=255)
                                                                                                                                            |                                                                                                                                           
GATTATTTCCGCCAATGAGATCTTCTGGATGTCAGCCGGGATATTCCTCGTCCTGCTGGGGCTGGTGTGGTTTGCTAAACCGCCATTTGGCGCAGGTGGCGGCGGAGGCGGTGCGCACTAAGTACAACTAAGCCAGTTCATTTGAACTGGCTTTTTTCAATTAATTGTGAAGATAGTTTACTGACTAGATGTGCAGTTCCTGCAACTTCTCTTTCGGCAGTGCCAGTTCTTCGTTGCTGTTGATGCGTACGTCACGTTCCAGAATGCTACGCGCAATATC  >  NC_000913/2814034‑2814313

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