Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,352,966 T→C 7.1% N100S (AAC→AGC)  PP_2067 ← EmrB/QacA family drug resistance transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,352,9660TC 7.1% 160.2 / 2.5 56N100S (AAC→AGC) PP_2067EmrB/QacA family drug resistance transporter
Reads supporting (aligned to +/- strand):  ref base T (24/28);  new base C (2/2);  total (26/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.17e-01

CACGCTGTTGGCCAGGCTGAACAGCGCCATGCCCATGGCCATCTTGCTGGCCGGCAGTTCGGTGATGATCAGCTGGAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGGCGCCCATGAGCAGGCCACCGAGGCGATCAGAAAGATGAGCGAGCCGGTCCACATGACCCGGCGCAGCGAGAACACGTCCACCAACCAGGCCGTCAGCGGAATCATCACGATTTCCGCCACCAGGTAGGCGGTGGA  >  NC_002947/2352821‑2353105
                                                                                                                                                 |                                                                                                                                            
cACGCTGTTGGCCAGGCTGAACAGCGCCATGCCCATGGCCATCTTGCTGGCCGGCAGTTCGGCGATGATCAGCTGGAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGTTgg                                                                                                                                           <  3:34057/149‑1 (MQ=255)
  cGCTGTTGGCCAGGCTGAACAGCGCCATGCCCATGGCCATCTTGCTGGCCGGCAGTTCGGTGATGATCAGCTGGAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGgc                                                                                                                                         <  8:176080/149‑1 (MQ=255)
    cTGTTGGCCAGGCTGAACAGCGCCATGCCCATGGCCATCTTGCTGGCCGGCAGTTCGGTGATGATCAGCTGGAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCCCCCGGATTACGATCATCGCCTCCAGGTTGGgcgc                                                                                                                                       >  2:270755/1‑149 (MQ=255)
    cTGTTGGCCAGGCTGAACAGCGCCATGCCCATGGCCATCTTGCTGGCCGGCAGTTCGGTGATGATCAGCTGGAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGgcgc                                                                                                                                       <  6:130816/149‑1 (MQ=255)
       ttGGCCAGGCTGAACAGCGCCATGCCCATGGCCATCTTGCTGGCCGGCAGTTCGGTGATGATCAGCTGGAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGGCGccc                                                                                                                                     <  6:300441/148‑1 (MQ=255)
         ggCCAGGCTGAACAGCGCCATGCCCATGGCCATCTTGCTGGCCGGCAGTTCGGTGATGATCAGCTGGAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGGCGCCCAt                                                                                                                                   <  7:298072/148‑1 (MQ=255)
          gCCAGGCTGAACAGCGCCATGCCCATGGCCATCTTGCTGGCCGGCAGTTCGGTGATGATCAGCTGGAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGGCGCCCAt                                                                                                                                   >  1:295936/1‑147 (MQ=255)
          gCCAGGCTGAACAGCGCCATGCCCATGGCCATCTTGCTGGCCGGCAGTTCGGTGATGATCAGCTGGAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGGCGCCCATg                                                                                                                                  <  4:61108/148‑1 (MQ=255)
          gCCAGGCTGAACAGCGCCATGCCCATGGCCATCTTGCTGGCCGGCAGTTCGGTGATGATCAGCTGGAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGGCGCCCATGa                                                                                                                                 <  7:297283/149‑1 (MQ=255)
          gCCAGGCTGAACAGCGCCATGCCCATGGCCATCTTGCTGGCCGGCAGTTCGGTGATGATCAGCTGGAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGGCGCCCATGa                                                                                                                                 <  5:350750/149‑1 (MQ=255)
          cgcAGGCTGAACAGCGCCATGCCCATGGCCATCTTGCTGGCCGGCAGTTCGGTGATGATCAGCTGGAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGGCGCCCATGa                                                                                                                                 <  3:300696/147‑1 (MQ=255)
           ccAGGCTGAACAGCGCCATGCCCATGGCCATCTTGCTGGCCGGCAGTTCGGTGATGATCAGCTGGAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGGCGCCCATGAg                                                                                                                                >  2:361621/1‑149 (MQ=255)
              ggCTGAACAGCGCCATGCCCATGGCCATCTTGCTGGCCGGCAGTTCGGTGATGATCAGCTGGAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGCTGGGCGCCCATGAGCAg                                                                                                                             >  8:49230/1‑149 (MQ=255)
              ggCTGAACAGCGCCATGCCCATGGCCATCTTGCTGGCCGGCAGTTCGGTGATGATCAGCTGGAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGCTGGGCGCCCATGAGCAg                                                                                                                             >  8:43390/1‑149 (MQ=255)
               gCTGAACAGCGCCATGCCCATGGCCATCTTGCTGGCCGGCAGTTCGGTGATGATCAGCTGGAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGGCGCCCATGAGCAgg                                                                                                                            <  1:361621/149‑1 (MQ=255)
                        cgcCATGCCCATGGCCATCTTGCTGGCCGGCAGTTCGGTGATGATCAGCTGGAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGGCGCCCATGAGCAGGCCACCGAgg                                                                                                                   <  1:103723/149‑1 (MQ=255)
                              gCCCATGGCCATCTTGCTGGCCGGCAGTTCGGTGATGATCAGCTGGAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGGCGCCCAGGCGCAGGCCACCGAGGCGag                                                                                                               >  8:194576/1‑146 (MQ=255)
                               cccATGGCCATCTTGCTGGCCGGCAGTTCGGTGATGATCAGCTGGAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGGCGCCCATGAGCAGGCCACCGAGGCGATCAg                                                                                                            <  6:159294/149‑1 (MQ=255)
                                ccATGGCCATCTTGCTGGCCGGCAGTTCGGTGATGATCAGCTGGAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTACAGGTTGGCCGCCCATGAGAAGGCCACCGAGGCGATCAg                                                                                                            >  6:33967/1‑148 (MQ=255)
                                  aTGGCCATCTTGCTGGCCGGCAGTTCGGTGATGATCAGCTGGAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGCTGGGCGCCCATGAGCAGGCCACCGAGGCGATCAGaaa                                                                                                         <  7:49230/149‑1 (MQ=255)
                                  aTGGCCATCTTGCTGGCCGGCAGTTCGGTGATGATCAGCTGGAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGCTGGGCGCCCATGAGCAGGCCACCGAGGCGATCAGaaa                                                                                                         <  7:43390/149‑1 (MQ=255)
                                     gCCATCTTGCTGGCCGGCAGTTCGGTGATGATCAGCTGGAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGGCGCCCATGAGCAGGCCACCGAGGCGATCAGAAAGAt                                                                                                      <  5:88701/149‑1 (MQ=255)
                                      ccATCTTGCTGGCCGGCAGTTCGGTGATGATCAGCTGGAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGGCGCCCATGAGCAGGCCACCGAGGCGATCAGAAAGATg                                                                                                     <  2:295936/149‑1 (MQ=255)
                                       cATCTTGCTGGCCGGCAGTTCGGTGATGATCAGCTGGAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGGCGCCCATGAGCAGGCCACCGAGGCGATCAGAAAGAt                                                                                                      >  6:37886/1‑147 (MQ=255)
                                       cATCTTGCTGGCCGGCAGTTCGGTGATGATCAGCTGGAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGGCGCCCATGAGCAGGCCACCGAGGCGATCAGAAAGATga                                                                                                    <  5:124799/149‑1 (MQ=255)
                                        aTCTTGCTGGCCGGCAGTTCGGTGATGATCAGCTGGAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGGCGCCCATGAGCAGGCCACCGAGGCGATCAGAAAGATgag                                                                                                   <  5:51587/149‑1 (MQ=255)
                                        aTCTTGCTGGCCGGCAGTTCGGTGATGATCAGCTGGAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGGCGCCCATGAGCAGGCCACCGAGGCGATCAGAAAGATgag                                                                                                   <  7:117622/149‑1 (MQ=255)
                                             acTGGCCGGCAGTTCGGTGATGATCAGCTGGAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGGCGCCCATGAGCAGGCCACCGAGGCGATCAGAAAGATgagcgagc                                                                                              >  8:39121/2‑149 (MQ=255)
                                               tGGCCGGCAGTTCGGTGATGATCAGCTGGAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGGTCACCCGGATTACGATCATCGCCTCCAGGTTGGGCGCCCATGAGCAGGCCACCGAGGCGATCAGAAAGATGAGCGAGCCg                                                                                            >  3:269997/1‑149 (MQ=255)
                                                   cGGCAGTTCGGTGATGATCAGCTGGAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGGCGCCCATGAGCAGGCCACCGAGGCGATCAGAAAGATGAGCGAGCCgg                                                                                           <  1:286145/146‑1 (MQ=255)
                                                   cGGCAGTTCGGTGATGATCAGCTGGAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGGCGCCCATGAGCAGGCCACCGAGGCGATCAGAAAGATGAGCGAGCCgg                                                                                           >  2:286145/1‑146 (MQ=255)
                                                                  gatCAGCTGGAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGGCGCCCATGAGCAGGCCACCGAGGCGATCAGAAAGATGAGCGAGCCGGTCCACATGACCCGgcg                                                                          >  3:136177/1‑147 (MQ=255)
                                                                  gatCAGCTGGAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGGCGCCCATGAGCAGGCCACCGAGGCGATCAGAAAGATGAGCGAGCCGGTCCACATGACCCGgcg                                                                          >  4:286752/1‑147 (MQ=255)
                                                                       gCTGGAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGGCGCCCATGAGCAGGCCACCGAGGCGATCAGAAAGATGAGCGAGCCGGTCCACATGACCCGGCGCAGc                                                                      >  2:352957/1‑146 (MQ=255)
                                                                         tGGAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGGCGCCCATGAGCAGGCCACCGAGGCGATCAGAAAGATGAGCGAGCCGGTCCACATGACCCGGCGCAGCGAGaa                                                                 >  8:127042/1‑148 (MQ=255)
                                                                           gAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGGCGCCCATGAGCAGGCCACCGAGGCGATCAGAAAGATGAGCGAGCCGGTCCACATGACCCGGCGCAGCGAGAaca                                                               >  3:40605/1‑149 (MQ=255)
                                                                           gAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGGCGCCCATGAGCAGGCCACCGAGGCGATCAGAAAGATGAGCGAGCCGGTCCACATGACCCGGCGCAGCGAGAac                                                                <  7:39121/148‑1 (MQ=255)
                                                                               gACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGGCGCCCATGAGCAGGCCACCGAGGCGATCAGAAAGATGAGCGAGCCGGTCCACATGACCCGGCGCAGCGAGAACACGTc                                                           >  8:186862/1‑149 (MQ=255)
                                                                                 cAGGGGAATCAGCACAGCGCCCGCAGCCCCCTGGATCACCCGGAATACGATCATCGCCTCCAGGTTGGGCGCCCATGAGCAGGCCACCGAGGCGATCAGAAAGATGAGCGAGCCGGTCCACATGACCCGGCGCAGCGAGAACACGTcca                                                         >  2:6320/1‑149 (MQ=255)
                                                                                  aGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGGCGCCCATGAGCAGGCCACCGAGGCGATCAGAAAGATGAGCGAGCCGGTCCACATGACCCTGCGCAGCGAGAACACGTccac                                                        <  5:37886/149‑1 (MQ=255)
                                                                                    gggAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGGCGCCCATGAGCAGGCCACCGAGGCGATCAGAAAGATGAGCGAGCCGGTCCACATGACCCGGCGCAGCGAGAACACGTccacc                                                       >  8:222213/1‑148 (MQ=255)
                                                                                                cgcgCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGGCGCCCATGAGCAGGCCACCGAGGCGATCAGAAAGATGAGCGAGCCGGTCCACATGACCCGGCGCAGCGAGAACACGTCCACCAACCAGGCCGTCa                                          <  3:279941/149‑1 (MQ=255)
                                                                                                       gCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGGCGCCCATGAGCAGGCCACCGAGGCGATCAGAAAGATGAGCGAGCCGGTCCACATGACCCGGCGCAGCGAGAACACGTCCACCAACCAGGCCGTCAGCGGAat                                   >  1:138610/1‑149 (MQ=255)
                                                                                                            ccccTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGGCGCCCATGAGCAGCCCACCGAGGCGAGCAGAAAGATGAGCGCGCCGGTCCACATGCCCCGGCGCAGCGAGAACACGTCCACCAACCAGTCAGTACGCGGAAtcatca                              >  8:125665/1‑149 (MQ=255)
                                                                                                              ccTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGGCGCCCATGAGCAGGCCACCGTGGCGTTCATAAAGATGAGCGAGCCGGGCCACATGACCCGGGAGCAGCGAGAACACGTCCACCAACCAGGCCGTCAGCGGAATCATCAc                             <  1:270755/149‑1 (MQ=255)
                                                                                                              ccTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGGCGCCCATGAGCAGGCCACCGAGGCGATCAGAAAGATGAGCGAGCCGGTCCACATGACCCGGCGCAGCGAGAACACGTCCACCAACCAGGCCGTCAGCGGAATCATCAc                             >  1:146672/1‑148 (MQ=255)
                                                                                                                       cccGGATTACGATCATCGCCTCCAGGTTGGGCGCCCATGAGCAGGCCACCGAGGCGATCAGAAAGATGAGCGAGCCGGTCCACATGACCCGGCGCAGCGAGAACACGTCCACCACCCAGGCCGTCAGCGGAATCATCACGATTTCCGcc                   >  5:288242/1‑149 (MQ=255)
                                                                                                                       cccGGATTACGATCATCGCCTCCAGGTTGGGCGCCCATGAGCAGGCCACCGAGGCGATCAGAAAGATGAGCGAGCCGGTCCACATGACCCGGCGCAGCGAGAACACGTCCACCAACCAGGCCGTCAGCGGAATCATCACGATTTCCGc                    >  2:272197/1‑148 (MQ=255)
                                                                                                                       cccGGATTACGATCATCGCCTCCAGGTTGGGCGCCCATGAGCAGGCCACCGAGGCGATCAGAAAGATGAGCGAGCCGGTCCACATGACCCGGCGCAGCGAGAACACGTCCACCAACCAGGCCGTCAGCGGAATCATCACGATTTCCGc                    >  5:202681/1‑148 (MQ=255)
                                                                                                                         cGGATTACGATCTTCGCCTCCTGGTTGGGCGCCCATGAGCAGGCCACCGAGGCGATCAGAAAGATGAGCGTGACTGTCCACAAGACCCGGAGCAGCGAGAAGACGTCCTCCATCCAGGCCGTCAGAGGAATCATCACGATTTCCGccac                 <  5:33967/149‑1 (MQ=255)
                                                                                                                          ggATTACGATCATCGCCTCCAGGTTGGGCGCCCATGAGCAGGCCACCGAGGCGATCAGAAAGATGAGCGAGCCGGTCCACATGACCCGGCGCAGCGAGAACACGTCCACCAACCAGGCCGTCAGCGGAATCATCACGATTTCCGccacc                <  3:286752/149‑1 (MQ=255)
                                                                                                                            aTTACGATCATCGCCTCCAGGTTGGGCGCCCATGAGCAGGCCACCGAGGCGATCAGAAAGATGAGCGAGCCGGTCCACATGACCCGGCGCAGCGAGAACACGTCCACCAACCAGGCCGTCAGCGGAATCATCACGATTTCCGCCACCAg              <  1:108527/149‑1 (MQ=255)
                                                                                                                                   tcatcGCCTCCAGGTTGGGCGCCCATGAGCAGGCCACCGAGGCGATCAGAAAGATGAGCGAGCCGGTCCACATGACCCGGCGCAGCGAGAACACGTCCACCAACCAGGCCGTCAGCGGAATCATCACGATTTCCGCCACCAGGTAGGCg       >  7:182459/1‑149 (MQ=255)
                                                                                                                                      tcGCCTCCAGGTTGGGCGCCCATGAGCAGGCCACCGAGGCGATCAGAAAGATGAGCGAGCCGGTCCACATGACCCGGCGCAGCGAGAACACGTCCACCAACCAGGCCGTCAGCGGAATCATCACGATTTCCGCCACCAGGTAGGCGGTg    <  4:40605/149‑1 (MQ=255)
                                                                                                                                        gCCTCCAGGTTGGGCGCCCATGAGCAGGCCACCGAGGCGATCAGAAAGATGAGCGAGCCGGTCCACATGACCCGGCGCAGCGAGAACACGTCCACCAACCAGGCCGTCAGCGGAATCATCACGATTTCCGCCACCAGGTAGGCGGTGga  <  7:127042/149‑1 (MQ=255)
                                                                                                                                        gCCTCCAGGTTGGGCGCCCATGAGCAGGCCACCGAGGCGATCAGAAAGATGAGCGAGCCGGACCACATGACCCGGCGCAGCGAGAACACGTCCACCAACCAGGCCGTCAGCGGAATCATCACGATTTCCGCCACCAGGTAGGCGGTgg   <  7:222213/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
CACGCTGTTGGCCAGGCTGAACAGCGCCATGCCCATGGCCATCTTGCTGGCCGGCAGTTCGGTGATGATCAGCTGGAACGACAGGGGAATCAGCACCGCGCCCGCAGCCCCCTGGATCACCCGGATTACGATCATCGCCTCCAGGTTGGGCGCCCATGAGCAGGCCACCGAGGCGATCAGAAAGATGAGCGAGCCGGTCCACATGACCCGGCGCAGCGAGAACACGTCCACCAACCAGGCCGTCAGCGGAATCATCACGATTTCCGCCACCAGGTAGGCGGTGGA  >  NC_002947/2352821‑2353105

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