Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,024,150 G→A 8.1% P120L (CCG→CTG)  emrB ← multidrug transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,024,1500GA 8.1% 140.5 / 3.1 51P120L (CCG→CTG) emrBmultidrug transporter permease
Reads supporting (aligned to +/- strand):  ref base G (26/19);  new base A (2/2);  total (28/23)
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.86e-01

GCTATCGGTAATCCAGCCCCCAAGAATCGGCCCGGCAATCGGCGCCACCACCGTGACCATCGCCAGCAGCGCCAGTGCCATCCCCCGTTTTGCCGGGGGGTAGACAGCAATCAGCAGGGTCTGGGTCATCGGATAAAGCGGTCCGGCGACCACGCCCTGCAAGACACGAAAGCCCACCAGCTCGGGCATCGACTGGGCAATACCGCACAGGAACGACGCCAGCACGAACAGCAGCGTGGCCCAGATGAACAGCTTCACTTCACCAAAGCGCCGGCTCAGCCAGCCG  >  NC_002947/4024011‑4024296
                                                                                                                                           |                                                                                                                                                  
gCTATCGGTAATCCAGCCCCCAAGAATCGGCCCGGCAATCGGCGCCACCACCGTGACCATCGCCAGCAGCGCCAGTGCCATCCCCCGTTTTGCCGGGGGGTAGACAGCAATCAGCAGGGTCTGGGTCATCGGATAAAGCGGTCCGGCGa                                                                                                                                           >  6:129742/1‑149 (MQ=255)
    tCGGTAATCCAGCCCCCAAGAATCGGCCCGGCAATCGGCGCCACCACCGTGACCAACGCCAGCAGCGCCAGTGCCATCCCACGTTTTGCCGGGGGGTAGACAGCAATAAGCAGCGACTGGGTCATCGGATAAAGCGGTCCGGCTACCAc                                                                                                                                       >  3:69921/1‑149 (MQ=255)
     cGGTAATCCAGCCCCCAAGAATCGGCCCGGCAATCGGCGCCCCCACCGTGACCATCGCCAGCAGCGCCAGTGCCATCCCCCGGTTTGCCGGGGGGTAGACCGCACTCAGCAGGGGCGGGGTCATCGGATAAAGCGGTCCGGCGACCACg                                                                                                                                      >  1:246763/1‑149 (MQ=255)
          aTCCAGCCCCCAAGAATCGGCCCGGCAATCGGCGCCACCACCGTGACCATCGCCAGCAGCGCCAGTGCCATCCCCCGTTTTGCCGGGGGGTAGACAGCAATCAGCAGGGTCTGGGTCATCGGATAAAGCGGTCCGGCGACCACGCCCt                                                                                                                                  <  6:304318/148‑1 (MQ=255)
               gCCCCCAAGAATCGGCCCGGCAATCGGCGCCACCACCGTGACCATCGCCAGCAGCGCCAGTGCCATCCCCCGTTTTGCCGGGGGGTAGACAGCAATCAGCAGGGTCTGGGTCATCGGATAAAGCGGTCCGGCGACCACGCCCTGCAAGa                                                                                                                            <  5:129742/149‑1 (MQ=255)
                  cccAAGAATCGGCCCGGCAATCGGCGCCACCACCGTGACCATCGCCAGCAGCGCCAGTGCCATCCCCCGTTTTGCCGGGGGGTAGACAGCAATCAGCAGGGTCTGGGTCATCGGATAAAGCGGTCCGGCGACCACGCCCTGCAAGacac                                                                                                                         <  7:91276/149‑1 (MQ=255)
                      aGAATCGGCCCGGCAATCGGCGCCACCACCGTGACCATCGCCAGCAGCGCCAGTGCCATCCCCCGTTTTGCCGGGGGGTAGACAGCAATCAGCAGGGTCTGGGTCATCGGATAAAGCGGTCCGGCGACCACGCCCTGCAAGACACGaaa                                                                                                                     <  7:225801/149‑1 (MQ=255)
                              cccGGCAATCGGCGCCACCACCGTGACCATCGCCAGCAGCGCCAGTGCCATCCCCCGTTTTGCCGGGGGGTAGACAGCAATCAGCAGGGTCTGGGTCATCGGATAAAGCGGTCCGGCGACCACGCCCTGCAAGACACGAAAGCccac                                                                                                               >  4:239074/1‑147 (MQ=255)
                                         gcgcCACCAGCGTGACCATCGCCAGCAGCGCCAGTGCCATCCCCCGTTTTGCCGGGGGGTAGACAGCAATCAGCAGGGTCTGGGTCATCGGATAAAGCGGTCCGGCGACCACGCCCt                                                                                                                                  <  5:298099/117‑1 (MQ=255)
                                         gcgcCACCACCGTGACCATCGCCAGCAGCGCCAGTGCCATCCCCCGTTTTGCCGGGGGGTAGACAGCAATCAGCAGGGTCTGGGTCATCGGATAAAGCGGTCCGGCGACCACGCCCt                                                                                                                                  >  6:298099/1‑117 (MQ=255)
                                           gcCACCACCGTGACCATCGCCAGCAGCGCCAGTGCCATCCCCCGTTTTGCCGGGGGGTAGACAGCAATCAGCAGGGTCTGGGTCATCGGATAAAGCGGTCCGGCGACCACGCCCTGCAAGa                                                                                                                            <  3:270000/121‑1 (MQ=255)
                                           gcCACCACCGTGACCATCGCCAGCAGCGCCAGTGCCATCCCCCGTTTTGCCGGGGGGTAGACAGCAATCAGCAGGGTCTGGGTCATCGGATAAAGCGGTCCGGCGACCACGCCCTGCAAGa                                                                                                                            >  4:270000/1‑121 (MQ=255)
                                            ccaccaCCGTGACCATCGCCAGCAGCGCCAGTGCCATCCCCCGTTTTGCCGGGGGGTAGACAGCAATCAGCAGGGTCTGGGTCATCGGATAAAGCAGTCCGGCGACCACGCCCTGCAAGACACGAAAGCCCACCAGCTCgg                                                                                                       >  5:165328/1‑141 (MQ=255)
                                            ccaccaCCGTGACCATCGCCAGCAGCGCCAGTGCCATCCCCCGTTTTGCCGGGGGGTAGACAGCAATCAGCAGGGTCTGGGTCATCGGATAAAGCAGTCCGGCGACCACGCCCTGCAAGACACGAAAGCCCACCAGCTCgg                                                                                                       <  6:165328/141‑1 (MQ=255)
                                              accaccGTGACCATCGCCAGCAGCGCCAGTGCCATCCCCCGTTTTGCCGGGGGGTAGACAGCAATCAGCAGTGTCTGGGTCATCGGATAAAGCGGTCCGGCGACCACGCCCTGCAAGACACGAAAGCCCACCAGCTCGGGCATCGACTg                                                                                             <  6:336584/149‑1 (MQ=255)
                                              accaccGTGACCATCGCCAGCAGCGCCAGTGCCATCCCCCGTTTTGCCGGGGGGTAGACAGCAATCAGCAGGGTCTGGGTCATCGGATAAAGCGGTCCGGCGACCACGCCCTGCAAGACACGAAAGCCCACCAGCTCGGGCATCGACTg                                                                                             >  4:62634/1‑149 (MQ=255)
                                                caccGTGACCATCGCCAGCAGCGCCAGTGCCATCCCCCGTTTTGCCGGGGGGTAGACAGCAATCAGCAGGGTCTGGGTCATCGGATAAAGCGGTCCGGTGACCACGCCCTGCAAGACACGAAAGCCCACCAGCTCGGGCATCGACTggg                                                                                           >  2:288963/1‑149 (MQ=255)
                                                caccGTGACCATCGCCAGCAGCGCCAGTGCCATCCCCCGTTTTGCCGGGGGGTAGACAGCAATCAGCAGGGTCTGGGTCATCGGATAAAGCGGTCCGGCGACCACGCCCTGCAAGACACGaaa                                                                                                                     <  4:158029/123‑1 (MQ=255)
                                                caccGTGACCATCGCCAGCAGCGCCAGTGCCATCCCCCGTTTTGCCGGGGGGTAGACAGCAATCAGCAGGGTCTGGGTCATCGGATAAAGCGGTCCGGCGACCACGCCCTGCAAGACACGaaa                                                                                                                     >  3:158029/1‑123 (MQ=255)
                                                  ccGTGACCATCGCCAGCAGCGCCAGTGCCATCCCCCGTTTTGCCGGGGGGTAGACAGCAATCAGCAGGGTCTGGGTCATCGGATAAATCAGTCCGGCGACCACGCCCTGCaa                                                                                                                              <  7:215789/112‑1 (MQ=255)
                                                  ccGTGACCATCGCCAGCAGCGCCAGTGCCATCCCCCGTTTTGCCGGGGGGTAGACAGCAATCAGCAGGGTCTGGGTCATCGGATAAAGCAGTCCGGCGACCACGCCCTGCaa                                                                                                                              >  8:215789/1‑112 (MQ=255)
                                                    gTGACCATCGCCAGCAGCGCCAGTGCCATCCCCCGTTTTGCCGGGGGGTAGACAGCAATCAGCAGGGTCTGGGTCATCGGATAAAGCGGTCCGGCGACCACGCCCTGCAAGACACGAAAGCCCACCAGCTCGGGCATCGACTGGGCAAt                                                                                       <  5:209438/149‑1 (MQ=255)
                                                     tGACCATCGCCAGCAGCGCCAGTGCCATCCCCCGTTTTGCCGGGGGGTAGACAGCAATCAGCAGGGTCTGGGTCATCGGATAAAGCGGTCCGGCGACCACGCCCTGCAAGACACGAAAGCCCACCAGCTCGGGCATCGACTGGGCAATa                                                                                      <  3:95324/149‑1 (MQ=255)
                                                        ccATCGCCAGCAGCGCCAGTGCCATCCCCCGTTTTGCCGGGGGGTAGACAGCAATCAGCAGGGTCTGGGTCATCGGATAAAGCGGTCCGGCGACCACGCCCTGCAAGACACGAAAGCCCACCAGCTCGGGCATCGACTGGGCAATACCg                                                                                   <  3:239074/149‑1 (MQ=255)
                                                         cATCGCCAGCAGCGCCAGTGCCATCCCCCGTTTTGCCGGGGGGTAGACAGCAATCAGCAGGGTCTGGGTCATCGGATAAAGCGGTCCGGCGACCACGCCCTGCAAGACACGAAAGCCCACCAGCTCGGGCATCGACTGGGCAATACCg                                                                                   >  4:111008/1‑148 (MQ=255)
                                                                agcagcGCCAGTGCCATCCCCCGTTTTGCCGGCGGGTAGACAGCAATCAGCAGGGTCTGGGTCATCGGATAAAGCGGTCCGGCGACCACGCCCTGCAAGACACGAAAGCCCACCAGCTCGGGCATCGACTGGGCAATAACGCACAGGaa                                                                           >  1:174050/1‑149 (MQ=255)
                                                                    gcgcCAGTGCCATCCCCCGTTTTGCCGGGGGGTAGACAGCAATCAGCAGGGTCTGGGTCATCGGATAAAGCGGTCCGGCGACCACGCCCTGCAAGACACGAAAGCCCACCAGCTCGGGCATCGACTGGGCAATACCGCACAGGAacgac                                                                       <  4:290895/149‑1 (MQ=255)
                                                                    gcgcCAGTGCCATCCCCCGTTTTGCCGGGGGGTAGACAGCAATCAGCAGGGTCTGGGTCATCGGATAAAGCGGTCCGGCGACCACGCCCTGCAAGACACGAAAGCCCACCAGCTCGGGCATCGACTGGGCAATACCGCACAGGAacgac                                                                       <  8:46658/149‑1 (MQ=255)
                                                                           tGCCATCCCCCGTTTTGCCGGGGGGTAGACAGCAATCAGCAGGGTCTGGGTCATCGGATAAAGCGGTCCGGCGACCACGCCCTGCAAGACACGAAAGCCCACCCGCTCGGGCATCGACTGGGCAATACCGCACAGGAACGACGCCAGCa                                                                >  3:121990/1‑149 (MQ=255)
                                                                           tGCCATCCCCCGTTTTGCCGGGGGGTAGACAGCAATCAGCAGGGTCTGGGTCATCGGATAAAGCGGTCCGGCGACCACGCCCTGCAAGACACGAAAGCCCACCAGCTCGGGCATCGACTGGGCAATACCGCACAGGAACGACGCCCGCa                                                                >  5:72397/1‑149 (MQ=255)
                                                                            ggcATCCCCCGTTTTGCCGGGGGGTAGACATCAATCAGCAGGGTCTGGTTCATCGGATATATTTGTCCGGCGACCACGCCCTGCAAGACACGAAAGCCCACCAGCTCTGGCATCGACTGGGCAATACCGCACAGGAACGACGCCAGCa                                                                <  4:69921/146‑1 (MQ=255)
                                                                            gCCATCCCCCGTTTTGCCGGGGGGTAGACAGCAATCAGCAGGGTCTGGGTCATCGGATAAAGCGGTCCGGCGACCACGCCCTGCAAGACACGAAAGCCCACCAGCTCGGGCATCGACTGGGCAATACCGCACAGGATCGACGCCAGCa                                                                <  6:109685/148‑1 (MQ=255)
                                                                                tCCCCCGTTTTGCCGGGGGGTAGACAGCAATCAGCAGGGTCTGGGTCATCGGATAAAGCGGTCCGGCGACCACGCCCTGCAAGACACGAAAGCCCACCAGCTCGGGCATCGACTGGGCAATACCGCACAGGAACGACGCCAGCACGAAc                                                           <  2:30895/149‑1 (MQ=255)
                                                                                   cccGTTTTGCCGGGGGGTAGACAGCAATCAGCAGGGTCTGGGTCATCGGATAAAGCGGTCCGGCGACCACGCCCTGCAAGACACGAAAGCCCACCAGCTCgg                                                                                                       >  2:287950/1‑102 (MQ=255)
                                                                                   cccGTTTTGCCGGGGGGTAGACAGCAATCAGCAGGGTCTGGGTCATCGGATAAAGCGGTCCGGCGACCACGCCCTGCAAGACACGAAAGCCCACCAGCTCGGGCATCGACTGGGCAATACCGCACAGGaa                                                                           >  7:8880/1‑130 (MQ=255)
                                                                                   cccGTTTTGCCGGGGGGTAGACAGCAATCAGCAGGGTCTGGGTCATCGGATAAAGCGGTCCGGCGACCACGCCCTGCAAGACACGAAAGCCCACCAGCTCGGGCATCGACTGGGCAATACCGCACAGGaa                                                                           <  8:8880/130‑1 (MQ=255)
                                                                                             cGGGGGGTAGACAGCAATCAGCAGGGTCTGGGTCATCGGATAAAGCGGTCCGGCGACCACGCCCTGCAAGACACGAAAGCCCACCAGCTCGGGCATCGACTGGGCAATACCGCACAGGAACGACGCCAGCACGAACagcagc                                                     <  6:208333/142‑1 (MQ=255)
                                                                                             cGGGGGGTAGACAGCAATCAGCAGGGTCTGGGTCATCGGATAAAGCGGTCCGGCGACCACGCCCTGCAAGACACGAAAGCCCACCAGCTCGGGCATCGACTGGGCAATACCGCACAGGAACGACGCCAGCACGAACAGCAGCgagctcg                                              >  5:208333/1‑143 (MQ=255)
                                                                                              ggggggTAGACAGCAATCATCAGGGTCTGGTTCATCGGATAAAGCTGTCCGGCGACCACGCCCTGCAAGACACGAAAGCCCACCAGCTCGGGCATCGACTTGGCAATACCGCACAGGAACGACCCCAGCACGAACAGCAGCGTGGCCCa                                             <  7:79796/149‑1 (MQ=255)
                                                                                                ggggTAGACCGCAATCAGCAGGGTCTGGGTCATCGGATAAAGCGGTCCGGCGACCACGCCCTGCAAGACACGAAAGCCCACCAGCTCGGGCATCGACTGGGCAATACCGCACAGGAACGACGCCAGCACGAACAGCAGCGTGGCCCAGa                                           <  2:246763/149‑1 (MQ=255)
                                                                                                  ggTAGACAGCAATCAGCAGGGTCTGGGTCATCGGATAAAGCGGTCCGGCGACCACGCCCTGCAAGACACGAAAGCCCACCAGCTCGGGCATCGACTGGGCAATACCGCACAGGAACGACGCCAGCACGAACAGCAGCGTGGCCCAGAt                                          >  2:125988/1‑148 (MQ=255)
                                                                                                    tAGACAGCAATCAGCAGGGTCTGGGTCATCGGATAAAGCGGTCCGGCGACCACGCCCTGCAAGACACGAAAGCCCACCAGCTCGGGCATCGACTGGGCAATACCGCACAGGAACGACGCCAGCACGAACAGCAGCGTGGCCCAGATGaa                                       >  4:212491/1‑149 (MQ=255)
                                                                                                                 gcagGGTCTGGGTCATCGGATAAAGCGGTCCGGCGACCACGCCCTGCAAGACACGAAAGCCCACCAGCTCGGGCATCGACTGGGCAATACCGCACAGGAACGACGCCAGCACGAACAGCAGCGTGGCCCAGATGAACAGcttcacttca                          >  7:248709/1‑149 (MQ=255)
                                                                                                                       tCTGGGTCATCGGATAAAGCGGTCCGGCGACCACGCCCTGCAAGAGACGAAAGCCCACCAGCTCGGGCATCGACTGGGCAATACCGCACAGGAACGACGCCAGCACGAACAGCAGCGTGGCCCAGATGAACAGCTTCACTTCACCAAAg                    >  1:243687/1‑149 (MQ=255)
                                                                                                                         tGGGTCATCGGATAAAGCGGTCCGGTGACCACGCCCTGCAAGACACGAAAGCCCACCAGCTCGGGCATCGACTGGGCAATACCGCACAGGAACGACGCCAGCACGAACAGCAGCGTGGCCCAGATGAACAGCTTCACTTCACCAAAgcg                  >  6:77738/1‑149 (MQ=255)
                                                                                                                         tGGGTCATCGGATAAAGCGGTCCGGCGACCACGCCCTGCAAGACACGAAAGCCCACCAGCTCGGGCATCGACTGGGCAATACCGCACAGGAACGACGCCAGCACGAACAGCAGCGTGGCCCAGATCAAAAGCATCACTTCACCAAAgcg                  >  4:125264/1‑149 (MQ=255)
                                                                                                                                tCGGATAAAGCGGTCCGGCGACCACGCCCTGCAAGACACGAAAGCCCACCAGCTCGGGCATCGGCTGGGCAATACCGCACAGGAACGACGCCAGCACGAACAGCAGCGTGGCCCAGATGAACAGCTTCACTTCACCAAAGCGCCGGCt            >  6:41934/1‑148 (MQ=255)
                                                                                                                                      aaaGCGGTCCGGCGACCACGCCCTGCAAGACACGAAAGCCCACCAGCTCGGGCATCGACTGGGCAATACCGCACAGGAACGACGCCAGCACGAACAGCAGCGTGGCCCAGATGATCAGCTTCACTTCACCAAAGCGCCGGCTcagccag     >  6:110017/1‑149 (MQ=255)
                                                                                                                                      aaaGCGGTCCGGCGACCACGCCCTGCAAGACACGAAAGCCCACCAGCTCGGGCATCGACTGGGCAATACCGCACAGGAACGACGCCAGCACGAACAGCAGCGTGGCCCAGATGAACAGCTTCACTTCACCAAAGCGCCGGCTcagccag     >  1:246891/1‑149 (MQ=255)
                                                                                                                                        aGCGGTCCGGCGACCACGCCCTGCAAGACACGAAAGCCCACCAGCTCGGGCATCGACTGGGCAATACCGCACAGGAACGACGCCAGCACGAACAGCAGCGTGGCCCAGATGAACAGCTTCACTTCACCAAAGCGCCGGCTCagccagcc   >  3:250597/1‑149 (MQ=255)
                                                                                                                                         gCGGTCCGGCGACCACGCCCTGCAAGACACGAAAGCCCACCAGCTCGGGCATCGACTGGGCAATACCGCACAGGAACGACGCCAGCACGAACAGCAGCGTGGCCCAGATGAACAGCTTCACTTCACCAAAGCGCCGGCTCAGCCAGCCg  <  4:250597/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                  
GCTATCGGTAATCCAGCCCCCAAGAATCGGCCCGGCAATCGGCGCCACCACCGTGACCATCGCCAGCAGCGCCAGTGCCATCCCCCGTTTTGCCGGGGGGTAGACAGCAATCAGCAGGGTCTGGGTCATCGGATAAAGCGGTCCGGCGACCACGCCCTGCAAGACACGAAAGCCCACCAGCTCGGGCATCGACTGGGCAATACCGCACAGGAACGACGCCAGCACGAACAGCAGCGTGGCCCAGATGAACAGCTTCACTTCACCAAAGCGCCGGCTCAGCCAGCCG  >  NC_002947/4024011‑4024296

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