Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,237,467 G→A 9.7% Q373* (CAG→TAG)  PP_1975 ← EmrB/QacA family drug resistance transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,237,4670GA 9.7% 106.5 / 3.2 41Q373* (CAG→TAG) PP_1975EmrB/QacA family drug resistance transporter
Reads supporting (aligned to +/- strand):  ref base G (18/19);  new base A (2/2);  total (20/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

CAAGGCAATACCGATGGCCCCGCCCAGGTTACGCAGAATATTGAACAGGCTCGACGCGGACCCGGCATCCTGTTGTTGGATATACGCGGTGGCAATCAGCGAGATGGTCACCATGATCATCGGCTGGCCCAGGGCGCGGATGATCTGGATATGGTTGAACTGCGGCCCGGCAAAGTCCGGGTTGAGCACCCCCGAACCGAAACTGGCAGCGCCGAACAGGCAAAACCCCAGCGCGCATAGTACCTTGGGCGATATCACCTTCATCAGTTGCGGCACCAGCGGGATCAG  >  NC_002947/2237321‑2237608
                                                                                                                                                  |                                                                                                                                             
cAAGGCAATACCGATGGCCCCGCCCAGGTTACGCAGAATATTGAACAGGCTCGACGCGGACCCGGCATCCTGTTGTTGGATATACGCGGTGGCAATCAGCGAGATGGTCACCATGATCATCGGCTGGCCCAGGGCGCGGATGATCTGGa                                                                                                                                             <  8:181542/149‑1 (MQ=255)
cAAGGCAATACCGATGGCCCCGCCCAGGTTACGCAGAATATTGAACAGGCTCGACGCGGACCCGGCATCCTGTTGTTGGATATACGCGGTGGCAATCAGCGAGATGGTCACCATGATCATCGGCTGGCCCAGGGCGCGGATGATCTGGa                                                                                                                                             <  1:214854/149‑1 (MQ=255)
      aaTACCGATGGCCCCGCCCAGGTTACGCAGAATATTGAACAGGCTCGACGCGGACCCGGCATCCTGTTGTTGGATATACGCGGTGGCAATCAGCGAGATGGTCACCATGATCATCGGCTGGCCCAGGGCGCGGATGATCTGGATATGGt                                                                                                                                       >  2:287980/1‑149 (MQ=255)
      aaTACCGATGGCCCCGCCCAGGTTACGCAGAATATTGAACAGGCTCGACGCGGACCCGGCATCCTGTTGTTGGATATACGCGGTGGCAATCAGCGAGATGGTCACCATGATCATCGGCTGGCCCAGGGCGCGGATGATCTAGATATGGt                                                                                                                                       >  3:32760/1‑149 (MQ=255)
           cGATGGCCCCGCCCAGGTTACGCAGAATATTGAACAGGCTCGACGCGGACCCGGCATCCTGTTGTTGGATATACGCGGTGGCAATCAGCGAGATGGTCACCATGATCATCGGCTGGCCCAGGGCGCGGATGATCTGGATATGGtt                                                                                                                                      <  6:487/145‑1 (MQ=255)
           cGATGGCCCCGCCCAGGTTACGCAGAATATTGAACAGGCTCGACGCGGACCCGGCATCCTGTTGTTGGATATACGCGGTGGCAATCAGCGAGATGGTCACCATGATCATCGGCTGGCCCAGGGCGCGGATGATCTGGATATGGtt                                                                                                                                      >  5:487/1‑145 (MQ=255)
           cGATGGCCCCGCCCAGGTTACGCAGAATATTGAACAGGCTCGACGCGGACCCGGCATCCTGTTGTTGGATATACGCGCTGGCAATCAGCGAGATGGTCACCATGATCATCGGCTGGCCCAGGGCGCGGATGATCTGGATATGGTTGAAc                                                                                                                                  >  1:97670/1‑149 (MQ=255)
            gATGGCCCCGCCCAGGTTACGCAGAATATTGAACAGGCTCGACGCGGACCCGGCATCCTGTTGTTGGATATACGCGGTGGCAATCAGCGAGATGGTCACCATGATCATCGGCTGGCCCAGGGCGCGGATGATCTAGATATGGTTGAACt                                                                                                                                 <  4:32760/149‑1 (MQ=255)
             aTGGCCCCGCCCAGGTTACGCAGAATATTGAACAGGCTCGACGCGGACCCGGCATCCTGTTGTTGGATATACGCGGTGGCAATCAGCGAGATGGTCACCATGATCATCGGCTGGCCCAGGGCGCGGATGATCTGGATATGGTTGAACTg                                                                                                                                >  6:122411/1‑149 (MQ=255)
              tGGCCCCGCCCAGGTTACGCAGAATATTGAACAGGCTCGACGCGGACCCGGCATCCTGTTGTTGGATATACGCGGTGGCAATCAGCGAGATGGTCACCATGATCATCGGCTGGCCCAGGGCGCGGATGATCTGGATATGGTTGAACTGc                                                                                                                               <  7:279441/149‑1 (MQ=255)
                  cccGCCCAGGTTACGCAGAATATTGAACAGGCTCGACGCGGACCCGGCATCCTGTTGTTGGATATACGCGGTGGCAATCAGCGAGATGGTCACCATGATCATCGGCTGGCCCAGGGCGCGGATGATCTGGATATGGTTGAACTGCGGcc                                                                                                                           >  8:181394/1‑149 (MQ=255)
                    cGCCCAGGTTACGCAGAATATTGAACAGGCTCGACGCGGACCCGGCATCCTGTTGTTGGATATACGCGGTGGCAATCAGCGAGATGGTCACCATGATCATCGGCTGGCCCAGGGCGCGGATGATCTGGATATGGTTGAACTGCGGCCCg                                                                                                                         <  6:308977/149‑1 (MQ=255)
                               cGCAGAATATTGAACAGGCTCGACGCGGACCCGGCATCCTGTTGTTGGATATACGCGGTGGCAATCAGCGAGATGGTCACCATGATCATCGGCTGGCCCAGGGCGCGGATGATCTGGATATGGTTGAACTGCGGCCCGGCAAAGTCCgg                                                                                                              <  5:122411/149‑1 (MQ=255)
                               cGCAGAATATTGAACAGGCTCGACGCGGACCCGGCATCCTGTTGTTGAATATACGCGGTGGAAATCAGCGAGATTGTAACCATGATCATCGGCTGGCCCAGGGCGCGGATGATCTGGATATGGTTGAACTGCGGCCCGGCAAAGTCCg                                                                                                               <  5:53133/148‑1 (MQ=255)
                                      tatTGAACAGGCTCGACGCGGACCCGGCATCCTGTTGTTGGATATACGCGGTGGCAATCAGCGAGATGGTCACCATGATCATCGGCTGGCCCAGGGCGCGGATGATCTGGATATGGTTGAACTGCGGCCCGGCAAAGTCCg                                                                                                               >  8:39259/1‑141 (MQ=255)
                                      tatTGAACAGGCTCGACGCGGACCCGGCATCCTGTTGTTGGATATACGCGGTGGCAATCAGCGAGATGGTCACCATGATCATCGGCTGGCCCAGGGCGCGGATGATCTGGATATGGTTGAACTGCGGCCCGGCAAAGTCCg                                                                                                               <  7:39259/141‑1 (MQ=255)
                                            aCAGGCTCGACGCGGACCCGGCATCCTGTTGTTGGATATACGCGGTGGCAATCAGCGAGATGGTCACCATGATCATCGGCTGGCCCAGGGCGCGGATGATCTGGATATGGTTGAACTGCGGCCCGGCAAAGTCCGGGTTGAGCAccccc                                                                                                 >  3:245946/1‑149 (MQ=255)
                                                gCTCGACGCGGACCCGGCATCCTGTTGTTGGATATACGCGGTGGCAATCAGCGAGATGGTCACCATGATCATCGGCTGGCCCAGGGCGCGGATGATCTGGATATGGTTGAACTGCGGCCCGGCAAAGTCCGGGTTGAGCACCCccgaac                                                                                             <  2:298557/149‑1 (MQ=255)
                                                  tCGACGCGGACCCGGCATCCTGTTGTTGGATATACGCGGTGGCAATCAGCGAGATGGTCACCATGATCATCGGCTGGCCCAGGGCGCGGATGATCTGGATATGGTTGAACTGCGGCCCGGCAAAGTCCGGGTTGAGCACCCccgaacc                                                                                            >  2:205799/1‑148 (MQ=255)
                                                          gACCCGGCATCCTGTTGTTGGATATACGCGGTGGCAATCAGCGAGATGGTCACCATGATCATCGGCTGGCCCAGGGCGCGGATGATCTGGATATGGTTGAACTGCGGCCCGGCAAAGTCCGGGTTGAGCACCCCCGAACCGAAACTGGc                                                                                   <  1:287980/149‑1 (MQ=255)
                                                                      tgttgtTGGATATACGCGGTGGCAATCAGCGAGATGGTCACCATGATCATCGGCTGGCCCAGGGCGCGGATGATCTGGATATGGTTGAACTGCGGCCCGGCAAAGTCCGGGTTGAGCACCCCCGAACCGAAACTGGCAGCGCCGAACAg                                                                       >  1:287049/1‑149 (MQ=255)
                                                                           ttgGATATACGCGGTGGCAATCAGCGAGATGGTCACCATGATCATCGGCTGGCCCAGGGCGCGGATGATCTGGATATGGTTGAACTGCGGCCCGGCAAAGTCCGGGTTGAGCACCCCCGAACCGAAACTGGCAGCGCCGAACAGGCaa                                                                   >  6:125996/1‑148 (MQ=255)
                                                                              gATATACGCGGTGGCAATCAGCGAGATGGTCACCATGATCATCGGCTGGCCCAGGGCGCGGATGATCTGGATATGGTTGAACTGCGGCCCGGCAAAGTCCGGGTTGAGCACCCCCGAACCGAAACTGGCAGCGCCGAACAGGCAAAAc                                                                <  2:97670/148‑1 (MQ=255)
                                                                                        gTGGCAATCAGCGAGATGGTCACCATGATCATCGGCTGGCCCAGGGCGCGGATGATCTGGATATGGTTGAACTGCGGCCCGGCAAAGTCCGGGTTGAGCACCCCCGAACCGAAACTGGCAGCGCCGAACAGGCAAAACCCCAGCGCGCa                                                     >  3:218360/1‑149 (MQ=255)
                                                                                         tGGCAATCAGCGAGATGGTCACCATGATCATCGGCTGGCCCAGGGCGCGGATGATCTGGATATGGTTGAACTGCGGCCCGGCAAAGTCCGGGTTGAGCACCCCCGAACCGAAACTGGCAGCGCCGAACAGGCAAAACCCCAGCGCGCa                                                     >  7:152123/1‑148 (MQ=255)
                                                                                           gCAATCAGCGAGATGGTCACCATGATCATCGGCTGGCCCAGGGCGCGGATGATCTGGATATGGTTGAACTGCGGCCCGGCAAAGTCCGGGTTGAGCACCCCCGAACCGAAACTGGCAGCGCCGAACAGGCAAAACCCCAGCGCGCATa                                                   <  7:181394/148‑1 (MQ=255)
                                                                                               tCAGCGAGATGGTCACCATGATCATAGGCTGGCACAGGGCGCGGAAGAACTGGATATGGTTGAACTGCGGCCCGGCAAAGTCCGGCGTGAGCACCCCCGCACCGAAACTGGCAGCGCCGAACCGGCAAAGCCCCAGGGCGCATAGGAcc                                              >  2:15609/1‑149 (MQ=255)
                                                                                                  gCGAGATGGTCACCATGATCATCGGCTGGCCCAGGGCGCGGATGATCTGGATATGGTTGAACTGCGGCCCGGCAAAGTCCGGGTTGAGCACCCCCGAACCGAAACTGGCAGCGCCGAACAGGCAAAACCCCAGCGCGCATAGTACCTTg                                           >  3:263298/1‑149 (MQ=255)
                                                                                                     agaTGGTCACCATGATCATCGGCTGGCCCAGGGCGCGGATGATCTGGATATGGTTGAACTGCGGCCCGGCAAAGTCCGGGTTGAGCACCCCCGAACCGAAACTGGCAGCGCCGAACAGGCAAAACCCCAGCGCGCATAGTACCTTGGGc                                        <  2:287049/149‑1 (MQ=255)
                                                                                                              ccATGATCATCGGCTGGCCCAGGGCGCGGATGATCTGGATATGGTTGAACTGCGGCCCGGCAAAGTCCGGGTTGAGCACCCCCGAACCGAAACTGGCAGCGCCGAACAGGCAAAACCCCAgc                                                          <  8:270742/122‑1 (MQ=255)
                                                                                                               cATGATCATCGGCTGGCCCAGGGCGCGGATGATCTGGATATGGTTGAACTGCGGCCCGGCAAAGTCCGGGTTg                                                                                                          <  2:160381/73‑1 (MQ=255)
                                                                                                               cATGATCATCGGCTGGCCCAGGGCGCGGATGATCTGGATATGGTTGAACTGCGGCCCGGCAAAGTCCGGGTTg                                                                                                          >  1:160381/1‑73 (MQ=255)
                                                                                                               cATGATCATCGGCTGGCCCAGGGCGCGGATGATCTGGATATGGTTGAACTGCGGCCCGGCAAAGTCCGGGTTGAGCACCCCCGAACCGAAACTGGCAGCGCCGAACAGGCAAAACCCCAGCGCGCATAGTACCTTGGGCGATATCACCt                              >  3:143891/1‑149 (MQ=255)
                                                                                                                          gCTGTCCAAGGGCGCGGATGATCTGGATATGGTTGATCTGCGGCCAGGCAAAGACCGGCTTGAGAACCCCCGAACCGAAACTGGCAGCGCCGCACAGGCAAAACCCCAGCGCGCATAGTACCTTGGGCTCTATCCCCTTCCTCAGTTGc                   >  6:267052/1‑149 (MQ=255)
                                                                                                                          gCTGGCCCAGGGCGCGGATGATCTGGATATGGTTGAACTGCGGCCCGGCAAAGTACGGGTTGAGCACCCCCGAACCGAAACTGGCAGCGCCGAACAGGCAAAACCCCAGCGCGCATAGTACCTTGGGCGATATCACCTTCATCAGTTGc                   >  1:234039/1‑149 (MQ=255)
                                                                                                                              gCCCAGGGCGCGGATGATCTAGATATGGTTGAACTGCGGCCCGGCAAAGTCCGGGTTGAGCACCCCCGAACCGAAACTGGCAGCGCCGAACAGGCAAAACCCCAGCGCGCATAGTACCTTGGGCGATATCACCTTCATCAGTTGCg                  >  7:232255/1‑146 (MQ=255)
                                                                                                                              gCCCAGGGCGCGGATGATCTAGATATGGTTGAACTGCGGCCCGGCAAAGTCCGGGTTGAGCACCCCCGAACCGAAACTGGCAGCGCCGAACAGGCAAAACCCCAGCGCGCATAGTACCTTGGGCGATATCACCTTCATCAGTTGCg                  <  8:232255/146‑1 (MQ=255)
                                                                                                                                    ggCGCGGATGATCTGGATATGGTTGAACTGCGGCCCGGCAAAGTCCGGGTTGAGCACCCCCGAACCGAAACTGGCAGCGCCGAACAGGCAAAACCCCAGCGCGCATAGTACCTTGGGCGATATCACCTTCATCAGTTGCGGCACCAGCg         <  1:205799/149‑1 (MQ=255)
                                                                                                                                       gcgGATGATCTGGATATGGTTGAACTGCGGCCCGGCAAAGTCCGGGTTGAGCACCCCCGAACCGAAACTGGCAGCGCCGAACAGGCAAAACCCCAGCGCGCATAGTACCTTGGGCGATATCACCTTCATCAGTTGCGGCACCAGCGGGa      <  4:245946/149‑1 (MQ=255)
                                                                                                                                          gatgatCTGGATATGGTTGAACTGCGGCCCGGCAAAGTCCGGGTTGAGCACCCCCGAACCGAAACTGGCAGCGCCGAACAGGCAAAACCCCAGCGCGCATAGTACCTTGGGCGATATCACCTTCATCAGTTGCGGCACCAGCGGGATCa   <  5:125996/149‑1 (MQ=255)
                                                                                                                                           atgatCTGGATATGGTTGAACTGCGGCCCGGCAAAGTCCGGGTTGAGCACCCCCGAACCGAAACTGGCAGCGCCGAACAGGCAAAACCCCAGCGCGCATAGTACCTTGGGCGATATCACCTTCATCAGTTGCGGCACCAGCGGGATCAg  <  8:287707/149‑1 (MQ=255)
                                                                                                                                                  ggATATGGTTGAACTGCGGCCCGGCAAAGTCCGGGTTGAGCACCCCCGAACCGAAACTGGCAGCGCCGAACAGGCAAAACCCCAGCGCGCAt                                                    >  8:263322/1‑92 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
CAAGGCAATACCGATGGCCCCGCCCAGGTTACGCAGAATATTGAACAGGCTCGACGCGGACCCGGCATCCTGTTGTTGGATATACGCGGTGGCAATCAGCGAGATGGTCACCATGATCATCGGCTGGCCCAGGGCGCGGATGATCTGGATATGGTTGAACTGCGGCCCGGCAAAGTCCGGGTTGAGCACCCCCGAACCGAAACTGGCAGCGCCGAACAGGCAAAACCCCAGCGCGCATAGTACCTTGGGCGATATCACCTTCATCAGTTGCGGCACCAGCGGGATCAG  >  NC_002947/2237321‑2237608

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: ‑