Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,584,487 A→T 8.6% L176Q (CTG→CAG)  PP_1388 ← EmrB/QacA family drug resistance transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,584,4870AT 8.6% 164.6 / 5.7 58L176Q (CTG→CAG) PP_1388EmrB/QacA family drug resistance transporter
Reads supporting (aligned to +/- strand):  ref base A (24/29);  new base T (3/2);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.29e-01

CGTGGCCCTGGCCCACCAAGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCGGGGGC  >  NC_002947/1584341‑1584630
                                                                                                                                                  |                                                                                                                                               
cGTGGCCCTGGCCCACCAAGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCTAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGc                                                                                                                                               <  4:194402/149‑1 (MQ=255)
       ctggcccACCAAGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCg                                                                                                                                        >  5:43396/1‑149 (MQ=255)
          gcccACCAAGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCa                                                                                                                                     <  8:223545/149‑1 (MQ=255)
          gcccACCAAGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCa                                                                                                                                     <  6:214482/149‑1 (MQ=255)
                 aaGGTGATGCCCAGCAGCATGATGCCAATGCCGATAATCAACATCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCTGCCGCTGGGGCTGGTTgct                                                                                                                              <  2:126714/149‑3 (MQ=255)
                 aaGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGAt                                                                                                                              <  1:145275/149‑1 (MQ=255)
                  aGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATc                                                                                                                             <  3:133673/149‑1 (MQ=255)
                  aGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATc                                                                                                                             <  3:128872/149‑1 (MQ=255)
                  aGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATc                                                                                                                             <  5:180672/149‑1 (MQ=255)
                    gTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGTCCCAGCGGCAGGTTGATCCa                                                                                                                           <  2:17208/149‑1 (MQ=255)
                     tGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCg                                                                                                                                        <  6:132661/135‑1 (MQ=255)
                     tGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAg                                                                                                                          <  5:305006/149‑1 (MQ=255)
                     tGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGAATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCg                                                                                                                                        >  5:132661/1‑135 (MQ=255)
                      gATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAg                                                                                                                          >  8:41200/1‑148 (MQ=255)
                       aTGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGaa                                                                                                                        <  3:299011/149‑1 (MQ=255)
                            cagcagcagGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACAcc                                                                                                                    >  5:285148/1‑148 (MQ=255)
                                agcagGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAgcg                                                                                                               >  3:199738/1‑149 (MQ=255)
                                   agGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCa                                                                                                            <  4:199738/149‑1 (MQ=255)
                                          cAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCa                                                                                                                                     <  7:59243/117‑1 (MQ=255)
                                          cAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCa                                                                                                                                     >  8:59243/1‑117 (MQ=255)
                                             ggCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTAc                                                                                                  >  7:124038/1‑149 (MQ=255)
                                                cGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCg                                                                                               >  2:288364/1‑149 (MQ=255)
                                                        aaCAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAAc                                                                                         >  1:151708/1‑147 (MQ=255)
                                                         aCAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCTGCCGCTGGGGCTGGTTgc                                                                                                                               >  8:15043/1‑107 (MQ=255)
                                                         aCAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCTGCCGCTGGGGCTGGTTgc                                                                                                                               <  7:15043/108‑2 (MQ=255)
                                                         aCAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAg                                                                                      <  8:124038/149‑1 (MQ=255)
                                                              aTGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAgccgcc                                                                                 <  7:41200/149‑1 (MQ=255)
                                                                   gcCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGCCAGGAACTCGGTCAACCAGCCGCCCAGCa                                                                            >  8:235532/1‑149 (MQ=255)
                                                                            gTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCg                                                                    <  6:285148/148‑1 (MQ=255)
                                                                               cACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGACCCGGCCa                                                                >  3:207736/1‑149 (MQ=255)
                                                                                  cTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAAc                                                                                         >  5:295977/1‑121 (MQ=255)
                                                                                  cTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAAc                                                                                         <  6:295977/121‑1 (MQ=255)
                                                                                    gtgCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTg                                                           >  1:312417/1‑149 (MQ=255)
                                                                                     tgCCTCCCGCCGCCGCACTGGCATCGCACCGAGCGCACGGCGAATGGCCCACACAGCAACCAGCCCCAGCGGCAGGGTGATCCAGAACACCCCCCGCCACGACCGGTACTCGGGCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTgg                                                          <  6:230976/149‑1 (MQ=255)
                                                                                            cgacgCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGAGCCAGAACACCCAGCGCCTCGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCgcggc                                                   <  1:303512/149‑1 (MQ=255)
                                                                                             gacgCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGc                                                                       <  5:173617/128‑1 (MQ=255)
                                                                                             gacgCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGc                                                                       >  6:173617/1‑128 (MQ=255)
                                                                                             gacgCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCgcggcg                                                  <  2:151708/149‑1 (MQ=255)
                                                                                             gacgCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCgcggcg                                                  >  2:252451/1‑149 (MQ=255)
                                                                                                gCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTa                                                >  2:20602/1‑148 (MQ=255)
                                                                                                    cACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACAtg                                            <  6:282062/148‑1 (MQ=255)
                                                                                                            tGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAg                                   <  1:288364/149‑1 (MQ=255)
                                                                                                                agcgagcgCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGc                               >  1:5147/1‑149 (MQ=255)
                                                                                                                 gcgagcgCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGcc                              >  7:60273/1‑149 (MQ=255)
                                                                                                                        cACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTAc                                               <  8:277727/125‑1 (MQ=255)
                                                                                                                        cACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTAc                                               >  7:277727/1‑125 (MQ=255)
                                                                                                                                 tGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGt              >  2:99601/1‑149 (MQ=255)
                                                                                                                                   gCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTc            >  2:302645/1‑149 (MQ=255)
                                                                                                                                    cccACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTcg           >  8:54078/1‑149 (MQ=255)
                                                                                                                                     ccACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTcgc          <  1:20602/149‑1 (MQ=255)
                                                                                                                                      cacaAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCg                                                                    <  1:254884/90‑1 (MQ=255)
                                                                                                                                      cacaAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCg                                                                    >  2:254884/1‑90 (MQ=255)
                                                                                                                                         ggatcAACCTGCCGCTGGGGCTGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCgg      >  7:179095/5‑149 (MQ=255)
                                                                                                                                          gatcAACCTGCCGCTGAGGCTGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCggg     >  8:101860/4‑149 (MQ=255)
                                                                                                                                          aaGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGaa                                    <  4:107593/118‑1 (MQ=255)
                                                                                                                                          aaGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGaa                                    >  3:107593/1‑118 (MQ=255)
                                                                                                                                           aGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCgggg    <  1:214118/149‑1 (MQ=255)
                                                                                                                                             cAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCGGGGGc  <  1:99601/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
CGTGGCCCTGGCCCACCAAGGTGATGCCCAGCAGCAGGCTGCCAAGGCCGAGAATCAACAGCATGGCGCCAAGGTAGTCCACCTGTGCCTCACGACGCTGCACTGGCATGTCAGCGAGCGCACGGCGAATGGCCCACAAAGCAACCAGCCCCAGCGGCAGGTTGATCCAGAACACCCAGCGCCACGACAGGTACTCGGTCAACCAGCCGCCCAGCACTGGCCCGGCCACGCTGGCCGCGGCGTACATGCTGCTGAAGTAGCCCTGGTAGCGGCCGCGTTCGCGCGGGGGC  >  NC_002947/1584341‑1584630

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